Intégrale de Jacobi

En mécanique céleste, l'intégrale de Jacobi ou constante de Jacobi (d'après Charles Gustave Jacob Jacobi) est la seule quantité conservée dans le problème à trois corps restreint circulaire ; contrairement au cas du problème à deux corps, l'énergie et le moment cinétique du système ne sont pas conservés séparément et une solution analytique générale n'est pas possible. Cette intégrale a été utilisée pour fournir des solutions à des cas particuliers du problème à trois corps.

Définition

Système de coordonnées synodique

Système co-rotationnel

Le système de coordonnées utilisé est appelé synodique ou système co-rotationnel, placé au barycentre, la ligne connectant les deux masses μ₁, μ₂ étant choisie comme l'axe x et la longueur unité égale à leur distance. Comme le système tourne en même temps que les deux masses, celles-ci restent stationnaires et positionnées à (−μ₂, 0) et à (+μ₁, 0)¹.

Dans le système de coordonnées (x,  y), l'intégrale de Jacobi s'exprime comme :

$C_J=n^2(x^2+y^2)+2(\frac{{\mu }_1}{r_1}+\frac{{\mu }_2}{r_2})-({\dot{x}}^2+{\dot{y}}^2+{\dot{z}}^2)$

où :

• $n=\frac{2\pi }{T}$ est le mouvement moyen (T est la période orbitale)
• ${\mu }_1=G{m}_1,{\mu }_2=G{m}_2$, pour les masses m₁, m₂ et G est la constante gravitationnelle
• $r_1,r_2$ sont les distances de la particule test aux deux masses

On notera que l'intégrale de Jacobi vaut moins deux fois l'énergie totale par unité de masse dans le référentiel de référence tournant : le premier terme représente l'énergie potentielle centrifuge, le second représente l'énergie potentielle gravitationnelle et le troisième est l'énergie cinétique. Dans ce système de référence, les forces qui agissent sur la particule sont les deux attractions gravitationnelles, la force centrifuge et la force de Coriolis. Puisque les trois premières dérivent de potentiels et que la dernière est perpendiculaire à la trajectoire, elles sont toutes conservatives, et donc l'énergie calculée dans ce système de référence (et donc l'intégrale de Jacobi) est une constante du mouvement. Une preuve calculatoire directe est donnée ci-dessous.

Système de coordonnées sidéral

Système de coordonnées inertiel.

Dans le système de coordonnées sidéral inertiel (ξ,  η,  ζ), les masses orbitent autour du barycentre. Dans ce système de coordonnées, l'intégrale de Jacobi s'exprime par :

$C_J=2(\frac{{\mu }_1}{r_1}+\frac{{\mu }_2}{r_2})+2n(\xi \dot{\eta }-\eta \dot{\xi })-({\dot{\xi }}^2+{\dot{\eta }}^2+{\dot{\zeta }}^2).$

Démonstration

Dans le système co-rotationnel, les accélérations peuvent être exprimées comme les dérivées d'une simple fonction scalaire

$U(x,y,z)=\frac{n^2}{2}(x^2+y^2)+\frac{{\mu }_1}{r_1}+\frac{{\mu }_2}{r_2}$

En utilisant une représentation lagrangienne des équations du mouvement :

[Eq.1] $\ddot{x}-2n\dot{y}=\frac{\delta U}{\delta x}$

[Eq.2] $\ddot{y}+2n\dot{x}=\frac{\delta U}{\delta y}$

[Eq.3] $\ddot{z}=\frac{\delta U}{\delta z}$

En multipliant respectivement les [Eq.1], [Eq.2] et [Eq.3] par $\dot{x},\dot{y}$ et $\dot{z}$ et en les additionnant, on obtient

$\dot{x}\ddot{x}+\dot{y}\ddot{y}+\dot{z}\ddot{z}=\frac{\delta U}{\delta x}\dot{x}+\frac{\delta U}{\delta y}\dot{y}+\frac{\delta U}{\delta z}\dot{z}=\frac{dU}{dt}$

L'intégration fournit

${\dot{x}}^2+{\dot{y}}^2+{\dot{z}}^2=2U-C_J$

où C_J est la constante d'intégration.

Le terme de gauche représente le carré de la vitesse v de la particule test dans le système co-rotationnel.

¹Ce système de coordonnées est non inertiel, ce qui explique l'apparition de termes relatifs aux accélérations centrifuge et de Coriolis.

Voir aussi

• Référentiel tournant
• Critère de Tisserand

Références

• Carl D. Murray and Stanley F. Dermot Solar System Dynamics [Cambridge, England: Cambridge University Press, 1999], pages 68–71. Modèle:ISBN