Improved Directional State Transition Tensors for Accurate Aerocapture Performance Analysis

Technische Samenvatting: Verbeterde Richtingsgeoriënteerde Toestands-Transitietensoren voor Nauwkeurige Aerocapture-Prestatieanalyse

Probleemstelling
Aerocapture is een kritieke technologie voor missies naar de buitenplaneten, die brandstofmassabesparingen biedt ten opzichte van propulsieve baaninjectie. De dynamica van aerocapture wordt echter gedomineerd door niet-conservatieve aerodynamische krachten die aanzienlijk variëren gedurende de baan, waardoor inherente niet-lineariteiten ontstaan. Traditionele methoden voor lineaire onzekerheidspropagatie, zoals Lineaire Covariantieanalyse (LCA), slagen er niet in deze niet-lineaire effecten nauwkeurig te vangen, vooral wanneer de invoeronzekerheden groot zijn of complexe atmosferische modellen worden gebruikt. Hoewel er hogere-orde methoden bestaan, zoals Polynoom Chaos-expansie (PCE) en de Unscented Transform (UT), lijden deze onder de vloek van de dimensionaliteit of hoge rekenkosten die implementatie aan boord belemmeren.

Toestands-Transitietensoren (STTs), die hogere-orde Taylor-reeksontwikkelingen van de oplossingsstroom representeren, bieden een rekenkundig efficiënt semi-analytisch alternatief voor baanscenario's. Hun toepassing op aerocapture is echter beperkt door de exponentiële groei in het aantal differentiaalvergelijkingen die nodig zijn om hogere-orde STTs te berekenen naarmate de toestandsdimensionaliteit toeneemt. Richtingsgeoriënteerde Toestands-Transitietensoren (DSTTs) zijn ontwikkeld om deze kosten te beperken door de toestand te projecteren op een basis met gereduceerde dimensionaliteit langs een dynamisch belangrijke richting. Vorig werk heeft aangetoond dat voor baandynamica de richting van maximale lineaire rek (de dominante eigenvector van de Cauchy-Green-tensor van de tweede orde, of CGT) relatief constant is en uitlijnt met hogere-orde rekrichtingen. Dit artikel onderzoekt of deze aannames ook gelden voor aerocapture en streeft naar verbetering van de DSTT-prestaties door effectievere reductiebasissen te identificeren.

Methodologie
De auteurs analyseren het aerocapture-probleem met behulp van een toestandsvector met 7 vrijheidsgraden, waaronder radiale afstand, breedtegraad, lengtegraad, snelheid, vluchtpadhoek, koershoek en de natuurlijke logaritme van de atmosferische dichtheid. De studie maakt gebruik van twee scenario's gebaseerd op een voorgestelde Uranus-Orbiter-Probe-missie: een basisinval en een inval met hoge snelheid.

De kernmethodologie omvat drie hoofdstappen:

1. Dynamische Decompositie: De auteurs decomponeren de aerocapture-dynamica in conservatieve (zwaartekracht) en dissipatieve (aerodynamische) componenten. Door de CGT's van deze afzonderlijke componenten te analyseren, tonen ze aan dat de richting van maximale lineaire rek in aerocapture tijdvariabel is, in tegenstelling tot bij baandynamica, vanwege de dominantie van aerodynamische krachten tijdens de piek van de dynamische druk.
2. Identificatie van Nieuwe Bases: Gezien het feit dat de richting van lineaire rek (CGT van de tweede orde) niet uitlijnt met hogere-orde rekrichtingen in aerocapture, stellen de auteurs het gebruik van Cauchy-Green-tensoren van hogere orde (HOCGT's) voor. Ze introduceren twee uitgebreide formuleringen:
   • Selectieve CGT's (sCGT's): Deze identificeren de richting van maximale rek voor een specifieke subset van toestanden (bijvoorbeeld positie, snelheid, vluchtpadhoek).
   • Kwantiteit-van-Interesse CGT's (qCGT's): Deze identificeren de richting van maximale rek voor een specifieke niet-lineaire functie van de toestand, zoals specifieke baanenergie of apoapsisstraal.
     De dominante eigenvectoren van deze tensoren worden berekend met de Shifted-Symmetric Higher-Order Power Method (SS-HOPM).
3. DSTT-construktie en Validatie: De auteurs construeren DSTT's van de derde orde met deze nieuwe richtingen (hoDSTT, sDSTT, qDSTT) en vergelijken ze met traditionele DSTT's die zijn opgebouwd met eigenvectoren van de CGT van de tweede orde (1-DSTT, 3-DSTT, 6-DSTT) en volledige STT's. De prestaties worden geëvalueerd aan de hand van de genormaliseerde Frobenius-normfout, de nauwkeurigheid van de propagatie van perturbaties en de voorspellingsnauwkeurigheid voor de apoapsisstraal en specifieke energie aan het einde.

Belangrijkste Bijdragen

• Identificatie van Tijdvariabele Rek: Het artikel toont aan dat voor aerocapture de richting van maximale lineaire rek niet constant is en niet uitlijnt met hogere-orde rekrichtingen, waardoor eerdere aannames afgeleid uit baandynamica ongeldig worden.
• Uitgebreide HOCGT-formuleringen: De ontwikkeling van sCGT's en qCGT's maakt het mogelijk om gevoeligheidsrichtingen te identificeren die zijn toegesneden op specifieke subsets van de toestand of specifieke grootheden van belang, in plaats van alleen de volledige toestandsnorm.
• Verbeterde DSTT-nauwkeurigheid: De studie toont aan dat het construeren van DSTT's langs richtingen afgeleid van HOCGT's van de derde en vierde orde (en hun uitgebreide varianten) de fout in de Frobenius-norm aanzienlijk verlaagt in vergelijking met traditionele DSTT's, met name in het gebied van piekdynamische druk waar niet-lineariteiten het sterkst zijn.
• Trade-off-analyse: De auteurs bieden een gedetailleerde analyse die aantoont dat het gebruik van één latente dimensie langs een richting met hogere-orde gevoeligheid betere nauwkeurigheid kan opleveren dan het gebruik van meerdere dimensies (bijvoorbeeld 6-DSTT) uitsluitend gebaseerd op lineaire rek.

Resultaten

• Benaderingsnauwkeurigheid: DSTT's die zijn opgebouwd met de nieuwe hogere-orde richtingen (hoDSTT, sDSTT, qDSTT) vertonen lagere genormaliseerde fouten in de Frobenius-norm dan traditionele DSTT's tijdens de kritieke fase van piekdynamische druk.
• Propagatie van Perturbaties: Voor perturbaties die zijn uitgelijnd met de richting van maximale niet-lineaire rek behouden de nieuwe DSTT's een hoge nauwkeurigheid. Zelfs voor perturbaties loodrecht op deze richting presteert de eendimensionale hoDSTT vergelijkbaar met of beter dan de Toestands-Transitiematrix (STM).
• Voorspelling van Grootheden van Belang:
  • Specifieke Energie: Het gebruik van sCGT's (selectie van samenstellende toestanden) of qCGT's (gebruik van de energiefunctie) voor het construeren van DSTT's verlaagde de mediaanfout in de energievoorspelling met ongeveer een halve orde van grootte in vergelijking met traditionele methoden.
  • Apoapsisstraal: De qDSTT-aanpak die direct gebruikmaakt van de functie voor de apoapsisstraal, presteerde slecht vanwege numerieke problemen die voortvloeien uit de discontinuïteit in de halve grote as in de buurt van de opvanggrens. De sDSTT-aanpak, die de samenstellende toestanden ($r, V, \gamma$) selecteert die worden gebruikt in de berekening van de apoapsis, presteerde echter beter dan andere DSTT-methoden bij het voorspellen van de apoapsisstraal.
• Rekenkundige Efficiëntie: De nieuwe DSTT's vereisen aanzienlijk minder termen om te integreren en te propagëren dan volledige STT's (bijvoorbeeld 14 termen voor een eendimensionale DSTT van de derde orde versus 343 termen voor een STT van de tweede orde), wat een rekenkundig efficiënt alternatief biedt voor geleiding en navigatie aan boord.

Betekenis
Het artikel stelt vast dat STT's en DSTT's praktische hulpmiddelen zijn voor de analyse van aerocapture-prestaties, mits de reductiebasis zorgvuldig wordt gekozen. Door verder te gaan dan de aanname dat richtingen van lineaire rek voldoende zijn, tonen de auteurs aan dat het benutten van dynamica van hogere orde en het aanpassen van de basis aan specifieke grootheden van belang leidt tot robuustere en nauwkeurigere onzekerheidspropagatie. Deze resultaten maken het gebruik mogelijk van efficiënte semi-analytische methoden voor stochastische geleiding en navigatie aan boord in sterk niet-lineaire, niet-conservatieve omgevingen, waarmee een kritieke kloof in de huidige mogelijkheden voor het ontwerpen van aerocapture-missies wordt gedicht. De voorgestelde methodologie is generaliseerbaar naar elk niet-conservatief, niet-lineair dynamisch systeem.