¿Por qué el algoritmo Velocity Verlet es superior a Euler en astrodinámica?

Porque el integrador Velocity Verlet es simpléctico y reversible en el tiempo. A diferencia del método de Euler hacia delante (que inyecta energía artificialmente en cada ciclo expandiendo las órbitas en espiral hacia el infinito) o de métodos disipativos, Verlet conserva el volumen en el espacio de fases y acota el error del hamiltoniano dentro de oscilaciones periódicas infinitesimales sin deriva energética acumulada.

¿Qué es el problema de los N cuerpos y por qué genera caos determinista?

Para dos cuerpos (N=2), Newton y Kepler demostraron que las órbitas tienen soluciones analíticas cerradas en forma de cónicas (elipses, parábolas o hipérbolas). Sin embargo, a partir de tres cuerpos (N≥3), Henri Poincaré descubrió en 1889 que no existe solución analítica general en términos de integrales algebraicas estándar: las interacciones no lineales mutuas generan caos determinista, donde una alteración infinitesimal en las condiciones de partida produce trayectorias totalmente divergentes e impredecibles a medio plazo.

¿Qué son los puntos de Lagrange y por qué son estables L4 y L5?

Los puntos de Lagrange son las cinco posiciones de equilibrio gravitatorio donde las fuerzas de atracción de dos cuerpos masivos en órbita circular mutua se equilibran exactamente con la fuerza centrífuga en un sistema de referencia rotatorio. Mientras que L1, L2 y L3 son puntos colineales inestables tipo punto de silla, los puntos L4 y L5 forman triángulos equiláteros y son estables gracias a la fuerza de Coriolis, atrapando asteroides troyanos en órbitas cerradas de libración.

¿Cómo se resuelven las colisiones inelásticas en el simulador?

Cuando la distancia entre dos cuerpos es menor o igual a la suma de sus radios físicos, se produce una fusión inelástica perfecta. Se conserva estrictamente la masa total (M = m1 + m2) y el momento lineal total (V = (m1·v1 + m2·v2) / M), recalculando la posición del nuevo cuerpo en el baricentro del choque y modelando el radio volumétrico equivalente.

¿Para qué sirve el parámetro de suavizado gravitacional (epsilon)?

El potencial de Plummer añade un término ε² en el denominador de la ley de Newton: F ~ G·m1·m2 / (r² + ε²)^(3/2). En la realidad, las partículas puntuales no existen: los planetas y estrellas son esferas continuas de materia. El parámetro epsilon evita que la fuerza tienda a infinito cuando dos masas pasan a distancias cercanas a cero en una simulación discreta de paso temporal finito.

Sandbox orbital de gravedad y caos de N cuerpos • Fijate Tools

Laboratorio de mecánica celeste 2D

Arrastra para disparar cuerpos con vector de velocidad o selecciona presets orbitales.

Cuerpos: 0

Fusiones: 0

Masa total: 0

Energía: 0

Deriva H: 0.00%

60 FPS

Escenarios predefinidos

Lanzador de cuerpos

Órbita circular automática asistida

Física gravitatoria

Velocidad temporal: 1.0×

Constante gravitacional (G): 1.00

Suavizado de Plummer (ε): 10 px

Dinámica de colisión:

Cámara y visuales

Seguimiento de cámara:

Longitud de estelas: 140 pts

Estelas luminosas Vectores velocidad Nombres celestes Marcador baricentro Cuadrícula espacial

Fundamentos de la mecánica celeste y el problema de N cuerpos

La ley de gravitación universal formulada por Isaac Newton en 1687 postula que dos cuerpos cualesquiera en el universo se atraen con una fuerza directamente proporcional al producto de sus masas e inversamente proporcional al cuadrado de la distancia que los separa. Cuando el número de cuerpos es dos ($N=2$), el sistema admite solución matemática exacta cerrada en forma de órbitas cónicas (elipses de Kepler, parábolas o hipérbolas).

Sin embargo, al incorporar un tercer cuerpo ($N=3$) o una multitud ($N > 3$), el acoplamiento cruzado de campos gravitatorios no lineales elimina la existencia de integrales primeras analíticas generales. Como demostró Henri Poincaré a finales del siglo XIX, los sistemas gravitatorios de tres o más cuerpos exhiben caos determinista: su evolución es completamente estricta y matemática, pero una perturbación tan ínfima como la posición de un átomo provoca trayectorias radicalmente distintas con el paso del tiempo.

1. Integración Velocity Verlet

El simulador utiliza un algoritmo numérico simpléctico de segundo orden. A diferencia de los integradores de Euler tradicionales, Verlet conserva el volumen del espacio de fases y acota la energía hamiltoniana total ($E = E_k + E_p$) con una deriva cercana al 0,001 %, evitando que los planetas salgan disparados artificialmente al infinito.

2. Suavizado de Plummer (ε)

Para evitar divisiones por cero y aceleraciones infinitas irreales cuando dos masas puntuales se cruzan a distancias infinitesimales, se emplea el potencial de Plummer, sumando un término suavizador ε² en el denominador: F ∝ G · m₁ · m₂ / (r² + ε²)^(3/2).

3. Colisiones inelásticas

En encuentros frontales donde d ≤ R₁ + R₂, se ejecuta una fusión inelástica perfecta. La cantidad de movimiento lineal total P = m₁ · v₁ + m₂ · v₂ y la masa combinada se conservan al 100 %, calculando el nuevo radio volumétrico y desencadenando ondas de choque de plasma.

4. Puntos de Lagrange (L1 a L5)

En sistemas de dos cuerpos masivos en rotación mutua, existen cinco posiciones singulares de equilibrio relativo. Los puntos colineales L1, L2 y L3 son inestables, mientras que L4 y L5 (situados a 60° a la cabeza y cola orbital) son cuencas estables protegidas por la fuerza de Coriolis que albergan asteroides troyanos.

Preguntas frecuentes sobre la simulación gravitatoria

¿Cómo se garantiza que la simulación no inventa ni pierde energía orbital?

El simulador emplea integración Velocity Verlet dividida en 6 subpasos temporales por cada fotograma a 60 FPS. Al ser un método simpléctico simétrico en el tiempo, el error de truncamiento numérico no se acumula de forma monótona, garantizando que el sistema solar y las órbitas en figura de 8 puedan orbitar indefinidamente sin desintegrarse.

¿Qué representa la solución de la figura en ocho de Chenciner y Montgomery?

Durante más de tres siglos se pensó que las únicas órbitas periódicas de 3 cuerpos eran las de Lagrange y Euler. En el año 2000, los matemáticos Alain Chenciner y Richard Montgomery demostraron la existencia de una coreografía cerrada donde tres masas idénticas se persiguen indefinidamente trazando una curva en forma de 8 sin colisionar jamás. Puedes experimentarla en vivo seleccionando dicho preset.

¿Por qué las galaxias forman colas de marea al colisionar?

En una colisión galáctica, las estrellas raramente chocan físicamente entre sí debido a las inmensas distancias interestelares. Sin embargo, el gradiente del campo gravitatorio del núcleo invasor arranca las estrellas exteriores de sus órbitas circulares, estirándolas en filamentos curvados de cientos de miles de años luz conocidos como colas de marea.

¿Cómo disparar un cuerpo para lograr una órbita circular perfecta?

Puedes activar la casilla «Órbita circular automática asistida». Al arrastrar y soltar el cuerpo, el simulador calculará en vivo la velocidad tangencial exacta $v = \sqrt{G \cdot M / r}$ perpendicular al radio vector respecto al cuerpo masivo más próximo, permitiéndote poner satélites y planetas en órbita estable con un solo clic.