Sintesi Tecnica: Il Problema della Massa Iniziale Massima nella Progettazione di Traiettorie a Bassa Spinta

Definizione del Problema

Il documento affronta una lacuna nella teoria del controllo ottimo a bassa spinta: la mancanza di una formalizzazione esplicita e di un trattamento analitico del problema della Massa Iniziale Massima (MIM). Mentre il campo è dominato dalle formulazioni di propellente minimo (massa finale massima) e di trasferimento a tempo minimo, il problema MIM pone una domanda distinta: per una determinata geometria di trasferimento, uno stato finale fissato e un tempo di trasferimento prestabilito, qual è la massa iniziale più grande da cui una sonda a bassa spinta può completare con successo il trasferimento?

Concettualmente, questo inverte il paradigma standard "data una massa iniziale, minimizzare tempo/propellente" in "dato un tempo di volo, massimizzare la massa iniziale". Sebbene questa prospettiva sia stata utilizzata informalmente in competizioni (ad es. GTOC) e per studi di fattibilità per l'estrazione mineraria di asteroidi, non è stata rigorosamente stabilita come una formulazione di controllo ottimo autonoma dotata delle proprie condizioni necessarie e proprietà strutturali.

Metodologia

Gli autori utilizzano il Principio di Minimum di Pontryagin (PMP) per derivare le condizioni necessarie del primo ordine per il problema MIM e confrontarle strutturalmente con il classico problema a tempo minimo.

1. Formulazione:

   • Dinamica: Il sistema è modellato utilizzando la dinamica standard a bassa spinta dove il vettore di stato $\mathbf{x}$ (ad es. elementi equinoziali modificati) e la massa $m$ evolvono sotto un comando di accelerazione $u \in [0,1]$ e una direzione della spinta $\hat{\mathbf{i}}$.
   • Obiettivo: Massimizzare la massa iniziale $m(t_0)$ soggetto a un tempo di trasferimento fissato $t_f - t_0$ e condizioni al contorno $\mathbf{x}(t_0) = \mathbf{x}_0, \mathbf{x}(t_f) = \mathbf{x}_f$. Questo viene convertito in una minimizzazione di $J = -m(t_0)$.
   • Analisi dell'Hamiltoniana: Gli autori derivano l'Hamiltoniana e la relativa funzione di commutazione (switching function). Una scoperta chiave è che il comando ottimale è sempre $u^*=1$ (spinta massima) nel caso non degenere, e la direzione della spinta segue una legge del vettore primer ($\hat{\mathbf{i}}^* \propto -\mathbf{B}(\mathbf{x})^T\boldsymbol{\lambda}$). Ciò risulta in una struttura di controllo continua, evitando i profili bang-bang discontinui spesso visti in altre formulazioni, il che facilita la convergenza rapida sia per i metodi numerici diretti che indiretti.

2. Trasversalità e Condizioni di Gauge:

   • Il documento riesamina criticamente la condizione di trasversalità per i problemi a tempo libero. Sostiene che la comune assunzione secondo cui l'Hamiltoniana debba essere zero al tempo finale ($H(t_f) + \boldsymbol{\lambda}(t_f)^T\dot{\mathbf{x}}_f(t_f) = 0$) per i problemi a tempo minimo sia una scelta di gauge piuttosto che una condizione necessaria derivata dal PMP per gli ottimi di confine. La vera condizione è una disuguaglianza unidirezionale ($H(t_f) + \boldsymbol{\lambda}(t_f)^T\dot{\mathbf{x}}_f(t_f) \le 0$).
   • Questa intuizione permette la normalizzazione del vettore costato $\lambda_m(t_0) = 1$ senza violare le condizioni necessarie.

3. Corrispondenza degli Estremali:

   • Utilizzando le condizioni derivate, gli autori dimostrano il Teorema 1: Esiste una corrispondenza biunivoca tra gli estremi del problema a tempo minimo (con massa iniziale fissata) e il problema della massa iniziale massima (con tempo di volo fissato). Nello specifico, una traiettoria che è ottima per il tempo dato una certa massa, è anche ottima per la massa risultante per un dato tempo di volo, e viceversa.

4. Applicazione dell'Omotopia:

   • Gli autori applicano la formulazione MIM al problema a tempo minimo a più rivoluzioni. Costruiscono un percorso di omotopia variando il numero di rivoluzioni ($N_{rev}$) mentre tracciano la superficie della massa iniziale massima $M(t_f, N_{rev})$. Questo approccio è confrontato con le tradizionali omotopie di continuazione della spinta, che spesso soffrono di cattivo condizionamento nei casi multi-rivoluzione.

Risultati Chiave

• Proprietà Numeriche: Il problema MIM produce un problema ai valori al contorno ben posto con sette equazioni e sette incognite (per endpoint fissi), evitando i problemi di indeterminazione di scala comuni nei metodi di shooting indiretti a tempo minimo. Il controllo ottimo è strettamente a spinta massima, garantendo una convergenza numerica fluida.
• Chiarificazione Teorica: Il documento chiarisce che la condizione di uguaglianza dell'Hamiltoniana terminale non è una condizione necessaria per l'ottimalità temporale, ma una scelta di normalizzazione. Ciò risolve un "errore comune" nella letteratura e valida l'uso di condizioni di gauge alternative per migliorare la stabilità numerica.
• Visualizzazione del Paesaggio di Soluzione: In un trasferimento di riferimento da un'orbita di trasferimento geostazionaria (GTO) a un'orbita geostazionaria (GEO), l'omotopia MIM ha mappato con successo il paesaggio delle soluzioni. Il metodo ha rivelato una superficie di massa iniziale massima continua e differenziabile.
  • Ha identificato molteplici traiettorie estreme, inclusi minimi e massimi locali del tempo di volo lungo le curve iso-massa.
  • Ha individuato con successo la soluzione globale ottima per il tempo (massa iniziale di 1500 kg) e ha dimostrato che alcune soluzioni precedentemente riportate (ad es. la soluzione E del benchmark) corrispondono a massimi locali del tempo di volo, non a minimi.
• Efficienza dell'Omotopia: La via di continuazione nel framework MIM (variando $N_{rev}$ per tracciare le curve iso-massa) è stata osservata essere fluida e ben condizionata, offrendo un'alternativa più robusta ai metodi di continuazione della spinta per navigare in complessi paesaggi multi-rivoluzione.

Significato e Rivendicazioni

Il documento sostiene che il problema della Massa Iniziale Massima meriti di essere elevato a una formulazione di controllo ottimo di "prima classe" nella progettazione di traiettorie a bassa spinta. La sua importanza risiede in tre aree principali:

1. Primitivo di Progettazione Complementare: Sebbene strutturalmente correlato al problema a tempo minimo, la formulazione MIM offre una prospettiva distinta che fissa il tempo e libera la massa. Ciò fornisce uno strumento complementare per la progettazione di missioni, particolarmente durante l'analisi preliminare dove devono essere valutati grandi set di trasferimenti.
2. Rigore Analitico: Fornendo la derivazione formale delle condizioni necessarie e correggendo la comprensione delle condizioni di trasversalità nei problemi a tempo ottimale, il documento fornisce una base teorica più rigorosa per l'ottimizzazione a bassa spinta.
3. Utilità Algoritmica: La formulazione MIM consente la costruzione di percorsi di omotopia fluidi che sono meno soggetti al cattivo condizionamento tipico dei metodi di continuazione della spinta. Ciò permette l'esplorazione sistematica di spazi di soluzione complessi (come i trasferimenti multi-rivoluzione) e la generazione di ampi database di traiettorie ottime per la modellazione surrogata e l'analisi di raggiungibilità.

Gli autori concludono che l'adozione della formulazione MIM come standard di progettazione può migliorare l'efficienza e la robustezza della progettazione di missioni a bassa spinta, specialmente negli approcci basati sui dati e nell'analisi preliminare della missione.