Improved Directional State Transition Tensors for Accurate Aerocapture Performance Analysis

Sintesi Tecnica: Tensori di Transizione di Stato Direzionali Migliorati per un'Analisi Precisa delle Prestazioni dell'Aerocattura

Enunciato del Problema
L'aerocattura è una tecnologia critica per le missioni verso i pianeti esterni, che offre risparmi di massa del propellente rispetto all'inserimento orbitale propulsivo. Tuttavia, la dinamica dell'aerocattura è dominata da forze aerodinamiche non conservative che variano significativamente lungo la traiettoria, creando non linearità intrinseche. I metodi tradizionali di propagazione dell'incertezza lineare, come l'Analisi della Covarianza Lineare (LCA), non riescono a catturare accuratamente questi effetti non lineari, specialmente quando le incertezze di ingresso sono elevate o sono utilizzati modelli atmosferici complessi. Sebbene esistano metodi di ordine superiore come l'Espansione del Caos Polinomiale (PCE) e la Trasformata Unscented (UT), questi soffrono della maledizione della dimensionalità o di costi computazionali elevati che ne ostacolano l'implementazione a bordo.

I Tensori di Transizione di Stato (STT), che rappresentano espansioni in serie di Taylor di ordine superiore del flusso della soluzione, offrono un'alternativa semi-analitica computazionalmente efficiente per scenari orbitali. Tuttavia, la loro applicazione all'aerocattura è stata limitata dalla crescita esponenziale del numero di equazioni differenziali richieste per calcolare STT di ordine superiore all'aumentare della dimensione dello stato. I Tensori di Transizione di Stato Direzionali (DSTT) sono stati sviluppati per mitigare questo costo proiettando lo stato in una base a dimensione ridotta lungo una direzione dinamicamente importante. Lavori precedenti hanno stabilito che, per la dinamica orbitale, la direzione di massimo allungamento lineare (il vettore proprio dominante del tensore di Cauchy-Green del secondo ordine, o CGT) è relativamente costante e si allinea con le direzioni di allungamento di ordine superiore. Questo articolo indaga se queste ipotesi valgano per l'aerocattura e cerca di migliorare le prestazioni dei DSTT identificando basi di riduzione più efficaci.

Metodologia
Gli autori analizzano il problema dell'aerocattura utilizzando un vettore di stato a 7 gradi di libertà che include la distanza radiale, la latitudine, la longitudine, la velocità, l'angolo di traiettoria di volo, l'angolo di assetto e il logaritmo naturale della densità atmosferica. Lo studio impiega due scenari basati su una missione proposta per un Orbiter-sonda di Urano: un ingresso di base e un ingresso ad alta velocità.

La metodologia centrale coinvolge tre fasi principali:

1. Decomposizione della Dinamica: Gli autori decompongono la dinamica dell'aerocattura in componenti conservative (gravitazionali) e dissipative (aerodinamiche). Analizzando i CGT di queste componenti separate, dimostrano che la direzione di massimo allungamento lineare è variabile nel tempo nell'aerocattura, a differenza della dinamica orbitale, a causa del dominio delle forze aerodinamiche durante il picco di pressione dinamica.
2. Identificazione di una Nuova Base: Riconoscendo che la direzione di allungamento lineare (CGT del secondo ordine) non si allinea con le direzioni di allungamento di ordine superiore nell'aerocattura, gli autori propongono l'uso di Tensori di Cauchy-Green di Ordine Superiore (HOCGT). Introducono due formulazioni aumentate:
   • CGT Selettivi (sCGT): Questi identificano la direzione di massimo allungamento per un sottoinsieme specifico di stati (ad esempio, posizione, velocità, angolo di traiettoria di volo).
   • CGT per Quantità di Interesse (qCGT): Questi identificano la direzione di massimo allungamento per una specifica funzione non lineare dello stato, come l'energia orbitale specifica o il raggio dell'apoapside.
     I vettori propri dominanti di questi tensori sono calcolati utilizzando il Metodo di Potenza di Ordine Superiore Simmetrico Spostato (SS-HOPM).
3. Costruzione e Validazione dei DSTT: Gli autori costruiscono DSTT del terzo ordine utilizzando queste nuove direzioni (hoDSTT, sDSTT, qDSTT) e li confrontano con i DSTT tradizionali costruiti utilizzando i vettori propri del CGT del secondo ordine (1-DSTT, 3-DSTT, 6-DSTT) e gli STT completi. Le prestazioni sono valutate utilizzando l'errore della norma di Frobenius normalizzata, l'accuratezza della propagazione delle perturbazioni e l'accuratezza della previsione per il raggio finale dell'apoapside e l'energia specifica.

Contributi Chiave

• Identificazione dell'Allungamento Variabile nel Tempo: Il documento dimostra che per l'aerocattura la direzione di massimo allungamento lineare non è costante e non si allinea con le direzioni di allungamento di ordine superiore, invalidando le ipotesi precedenti derivate dalla dinamica orbitale.
• Formulazioni Aumentate HOCGT: Lo sviluppo di sCGT e qCGT permette l'identificazione di direzioni di sensibilità adattate a sottoinsiemi specifici dello stato o a specifiche quantità di interesse, piuttosto che alla sola norma dello stato completo.
• Migliorata Accuratezza dei DSTT: Lo studio mostra che la costruzione di DSTT lungo direzioni derivate da HOCGT del terzo e quarto ordine (e le loro varianti aumentate) riduce significativamente l'errore della norma di Frobenius rispetto ai DSTT tradizionali, in particolare durante la regione di picco di pressione dinamica dove le non linearità sono più forti.
• Analisi dei Compromessi: Gli autori forniscono un'analisi dettagliata che mostra come l'uso di una singola dimensione latente lungo una direzione sensibile di ordine superiore possa produrre una migliore accuratezza rispetto all'uso di più dimensioni (ad esempio, 6-DSTT) basate esclusivamente sull'allungamento lineare.

Risultati

• Accuratezza dell'Approssimazione: I DSTT costruiti utilizzando le nuove direzioni di ordine superiore (hoDSTT, sDSTT, qDSTT) presentano errori della norma di Frobenius normalizzata inferiori rispetto ai DSTT tradizionali durante la fase critica di picco di pressione dinamica.
• Propagazione delle Perturbazioni: Per le perturbazioni allineate con la direzione di massimo allungamento non lineare, i nuovi DSTT mantengono un'alta accuratezza. Anche per le perturbazioni ortogonali a questa direzione, il singolo-dimensionale hoDSTT performa in modo comparabile o migliore rispetto alla Matrice di Transizione di Stato (STM).
• Previsione della Quantità di Interesse:
  • Energia Specifica: L'uso di sCGT (selezionando stati costituenti) o qCGT (utilizzando la funzione di energia) per costruire DSTT ha ridotto l'errore mediano di previsione dell'energia di circa mezzo ordine di grandezza rispetto ai metodi tradizionali.
  • Raggio dell'Apoapside: L'approccio qDSTT che utilizza direttamente la funzione del raggio dell'apoapside ha performato male a causa di problemi numerici derivanti dalla discontinuità nel semiasse maggiore vicino al confine di cattura. Tuttavia, l'approccio sDSTT, che seleziona gli stati costituenti ($r, V, \gamma$) utilizzati nel calcolo dell'apoapside, ha superato gli altri metodi DSTT nella previsione del raggio dell'apoapside.
• Efficienza Computazionale: I nuovi DSTT richiedono termini significativamente inferiori da integrare e propagare rispetto agli STT completi (ad esempio, 14 termini per un DSTT a singola dimensione del terzo ordine contro 343 termini per un STT del secondo ordine), offrendo un'alternativa computazionalmente efficiente per la guida e la navigazione a bordo.

Significato
Il documento stabilisce che gli STT e i DSTT sono strumenti pratici per l'analisi delle prestazioni dell'aerocattura, a condizione che la base di riduzione sia scelta con cura. Superando l'assunzione che le direzioni di allungamento lineare siano sufficienti, gli autori dimostrano che sfruttare la dinamica di ordine superiore e adattare la base a specifiche quantità di interesse porta a una propagazione dell'incertezza più robusta e accurata. Questi risultati abilitano l'uso di metodi semi-analitici efficienti per la guida e la navigazione stocastica a bordo in ambienti altamente non lineari e non conservativi, colmando una lacuna critica nelle attuali capacità di progettazione delle missioni di aerocattura. La metodologia proposta è generalizzabile a qualsiasi sistema dinamico non conservativo e non lineare.