Dibuat oleh: roberto.c.alfredo di physics pada
1. El problema de siempre: ¿cómo avanzar en el tiempo?
Imagina que tienes un planeta girando en torno a una estrella. Tú conoces sus posiciones y velocidades en el presente, y la fuerza que actúa (gravedad). La pregunta eterna: ¿dónde estará un pasito después?
Aquí entra la integración numérica. En la compu, el tiempo no es continuo: vivimos de ticks discretos \(\Delta t\). Pero cuidado: si eliges mal el método, tu planeta puede terminar fugándose a otra galaxia solo porque el algoritmo le metió más energía de la cuenta. 😅
2. El primo bruto: Euler
El método de Euler dice:
\[ \vec x_{t+\Delta t} \approx \vec x_t + \vec v_t \,\Delta t, \qquad \vec v_{t+\Delta t} \approx \vec v_t + \vec a_t\,\Delta t. \]
Es directo y fácil, pero no conserva la energía. Resultado: órbitas que deberían ser elipses terminan como espirales abiertas o colapsos al centro. Un desastre.
Analogía musical: es como tocar batería sin metrónomo. Empiezas bien, pero al tercer compás ya vas a destiempo. 🥁
3. La idea de Verlet (1967)
Loup Verlet, trabajando en dinámica molecular, pensó: “¿Y si uso tanto la aceleración presente como la futura para equilibrar la velocidad?” Boom:
\[ \vec x_{t+\Delta t} = \vec x_t + \vec v_t \,\Delta t + \tfrac{1}{2}\vec a_t\,\Delta t^2, \]
\[ \vec v_{t+\Delta t} = \vec v_t + \tfrac{1}{2}\left(\vec a_t + \vec a_{t+\Delta t}\right)\Delta t. \]
Esto es el Velocity-Verlet, el que usamos en el bloque N-cuerpos. Lo mágico: es un integrador simpléctico. Palabra rimbombante que significa que respeta la geometría del espacio de fases: la energía no se “derrite” tan fácil.