Integración de Verlet
Es particularmente apropiado en las situaciones en que la expresión de la segunda derivada solo es función de las variables, dependiente o independiente, sin participar la primera derivada.
Este es el caso de numerosos problemas de la dinámica newtoniana, por lo que se emplea frecuentemente en astronomía y mecánica molecular.
Dado un sistema de coordenadas, la posición de una partícula a lo largo del tiempo puede ser expresada por una función de las coordenadas del vector y del tiempo.
En 1967, el matemático francés Loup Verlet presentó la primera versión de su modelo,[1] denominada Integración de Verlet, caracterizada por su simplicidad sin pérdida de exactitud y estabilidad.
Posteriormente, en 1985,[2] se propuso una ligera corrección a la Integración de Verlet, conocida como algoritmo de Verlet con velocidad, que mejora la precisión y estabilidad de las soluciones.
Desarrollando en serie de Taylor el valor de la función para dos incrementos de la variable independiente: dt y - dt se obtiene: Ec.
Sumando ambas ecuaciones, limitando la serie hasta el término de cuarto grado, y reagrupando términos se obtiene la ley de recursión para las posiciones: Ec.
Al despreciar el último término el error cometido es de cuarto grado, inferior a la mayoría de los algoritmos convencionales.
(4) proporciona la trayectoria sin intervenir la velocidad, lo que simplifica los cálculos.
Caso de necesitar conocer las velocidades, estas se obtienen directamente por simple diferencia de las Ec.
Al comenzar la serie, se halla el tercer término (
El primer término, es un dato del problema, pero el segundo es desconocido y la Ec.
Para obviar esta limitación, el segundo valor de la serie (
) se calcula a partir de los valores iniciales con ayuda de algún método explícito que permita esta operación, por ejemplo el de Euler, Ec.
Y ya se tienen todos los términos necesarios para proseguir el cálculo recursivo.
es suficientemente pequeño, suele corregirse automáticamente.
Como es conocido el centro de masas del conjunto puede ser tomado como origen de coordenadas y las órbitas se mueven en un plano.
Adoptando un sistema de coordenadas y posiciones iniciales como se muestra en el gráfico, el problema se reduce a un sistemas de cuatro variables dependientes, las coordenadas de posición de los dos cuerpos, una variable independiente, el tiempo, todo ello ligado por cuatro ecuaciones que se deben resolver secuencialmente en cada paso de la iteración.
La figura muestra los cálculos realizados en una Hoja de Cálculo y la salida de un programa gráfico escrito en javascript El dibujo se ha obtenido trazando varios periodos con el propósito de mostrar la falta de coincidencia entre órbitas consecutivas, motivado por haber seleccionado un
Resulta sencillo mejorar la precisión reduciendo el intervalo o, como se muestra a continuación, aplicando el siguiente algoritmo.
Como mejora del original algoritmo, se ha propuesto una fórmula aplicable desde el principio donde interviene explicitamente la velocidad, que es calculada por un procedimiento algo más refinado.
En la expresión del desarrollo en serie, Ec.
Limitando la serie al término de segundo grado,[3] el error en la posición es de orden tercero, se encuentra la forma iterativa para determinar la posición en función de velocidad (d
Para calcular las velocidades, operando sobre las expresiones obtenidas en las Ec.
(6)[4] se llega a la siguiente expresión: Ec.
(8) presentan gran analogía con las utilizadas en el procedimiento de Integración del salto de rana; en Verlet, posición y velocidad se obtienen en el mismo punto de la trayectoria, mientras que en salto de rana se calculan intercalados.
se necesita el valor previo de
Pero a diferencia de la Integración de Verlet, no es necesario recurrir a ninguna suposición externa al propio procedimiento Llamando t0, v0 y a0 a los valores de la respectivas variables en el instante inicial y t1 = t0 +
, v1 y a1 para el siguiente término: A partir de la Ec.
el perfil de las órbitas parece nítido por completo, lo que significa que en las sucesivas revoluciones se ha utilizado idéntico camino(hasta donde alcanza la precisión observable), siendo preferible esta vía para mejorar los resultados.
