EDO rígidas y estabilidad

Las ecuaciones diferenciales rígidas contienen procesos que evolucionan en escalas de tiempo muy separadas, por lo que los métodos explícitos se ven obligados a dar pasos impracticamente pequeños para la estabilidad; su solución eficiente requiere métodos implícitos con fuertes propiedades de estabilidad.

Definition

Una ecuación diferencial se denomina rígida cuando admite componentes de solución que decaen en escalas de tiempo muy diferentes, de modo que la estabilidad numérica, en lugar de la precisión, dicta el tamaño del paso; la teoría de la estabilidad analiza qué métodos pueden dar pasos grandes sin crecimiento de errores.

Scope

Este tema cubre el fenómeno y la definición informal de rigidez, la ecuación de prueba lineal y la región de estabilidad absoluta, los conceptos de estabilidad A, estabilidad A(alfa) y estabilidad L, por qué los métodos explícitos fallan en problemas rígidos, y los métodos implícitos —Runge-Kutta implícito y fórmulas de diferenciación hacia atrás— que los resuelven.

Core questions

• ¿Qué hace que un problema sea rígido y por qué derrota a los métodos explícitos?
• ¿Cómo se define la región de estabilidad absoluta a través de la ecuación de prueba lineal?
• ¿Qué requieren la estabilidad A y la estabilidad L, y por qué son importantes para los problemas rígidos?
• ¿Qué métodos proporcionan la estabilidad necesaria para sistemas rígidos y diferenciales-algebraicos?

Key theories

Estabilidad absoluta y la ecuación de prueba

La aplicación de un método a la ecuación de prueba lineal escalar produce un factor de amplificación; el conjunto de productos de tamaño de paso por valor propio para los cuales este factor tiene una magnitud de como máximo uno es la región de estabilidad absoluta del método, que debe contener los valores propios rígidos del problema para permitir pasos grandes.

Estabilidad A y estabilidad L

Un método es A-estable si su región de estabilidad contiene todo el semiplano izquierdo, por lo que es estable para todos los modos de decaimiento independientemente del tamaño del paso, y L-estable si además amortigua completamente los modos muy rígidos; estas propiedades distinguen los métodos implícitos adecuados para problemas rígidos.

Mechanisms

En un problema rígido, el modo de decaimiento más rápido tiene un valor propio negativo grande; la región de estabilidad acotada de un método explícito obliga a que el tamaño del paso resuelva ese modo incluso mucho después de que haya desaparecido físicamente, lo que hace que el cálculo sea desesperadamente lento. Los métodos implícitos, como el método de Euler hacia atrás, los esquemas de Runge-Kutta implícitos y las fórmulas de diferenciación hacia atrás, tienen regiones de estabilidad que cubren el semiplano izquierdo (o la mayor parte de él), por lo que permanecen estables en pasos grandes y permiten que el tamaño del paso se elija solo por la precisión. Cada paso requiere entonces resolver un sistema algebraico (generalmente no lineal), típicamente mediante una iteración de Newton utilizando el jacobiano.

Clinical relevance

La rigidez es omnipresente en las redes de reacciones químicas, la combustión, los circuitos eléctricos, los sistemas de control y las discretizaciones de ecuaciones diferenciales parciales parabólicas mediante el método de las líneas; reconocer la rigidez y seleccionar un solucionador implícito apropiadamente estable es esencial para obtener resultados en un tiempo factible, y la mayoría del software de producción de EDO incluye detección automática de rigidez y conmutación.

History

La noción de rigidez fue identificada por Curtiss y Hirschfelder en 1952, y la teoría de estabilidad de apoyo —estabilidad A y las barreras de orden— fue desarrollada por Dahlquist; los códigos de fórmulas de diferenciación hacia atrás de Gear y, posteriormente, los métodos de Runge-Kutta implícitos de alto orden establecieron el conjunto de herramientas prácticas para problemas rígidos y diferenciales-algebraicos.

Key figures

• Germund Dahlquist
• C. William Gear
• Ernst Hairer
• Gerhard Wanner

Seminal works

• hairer1996
• iserles2008

Frequently asked questions

¿Qué hace exactamente que una EDO sea rígida?

La rigidez surge cuando el sistema tiene componentes que decaen mucho más rápido de lo que evoluciona la solución de interés. No existe una definición única y precisa, pero la característica práctica es que los métodos explícitos se ven obligados a usar pasos muy pequeños para la estabilidad, incluso cuando la precisión permitiría pasos grandes.

¿Por qué los problemas rígidos requieren métodos implícitos?

Los métodos implícitos pueden tener regiones de estabilidad que cubren todo el semiplano izquierdo (estabilidad A), por lo que permanecen estables con tamaños de paso grandes para modos de decaimiento rápido. Los métodos explícitos tienen regiones de estabilidad acotadas, lo que obliga a pasos diminutos y los hace poco prácticos para problemas rígidos.