En matemáticas, el método de Euler semi-implícita , también llamado simpléctico Euler , Euler semi-explícita , Euler-Cromer , y Newton-Stormer-Verlet (NSV) , es una modificación del método de Euler para resolver ecuaciones de Hamilton , un sistema de ordinario ecuaciones diferenciales que surge en la mecánica clásica . Es un integrador simpléctico y, por lo tanto, produce mejores resultados que el método de Euler estándar.
Configuración
El método de Euler semi-implícito se puede aplicar a un par de ecuaciones diferenciales de la forma
- [ cita requerida ]
donde f y g son funciones dadas. Aquí, x y v pueden ser escalares o vectores. Las ecuaciones de movimiento en la mecánica hamiltoniana toman esta forma si el hamiltoniano es de la forma
Las ecuaciones diferenciales deben resolverse con la condición inicial
El método
El método de Euler semi-implícito produce una solución discreta aproximada iterando
donde Δ t es el paso de tiempo y t n = t 0 + n Δ t es el tiempo después de n pasos.
La diferencia con el método de Euler estándar es que el método de Euler semi-implícito usa v n +1 en la ecuación para x n +1 , mientras que el método de Euler usa v n .
Aplicando el método con paso de tiempo negativo al cálculo de de y el reordenamiento conduce a la segunda variante del método de Euler semi-implícito
que tiene propiedades similares.
El Euler semi-implícito es un integrador de primer orden , al igual que el método de Euler estándar. Esto significa que comete un error global del orden de Δt. Sin embargo, el método de Euler semi-implícito es un integrador simpléctico , a diferencia del método estándar. Como consecuencia, el método de Euler semi-implícito casi conserva la energía (cuando el hamiltoniano es independiente del tiempo). A menudo, la energía aumenta de manera constante cuando se aplica el método estándar de Euler, lo que lo hace mucho menos preciso.
Alternar entre las dos variantes del método de Euler semi-implícito conduce en una simplificación a la integración Störmer- Verlet y en una simplificación ligeramente diferente a la integración de salto , aumentando tanto el orden del error como el orden de conservación de la energía. [1]
La región de estabilidad del método semi-implícito fue presentada por Niiranen [2] aunque el Euler semi-implícito fue llamado erróneamente Euler simétrico en su artículo. El método semi-implícito modela el sistema simulado correctamente si las raíces complejas de la ecuación característica están dentro del círculo que se muestra a continuación. Para las raíces reales, la región de estabilidad se extiende fuera del círculo para el cual el criterio es
Como puede verse, el método semi-implícito puede simular correctamente tanto sistemas estables que tienen sus raíces en el semiplano izquierdo como sistemas inestables que tienen sus raíces en el semiplano derecho. Esta es una clara ventaja sobre Euler hacia adelante (estándar) y Euler hacia atrás. El Euler hacia adelante tiende a tener menos amortiguamiento que el sistema real cuando las partes reales negativas de las raíces se acercan al eje imaginario y Euler hacia atrás puede mostrar que el sistema es estable incluso cuando las raíces están en el semiplano derecho.
Ejemplo
El movimiento de un resorte que satisface la ley de Hooke viene dado por
El Euler semi-implícito para esta ecuación es
Sustituyendo en la segunda ecuación con la expresión dada por la primera ecuación, la iteración se puede expresar en la siguiente forma matricial
y dado que el determinante de la matriz es 1, la transformación conserva el área.
La iteración conserva la energía modificada funcional exactamente, lo que lleva a órbitas periódicas estables (para un tamaño de paso suficientemente pequeño) que se desvían por de las órbitas exactas. La frecuencia circular exacta aumenta en la aproximación numérica por un factor de .
Referencias
- ^ Hairer, Ernst; Lubich, Christian; Wanner, Gerhard (2003). "Integración numérica geométrica ilustrada por el método Störmer / Verlet". Acta Numerica . 12 : 399–450. CiteSeerX 10.1.1.7.7106 . doi : 10.1017 / S0962492902000144 .
- ^ Niiranen, Jouko: algoritmo de Euler simétrico rápido y preciso para simulaciones electromecánicas Proceedings of the Electrimacs'99, 14-16 de septiembre de 1999 Lisboa, Portugal, vol. 1, páginas 71 - 78.
- Giordano, Nicholas J .; Hisao Nakanishi (julio de 2005). Física Computacional (2ª ed.). Benjamin Cummings. ISBN 0-13-146990-8.
- MacDonald, James. "El método Euler-Cromer" . Universidad de Delaware . Consultado el 11 de abril de 2013 .
- Vesely, Franz J. (2001). Física Computacional: Introducción (2ª ed.). Saltador. págs. 117 . ISBN 978-0-306-46631-1.