Quasi-frozen polar orbits around the Moon

Sintesi Tecnica: Orbite Polari Quasi-Congelate attorno alla Luna

Problematica
Lo studio affronta la sfida di mantenere traiettorie stabili a lungo termine per satelliti artificiali in orbite lunari a bassa quota. A causa della distribuzione di massa altamente non uniforme della Luna (mascon), le orbite non gestite sono soggette a severe perturbazioni che tipicamente causano il collasso delle traiettorie sulla superficie lunare in tempi brevi. L'obiettivo primario è identificare specifiche configurazioni orbitali che evitino eventi di impatto prematuri bilanciando i drift secolari degli elementi orbitali.

Metodologia
Gli autori impiegano un approccio a doppio framework che combina la derivazione analitica e la validazione numerica:

1. Modellazione Dinamica:

   • Perturbazioni: Il modello tiene conto del campo gravitazionale non uniforme della Luna e delle perturbazioni gravitazionali di terzo corpo da parte della Terra. La pressione di radiazione solare (SRP) è considerata trascurabile per i tipici satelliti ed è esclusa.
   • Modello 18 × 0 (Analitico): Un modello analitico semplificato che incorpora gli armonici zonali fino all'ordine 18 ($J_2$ a $J_{18}$). Questo modello utilizza un metodo a media singola per eliminare i termini a breve periodo.
   • Modello 18 × 3 (Analitico/Numerico Potenziato): Un modello esteso che include gli armonici zonali fino a $J_{18}$ più specifici termini settoriali e tesserali di ordine e grado tre ($C_{22}, S_{22}, C_{31}, S_{31}, C_{32}, S_{32}, C_{33}, S_{33}$).
   • Trattamento del Terzo Corpo: L'influenza della Terra è modellata utilizzando una formulazione a media singola per un'orbita ellittica e inclinata.

2. Analisi delle Orbite Congelate (Frozen Orbits):

   • Lo studio definisce "orbite congelate" come configurazioni in cui le variazioni secolari dell'eccentricità ($\dot{e}$) e dell'inclinazione ($\dot{i}$) sono nulle.
   • L'analisi si concentra sul fissare l'argomento del perilune ($g$) a $90^\circ$ e $270^\circ$.
   • Poiché il modello analitico con armonici zonali di alto grado rende matematicamente intrattabile la derivazione di un'espressione esplicita $e = f(i, a)$, gli autori utilizzano un metodo numerico di ricerca delle radici semi-analitico per risolvere $\dot{g} = 0$ e mappare i contorni di equilibrio nel piano di fase $(e, i)$.

3. Validazione:

   • Integrazione Numerica: Le equazioni del moto derivate dal Modello 18 × 3 sono integrate utilizzando il metodo Runge-Kutta-Fehlberg (rkf45) in Maple.
   • Benchmarking ad Alta Fedeltà: I risultati sono validati incrociati utilizzando il General Mission Analysis Tool (GMAT) della NASA. La configurazione GMAT utilizza un modello di potenziale selenico ad alto ordine e grado 50 × 50 e le effemeridi planetarie DE405 per la perturbazione terrestre.

Risultati Chiave

• Condizioni di Orbita Congelata: L'analisi identifica che esistono configurazioni polari stabili dove i coefficienti della gravità zonale lunare bilanciano i drift secolari. Lo studio analizza sia le configurazioni $g = 90^\circ$ che $g = 270^\circ$, osservando che, sebbene entrambe forniscano soluzioni congelate nel modello analitico, la configurazione $g = 270^\circ$ è specificamente evidenziata per generare orbite polari congelate e stabili nell'ambiente lunare quando si considera il contesto dinamico completo.
• Dipendenza dal Semiasse Maggiore:
  • A un semiasse maggiore di 1838 km (circa 100 km di altitudine), l'eccentricità congelata richiesta è elevata ($e \approx 0,049$). La propagazione numerica sotto il Modello 18 × 3 rivela che, anche con queste condizioni iniziali, l'alta eccentricità rende l'orbita suscettibile a un impatto rapido (entro circa 250 giorni) a causa delle perturbazioni non zonali.
  • A un semiasse maggiore di 1938 km (circa 200 km di altitudine), l'eccentricità richiesta si riduce a $e \approx 0,0241$. Questa configurazione dimostra una stabilità significativamente maggiore.
• Comportamento Quasi-Congelato: Propagando le condizioni iniziali derivate dal Modello analitico 18 × 0 utilizzando il più complesso Modello 18 × 3, le orbite non rimangono perfettamente congelate ma esibiscono oscillazioni limitate e a bassa ampiezza in eccentricità e argomento del perilune. Queste sono classificate come orbite quasi-congelate.
• Verifica GMAT: Le simulazioni GMAT (modello 50 × 50 + Terra) confermano che la configurazione quasi-congelata identificata a $a = 1938$ km mantiene una traiettoria limitata per lunghe durate (1000 giorni testati). Il framework semi-analitico semplificato (Modello 18 × 3 + Terra) corrisponde strettamente ai risultati GMAT ad alta fedeltà, mostrando solo deviazioni minori nell'eccentricità e nell'argomento del perilune.

Significatività e Rivendicazioni
Il documento sostiene di fornire una metodologia semi-analitica robusta per l'identificazione di orbite polari congelate a bassa quota attorno alla Luna. Utilizzando formulazioni zonali di alto ordine (fino a $J_{18}$) e validandole rispetto a modelli numerici ad alta fedeltà, gli autori stabiliscono una linea di base affidabile per la progettazione di missioni.

Il contributo primario è l'identificazione di specifiche condizioni iniziali ($a = 1938$ km, $e = 0,0241$, $i = 90^\circ$, $g = 270^\circ$, $h = 90^\circ$) che minimizzano la necessità di manovre di mantenimento della stazione (station-keeping). Gli autori affermano che questi risultati offrono una base pratica per la selezione di orbite operative per costellazioni lunari permanenti e per la previsione di traiettorie "graveyard", riducendo così il consumo di carburante rispetto agli stati orbitali vicini. Lo studio sottolinea che, sebbene il modello analitico sia una semplificazione, esso cattura la dinamica primaria con alta fedeltà, rendendolo adatto per la progettazione preliminare di missioni.