Développements limités

Auteurs: Marc Fouchard, Jérôme Thiébaut, S. Renner, Stéphane Erard

Lentille gravitationnelle

Auteur : Jérôme Thiébaut

L'idée de base de la relativité générale est que la matière, par sa masse, courbe l'espace. Ainsi, une planète orbitant autour d'une étoile n'est pas soumise à une force de gravitation mais circule librement sur un espace courbé par l'étoile. Il s'ensuit que la lumière, bien que dépourvue de masse, est également déviée par la présence d'un corps massif. Si un corps massif se situe entre une galaxie lointaine et un observateur, celui ci va donc dévier la lumière de la galaxie et déformer son image. C'est ce qu'on appelle une lentille gravitationnelle. Dans le cas où les trois objets sont parfaitement alignés, l'image de la galaxie se déforme pour former un anneau autour de la lentille appelé anneau d'Einstein. Le but de cet exercice est de déterminer le diamètre angulaire de cet anneau en fonction des caractéristiques du système (masse et distances).

Lentille Gravitationnelle
gravlens.jpg
Lentille gravitationnelle : la déviation de la lumière par un fort potentiel gravitationnel (l'amas de galaxies 0024+1654) conduit à de multiples images d'un objet situé derrière le centre de masse du déflecteur.
Crédit : HST

Ex : Lentille gravitationnelle

Auteur: Jérôme Thiébaut

exerciceLentille gravitationnelle

Difficulté :    Temps : 30 min

La figure montre le schéma du principe d'une lentille gravitationnelle. La source, S, est déviée par la lentille, L, et son image, I, est donc vue par l'observateur, O, selon un angle theta_I au lieu de theta_S.

Principe d'une lentille gravitationnelle
lg.png
Crédit : Jérôme Thiébaut
Question 1)

Exprimer les distances AS et AI en fonction des angles theta_S et theta_I et de la distance OS.

Question 2)

Les calculs relativistes montrent que la distance SI vaut alpha*LSoù l'angle de déviation, alpha, vaut 4GM/c^2*theta_I*OL, avec G la constante de gravitation, M la masse de la lentille et c la vitesse de la lumière. Dans l'approximation des petits angles, déterminer l'équation du second degré à laquelle obéit theta_I en fonction de theta_S, et des grandeurs caractéristiques du système, LS, OS, OL et M.

Question 3)

Dans le cas où la source, la lentille et l'observateur sont parfaitement alignés, déterminer le rayon angulaire (rayon d'Einstein) sous lequel sera vu la source toujours en fonction des grandeurs caractéristiques du système.


Repérage des astres

Auteur : S. Renner

Date de création: 30 novembre 2009

On propose ici un exercice en lien avec le passage des coordonnées locales (ou horizontales) aux coordonnées horaires (angle horaire et déclinaison).

Il est donc préférable de se familiariser avec les systèmes de coordonnées utilisés en astronomie pour le repérage des étoiles dans le ciel, ainsi qu'avec les formules de base de la trigonométrie sphérique.


Ex: Repérage des astres

Auteur: S. Renner

exerciceRepérage des astres

Difficulté : ☆☆   Temps : 1h

Dans le triangle sphérique PZM (voir la fiche de résolution du triangle sphérique), on obtient des relations qui permettent de passer des coordonnées horizontales aux coordonnées horaires, et inversement.

Question 1)

La distance zénithale z dépendant seulement de la variable H (angle horaire), démontrer que {dz \over dH} = \cos \phi \sin a , puis {da \over dH} = {\cos \delta \cos S \over \sin z}.

Question 2)

En déduire qu'au voisinage du passage au méridien d'un astre, sa distance zénithale z varie comme le carré de l'angle horaire (en se limitant au 2ème ordre).


Equation de Kepler

Auteur : Marc Fouchard

Dans le problème de deux corps (voir cet exercice), on sait que le déplacement d'un corps par rapport à l'autre se fait sur une conique dont le deuxième corps occupe l'un des foyers (voir aussi cet . Une fois la conique fixée il ne reste alors qu'à positionner le corps sur son orbite. Pour cela on utilise une quantité qu'on appelle anomalie. On définie trois types d'anomalie: l'anomalie moyenne, l'anomalie vraie et l'anomalie excentrique.

Les trois anomalies
kepler.gif
Crédit : Astrophysique sur Mesure / Bessou Fouchard

L'animation ci-dessus montre le lien entre les 3 anomalies. Comme on peut le voir, l'anomalie moyenne correspond en fait à un temps. Il n'existe pas de relation géométrique entre l'anomalie moyenne et les autres anomalies. En revanche, il existe une relation (voir cet exercice), appelée équation de Kepler, qui relie l'anomalie moyenne à l'anomalie excentrique. Cette relation est:

M=E-e\sin E,

M est l'anomalie moyenne, E l'anomalie excentrique et e l'excentricité de la trajectoire.

On voit que connaissant E il est facile d'avoir M, mais en revanche connaissant M il n'est pas possible d'avoir E sous forme analytique. L'objet de cette exercice est justement de déterminer un algorithme puissant d'inversion de cette équation.


Ex: Equation de Kepler

Auteur: Marc Fouchard

exerciceEquation de Kepler

Difficulté : ☆☆   Temps : 1h30

Question 1)

On souhaite résoudre une équation du type f(x)=0. Soit a une solution à l'équation. On supposera dans la suite que, sur un intervalle I contenant a, f' ne s'annulle jamais et que f est de classe \mathcal{C}^{\infty} sur cet intervalle. Soit x_n, la suite définie par :

x_{n+1}=x_n+\delta_n,

avec \delta_n, la solution de f(x_n)+\delta_n f'(x_n)+\cdots+\frac{\delta^n}{k!}f^{(k)}(x_n)=0 (k étant fixé).

On souhaite montrer que si la suite x_n converge vers a, alors elle converge au moins avec l'ordre k+1.

Calculer f(x_{n+1}) en fonction de \delta_n et un nombre \zeta_1 \in \lbrack x_n,x_{n+1}\rbrack.

Question 2)

Soit \epsilon_{n+1} l'erreur de x_{n+1} définie par \epsilon_{n+1}=a-x_{n+1}. Cacluler \epsilon_{n+1} en fonction de \delta_n, \zeta_1 et un nombre \zeta_2 \in \lbrack x_{n+1},a\rbrack.

Question 3)

Montrer que:

\frac{\epsilon_n}{\delta_n}=\frac{f'(x_n)}{f'(\zeta_3)}+\frac{\delta_n}{2!}\frac{f''(x_n)}{f'(\zeta_3)}+\cdots+\frac{\delta_n^{k-1}}{k!}\frac{f^{(k)}(x_n)}{f'(\zeta_3)}, avec \zeta_3 un nombre entre a et x_n.

Question 4)

En déduire que \lim_{n\to +\infty}\frac{\epsilon_n}{\delta_n}=1.

Question 5)

En déduire que \lim_{n\to+\inty}\frac{\epsilon_{n+1}}{\epsilon_n^{k+1}}=-\frac{1}{(k+1)!}\frac{f^{(k+1)}(a)}{f'(a)}, et que la suite x_n converge vers a au moins à l'ordre k+1.

Question 6)

Il nous faut maintenant calculer la solution \delta_n de l'équation qui le défini. Avant tout, jusitfier pourquoi il suffit d'avoir une solution approxée à \delta_n^{k+1} près.

Question 7)

Montrer que pour k=3, on a:

\delta_n=-\frac{f(x_n)}{f'(x_n)+\frac{1}{2}\delta_n f''(x_n)+\frac{1}{6}\delta_n^2f'''(x_n)}

Question 8)

Montrer qu'en utilisant la séquance suivante:

\delta_{n1}=-\frac{f(x_n)}{f'(x_n)},

\delta_{n2}=-\frac{f(x_n)}{f'(x_n)+\frac{1}{2}\delta_{n1}f''(x_n)},

\delta_{n3}=-\frac{f(x_n)}{f'(x_n)+\frac{1}{2}\delta_{n2}f''(x_n)+\frac{1}{6}\delta_{n2}f'''(x_n)},

on obtient une approximation de \delta_n à \delta_n^4 près. On pourra procéder par étape en montrant que \delta_n-\delta_{n1}=\mathcal{O}(\delta_n^2), puis que \delta_n-\delta_{n2}=\mathcal{O}(\delta_n^3) et enfin que \delta_n-\delta_{n3}=\mathcal{O}(\delta_n^4).

Question 9)

Montrer que l'on peut appliquer l'algorithme précédent pour inverser l'équation de Kepler sur l'intervalle [0,2\pi]. Aplliquer l'algorithme à quelque exemple et remarquer qu'on obtient une erreur inférieure à 10^{-12} en moins de 4 itérations.


Effet de marée et limite de Roche

Auteur : S. Renner

Date de création: 31 janvier 2011

On cherche à déterminer la distance, que l'on notera d_R, en dessous de laquelle un satellite commence à se disloquer sous l'action des forces de marée causées par la planète autour de laquelle il orbite.

Cette distance théorique s'appelle la limite de Roche. Elle tire son nom de l'astronome français qui l'a formulée en 1850.


Ex: Effet de marée et limite de Roche

Auteur: S. Renner

exerciceEffet de marée et limite de Roche

Difficulté : ☆☆   Temps : 1h

On se place dans un référentiel galiléen centré sur la planète, supposée sphérique, homogène, de masse M, rayon R, et masse volumique \rho_P.

On suppose que le satellite est constitué de deux sphères homogènes identiques (de masse m, de rayon r et de masse volumique \rho), et qu'il est en orbite circulaire de rayon d autour de la planète.

roche_limite.jpeg
Bilan des forces agissant sur une masse m dans le repère lié à la planète (supposé galiléen) : {\bf F_{P1}} est la force d'attraction gravitationnelle de la planète, {\bf F_m} celle de la deuxième sphère de masse m, {\bf C_1} la force de contact entre les deux sphères.
Crédit : Astrophysique sur Mesure / Renner
Question 1)

En appliquant le principe fondamental de la dynamique (PFD) au système des deux masses m, montrer qu'au premier ordre en r/d la vitesse angulaire de la comète est \displaystyle \omega = \sqrt{\frac{G M}{d^3}}.

Question 2)

Appliquer le PFD à l'une des masses m, et trouver un critère de fragmentation du satellite (contact rompu entre les deux sphères) au premier ordre en r/d.

Question 3)

En déduire l'expression de la limite de Roche d_R en fonction de M, m, r, puis en fonction de \rho_P, \rho, R.

Question 4)

Retrouver l'expression de la limite de Roche en écrivant que la différence de force entre les deux masses m due à l'attraction gravitationnelle de la planète est supérieure à la force de gravitation mutuelle entre les deux sphères.


Générateur de coefficients de Taylor

Auteur : Marc Fouchard

Date de création: 8 Mai 2013

En astronomie on a souvent besoin de reproduire le mouvement des objets en effectuant des intégraitons numériques. Quelque soit la méthode, elle repose toujours sur des approximations et elle nécessite un certain temps de calcul. Ainsi l'objectif d'une méthode est de trouver le meilleur compromis entre temps de calcul et précision.

Une méthode très performante est la méthode dite de Taylor.

Soit l'équation différentielle suivante :

\mathbf{y}' = f(t,\mathbf{y}), \quad \mathbf{y}(t_0)=\mathbf{y}_0.

Soit h, un pas d'intégration. Le but est de trouver \mathbf{y}(t_0+h) avec la meilleure précision possible.

La méthode de Taylor repose sur le développement de \mathbf{y} en t_0 :

\mathbf{y}(t_0+h)=\mathbf{y}(t_0) + \mathbf{y}'(t_0) h + \frac{1}{2!}\mathbf{y}''(t_0) h ^2 + \cdots + \frac{1}{k!}\mathbf{y}^{(k)}(t_0) h^k + \cdots

Le problème principal est donc d'estimer les coefficients du développement de Taylor de la solution.


Ex: générateur de coefficients de Taylor

Auteur: Marc Fouchard

exerciceGénérateur de coefficients de Taylor

Difficulté : ☆☆   Temps : 30mn

Question 1)

Soit f, une fonction analytique. On note (f)_i, le ième coefficient de Taylor de la fonction f en t_0 défini par :

(f)_i = \frac{f^{(i)}(t_0)}{i!}.

Soit p et q, deux fonctions analytiques, déterminer (p \pm q )_i et (pq)_i

Question 2)

Déterminer \left(\frac{p}{q}\right)_i.

Question 3)

Montrer que :

(p^a)_i=\frac{1}{p}\sum_{r=0}^{i-1} \left ( a -\frac{r(a+1)}{i}\right) (p)_{i-r}(p^a)_r


Limites de la loi de Planck

Auteur : Stéphane Erard

Date de création: 30 Mai 2013

La loi de Planck donne l'expression exacte du spectre du corps noir. Historiquement elle a été dérivée après des approximations valables aux grandes et courtes longueurs d'onde, qui sont toujours utilisées dans certaines situations. L'exercice propose de retrouver ces approximations à partir de la loi complète.


Ex: Limites de la loi de Planck

Auteur: Stéphane Erard

exerciceLimites de la loi de Plack

Difficulté :    Temps : 30 min

On connaît la luminance du corps noir en fonction de la longueur d'onde, donnée par la loi de Planck (voir par exemple l'exercice sur la loi de Wien) : B_{\lambda} = \frac{2hc^2}{\lambda^5 (e^{hc/kT\lambda} -1) }

c est la vitesse de la lumière dans le vide, h la constante de Planck, k la constante de Boltzmann, \lambda la longueur d'onde et T la température du corps noir.

Cette expression donne la luminance directionnelle, en W\,m^{-2}\,sr^{-1}\,\mu m^{-1}, du corps noir à toutes les longueurs d'onde.

Question 1)

Donner une expression de cette loi à courtes longueurs d'onde. Cette expression est connue sous le nom de loi ou distribution de Wien.

Question 2)

Donner une expression de cette loi à grandes longueurs d'onde. Cette expression est connue sous le nom de loi de Rayleigh-Jeans. Commentaire ?

Question 3)

Tracer les graphiques de ces deux expressions et de la loi de Planck en échelle log/log, comparer. Quel problème pose l'approximation de Rayleigh-Jeans ?


Réponses aux exercices

pages_dl/lensing-ex.html

Exercice 'Lentille gravitationnelle'


pages_dl/exo-reperage-astre.html

Exercice 'Repérage des astres'


pages_dl/eq-kepler-exo.html

Exercice 'Equation de Kepler'


pages_dl/exo-limite-roche.html

Exercice 'Effet de marée et limite de Roche'


pages_dl/exo-coef-taylor.html

Exercice 'Générateur de coefficients de Taylor'


pages_dl/exo-dl-planck.html

Exercice 'Limites de la loi de Plack'