Équation de Tsiolkovski

L'équation de Tsiolkovski est l'équation fondamentale de l'astronautique, reliant l'accroissement de vitesse au cours d'une phase de propulsion d'un astronef doté d'un moteur à réaction au rapport de sa masse initiale à sa masse finale.

On la doit à Constantin Tsiolkovski et, indépendamment, à Hermann Oberth.

Histoire

L'équation de Tsiolkovski est considérée comme l'équation fondamentale de l'astronautique. Son éponyme est Constantin Tsiolkovski (1857-1935), qui l'a déduite puis publiée en 1903.

Bien que cette équation soit souvent attribuée au grand penseur autodidacte Tsiolkovski, une forme de l'équation figure déjà dans un traité du mathématicien britannique William Moore (en) (fl. c.1806-1823) paru en 1813 puis dans un article du général-major belge Casimir-Érasme Coquilhat (1811-1890) paru en 1873 .

L'expérience de la barque, de Tsiolkovski

Expérience de la barque, de Tsiolkovski.

Afin de faire comprendre le principe de la propulsion à réaction, Constantin Tsiolkovski a proposé une fameuse expérience « de la barque ». Une personne se trouve démunie d'avirons dans une barque à l'écart de la rive. Elle veut rejoindre cette rive. Elle remarque que la barque est chargée d'une certaine quantité de pierres et a l'idée de lancer, une à une et le plus vivement possible, ces pierres dans la direction opposée à la rive. Effectivement, à la quantité de mouvement des pierres jetées dans un sens correspond une quantité de mouvement égale pour la barque dans l'autre sens.

Énoncé

L'équation de Tsiolkovski s'écrit :

$\Delta {\vec {v}}={\vec {v}}_{\mathrm {f} }-{\vec {v}}_{\mathrm {i} }=-{\vec {v}}_{\mathrm {e} }\,\ln {\frac {m_{\mathrm {i} }}{m_{\mathrm {f} }}}$, vectoriellement, ou $\Delta v=v_{\mathrm {e} }\,\ln {\frac {m_{\mathrm {i} }}{m_{\mathrm {f} }}}$, numériquement

où :

• $\Delta v$ est la variation de vitesse de l'astronef entre la fin et le début de la phase propulsée considérée ;
• $v_{\mathrm {e} }$ est la vitesse d'éjection des gaz ;
• $m_{\mathrm {i} }$ est la masse totale de l'astronef au début de la phase propulsée (indice i pour initial) ;
• $m_{\mathrm {f} }$ est la masse totale de l'astronef à l'issue de la phase propulsée (indice f pour final), exprimée dans la même unité que $m_{\mathrm {i} }$ ;
• $\ln$ est la fonction logarithme népérien.

Établissement

Démonstration

Cette équation est établie en intégrant l'équation de conservation de la quantité de mouvement entre le début et la fin de la phase propulsée sous les hypothèses suivantes :

• l'étude du mouvement est faite dans un référentiel d'inertie ;
• l'astronef n'est soumis qu'à la force de poussée fournie par ses moteurs, aucune autre action extérieure (gravité, efforts aérodynamiques) n'est prise en compte (voir à la fin de l'article pour la prise en compte de la gravité) ;
• la vitesse d'éjection des gaz est constante.

À un instant donné, lorsque le vaisseau de masse $m$ se déplaçant à la vitesse ${\vec {v}}$ éjecte une petite quantité d'ergol à la vitesse ${\vec {v_{\mathrm {e} }}}$, on note $\mathrm {d} m$ sa variation de masse et $\mathrm {d} {\vec {v}}$ sa variation de vitesse. La variation de quantité de mouvement du système isolé (vaisseau + ergol éjecté) est nécessairement nulle, il vient donc :

$m\cdot \mathrm {d} {\vec {v}}+(-\mathrm {d} m)\cdot {\vec {v_{\mathrm {e} }}}={\vec {0}}$.

Pour obtenir la variation de vitesse $\Delta {\vec {v}}$ de l'astronef quand sa masse passe de $m_{\mathrm {i} }$ à $m_{\mathrm {f} }$, on peut intégrer cette petite variation de vitesse :

$\Delta {\vec {v}}=\int _{{\vec {v}}={\vec {v_{\mathrm {i} }}}}^{\vec {v_{\mathrm {f} }}}\mathrm {d} {\vec {v}}=\int _{m=m_{\mathrm {i} }}^{m_{\mathrm {f} }}{\frac {\vec {v_{\mathrm {e} }}}{m}}\mathrm {d} m={\vec {v_{\mathrm {e} }}}\int _{m_{\mathrm {i} }}^{m_{\mathrm {f} }}{\frac {\mathrm {d} m}{m}}={\vec {v_{\mathrm {e} }}}\ln {\frac {m_{\mathrm {f} }}{m_{\mathrm {i} }}}=-{\vec {v_{\mathrm {e} }}}\ln {\frac {m_{\mathrm {i} }}{m_{\mathrm {f} }}}$.

Puisque $m_{\mathrm {i} }>m_{\mathrm {f} }$, $\ln(m_{\mathrm {i} }/m_{\mathrm {f} })>0$, la variation de vitesse du vaisseau a donc comme prévu la même direction que la projection des ergols et en sens opposé.

Commentaires

Il est souvent dit que pour trouver cette équation, il faut que le débit massique d'ergol soit constant pendant la phase de propulsion, mais ce n'est pas obligatoire, même si cela simplifie le travail d'intégration dans un premier temps.^{[réf. nécessaire]}

L'équation est valable aussi bien lors d'une phase d'accélération (la poussée est dans la direction de la vitesse, $\Delta v$ est positif : c'est un accroissement de vitesse) que de décélération (la poussée est de direction opposée à la vitesse, $\Delta v$ est négatif : c'est une réduction de vitesse).

La différence entre la masse initiale $m_{\mathrm {i} }$ et la masse finale $m_{\mathrm {f} }$ correspond à la masse que la fusée a éjectée durant sa propulsion ; on l'appelle masse d'appui^{[réf. nécessaire]} (« d'appui » parce que c'est la masse sur laquelle la fusée s'est appuyée pour se propulser).

Pour les fusées thermo-chimiques (Ariane, Soyouz, navette, etc.), la masse d'appui est la masse des ergols (poudre, ou dioxygène et dihydrogène), laquelle est également source d'énergie chimique : c'est donc la masse d'appui elle-même qui contient l'énergie qui servira à sa propre éjection. Ce n'est plus le cas pour les moteurs ioniques (qui représentent sans doute l'avenir de la conquête spatiale). Ceux-ci sont régis tout pareillement par l'équation de Tsiolkovski, mais leur masse d'appui est constituée d'un gaz neutre (du xénon) ; c'est la très forte vitesse d'éjection de cette masse d'appui qui rend ces moteurs très économes en masse d'appui (il leur faut cependant une source d'énergie pour réaliser l'éjection). À ce titre, le fonctionnement des moteurs ioniques est comparable à celui des fusées à eau, dans lesquelles l'eau n'est utilisée que pour sa masse (l'énergie résidant dans l'air comprimé).

Dans le cas où la phase propulsée est réalisée au moyen de plusieurs étages fonctionnant successivement, la même équation de Tsiolkovski peut être utilisée pour le vol de chaque étage. On peut ainsi montrer l'intérêt de telles fusées à plusieurs étages. Voir l'exemple dans la section suivante.

Malgré l'apparente simplicité de cette équation et des hypothèses qui la sous-tendent, elle constitue une approximation utile au calcul des manœuvres de changement d'orbite, ces manœuvres étant qualifiées d'impulsionnelles, c’est-à-dire effectuées en un temps suffisamment bref pour que les hypothèses de l'équation de Tsiolkovski restent approximativement valables.

Temps nécessaire

Si le vaisseau utilise un débit massique d'ergol constant $q$ on peut écrire :

$\Delta t={\frac {\Delta m}{q}}$.

Or l'équation de Tsiolkovski peut s'écrire :

${\frac {\Delta v}{v_{\mathrm {e} }}}=\ln {\frac {m_{\mathrm {i} }}{m_{\mathrm {i} }-\Delta m}}$.

C'est-à-dire, en changeant de signe et en passant à l'exponentielle :

${\frac {m_{\mathrm {i} }-\Delta m}{m_{\mathrm {i} }}}=\mathrm {e} ^{-\Delta v/v_{\mathrm {e} }}$.

On en tire l'expression de $\Delta m$ qu'on reporte dans celle de $\Delta t$ :

$\Delta t={\frac {m_{\mathrm {i} }}{q}}\,\left(1-\mathrm {e} ^{-\Delta v/v_{\mathrm {e} }}\right)$.

Exemple

L'exemple qui suit a pour objet de montrer l'intérêt des fusées à plusieurs étages.

Soit une fusée à deux étages ayant les caractéristiques suivantes :

• la masse d'ergols embarqués par chaque étage (premier étage : 100 tonnes ; deuxième étage : 20 tonnes) représente 10 fois sa masse à vide ;
• la vitesse d'éjection des gaz est de 4 000 m/s ;

et supposons qu'elle emporte une charge utile de 2 t. Résumons ces données dans un tableau :

Étage	Masse d'ergols (t)	Masse à vide (t)	Masse totale (t)	Vitesse d'éjection des gaz (m/s)
Premier étage	$m_{\mathrm {e} {\mathfrak {1}}}=100$	$m_{\mathrm {v} {\mathfrak {1}}}=10$	$m_{\mathrm {t} {\mathfrak {1}}}=110$	$v_{\mathrm {e} }=4\,000$
Deuxième étage	$m_{\mathrm {e} {\mathfrak {2}}}=20$	$m_{\mathrm {v} {\mathfrak {2}}}=2$	$m_{\mathrm {t} {\mathfrak {2}}}=22$	$v_{\mathrm {e} }=4\,000$
Charge utile			$m_{\mathrm {cu} }=2$
Total fusée	$m_{\mathrm {e} {\mathfrak {}}}=120$	$m_{\mathrm {v} {\mathfrak {}}}=12$	$m_{\mathrm {t} {\mathfrak {}}}=134$

On peut alors mener les calculs d'incréments de vitesse, comme suit, en employant deux fois l'équation de Tsiolkovski, aux étapes 3 et 6 :

Étape de calcul		Formule	Masse (t)	Vitesse (m/s)
1	Masse à l'allumage du premier étage	$m_{\mathrm {i} {\mathfrak {1}}}=m_{\mathrm {t} {\mathfrak {}}}$	$134$
2	Masse à l'extinction du premier étage	$m_{\mathrm {f} {\mathfrak {1}}}=m_{\mathrm {i} {\mathfrak {1}}}-m_{\mathrm {e} {\mathfrak {1}}}$	$34$
3	Incrément de vitesse du premier étage	$\Delta v_{1}=v_{\mathrm {e} }\ln {\tfrac {m_{\mathrm {i} {\mathfrak {1}}}}{m_{\mathrm {f} {\mathfrak {1}}}}}$		$5\,486$
4	Masse à l'allumage du second étage	$m_{\mathrm {i} {\mathfrak {2}}}=m_{\mathrm {f} {\mathfrak {1}}}-m_{v{\mathfrak {1}}}$	$24$
5	Masse à l'extinction du second étage	$m_{\mathrm {f} {\mathfrak {2}}}=m_{\mathrm {i} {\mathfrak {2}}}-m_{e{\mathfrak {2}}}$	$4$
6	Incrément de vitesse du deuxième étage	$\Delta v_{2}=v_{\mathrm {e} }\ln {\tfrac {m_{\mathrm {i} {\mathfrak {2}}}}{m_{\mathrm {f} {\mathfrak {2}}}}}$		$7\,167$
7	Vitesse finale	$\!\,\Delta v=\Delta v_{1}+\Delta v_{2}$		$12\,653$

Par comparaison, une fusée comportant un seul étage avec la même quantité totale d'ergols (120 t) et la même masse à vide totale (12 t) imprimerait à une charge utile de même masse (2 t) une vitesse environ 30 % inférieure :

Étape de calcul		Formule	Masse (t)	Vitesse (m/s)
1	Masse à l'allumage de l'étage (unique)	$m_{\mathrm {i} {\mathfrak {}}}=m_{\mathrm {t} {\mathfrak {}}}$	$134$
2	Masse à l'extinction de l'étage	$m_{\mathrm {f} {\mathfrak {}}}=m_{\mathrm {i} {\mathfrak {}}}-m_{\mathrm {e} {\mathfrak {}}}=m_{\mathrm {v} {\mathfrak {}}}+m_{cu{\mathfrak {}}}$	$14$
3	Vitesse finale	$\Delta v=v_{\mathrm {e} }\ln {\tfrac {m_{\mathrm {i} }}{m_{\mathrm {f} }}}$		$9\,034$

Pertes par pesanteur

Les calculs ont été effectués dans l'hypothèse d'une absence de pesanteur (manœuvres en orbite). Lorsque cette gravité agit, un terme simple doit être ajouté à l'équation de Tsiolkovski. Celle-ci devient :

$\Delta {\vec {v}}=-{\vec {v}}_{\mathrm {e} }\,\ln {\frac {m_{\mathrm {i} }}{m_{\mathrm {f} }}}+\surd {\bigl (}2{\vec {g}}\,\Delta h{\bigr )}$, vectoriellement, ou $\Delta v=\ v_{\mathrm {e} }\,\ln {\frac {m_{\mathrm {i} }}{m_{\mathrm {f} }}}-\surd {\bigl (}2g\Delta h{\bigr )}$, si l'on projette l'équation radialement,

${\vec {g}}$ étant l'accélération locale de la pesanteur et $\Delta h$ la différence entre l'altitude de satellisation et celle du point de départ de la fusée. On fait l'hypothèse que $g$ est constant pendant la propulsion (alors que $g$ diminue très légèrement avec l'altitude).

Voir aussi

Articles connexes

• Moteur à réaction
• Propulsion électrique et ionique

Publications originales

• [Moore 1813] (en) Williams Moore, A treatise on the motion of rockets : to which is added an essay on navel gunnery, in theory and practice ; designed for the use of the army and navy, and in all places of military, naval and scientific instruction, Londres, 1813.
• [Coquilhat 1873] Casimir-Érasme Coquilhat, « Trajectoires des fusées volantes dans le vide », Mémoires de la Société royale des sciences de Liège, 2^e série, t. V,‎ nov. 1873, art. n^o 5, [2]−33 p. (OCLC 248870034, lire en ligne), réimpr. :
  • [Coquilhat 2018] Casimir-Érasme Coquilhat, Trajectoires des fusées volantes dans le vide : 1873, Dijon, Nielrow, janv. 2018, 1^re éd., 1 vol., 51 (ISBN 978-2-9559619-5-7, EAN 9782955961957, OCLC 1043726396, BNF 45456511).
• [Tsiolkovski 1903] (ru) Konstantin Èduardovič Ciolkovskij, « Исследовние мировых пространств реактивными приборами » [« L'exploration de l'espace cosmique par des engins à réaction »], Научное обозрение, vol. 7, n^o 5,‎ 1903, p. 45-75 (résumé) — publication originale de l'article de Tsiolkovski énonçant l'équation fondamentale de l'astronautique.

Bibliographie

 : document utilisé comme source pour la rédaction de cet article.

• [Bouchez 2010] Marc Bouchez, « Propulsion aérospatiale : lois simplifiées pour le dimensionnement et exercices d'application », Techniques de l'ingénieur,‎ 10 juillet 2010, art. BM 3 003, § 1.2 (« Formule de Tsiolkovski : cas d'une accélération constante en palier »), p. 3, col. 2 (résumé, lire en ligne).
• [Feodosiev et Siniarev 1959] (en) V. I. Feodosiev et G. B. Siniarev (trad. du russe par S. N. Samburof), Introduction to rocket technology [« Introduction à la technologie des fusées »], New York et Londres, Academic Press, janvier 1959 (réimpr. mai 2014), 1^re éd., X-344 p., in-8^o (15,2 × 23,5 cm) (ISBN 978-1-4832-3201-0, EAN 9781483232010, OCLC 840470516, BNF 44667202, DOI 10.1016/C2013-0-12448-8, présentation en ligne, lire en ligne), chap. I^er (OCLC 7333041790, DOI 10.1016/B978-1-4832-3201-0.50006-1, résumé), section B (« Tsiolkovski equation for the ideal velocity of a rocket ») [« Équation de Tsiolkovski pour la vitesse idéale d'une fusée »], p. 16-22.
• [Macdonald et Badescu 2014] (en) Malcolm Macdonald et Viorel Badescu (éditeur), The international handbook of space technology [« Le guide / manuel international de la technologie spatiale »], Heidelberg et Chichester, Springer et Praxis, coll. « Praxis books / Astronautical engineering », février 2014, 1^re éd., XV-731 p., 21,6 × 27,9 cm (ISBN 978-3-642-41100-7 et 978-3-662-50608-0, OCLC 958469721, BNF 44707981, DOI 10.1007/978-3-642-41101-4, SUDOC 199104573, lire en ligne) :
  • [Macdonald, Norris et Spencer 2014] (en) Malcolm Macdonald, Pat Norris et David B. Spencer, « Introduction », dans op. cit., chap. 1^er [« Introduction (générale) »], p. 1-24 (OCLC 7327569974, DOI 10.1007/978-3-642-41101-4_1). 
  • [Bonnal et al. 2014] (en) Christophe Bonnal, Alessandro Ciucci, Michael H. Obersteiner et Oskar Haidn, « Launch systems », dans op. cit., chap. 8, p. 165-195 [« Systèmes de lancement »] (OCLC 7326880462, DOI 10.1007/978-3-642-41101-4_8. 
• [Rax 2005] Jean-Marcel Rax (préf. de Bernard Bigot), Physique des plasmas : cours et applications, Paris, Dunod, coll. « Sciences sup », août 2005, 1^re éd., XII-426 p., 17 × 24 cm (ISBN 978-2-10-007250-7, EAN 9782100072507, OCLC 77230441, BNF 40025642, SUDOC 090320409, présentation en ligne, lire en ligne).
• [Serra, Jung et Picard 2013] (en) Jean-Jacques Serra, Philippe Jung et Théo Pirard, « Casimir Coquilhat's theory on rocket motion : the rocket equation established in 1871 ! », dans John Harlow (éditeur), History of rocketry and astronautics (actes du 42^e symposium d'histoire de l'Académie internationale d'astronautique, tenu à Glasgow en 2008), San Diego, Univelt, coll. « AAS history series » (n^o 39) et coll. « IAA history symposia » (n^o 28), 2013, 1^re éd., XII-345 p., 25 cm (ISBN 978-0-87703-589-3 et 978-0-87703-590-9, EAN 9780877035893, OCLC 840430825, présentation en ligne), 2^e partie, chap. 7, p. 121-134 (résumé).