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

Technische Samenvatting: Spectrale Koopman-benadering voor de Reconstructie van de Toestandsruimte-geometrie van het Cislunaire Beperkte Drie-Lichamenprobleem

Probleemstelling
De cislunaire omgeving, die wordt beheerst door de gravitationele wisselwerking tussen de Aarde, de Maan en de Zon, vormt een hoogst nietlineair dynamisch systeem dat gekenmerkt wordt door libratiepunten, invariante manifolden en chaotische transport. Traditioneel missieontwerp en trajectanalyse voor deze regio vertrouwen op het Circulaire Beperkte Drie-Lichamenprobleem (CR3BP) en numerieke optimalisatietechnieken (bijv. differentiële correctie, multiple shooting). Echter, deze bestaande methoden kampen met drie hardnekkige beperkingen:

1. Fidelity Mismatches: Discrepanties tussen vereenvoudigde modellen en de hoogwaardige realiteit.
2. Computationele Kosten: De noodzaak voor kostbare herberekening wanneer navigatie-updates plaatsvinden.
3. Schaalbaarheidsproblemen: Slechte schaalbaarheid van Monte Carlo-methoden voor langetermijn onzekerheidsvoortplanting.

Bovendien falen standaard lokale linearisatietechnieken er vaak in om de globale fase-ruimtegeometrie te vangen die georganiseerd wordt door sterke nietlineaire structuren zoals resonanties en transportkanalen. Hoewel operator-theoretische benaderingen (specifiek de Koopman-operator) een globale lineaire representatie bieden van nietlineaire systemen, hebben bestaande data-gedreven benaderingen (zoals Extended Dynamic Mode Decomposition) vaak zorgvuldige selectie van de basis nodig en hebben ze moeite om het intrinsieke globale spectrum nauwkeurig te herstellen.

Methodologie
De auteurs stellen een nieuwe benadering voor op basis van de padintegraalformulering van het Koopman-spectrum om de globale fase-ruimtegeometrie van het planaire CR3BP te reconstrueren zonder afhankelijk te zijn van ad hoc dictionary-selectie.

1. Koopman Operator Framework: Het systeemdynamica $\dot{x} = f(x)$ wordt geanalyseerd via de evolutie van observables. De Koopman-operator $U_t$ biedt een lineaire beschrijving van de nietlineaire flow. De primaire eigenfuncties $\phi_\lambda$ van deze operator voldoen aan de vergelijking $\frac{\partial \phi_\lambda}{\partial x}f = \lambda \phi_\lambda$.
2. Padintegraalformulering: In plaats van de operator te benaderen in een einddimensionale subspace, berekenen de auteurs de primaire eigenfuncties direct. Door het vectorveld te ontleden in lineaire ($Ax$) en nietlineaire ($F_n(x)$) componenten, wordt het nietlineaire deel van de eigenfunctie $h_\lambda(x)$ uitgedrukt als een padintegraal:
   $$h_\lambda(x) = \int_0^\infty e^{-\lambda \tau} w_\lambda^\top F_n(s_\tau(x)) d\tau$$
   waarbij $s_\tau(x)$ de systeemtraject representeert.
3. Numerieke Implementatie: De eigenfuncties worden benaderd door de integraal over een eindig tijdvenster $T$ te evalueren. Onstabiele eigenfuncties (geassocieerd met $\lambda_u > 0$) worden berekend met behulp van voorwaartse trajecten, terwijl stabiele eigenfuncties (geassocieerd met $\lambda_s < 0$) worden berekend met behulp van achterwaartse trajecten.
4. Geometrie Reconstructie: De auteurs maken gebruik van het eigenschap dat de gezamenlijke nul-niveau curves van deze primaire eigenfuncties overeenkomen met de invariante manifolden (stabiel, onstabiel en centrum) van de evenwichtspunten (Lagrangepunten) van het systeem.

Belangrijkste Bijdragen

• Globale Lineaire Representatie: Het artikel toont aan dat het Koopman-framework een globale lineaire representatie van de CR3BP kan genereren, wat de beperkingen van lokale linearisatie overwint.
• Directe Berekening van Eigenfuncties: Door gebruik te maken van de padintegraal-benadering berekent deze methode de primaire eigenfuncties direct uit tijdreeksdata, waardoor de noodzaak voor handmatige dictionary-ontwerp die vereist is bij methoden zoals EDMD wordt omzeild.
• Identificatie van Invariante Manifolds: Het werk stelt vast dat de nul-niveau verzamelingen van de berekende eigenfuncties de fase-ruimte kenmerken coderen, specif gezien de identificatie van stabiele, onstabiele en centrum manifolden geassocieerd met de vijf Lagrangepunten ($L_1$ tot en met $L_5$).
• Visualisatie van Transportkanalen: De benadering visualiseert succesvol de "onzichtbare snelwegen" van de cislunaire omgeving, waarbij getoond wordt hoe de onstabiele manifolden van het ene evenwichtspunt kunnen uitlijnen met de stabiele manifolden van een ander om transportkanalen te vormen.

Resultaten
De voorgestelde methode werd toegepast op het planaire 4D CR3BP-model (Aarde-Maan systeem). De simulatieresultaten bevatten:

• Reconstructie van Zero-Velocity Curves (ZVCs): De methode herstelt de energetisch toegestane regio's en verboden zones gedefinieerd door de Jacobi-constante.
• Lokale Geometrie bij Evenwichtspunten:
  • $L_1$ en $L_2$: De eigenfuncties schetsen duidelijk de transportpoorten. De onstabiele manifolden strekken zich naar buiten uit door de "halzen" (wat ontsnapping van de Aarde of toegang tot de buitenwereld mogelijk maakt), terwijl de stabiele manifolden naar binnen buigen (wat trajecten vangt).
  • $L_3$: Vertoont fragiele manifolden die naar externe regio's lekken in plaats van sterke transportcorridors te vormen.
  • $L_4$ en $L_5$: De centrum-eigenfuncties domineren en produceren gesloten lussen die overeenkomen met stabiele tadpole- en Trojan-banen, consistent met de lineaire stabiliteit van deze driehoekige punten.
• Globale Fase-Ruimte Partitionering: De nul-niveau curves van de eigenfuncties partitioneren de fase-ruimte spectraal. In tegenstelling tot energie-gebaseerde Jacobi-oppervlakken die definiëren waar beweging mogelijk is, onthullen de Koopman-curves hoe trajecten asymptotisch georganiseerd worden binnen die regio's, waarbij zij separatrices en invariante grenzen identificeren.
• Transportmechanisme: De resultaten illustreren een transportkanaal tussen $L_1$ en $L_2$, waarbij de onstabiele manifold van $L_1$ en de stabiele manifold van $L_2$ een invariante buis begrenzen, wat transfer tussen de Aarde- en Maanregio faciliteert.

Betekenis en Claims
Het artikel claimt dat de spectrale Koopman-benadering een krachtig instrument biedt voor de cislunaire dynamica door de focus te verschuiven van individuele trajecten naar de evolutie van observables. De primaire betekenis ligt in het vermogen om:

1. Globale Geometrie te Herstellen: De globale fase-ruimtegeometrie van de CR3BP te reconstrueren uit tijdreeksdata, waarbij nietlineaire structuren worden gevangen die lokale methoden missen.
2. Downstream Taken Mogelijk te Maken: Een fundament te leggen voor onzekerheidsvoortplanting, bereikbaarheidsanalyse en controleontwerp. Omdat de Koopman-operator lineair is, postuleren de auteurs dat onzekerheidsanalyse en controlestrategie-ontwerp voor het nietlineaire CR3BP vereenvoudigd en schaalbaarder kunnen worden gemaakt vergeleken met traditionele nietlineaire methoden.
3. Data-gedreven Ontdekking: Een nieuwe richting te openen voor data-gedreven ontdekking in de ruimtedynamica, waarbij de lineariteit van de operator de analyse en het ontwerp van controlestrategieën voor inherent nietlineaire systemen mogelijk maakt.

De auteurs concluderen dat hoewel het huidige werk zich richt op het herstellen van de geometrie van de CR3BP, het framework veelbelovend is voor toekomstige toepassingen in real-time trajectoptimalisatie en robuuste controle in complexe multi-body omgevingen.