Quasi-frozen polar orbits around the Moon

Resumen Técnico: Órbitas Polares Cuasi-congeladas alrededor de la Luna

Planteamiento del Problema
El estudio aborda el desafío de mantener trayectorias estables a largo plazo para satélites artificiales en órbitas polares de baja altitud lunar. Debido a la distribución de masa altamente no uniforme de la Luna (mascones), las órbitas no gestionadas están sujetas a severas perturbaciones que típicamente causan el colapso de las trayectorias hacia la superficie lunar en periodos de tiempo cortos. El objetivo principal es identificar configuraciones orbitales específicas que prevengan eventos de impacto prematuros mediante el equilibrio de las derivas seculares en los elementos orbitales.

Metodología
Los autores emplean un enfoque de marco dual que combina la derivación analítica y la validación numérica:

1. Modelado Dinámico:

   • Perturbaciones: El modelo tiene en cuenta el campo gravitatorio no uniforme de la Luna y las perturbaciones gravitatorias de tercer cuerpo de la Tierra. La Presión de Radiación Solar (SRP) se considera insignificante para satélites típicos y se excluye.
   • Modelo 18 × 0 (Analítico): Un modelo analítico simplificado que incorpora armónicos zonales hasta el orden 18 ($J_2$ a $J_{18}$). Este modelo utiliza un método de promedio único para eliminar los términos de periodo corto.
   • Modelo 18 × 3 (Analítico Mejorado/Numérico): Un modelo expandido que incluye los armónicos zonales hasta $J_{18}$ más términos sectoriales y teserales específicos de orden y grado tres ($C_{22}, S_{22}, C_{31}, S_{31}, C_{32}, S_{32}, C_{33}, S_{33}$).
   • Tratamiento de Tercer Cuerpo: La influencia de la Tierra se modela utilizando una formulación de promedio único para una órbita elíptica e inclinada.

2. Análisis de Órbita Congelada:

   • El estudio define las "órbitas congeladas" como configuraciones donde las variaciones seculares en la excentricidad ($\dot{e}$) y la inclinación ($\dot{i}$) son cero.
   • El análisis se centra en fijar el argumento del perilo ($g$) en $90^\circ$ y $270^\circ$.
   • Debido a que el modelo analítico con armónicos zonales de alto grado hace que la derivación de una expresión explícita $e = f(i, a)$ sea matemáticamente intratable, los autores utilizan un barrido numérico de búsqueda de raíces semi-analítico para resolver $\dot{g} = 0$ y mapear los contornos de equilibrio en el plano de fase $(e, i)$.

3. Validación:

   • Integración Numérica: Las ecuaciones de movimiento derivadas del Modelo 18 × 3 se integran utilizando el método Runge-Kutta-Fehlberg (rkf45) en Maple.
   • Benchmarking de Alta Fidelidad: Los resultados se validan de forma cruzada utilizando la herramienta General Mission Analysis Tool (GMAT) de la NASA. La configuración de GMAT utiliza un modelo de potencial selenográfico de orden y grado 50 × 50 de alta fidelidad y efemérides planetarias DE405 para la perturbación de la Tierra.

Resultados Clave

• Condiciones de Órbita Congelada: El análisis identifica que existen configuraciones polares estables donde la gravedad zonal lunar equilibra las derivas seculares. El estudio analiza tanto las configuraciones de $g = 90^\circ$ como las de $g = 270^\circ$, señalando que, si bien ambas producen soluciones congeladas en el modelo analítico, la configuración de $g = 270^\circ$ se destaca específicamente por generar órbitas polares estables y congeladas dentro del entorno lunar al considerar el contexto dinámico completo.
• Dependencia del Semieje Mayor:
  • En un semieje mayor de 1838 km (aprox. 100 km de altitud), la excentricidad congelada requerida es alta ($e \approx 0.049$). La propagación numérica bajo el Modelo 18 × 3 revela que, incluso con estas condiciones iniciales, la alta excentricidad hace que la órbita sea susceptible a un impacto rápido (dentro de ~250 días) debido a perturbaciones no zonales.
  • En un semieje mayor de 1938 km (aprox. 200 km de altitud), la excentricidad requerida se relaja a $e \approx 0.0241$. Esta configuración demuestra una estabilidad significativamente mayor.
• Comportamiento Cuasi-Congelado: Al propagar las condiciones iniciales derivadas del Modelo analítico 18 × 0 utilizando el modelo más complejo 18 × 3, las órbitas no permanecen perfectamente congeladas, sino que exhiben oscilaciones acotadas de baja amplitud en la excentricidad y el argumento del perilo. Estas se clasifican como órbitas cuasi-congeladas.
• Verificación GMAT: Las simulaciones de GMAT (modelo 50 × 50 + Tierra) confirman que la configuración cuasi-congelada identificada en $a = 1938$ km mantiene una trayectoria acotada durante periodos prolongados (1000 días probados). El marco semi-analítico simplificado (Modelo 18 × 3 + Tierra) coincide estrechamente con los resultados de alta fidelidad de GMAT, mostrando solo desviaciones menores en la excentricidad y el argumento del perilo.

Significancia y Reivindicaciones
El artículo afirma proporcionar una metodología semi-analítica robusta para identificar órbitas polares congeladas de baja altitud alrededor de la Luna. Al utilizar formulaciones zonales de alto orden (hasta $J_{18}$) y validarlas contra modelos numéricos de alta fidelidad, los autores establecen una base fiable para el diseño de misiones.

La contribución principal es la identificación de condiciones iniciales específicas ($a = 1938$ km, $e = 0.0241$, $i = 90^\circ$, $g = 270^\circ$, $h = 90^\circ$) que minimizan la necesidad de maniobras de mantenimiento de estación (station-keeping). Los autores sostienen que estos resultados ofrecen una base práctica para la selección de órbitas operativas para constelaciones lunares permanentes y para la predicción de trayectorias de "cementerio", reduciendo así el consumo de combustible en comparación con estados orbitales vecinos. El estudio enfatiza que, si bien el modelo analítico es una simplificación, captura la dinámica primaria con alta fidelidad, lo que lo hace adecuado para el diseño preliminar de misiones.