Los métodos espectrales son una clase de técnicas usadas en matemáticas aplicadas para la solución numérica de ciertas ecuaciones diferenciales , a veces usando la Transformada Rápida de Fourier . La idea es representar la solución de ecuaciones diferenciales como una suma de algunas " funciones base " (como cómo las series de Fourier son la suma de sinusoides ), y luego elegir los coeficientes en la suma que mejor satisfagan las ecuaciones dadas.
Los métodos espectrales y los métodos de elementos finitos están estrechamente relacionados y se basan en las mismas ideas. La principal diferencia es que los métodos espectrales utilizan funciones de base distintas de cero en todo el dominio de definición, mientras que los métodos de elementos finitos utilizan funciones de base distintas de cero solo en pequeños subdominios. En otras palabras, los métodos espectrales adoptan un enfoque global , mientras que los métodos de elementos finitos utilizan un enfoque local . En parte por esta razón, los métodos espectrales tienen excelentes propiedades de la llamada "convergencia exponencial", que es la más rápida posible si la solución es suave . Sin embargo, no se conoce ningún método espectral tridimensional de un dominio de conteo completo (la onda de choque no es suave) [1] . El método de elementos finitos, en el que la potencia del elemento es muy alta o aumenta a medida que disminuye el parámetro de red h , a veces se denomina método de elementos espectrales .
Los métodos espectrales se pueden usar para resolver ecuaciones diferenciales ordinarias (EDO), ecuaciones diferenciales parciales y problemas de valores propios que involucran ecuaciones diferenciales. Cuando se aplican métodos espectrales a ecuaciones diferenciales parciales dependientes del tiempo, la solución generalmente se escribe como la suma de funciones de base con coeficientes dependientes del tiempo. Sustituyendo tal suma en una ecuación diferencial parcial se obtiene un sistema de ecuaciones diferenciales ordinarias en los coeficientes, que se puede resolver usando cualquiera. método numérico de ecuaciones diferenciales ordinarias . El problema de encontrar valores propios para ecuaciones diferenciales ordinarias se reduce de manera similar al problema de encontrar los valores propios de una matriz.
Steven Orsaga ha desarrollado métodos espectrales en una serie de artículos. Desde 1969, se han desarrollado para métodos de Fourier para problemas geométricos periódicos, métodos espectrales polinómicos para problemas geométricos finitos e ilimitados, métodos pseudoespectrales para problemas fuertemente no lineales, métodos espectrales iterativos para resolver problemas de estado estacionario y otros problemas. La implementación del método espectral generalmente termina con una colocación o el método Galerkin , o el enfoque Tau[ aclarar ] .
Los métodos espectrales son computacionalmente menos costosos que los métodos de elementos finitos, pero se vuelven menos precisos para problemas con geometrías complejas y coeficientes discontinuos. Este aumento en el error es una consecuencia. Fenómenos de Gibbs
Aquí asumimos una comprensión del cálculo multivariante básico y las series de Fourier . Si g(x,y) es una función compleja conocida de dos variables reales y g es periódica en x e y ( g(x,y)=g(x+2π,y)=g(x,y+2π)) , entonces nos interesa encontrar una función f(x,y) tal que
para todo x,ydonde la expresión de la izquierda es la segunda derivada parcial de f con respecto a xey, respectivamente. Esta es la ecuación de Poisson y puede interpretarse físicamente como algún tipo de problema de transferencia de calor o problema en la teoría del potencial, entre otras posibilidades.
Si escribimos f y g como series de Fourier
Y sustituimos en la ecuación diferencial, obtenemos la ecuación:
Hemos cambiado la derivación parcial por la sumatoria, lo cual es legal si suponemos, por ejemplo, que f tiene una segunda derivada continua. De acuerdo con el teorema de unicidad para la expansión de Fourier, entonces debemos igualar los coeficientes de Fourier elemento por elemento, lo que da
(*)que es una fórmula explícita para los coeficientes de Fourier a j , k .
Con condiciones de frontera periódicas , la ecuación de Poisson tiene solución sólo si b 0 , 0 = 0 . Así, podemos elegir libremente un 0 , 0 . Esto corresponde a la elección de la constante de integración.
Para convertir esto en un algoritmo, solo se calcula un número finito de frecuencias. Esto da un error que puede demostrarse que es proporcional a , donde y es la frecuencia procesada más alta.
AlgoritmoDado que estamos interesados en una ventana finita de frecuencias (de tamaño n , digamos ), esto se puede hacer usando el algoritmo Fast Fourier Transform . Por lo tanto, globalmente, el algoritmo se ejecuta en tiempo O ( n log n ).
Deseamos resolver la ecuación transitoria no lineal de Burgers utilizando un enfoque especial.
Si se da en el dominio periódico , encontramos , tal que
donde ρ es el coeficiente de viscosidad . esto se convierte en
donde corresponde al producto escalar . La integración por partes y el uso de la periodicidad da
Para aplicar el método de Fourier- Galerkin , elegimos
y
donde _ Esto reduce el problema a encontrar , tal que
Usando la relación de ortogonalidad , donde es el delta de Kronecker , simplificamos los tres elementos anteriores para cada
Recopilamos tres términos para cada uno y obtenemos
Dividir por y finalmente obtener
Con las condiciones iniciales de la transformada de Fourier y la definición de , este par de ecuaciones diferenciales ordinarias se pueden integrar a lo largo del tiempo (usando, por ejemplo, la técnica de Runge-Kutta ) para encontrar una solución. El término no lineal es una convolución y existen varias técnicas basadas en transformaciones para calcularlo de manera eficiente.
Se puede demostrar que si es infinitamente derivable, entonces un algoritmo numérico que usa la transformada rápida de Fourier converge más rápido que cualquier polinomio en una red de tamaño h. Es decir, para cualquier n>0 existe , tal que el error es menor para todos los valores suficientemente pequeños de . Decimos que el método espectral tiene un orden para cualquier n>0.
Dado que el método de elementos espectrales es un método de elementos finitos de muy alto orden , existe una similitud en las propiedades de convergencia. Sin embargo, el método espectral se basa en la expansión de valores propios de un problema de valor límite específico, mientras que el método de elementos finitos no utiliza esta información y funciona para problemas de valor límite elípticos arbitrarios .
ecuaciones diferenciales | Métodos para resolver|||||
---|---|---|---|---|---|
Métodos de cuadrícula |
| ||||
Métodos sin cuadrícula |