Simuladores / N cuerpos: integradores y caos

N cuerpos: integradores y caos

SecundariaUniversitaria

RK4 frente a leapfrog simpléctico en dos cuerpos excéntricos, la figura en ocho y el problema pitagórico; error de energía en escala log y divergencia de trayectorias vecinas (Lyapunov).

\ddot{\mathbf r}_i = \sum_{j\ne i} \frac{G m_j(\mathbf r_j-\mathbf r_i)}{|\mathbf r_j-\mathbf r_i|^3},\quad \tilde H = H + O(h^2),\quad \lambda = \lim_{t\to\infty}\tfrac1t\ln\tfrac{\delta(t)}{\delta_0}

EscenarioIntegradorPaso h = 1.0e-2Velocidad = 1.5 unidades de t / s Perturbar x₁ en δ₀ = 10⁻⁶ (copia vecina)

tres masas iguales, T ≈ 6.3259 · t = 0.00

Librería: src/lib/math/nbody.tssrc/lib/math/integrators.ts