Calcolo quantità propellente per raggiungere la Luna: ottimizza le missioni spaziali

www.matematicamente.it M. Giancola Calcolo della quantità minima…

Calcolo della quantità minima di propellente

necessaria a raggiungere la Luna

di Marco Giancola

Un aspetto molto importante della progettazione di una missione spaziale è rappresentato dalla

minimizzazione della quantità di propellente di cui deve essere dotato il veicolo spaziale. Infatti la

massima massa che può essere trasportata dal veicolo costituisce uno dei principali vincoli della

missione.

Un veicolo spaziale, per modificare la sua traiettoria (ad esempio per passare da un’orbita circolare

ad una ellittica), deve eseguire delle manovre impulsive, ossia deve azionare il proprio propulsore,

il quale genera una spinta che modifica bruscamente la velocità del veicolo e lo porta sulla nuova

traiettoria desiderata. Ad ogni manovra impulsiva corrisponde quindi una certa variazione di

velocità, detta anche impulso, che è pari a r r

∆ = −

V V V

2 1

r r

dove e sono le velocità dell’astronave rispettivamente prima e dopo la manovra (con V, in

V V

1 2

genere, si indica la velocità del veicolo rispetto ad un sistema di riferimento concentrico e solidale

con il corpo celeste che, col proprio campo gravitazionale, influenza maggiormente il moto del

veicolo). ∆V

La somma dei associati alle varie manovre impulsive, che il veicolo deve compiere, costituisce

∆V

il totale di cui necessita il veicolo per poter percorrere la sua traiettoria. E’ ovvio che il consumo

∆V ∆V

di propellente è direttamente proporzionale al totale, quindi minimizzare il totale equivale a

minimizzare il consumo di propellente.

∆V,

Noto il la quantità di propellente di cui necessita il veicolo si ricava dall’equazione di

Tsiolkovsky: ∆

 

V

 

= −

v

m m e 1

 

p  

Dove è la massa del propellente, m quella del veicolo privo di propellente, mentre v è la

m p

velocità (rispetto al veicolo) con la quale viene espulso il propellente combusto.

Nel 1991, il matematico Theodore H. Sweetser, ricercatore del Jet Propulsion Laboratory, calcolò il

∆V totale minimo teorico necessario per realizzare un trasferimento lunare. Il suo calcolo è basato

sull’analisi della costante di Jacobi. La formula che esprime tale costante è uno dei risultati che ha

ottenuto Jacobi, nella prima metà del 19° secolo, studiando il problema ristretto dei tre corpi,

problema che approssima quello del trasferimento Terra-Luna. Il lavoro di Sweetser faceva parte

della progettazione di una missione lunare della NASA, che consisteva nel trasferire una sonda,

chiamata Lunar Observer, da un’orbita circolare terrestre di 167 km di altitudine ad un’orbita

circolare lunare di 100 km di altitudine.

Consideriamo il moto di un veicolo spaziale diretto verso la Luna. Supponiamo che tale moto si

svolga esclusivamente nel piano dell’orbita lunare e che le uniche forze agenti sul veicolo siano

l’attrazione gravitazionale terrestre e quella lunare. Ipotizziamo inoltre che l’orbita della Luna sia

perfettamente circolare. Grazie a queste ipotesi, il calcolo della traiettoria del veicolo si traduce in

un problema ristretto dei tre corpi.

Utilizziamo un sistema di riferimento cartesiano xy individuante il piano nel quale si svolge il moto

dei tre corpi, la cui origine coincide col baricentro del sistema Terra-Luna, con l’asse x passante per

i centri dei due corpi celesti e orientato verso la Luna. Tale sistema di riferimento viene denominato

sinodico.

Luglio 2010 1

www.matematicamente.it M. Giancola Calcolo della quantità minima…

(x(t),y(t)) = posizione del veicolo nel sistema sinodico all’istante t;

r ( ( ) ( )

)

= & &

v x t , y t = velocità rotazionale = velocità del veicolo nel sistema sinodico;

r

= = +

2 2

& &

v v x y = modulo della velocità rotazionale del veicolo;

ω = velocità angolare con cui il sistema Terra-Luna ruota intorno al proprio baricentro

−

⋅ 6

= 2

,

66 10 rad/s;

ρ = +

2 2 = distanza del veicolo dal baricentro del sistema Terra-Luna;

x y

r = distanza del veicolo dal centro della Terra;

T

r = distanza del veicolo dal centro della Luna;

L

µ 3 2

= costante gravitazionale terrestre = 398600,49 ;

km / s

T

µ 3 2

= costante gravitazionale lunare = 4902,79 .

km / s

L

Se l’astronave non esegue nessuna manovra propulsiva, e quindi le uniche accelerazioni che subisce

sono quelle causate dall’attrazione gravitazionale del sistema Terra-Luna, la quantità

µ µ

2 2

ω ρ

− + + +

2 2 2 T L

v ,

r r

T L

che prende il nome di integrale di Jacobi, rimane costante lungo tutta la sua traiettoria, ovvero

µ µ µ µ

2 2 2 2

( ) ( ) ( ) ( )

ω ρ ω ρ

− + + + = − + + + ∀ ≤ ≤

2 2 2 2 2 2

T L T L

v t t v t t t t t t

( ) ( )

( ) ( )

i i i f

r t r t r t r t

T L T i L i

o, più semplicemente, µ µ

2 2

ω ρ

− + + + =

2 2 2 T L

v C (1)

r r

T L

dove t è l’istante in cui l’astronave comincia a percorrere la sua traiettoria, è l’istante in cui

t

i f

termina il percorso e C è la costante di Jacobi. Ponendo:

µ µ

ω ρ

2 2

= + +

T L

J 2 r r

T L

la (1) diventa: − =

2 (2)

2 J v C

che, per v = 0, si riduce a =

J C 2

Quest’ultima è l’equazione della curva di Hill, o curva a velocità zero, corrispondente al valore C

della costante di Jacobi.

Le condizioni iniziali del moto del veicolo (posizione e velocità iniziali) determinano il valore di C

tramite la (1). Di conseguenza, rimane anche individuata la curva di Hill corrispondente e la zona

>

2

dello spazio nella quale, essendo , può avere luogo il moto del veicolo e dalla quale esso non

v 0 ∆C

può uscire. Ogni manovra impulsiva causerà una variazione della costante di Jacobi, facendole