Obtention pragmatique de développements

Auteurs: Alain Vienne, S. Renner

Séries du problème des 2-corps

Auteur: Alain Vienne

Souvent, lorsque l'on considère le problème des 2 corps et quand l'excentricité e de l'orbite est petite, on a besoin de l'approximation suivante:

V=M+2e \sin M + \frac{5}{4}e^2 \sin 2M  +O (e^3)

\frac{r}{a}=(1+\frac{e^2}{2})-e\cos M - \frac{e^2}{2} \cos 2M  +O (e^3)

r et V sont les coordonnées polaires du corps, V étant compté à partir du péricentre (anomalie vraie). M est le temps ou plus précisément c'est l'anomalie moyenne M=\frac{2\pi}{T} (t-t_0) avec T la période, t le temps et t_0 l'instant de passage au péricentre. a et e sont respectivement le demi-grand axe et l'excentricité de l'orbite.

Dans l'exercice qui est proposé ci-après, on utilisera en plus des anomalies moyennes et vraies, l'anomalie excentrique. L'animation qui est donnée ici visualise leur évolution dans le cas d'une excentricité de 0,7.

Evolution des 3 anomalies (respectivement vraie, excentrique et moyenne)
figures/kepler.gif
Crédit : Astrophysique sur Mesure / Bessou

Si l'existence des développements ci-dessus est admise, l'astronome s'autrorise alors une démarche pragmatique pour les obtenir. Sa démarche n'a pas la rigueur du mathématicien. Dans le cas présenté ici, elle n'a pas non plus une grande efficacité si on souhaite "pousser" le développement plus loin en ordre.Elle a le seul avantage de pouvoir "se tirer d'affaire" dans le cas qui le préoccupe.

Dans l'exercice qui est proposé, on utilisera au besoin le développement de Taylor, l'intégration d'un développement, la substitution de développements.

complementExercices reliés

Pour la recherche de la solution du problème des 2-corps, on peut voir cet exercice. Il y aussi celui qui utilse le théorème de Lagrange. D'autres exercices sur le problème de 2 corps existent sur ce site. On en trouvera, entre autres, sur l'équation de Kepler et son inversion, sur les solutions géométriques du problème de 2 corps, sur le problème de 2 corps perturbé et sur l'excentricité limite dans les développements du problème de 2 corps.


Ex : Séries du problème des 2-corps

Auteur: Alain Vienne

exercice Séries du problème des 2-corps

Difficulté : ☆☆   Temps : 1h

introductionIntroduction

Quand on intègre le problème des 2-corps, la loi des aires permet d'écrire:

M=(1-e^2)^{3/2} \int_0^V \frac{dv}{(1+e \cos v)^2}

Question 1)

On utilisant cette relation et en négligeant les termes d'ordre supérieur ou égal à 3 en excentricité, montrer que l'on a:

M=V-2e \sin V +\frac{3}{4}e^2 \sin 2V + O(e^3)

Question 2)

On sait que \frac{r}{a} = 1 - e \cos E et M=E-e\sin E (E est l'anomalie excentrique). En déduire le développement de \frac{r}{a} en puissance de e et en fonction de M (limité à l'ordre 2)

Question 3)

A partir du développement obtenu à la première question, déduire celui qui donne V en fonction de M (limité à l'ordre 2 en e).

remarqueRemarque

Dans cette dernière question, le calcul est fait "en crabe", il faut donc veiller à la discussion sur l'ordre en excentricité. Plus généralement, tous ces calculs supposent l'existence des développements recherchés. Cette supposition et l'unicité ont permis d'éviter de se soucier des conditions d'application des théorèmes utilisés.


Séries du problème des 2-corps (Lagrange)

Auteur: S. Renner

Date de création: 10 avril 2013

Dans le problème des 2 corps, lorsque l'excentricité e de l'orbite est petite, on peut écrire :

V=M+2e \sin M + \frac{5}{4}e^2 \sin 2M  +O (e^3),

V est l'angle entre la direction du péricentre et la position du corps sur son orbite (anomalie vraie), M est le temps ou plus précisément l'anomalie moyenne M=\frac{2\pi}{T} (t-t_0), avec T la période, t le temps et t_0 l'instant de passage au péricentre.

L'animation donnée ci-après montre l'évolution des anomalies vraie et moyenne (et excentrique) dans le cas d'une excentricité e=0.7.

Evolution des 3 anomalies (respectivement vraie, excentrique et moyenne)
figures/kepler.gif
Crédit : Astrophysique sur Mesure / Bessou

Dans l'exercice proposé, on établit le développement ci-dessus à l'aide du théorème d'inversion de Lagrange.


Ex : Séries du problème des 2-corps (Lagrange)

Auteur: S. Renner

exercice Séries du problème des 2-corps (Lagrange)

Difficulté : ☆☆   Temps : 1h

introductionIntroduction

Quand on intègre le problème des 2-corps, la loi des aires permet d'écrire:

M=(1-e^2)^{3/2} \int_0^V \frac{dv}{(1+e \cos v)^2}

Question 1)

On utilisant cette relation et en négligeant les termes d'ordre supérieur ou égal à 3 en excentricité, montrer que l'on a:

M=V-2e \sin V +\frac{3}{4}e^2 \sin 2V + O(e^3)

Question 2)

Utiliser le théorème d'inversion de Lagrange pour montrer que V=M+2e \sin M + \frac{5}{4}e^2 \sin 2M  +O (e^3).


Description épicyclique du mouvement keplerien

Auteur: S. Renner

Date de création: 3 février 2010

Souvent en dynamique du système solaire, l'excentricité est très faible, et il est donc utile de considérer des approximations au premier ordre en excentricité, en particulier pour des systèmes vus dans des repères tournants. Cette approche est par exemple intéressante pour décrire la dynamique des anneaux planétaires, ou les effets de l'aplatissement d'une planète sur les orbites de satellites.

On propose ici un exercice qui porte sur la description du mouvement keplerien par des épicycles : le mouvement elliptique d'une particule P autour d'un foyer F est vu dans un repère centré sur un point G (le centre guide) en orbite circulaire uniforme autour de F (de rayon a égal au demi-grand axe de la particule, et de vitesse angulaire égale au moyen mouvement n=2 \pi/T, où T est la période orbitale de P).

L'exercice est largement inspiré du livre Solar System Dynamics (C.D. Murray & S.F. Dermott, 1999).


Ex: Description épicyclique du mouvement keplerien

Auteur: S. Renner

exerciceDescription épicyclique du mouvement keplerien

Difficulté : ☆☆   Temps : 2h

introductionIntroduction

Dans le problème à deux corps, on considère le mouvement d'une particule P sur une ellipse de foyer F et de demi-grand axe a. On note n le moyen mouvement de la particule, et f son anomalie vraie.

Soit G un point fictif tournant autour du foyer F sur une orbite circulaire de rayon a égal au demi-grand axe de P, avec une vitesse angulaire égale au moyen mouvement n de la particule. On note M l'angle entre la ligne FG et la direction du péricentre de la particule (c'est donc l'anomalie moyenne de P et on a n = dM/dt).

remarqueRemarque

Les questions 5 à 8 correspondent exactement à la représentation de Ptolémée du mouvement du Soleil autour de la Terre : un mouvement circulaire de rayon R=a avec la Terre au foyer F, et uniforme par rapport à F' (appelé l'équant). Le modèle de Ptolémée était donc du premier ordre en excentricité, et le triomphe de Kepler fut d'élaborer une théorie à l'ordre 2.

Question 1)

Ecrire les coordonnées x, y de P dans le repère orthonormé ($G$,{\bf \hat{x}},{\bf \hat{y}}) tel que {\bf \hat{x}}={\bf FG}/a.

Question 2)

A partir de la loi des aires, montrer que: \displaystyle dM= \frac{(1-e^2)^{3/2}}{(1+e \cos f)^2} df.

Question 3)

Ainsi en intégrant on obtient: \displaystyle M= (1-e^2)^{3/2} \int_0^{f} \frac{du}{(1+e \cos u)^2} . Montrer à l'aide de cette relation que \displaystyle f-M \simeq 2 e \sin M + O(e^2).

Question 4)

En déduire que par rapport à G, la particule P suit une orbite elliptique de demi-grand axe 2ae et de demi-petit axe ae. Dans quel sens s'effectue ce mouvement et avec quelle période?

Question 5)

Soit R la distance entre P et le centre O de l'ellipse. Montrer que \displaystyle R \simeq a (1 - \frac{1}{2}  e^2 \sin ^2 M).

Question 6)

Montrer qu'au premier ordre en excentricité, la trajectoire de P est alors un cercle centré sur O, et que l'angle \widehat{POF} est confondu avec l'anomalie excentrique E.

Question 7)

On note g l'angle \widehat{FF'P}, où F' est le second foyer (vide) de l'ellipse. Ecrire \cos g et en déduire que cos g \simeq \cos M + O(e^2). On rappelle que FP+F'P=2a, et on utilisera le développement: r/a \simeq 1 - e \cos M + (e^2/2)(1 \cos 2 M) (voir cet exercice).

Question 8)

Ainsi, le mouvement de la particule P, vu du foyer F', est uniforme avec une vitesse angulaire égale au moyen mouvement n. Que peut-on dire des droites F'P et FG?


Sphère d'influence

Auteur: S. Renner et A. Vienne

Date de création: 3 mars 2010

Nous allons définir la notion de sphère d'influence en s'intéressant au problème à trois corps Soleil + Jupiter + satellite, de masses respectives M>>m>>m' (qui constitue ce que l'on appelle un problème de Kepler perturbé).


Ex: Sphère d'influence

Auteurs: S. Renner, A. Vienne

exerciceSphère d'influence

Difficulté : ☆☆   Temps : 3h

introductionIntroduction

On considère le système gravitationnel formé du Soleil $S$ de masse M_\odot \equiv 1, de Jupiter J de masse $m$ petite devant 1 (m=1/1047) et d'un troisième corps A de masse m_{A} petite devant $m$ (par exemple un astéroïde ou une sonde spatiale). On note K la constante de gravitation universelle, \overrightarrow{r}=\overrightarrow{SJ}, \overrightarrow{p}=\overrightarrow{SA}, et \overrightarrow{\Delta}=\overrightarrow{p}-\overrightarrow{r}.

Question 1)

Ecrire les équations du mouvement héliocentrique de J (c'est-à-dire de \overrightarrow{r}).

Question 2)

Quelle est la nature du mouvement de J si la quantité m_{A} est négligeable?

Question 3)

On suppose désormais que le mouvement de J est circulaire (c'est-à-dire que m_A est négligeable et que r est constant). Montrer que l'équation du mouvement héliocentrique de A (c'est-à-dire de \overrightarrow{p}) est \frac{d^{2}\overrightarrow{p}}{dt^{2}}=\overrightarrow{A_{0}} + \overrightarrow{A_{1}}  \textrm{, avec }  \overrightarrow{A_{0}}=-K\frac{\overrightarrow{p}}{p^{3}} \textrm{ et }  \overrightarrow{A_{1}}=-Km \Big{(} \frac{\overrightarrow{\Delta}}{\Delta^{3}}+ \frac{\overrightarrow{r}}{r^3} \Big{)}.

Question 4)

Montrer que l'équation du mouvement jovicentique de A (c'est-à-dire de \overrightarrow{\Delta}) est \frac{d^{2}\overrightarrow{\Delta}}{dt^{2}}= \overrightarrow{B_{0}}+ \overrightarrow{B_{1}} \textrm{, avec } \overrightarrow{B_{0}}=-Km\frac{\overrightarrow{\Delta}} {\Delta^{3}} \textrm{ et } \overrightarrow{B_{1}}=-K \Big{(} \frac{\overrightarrow{p}}{p^{3}}-\frac{\overrightarrow{r}}{r^3} \Big{)}.

Question 5)

Dans la suite, on cherche à étudier la surface $(S)$ définie par \frac{\left\Vert \overrightarrow{A_{1}}\right\Vert } {\left\Vert \overrightarrow{A_{0}}\right\Vert } =\frac{\left\Vert \overrightarrow{B_{1}}\right\Vert } {\left\Vert \overrightarrow{B_{0}}\right\Vert }. On va voir notamment qu'elle a presque la forme d'une sphère, appelée sphère d'influence de Jupiter. Justifier cette appellation en expliquant ce qui se passe pour le mouvement de A lorsque celui-ci est respectivement très à l'intérieur ou très à l'extérieur de cette sphère.

Question 6)

On pose \displaystyle u=\frac{\Delta}{r} et \varphi=\textrm{angle}(\overrightarrow{JS},\overrightarrow{JA}). Montrer que l'on a A_{1}=\frac{Km}{\Delta^{2}}\sqrt{1-2u^{2}\cos\varphi+u^{4}}.

Question 7)

Montrer que l'on a B_{1}=\frac{K}{r^{2}}\sqrt{1+\Big{(}\frac{r}{p}\Big{)}^{4}-2\Big{(} \frac{r}{p}\Big{)}^{3} \frac{\overrightarrow{p}.\overrightarrow{r}}{r^{2}}} .

Question 8)

En exprimant \displaystyle \Big{(} \frac{p}{r} \Big{)}^{2} et \displaystyle  \frac{\overrightarrow{p}.\overrightarrow{r}}{r^{2}} en fonction de \varphi et u, montrer que B_{1}=\frac{K}{r^{2}}\frac{\sqrt{1+(1-2u\cos\varphi+u^{2})^{2} -2(1-u\cos\varphi)\sqrt{1-2u\cos\varphi+u^{2}}}}{1-2u\cos\varphi+u^{2}}.

Question 9)

En déduire que l'équation de la surface (S) peut s'écrire m^{2}=f(u,\varphi), et donner l'expression de f(u,\varphi).

Question 10)

Puisque m est petit devant 1, on peut considérer que, sur la surface (S), u est également petit. Montrer que le développement de f(u,\varphi) suivant les puissances de u, limité à son terme de plus bas degré vaut \displaystyle u^{5}\sqrt{1+3\cos^2\varphi}. Justifier alors le mot sphère dans l'expression sphère d'influence qui désigne la surface (S).

Question 11)

Calculer le rayon de la sphère d'influence de Jupiter sachant que $r=5.2$ UA.


Réponses aux exercices

pages_develop-appli/exo-series-2corps.html

Exercice ' Séries du problème des 2-corps'


pages_develop-appli/exo-series-2corps-lagrange.html

Exercice ' Séries du problème des 2-corps (Lagrange)'


pages_develop-appli/exo-epicycle-kepler.html

Exercice 'Description épicyclique du mouvement keplerien'


pages_develop-appli/exo-sphere-influence.html

Exercice 'Sphère d'influence'