Spectral Koopman Approach for Reconstructing State-space Geometry of Cislunar Restricted 3-Body Problem

Sintesi Tecnica: Approccio Spettrale di Koopman per la Ricostruzione della Geometria dello Spazio degli Stati del Problema dei Tre Corpi Ristretto Cislunare

Definizione del Problema
L'ambiente cislunare, governato dall'interazione gravitazionale tra Terra, Luna e Sole, presenta un sistema dinamico altamente non lineare caratterizzato da punti di librazione, varietà invarianti (manifolds) e trasporto caotico. La progettazione di missioni tradizionali e l'analisi delle traiettorie in questa regione si affidano al Problema dei Tre Corpi Ristretto Circolare (CR3BP) e a tecniche di ottimizzazione numerica (ad esempio, correzione differenziale, multiple shooting). Tuttavia, questi metodi esistenti affrontano tre persistenti limitazioni:

1. Discrepanze di Fedeltà: Discrepanze tra modelli semplificati e la realtà ad alta fedeltà.
2. Costo Computazionale: La necessità di ricalcoli costosi quando avvengono aggiornamenti di navigazione.
3. Problemi di Scalabilità: Scarsa scalabilità dei metodi Monte Carlo per la propagazione dell'incertezza su lunghi orizzonti temporali.

Inoltre, le tecniche di linearizzazione locale standard falliscono nel catturare la geometria globale dello spazio delle fasi organizzata da strutture non lineari forti come risonanze e canali di trasporto. Tuttavia, gli approcci basati sulla teoria degli operatori (specificamente l'operatore di Koopman) offrono una rappresentazione lineare globale dei sistemi non lineari, ma le attuali approssimazioni basate sui dati (ad esempio, l'Estesa Decomposizione dei Modi Dinamici - EDMD) richiedono spesso una selezione accurata del dizionario di funzioni base e faticano a recuperare accuratamente lo spettro globale intrinseco.

Metodologia
Gli autori propongono un nuovo approccio basato sulla formulazione dell'integrale di cammino dello spettro di Koopman per ricostruire la geometria globale dello spazio delle fasi del CR3BP planare senza fare affidamento su una selezione ad hoc di dizionari.

1. Framework dell'Operatore di Koopman: Il sistema dinamico $\dot{x} = f(x)$ viene analizzato attraverso l'evoluzione degli osservabili. L'operatore di Koopman $U_t$ fornisce una descrizione lineare del flusso non lineare. Le autofunzioni principali $\phi_\lambda$ di questo operatore soddisfano l'equazione $\frac{\partial \phi_\lambda}{\partial x}f = \lambda \phi_\lambda$.
2. Formulazione dell'Integrale di Cammino: Invece di approssimare l'operatore in un sottospazio a dimensione finita, gli autori calcolano direttamente le autofunzioni principali. Decomponendo il campo vettoriale in componenti lineari ($Ax$) e non lineari ($F_n(x)$), la parte non lineare dell'autofunzione $h_\lambda(x)$ è espressa come un integrale di cammino:
   $$h_\lambda(x) = \int_0^\infty e^{-\lambda \tau} w_\lambda^\top F_n(s_\tau(x)) d\tau$$
   dove $s_\tau(x)$ rappresenta la traiettoria del sistema.
3. Implementazione Numerica: Le autofunzioni sono approssimate valutando l'integrale su un orizzonte temporale finito $T$. Le autofunzioni instabili (associate a $\lambda_u > 0$) sono calcolate utilizzando traiettorie in avanti, mentre le autofunzioni stabili (associate a $\lambda_s < 0$) sono calcolate utilizzando traiettorie all'indietro.
4. Ricostruzione della Geometria: Gli autori utilizzano la proprietà per cui le curve di livello zero congiunte di queste autofunzioni principali corrispondono alle varietà invarianti (stabili, instabili e centrali) dei punti di equilibrio (punti di Lagrange) del sistema.

Contributi Chiave

• Rappresentazione Lineare Globale: Il documento dimostra che il framework di Koopman può generare una rappresentazione lineare globale del CR3BP, superando i limiti della linearizzazione locale.
• Calcolo Diretto delle Autofunzioni: Sfruttando l'approccio dell'integrale di cammino, il metodo calcola direttamente le autofunzioni principali dai dati di serie temporali, evitando la necessità di una progettazione manuale del dizionario richiesta da metodi come l'EDMD.
• Identificazione delle Varietà Invarianti: Il lavoro stabilisce che gli insiemi di livello zero delle autofunzioni calcolate codificano le caratteristiche dello spazio delle fasi, identificando specificamente le varietà stabili, instabili e centrali associate ai cinque punti di Lagrange ($L_1$ attraverso $L_5$).
• Visualizzazione dei Canali di Trasporto: L'approccio visualizza con successo le "autostrade invisibili" dell'ambiente cislunare, mostrando come le varietà instabili di un equilibrio possano allinearsi con le varietà stabili di un altro per formare canali di trasporto.

Risultati
Il metodo proposto è stato applicato al modello CR3BP 4D planare (sistema Terra-Luna). I risultati della simulazione includono:

• Ricostruzione delle Curve di Velocità Zero (ZVCs): Il metodo recupera le regioni ammesse energeticamente e le zone proibite definite dalla costante di Jacobi.
• Geometria Locale agli Equilibri:
  • $L_1$ e $L_2$: Le autofunzioni delineano chiaramente i gateway di trasporto. Le varietà instabili si estendono verso l'esterno attraverso i "colli" (permettendo l'uscita dalla Terra o l'ingresso verso l'esterno), mentre le varietà stabili si curvano verso l'interno (catturando le traiettorie).
  • $L_3$: Esibisce varietà fragili che perdono verso le regioni esterne piuttosto che formare forti corridoi di trasporto.
  • $L_4$ e $L_5$: Le autofunzioni centrali dominano, producendo loop chiusi corrispondenti a orbite "tadpole" e "Trojan" stabili, coerentemente con la stabilità lineare di questi punti triangolari.
• Partizionamento Globale dello Spazio delle Fasi: Le curve di livello zero delle autofunzioni partizionano lo spazio delle fasi spettralmente. A differenza delle superfici di Jacobi basate sull'energia, che definiscono dove il moto è possibile, le curve di Koopman rivelano come le traiettorie si organizzano asintoticamente all'interno di quelle regioni, identificando separatrici e confini invarianti.
• Meccanismo di Trasporto: I risultati illustrano un canale di trasporto tra $L_1$ e $L_2$, dove la varietà instabile di $L_1$ e la varietà stabile di $L_2$ delimitano un tubo invariante, facilitando il trasferimento tra le regioni Terra e Luna.

Significato e Rivendicazioni
Il documento afferma che l'approccio spettrale di Koopman offre uno strumento potente per la dinamica cislunare, spostando l'attenzione dalle singole traiettorie all'evoluzione degli osservabili. Il significato primario risiede nella capacità di:

1. Recuperare la Geometria Globale: Ricostruire la geometria globale dello spazio delle fasi del CR3BP dai dati di serie temporali, catturando strutture non lineari che i metodi locali perdono.
2. Abilitare Compiti Successivi: Fornire una base per la propagazione dell'incertezza, l'analisi di raggiungibilità e la progettazione del controllo. Poiché l'operatore di Koopman è lineare, gli autori postulano che l'analisi dell'incertezza e la progettazione di strategie di controllo per il CR3BP non lineare possano essere semplificate e rese più scalabili rispetto ai metodi non lineari tradizionali.
3. Scoperta Basata sui Dati: Aprire una nuova direzione per la scoperta basata sui dati nella dinamica spaziale, dove la linearità dell'operatore permette l'analisi e la progettazione di strategie di controllo per sistemi intrinsecamente non lineari.

Gli autori concludono che, sebbene l'attuale lavoro si concentri sul recupero della geometria del CR3BP, il framework promette applicazioni future nella ottimizzazione delle traiettorie in tempo reale e nel controllo robusto in ambienti complessi multi-corpo.