Métodos de Diferencias Finitas

Los métodos de diferencias finitas aproximan las derivadas mediante cocientes de diferencias en una cuadrícula, transformando una ecuación diferencial en un sistema de ecuaciones algebraicas para los valores de la solución en los puntos de la cuadrícula.

Definition

Un método de diferencias finitas es una discretización de una ecuación diferencial en la que las derivadas se reemplazan por cocientes de diferencias de la incógnita evaluados en una cuadrícula estructurada, produciendo ecuaciones algebraicas cuya solución aproxima la solución de la ecuación diferencial en los puntos de la cuadrícula.

Scope

Este tema abarca la construcción de aproximaciones de diferencias a partir de expansiones de Taylor, la discretización de EDP elípticas, parabólicas e hiperbólicas, esquemas de avance temporal explícitos e implícitos (como Euler hacia adelante, Euler hacia atrás y Crank-Nicolson), el análisis de estabilidad de von Neumann y el marco de consistencia-estabilidad-convergencia especializado en esquemas de diferencias.

Core questions

• ¿Cómo se derivan las aproximaciones de diferencias precisas a las derivadas y se cuantifica su error de truncamiento?
• ¿En qué se diferencian los esquemas de avance temporal explícitos e implícitos en cuanto a estabilidad y costo?
• ¿Cómo determina el análisis de von Neumann la estabilidad de un esquema de diferencias?
• ¿Cómo establece el tipo de ecuación el esquema apropiado y las posibles restricciones en el tamaño del paso?

Key theories

Consistencia, estabilidad y convergencia

Un esquema de diferencias es consistente si su error de truncamiento se anula a medida que se refina la cuadrícula y es estable si los errores no crecen de forma ilimitada; según el teorema de equivalencia de Lax, estos dos aspectos garantizan la convergencia a la solución verdadera para problemas lineales bien planteados.

Análisis de estabilidad de von Neumann

Descomponer el error en modos de Fourier en una cuadrícula uniforme reduce la estabilidad a acotar un factor de amplificación para cada modo; el esquema es estable cuando ningún modo se amplifica, lo que produce condiciones explícitas de tamaño de paso como los límites de difusión y CFL.

Mechanisms

Las fórmulas de diferencias se construyen combinando expansiones de Taylor en puntos de cuadrícula vecinos para cancelar términos de orden bajo y aislar una derivada, con el término principal restante dando el error de truncamiento y el orden del método. Para problemas dependientes del tiempo, los esquemas explícitos actualizan cada nuevo valor directamente a partir de los antiguos, pero deben respetar un límite de estabilidad (un límite de número de difusión para ecuaciones parabólicas, la condición CFL para las hiperbólicas), mientras que los esquemas implícitos como Crank-Nicolson acoplan los nuevos valores en un sistema lineal que es incondicionalmente estable pero requiere una resolución en cada paso. El análisis de von Neumann sustituye los modos de Fourier para probar la estabilidad y determinar estos límites.

Clinical relevance

Los métodos de diferencias finitas se utilizan ampliamente para problemas en dominios regulares y cuadrículas estructuradas: conducción de calor y difusión, propagación de ondas y modelado sísmico, electromagnetismo computacional (el método de diferencias finitas en el dominio del tiempo) y fijación de precios de opciones a través de la ecuación de Black-Scholes; su simplicidad y facilidad de extensión a órdenes superiores los convierten en una primera opción cuando la geometría es simple.

History

La base matemática fue establecida por el artículo de Courant-Friedrichs-Lewy de 1928 sobre ecuaciones de diferencias para EDP; el análisis de estabilidad de von Neumann en tiempos de guerra y el teorema de equivalencia de Lax de la década de 1950 establecieron la teoría moderna, y los métodos de diferencias siguen siendo un pilar de la física y la ingeniería computacionales.

Key figures

• Richard Courant
• Kurt Friedrichs
• Hans Lewy
• John von Neumann
• Randall J. LeVeque

Seminal works

• leveque2007
• morton2005

Frequently asked questions

¿Cuál es la diferencia entre un esquema explícito y uno implícito?

Un esquema explícito calcula cada nuevo valor de la cuadrícula directamente a partir de valores conocidos, pero es estable solo para pasos de tiempo pequeños, mientras que un esquema implícito resuelve un sistema acoplado para todos los valores nuevos a la vez, lo que permite pasos de tiempo estables mucho mayores a costa de una resolución lineal por paso.

¿Por qué podrían preferirse las diferencias finitas a los elementos finitos?

En geometrías simples y regulares, las diferencias finitas son fáciles de implementar, económicas y sencillas de hacer de orden superior. Los elementos finitos se vuelven ventajosos principalmente cuando el dominio tiene una forma compleja o el problema tiene una formulación variacional natural.