Tag: free return trajectory

Conferenza Traiettorie di ritorno libero

By articolidiastronomia on 28/10/2023

Sapere quando è possibile andare verso la Luna è solo una prima parte del problema: per andarci occorre progettare anche una traiettoria che sia compatibile con i requisiti di missione, conceda un alto grado di affidabilità e minimizzi il rischio.

Tra le famiglie di traiettorie che rispondevano ai requisiti minimi di missione espressi sopra, la NASA scelse quelle chiamate Traiettorie di Ritorno Libero (Free Return Trajectory). Queste traiettorie vennero sfruttate dalle prime missioni Apollo: per ultime missioni invece, quando maggiore era la confidenza con il viaggio, vennero utilizzate traiettorie differenti.

Interessante notare che queste traiettorie erano note già alla fine del XIX secolo grazie al lavoro di alcuni scienziati russi e vennero in seguito riadattate al progetto Apollo. La conferenza intende approfondire questo argomento facendo uso di modelli 3D, molta lavagna e meccanica celeste.

La conferenza si terrà 10 Novembre presso la sede del GAV – Gruppo Astrofili Villasanta in via Bestetti 8 – Villasanta con inizio alle ore 21:00 circa.

Piazza Papa Giovanni Ventitreesimo

Alessandro Fumagalli

Traiettoria di ritorno libero – Conclusione

By articolidiastronomia on 19/09/2021 • ( 1 commento )

Riprendiamo l’analisi degli ultimi due metodi risolutivi per il tracciamento della traiettoria di ritorno libero.

Analisi R3B (Restricted 3 Body). Questo metodo sfrutta l’approssimazione del problema dei tre corpi: la Terra, la Luna e l’Apollo. Dato che la massa dell’Apollo è molto minore rispetto alle masse degli altri due oggetti celesti (Terra e Luna), i suoi effetti sugli altri oggetti possono essere ignorati. Consideriamo quindi il moto del sistema Terra-Luna intorno al comune centro di massa (CoM).

Si può dimostrare (qui) che le equazioni del moto della massa più piccola (la capsula m di massa trascurabile rispetto a m₁ e m₂) sono le seguenti:

l sistema di equazioni riportato sopra rappresenta un sistema di equazioni differenziali ordinarie non lineari in (x, y, z) con i seguenti parametri:

\begin{cases} \mu_{1} = G m_{1}\\ \mu_{2} = G m_{2}\\ \end{cases} \\[5pt] \begin{cases} \pi_{1} = \frac{m_{1}}{m_{1} + m_{2}}\\ \pi_{2} = \frac{m_{2}}{m_{1} + m_{2}} \end{cases} \\[5pt] \Omega= \sqrt{\frac{\mu_{2}}{r_{12}^3}} \\[5pt] r_{12}= \sqrt{(x_{1}-x_{2})^2} \\[5pt] r_{1}= \sqrt{(x+\pi_{2}r_{12})^2+y^2} \\[5pt] r_{2}= \sqrt{(x+\pi_{1}r_{12})^2+y^2}

Nel sistema di equazioni sopra riportato:

x₁ è la posizione di m₁ rispetto al CoM
x₂ è la posizione di m₂ rispetto al CoM
$\Omega$ è la velocità angolare riferita al sistema corotante Terra – Luna rispetto al CoM dei due corpi.
r₁₂ è la distanza fra CoM della Terra e la Luna
r₁ rappresenta la distanza di m₁ dal centro di gravità
r₂ rappresenta la distanza di m₂ dal centro di gravità
La derivata seconda ( $\ddot{\boldsymbol{x}}, \ddot{\boldsymbol{y}}, \ddot{\boldsymbol{z}}$ ) rappresenta l’accelerazione.

Risolvendo il sistema in si ricava la posizione (x, y, z), velocità ( $\dot{\boldsymbol{x}}, \dot{\boldsymbol{y}}, \dot{\boldsymbol{z}}$ ) e accelerazione $\ddot{\boldsymbol{x}}, \ddot{\boldsymbol{y}}, \ddot{\boldsymbol{z}}$ ) della capsula.

La soluzione si ottiene per integrazione numerica, come ad esempio il metodo di iterativo Runge-Kutta che si ottiene applicando i seguenti passi:

- 1 Si parte dalla condizione iniziale (x₀, y₀, z₀), per esempio la posizione del CSM al momento della  TLI per t = t₀

- 2 Si calcola l'accelerazione della sonda sotto l'azione gravitazionale di Luna e della Terra per un piccolo quanto di tempo t (dipende dalla precisione).

- 3 Si ricalcolano l'equazioni del moto per trovare i nuovi valori delle variabili di stato di posizione (x, y, z) e velocità (v_x, v_y, v_z).  Si annotano i nuovi valori: essi rappresentano la nuova posizione della capsula nello spazio. Si dice che si 'propaga' la soluzione.

- 4 Si passa al quanto successivo di tempo t = t +  Δt e si torna al punto 2. Una volta raccolto un set di punti (x, y, z) sufficienti per l'analisi si esce

- 5 Si disegna la traiettoria di (x, y, z) in un sistema di coordinate cartesiane.

La soluzione del sistema che si ottiene è molto sensibile alle condizioni iniziali di ( $\overrightarrow{x}, \overrightarrow{v}$ ).

In generale la traiettoria dell’Apollo nello spazio 3D non è detto che sia posizionata sullo stesso piano dell’orbita lunare, Intuitivamente possiamo immaginare di ‘stirare’ la traiettoria disegnata nel piano (vedi analisi RC3B) verso la direzione positiva e negativa dell’asse z. Di quanto viene modificata dipende dall’angolo formato fra i due piani:

il piano orbitale lunare rispetto all’eclittica
l’azumuth di lancio del Saturn V, e quindi dell’inclinazione del piano orbitale del CSM in orbita LEO.

Come già visto, i due piani hanno in comune il punto antipodale, dove avviene la TLI.

Analisi RC3B (Restricted Circular 3 Body). Se per ipotesi il moto della capsula è planare ( $\ddot{\boldsymbol{z}}=0$ ) (ovvero la traiettoria giace sullo stesso piano del moto dei 2 corpi celesti di massa più grande) possiamo trascurare l’ultima equazione.

Sistema di riferimento usato per l’equazione del moto. Nel caso di analisi R3B il punto (x, y, z) è libero di muoversi in 3 dimensioni.
Nel caso RC3B per ipotesi il movimento è limitato nel piano del moto di m₁ e m₂

Fonte: *Design of Trajectory and Perturbation Analysis for Satellite Orbital Parameters – Saurabh PandeyRishabh KumarAman Dalmia, Springer*

Il fatto di rimuovere la componente z significa che nei grafici seguenti viene disegnata la proiezione lungo il piano orbitale della Luna della traiettoria del CSM.

Bisogna calcolare in maniera molto precisa posizione e velocità iniziali affinché la traiettoria dell’Apollo circumnavighi la Luna e ritorni in prossimità della Terra con una condizione compatibile con la manovra di rientro. Condizioni anche di poco fuori dalla condizione ottimale fanno percorrere alla sonda traiettorie che andranno in direzione della Terra, ma non necessariamente nel punto di ritorno previsto.

La traiettoria della sonda potrebbe per esempio descrivere anche un moto oscillatorio tra i due corpi, oppure portarsi nella traiettoria giusta per il rientro ma dopo parecchi oscillazioni, quindi un tempo troppo lungo affinché si possa garantire la sopravvivenza dell’equipaggio.

Le figure seguenti mostrano i risultati di diverse simulazioni fatte dall’autore grazie al software sviluppato da Adrien (blog.nodraak.fr) in python. In entrambi i diagrammi viene tracciata una famiglia di traiettorie FRT considerando un periodo di dieci giorni di viaggio. Nel primo diagramma viene mostrato come cambia la traiettoria della FRT al variare del Δv a parità di angolo α di inserimento della TLI; viceversa, il secondo grafico mostra come cambia la FRT al variare dell’angolo α a parità di Δv.

In tutte le simulazioni l’Apollo si trova in un’orbita LEO ad un’altezza di cento miglia nautiche.

I due ingrandimenti a sinistra mostrano come la traiettoria cambi di molto al variare delle condizioni iniziali.

Variare di poco l’angolo α fa la differenza tra schiantarsi sulla Terra in fase di rientro o avvicinarsi entrare nella sfera di influenza terrestre. Anche se in ogni caso servono correzioni di rotta per garantire l’inserimento nel corridoio di rientro, la simulazione permette di escludere a priori alcuni parametri orbitali.

Durante lo sviluppo del progetto Apollo, la responsabile della pianificazione e verifica delle traiettorie FRT era Frances “Poppy” Northcutt: la prima donna ad essere ammessa al Mission Control Room (MCC), all’epoca un ambiente maschile. Poppy si era laureata in matematica in Texas perché voleva evitare lavori da donne ed avere una marcia in più nella ricerca di un lavoro.

Fonte: https://www.science101.com/woman-mission-control-story-apollo-missions/

Iniziò a lavorare alla TRW (un’azienda del settore avionico) prima di collaborare con NASA per il progetto Apollo come “donna computer” (computress), poi passò al Mission Planning & Analysis Room diventando la prima donna a lavorare al MCC. Lei divenne la responsabile sia della definizione del profilo di volo delle FRT per l’Apollo 8 che per riprogettare la FRT in tempo reale per la missione di salvataggio dell’Apollo 13: per i suoi contributi ha ricevuto La Medaglia Presidenziale delle Libertà.

Da quando il progetto Apollo è stato chiuso, Poppy ha lasciato la NASA per occuparsi della sua vera passione: la legge. Dal 1984 esercita come avvocato.

Bibliografia

Digital Apollo Human and Machine in Spaceflight, MIT Press
How Apollo flew to the Moon, David Woods Springer

Traiettoria di ritorno libero – Parte III

By articolidiastronomia on 12/09/2021 • ( 1 commento )

L’analisi di una traiettoria FRT è un esempio di applicazione del problema degli n corpi, che non ammette soluzioni in forma chiusa, ma si può modellare implementando tecniche di integrazione numerica.

Nel nostro caso possiamo identificare tre diversi metodi risolutivi:

Analisi per approssimazione conica (PCA)
Analisi del problema dei tre corpi ristretto (R3B)
Analisi del problema dei tre corpi ristretto su orbite circolari (RC3B)

Il primo è il modo più semplice ma anche più approssimativo: si ipotizza che la capsula si muova sotto l’influenza gravitazionale di un corpo a seconda della sfera di influenza ed è usato per avere un’idea iniziale della traiettoria che compirà la sonda. Le altre due soluzioni si basano sul problema ristretto dei 3 corpi. L’ipotesi di fondo è che la massa della capsula è molto minore della massa degli altri 2 corpi (Terra e Luna) e quindi sarà soggetta all’accelerazione di entrambi i corpi a seconda della loro distanza: la traiettoria è più precisa del metodo PCA.

Analisi PCA: questo metodo considera solo gli effetti del corpo gravitazionalmente dominante in cui si trova in quel momento la traiettoria della capsula. Il corpo dominante è definito dalla sfera di influenza R_SOI (rispetto al Sole) secondo la relazione:

Dove R è la distanza Terra – Sole, M_x è all’occorrenza la massa della Terra o della Luna, M_S è la massa del Sole ed R è la distanza Terra – Sole.

*Pianeta*	*SOI*
Terra	9,24643 * 10⁵ Km
Luna	0,66168 * 10⁵ Km

Valori di R_SOI per la Terra e la Luna. Notare come il valore della SOI della Terra è circa un’ordine di grandezza maggiore rispetto alla SOI lunare

Dato che la R_SOI terrestre include la R_SOI lunare, questo approccio esclude a priori l’influenza terrestre quando l’Apollo si trova all’interno della sfera gravitazionale lunare.

Per questo motivo l’analisi PCA è impiegata quando ancora si sta affrontando la fase di pianificazione iniziale del viaggio.

Il procedimento prevede di studiare separatamente il tratto del percorso dell’Apollo in regime della dinamica dei due corpi (uno alla volta) all’interno della propria SOI.

La traiettoria PCA si sviluppa in tre fasi:

traiettoria geocentrica: l’Apollo subisce solo l’influenza gravitazionale terrestre
traiettoria selenocentrica: l’Apollo subisce solo l’influenza gravitazionale lunare
la traiettoria finale è composta dall’unione delle due traiettorie

Grazie alla TLI, nella traiettoria geocentrica il CSM percorre un’orbita ellittica fino a trovarsi a distanza r₁ dalla Terra con una velocità $\overrightarrow{v_{1}}$ (dell’ordine di qualche centinaio m/s) cioè al raggiungimento della SOI nel punto S intersezione con la sfera r₂=R_SOI.

Fonte: *SEEDS – Astrodynamics, Guido Colasurdo e Giulio Avanzini – Politecnico di Torino*

Il punto S marca il passaggio nella seconda fase: la traiettoria selenocentrica: consideriamo ora come sistema di riferimento il piano selenocentrico non rotante e osserviamo cosa succede fino all’uscita dalla SOI al punto T.

Dal punto di vista di un osservatore solidale con la Luna essa si muove lungo la traiettoria intorno alla Terra con velocità $\overrightarrow{v_{m}}$ e vede arrivare l’Apollo con velocità $\overrightarrow{v_{2}} = (\overrightarrow{v_{1}} - \overrightarrow{v_{m}})$ composizione vettoriale delle due velocità. La velocità risultante d’ingresso $\overrightarrow{v_{2}}$ deve essere pianificata in fase di progetto tale da raggiungere la SOI con una velocità maggiore della velocità di fuga lunare al confine della SOI, in tal caso l’Apollo arriva intorno alla Luna con un’orbita iperbolica (con fuoco nella Luna).

A mano a mano che ci avvicina al pericinzio, la velocità $\overrightarrow{v_{2}}$ fino a $\overrightarrow{v_{max}}$ quindi si allontana e raggiunge il confine della SOI con velocità di uscita $| \overrightarrow{v_{4}}| = |\overrightarrow{v_{2}}|$ ma direzione differente, ovvero lascia la SOI con velocità geocentrica $\overrightarrow{v_{5}} = (\overrightarrow{v_{4}} + \overrightarrow{v_{m}})$ in direzione della Terra.

Esempio di traiettoria per la Luna ove sono evidenziati nel ramo di andata la combinazione della traiettoria geocentrica e selenocentrica.

Fonte: *Design of Trajectory and Perturbation Analysis for Satellite Orbital Parameters – Saurabh PandeyRishabh KumarAman Dalmia, Springer*

L’Apollo effettua un fly-by lunare: il cambiamento della direzione del vettore velocità dipende da quanto siamo vicini alla Luna: se il fly-by viene effettuato più distante, allora il cambiamento di direzione sarà minore, se il fly-by viene effettuato più vicino al perilunio r₃ più stretta sarà il cambiamento di direzione del moto, e quindi un vantaggio in prospettiva di un allunaggio.

In ogni caso sono possibili due casi estremi:

r₃ < 1738 Km (raggio lunare) la sonda si schianta sulla Luna.
r₃ > 1738 Km (raggio lunare) l’Apollo può entrare in orbita lunare (Lunar Orbit Insertion, LOI) o effettuare un fly-by.

Verso la Terra⁵. Una volta abbandonata la SOI, l’Apollo si inserisce nel ramo di ritorno della FRT: quest’ultimo è ancora un arco di ellisse. Per compensare gli errori sempre presenti, alcune correzioni di rotta lungo la traiettoria distribuiti fra ramo di andata e ritorno sono possibili, anche se non necessarie. Le correzioni orbitali sono necessarie per includere:

errori di approssimazione
per accelerare il rientro del CSM (come nel caso dell’Apollo 13)
per sincronizzarsi con il punto esatto di rientro a Terra della capsula (la Terra ruota su sé stessa mentre la sonda è in navigazione e la procedura di rientro in atmosfera deve avvenire all’interno di un corridoio preciso).

(continua)

Nota

⁵Nel presente articolo non viene considera alcuna manovra orbitale di inserimento lunare (LOI) o reinserimento in traiettoria terrestre (TEI) dopo il docking fra LM e CSM. Si considera solo la traiettoria in caduta libera in maniera simile a quanto successo con la missione Apollo 13 (a parte le correzioni di rotta per l’immissione in FRT sul ramo di andata e la correzione di rotta sul ramo di ritorno per velocizzarne il rientro).

Bibliografia