LES PASSAGES DE MERCURE ET DE VÉNUS DEVANT LE SOLEIL
FICHES PÉDAGOGIQUES: FICHE N° 02

Afficher cette fiche au format Acrobat Reader (204Ko)

Fiche n°02 : Comment calculer les circonstances locales dun passage de Vénus devant le Soleil.

I. Positions des astres par rapport au centre de la Terre.

La position dun corps céleste dans lespace est connue lorsque lon connaît ses coordonnées dans un repère tridimensionnel orthonormé que lon sest fixé. La première chose à définir est donc ce repère tridimensionnel. Pour cela il faut définir trois directions orthogonales particulières définissant trois axes Ox, Oy, Oz et une origine O. On pourrait prendre des directions quelconques, mais il est judicieux de choisir des directions particulières qui varient le moins possible par rapport au globe terrestre. Cest pourquoi on prend comme direction de laxe Oz la direction de laxe de rotation de la Terre sur elle-même (direction du pôle céleste). Le plan normal à cette direction (plan Oxy) est un plan parallèle à léquateur terrestre. Dans ce plan on choisit également une direction Ox particulière qui bouge le moins possible au cours du temps. On choisit la direction de léquinoxe de printemps intersection du plan de léquateur terrestre et du plan orbital de la Terre (plan de lécliptique). En fait, ce trièdre nest pas totalement fixe car laxe de rotation de la Terre est animé de petits mouvements. Mais à chaque instant on peut définir un repère instantané qui porte le nom de repère équatorial vrai. Ce repère peut être centré sur le centre de notre choix, si cest le Soleil, le repère est dit héliocentrique, si cest le centre de la Terre le repère est dit géocentrique et si cest un lieu quelconque sur la Terre le repère est dit topocentrique.

Dans un repère cartésien, la position de chaque corps du système solaire est connu grâce à ses coordonnées cartésiennes (x,y,z). En réalité les astronomes mesurent plutôt des angles que des distances ; cest pourquoi, parallèlement aux coordonnées cartésiennes ils utilisent également un système de coordonnées polaires construit à partir du même repère. Ce système de coordonnées comportent trois coordonnées, une longueur (D) et deux angles (a et d). La longueur ou rayon vecteur est la distance entre lorigine O du repère et le corps. Langle a sappelle lascension droite et est langle entre la projection orthogonale de la direction du corps dans le plan de léquateur (Oxy) et laxe Ox. Le second angle d sappelle la déclinaison, cest langle entre la direction de lastre et sa projection orthogonale dans le plan Oxy, voir figure 1.

Lorsque lon utilise uniquement les deux angles, tous les corps sont supposés être à la même distance ( D = 1), ce système de coordonnées polaires porte alors le nom de sphère céleste.



Figure 1. – La sphère céleste équatoriale.

À chaque instant t les théories planétaires permettent de calculer les coordonnées cartésiennes apparentes géocentriques des planètes qui donnent la position apparente de chaque planète telle quelle est vue depuis le centre de la Terre. Ces coordonnées ne sont pas les positions géométriques des corps, elles tiennent compte de toutes les corrections et aberrations liées aux mouvements des corps et à la vitesse finie de la propagation de la lumière (nutation, précession, aberration planétaire et temps de lumière), seule la réfraction atmosphérique nest pas prise en compte.

On passe dun système de coordonnées à lautre à laide des formules suivantes :

(1) (2)

Soient a, d et D les coordonnées équatoriales apparentes géocentriques du Soleil, av, dv et Dv les coordonnées équatoriales apparentes géocentriques de Vénus, alors Da = av - a et Dd = dv - d. sont les coordonnées différentielles équatoriales sphériques de la planète par rapport au Soleil.



Figure 2. - Coordonnées différentielles.

Sur la sphère céleste équatoriale géocentrique, soit S le centre du Soleil et V le centre de Vénus. Si lon appelle D langle entre S et V et w langle entre la direction du pôle et la direction SV compté dans le sens trigonométrique pour un observateur situé au centre de la sphère céleste géocentrique. Comme le système de coordonnées sphériques est direct, les axes rectangulaires, vus depuis le centre de la sphère sont orientés comme sur la figure 2. Le plan xSy est le plan tangent à la sphère céleste en S. V est la projection de V sur ce plan tangent, alors les coordonnées différentielles planes (X,Y) de SV sont données en radians par :

(3)

La distance SV entre le centre du Soleil et le centre de Vénus dans le plan tangent est alors égale à : .

Si lon note d et d les demi-diamètres apparents du Soleil et de Vénus. Les instants t1 et t4 des contacts extérieurs géocentriques sont solutions de léquation SV=d + d et les instants t2 et t3 des contacts intérieurs géocentriques sont les solutions de léquation SV=d – d. Ces calculs se font par approximations successives.

Les projections des demi-diamètres apparents de Vénus et du Soleil se calculent à partir des rayons rv et r de Vénus et du Soleil exprimés en kilomètres et des distances Dv et D de Vénus et du Soleil à la Terre par les formules suivantes :

    
et

Linstant tm du minimum de distance entre Vénus et le Soleil se calcule de même en calculant le zéro de la fonction et la distance minimale est donnée par la valeur de SV à cet instant.

II. Le calcul des circonstances locales.

Figure 3. - Coordonnées géographiques.

La position dun point M à la surface du globe terrestre est définie par ses coordonnées géographiques (figure 3) : sa latitude j par rapport à léquateur et sa longitude l par rapport à un méridien origine (méridien de Greenwich). La Terre tourne autour de son axe dun mouvement quasi uniforme en raison de un tour en 23h 56m 4s de temps universel (révolution sidérale de la Terre). Pour connaître à un instant t la position dun observateur par rapport au repère Oxyz, il suffit de connaître à cet instant langle que fait le méridien de Greenwich avec la direction de laxe Ox. Cet angle TG sappelle le temps sidéral vrai à Greenwich à linstant t. Langle TL que fait le méridien du lieu dobservation M avec laxe Ox au même instant, et qui porte le nom de temps sidéral local, est alors égal au temps sidéral vrai à Greenwich plus ou moins la longitude du lieu (en fonction de la position est ou ouest du lieu par rapport à Greenwich). Ainsi si les longitudes sont comptées positivement vers louest et négativement vers lest on a :

Alors les coordonnées cartésiennes du vecteur OM ( x,h,z) du lieu dans le repère équatorial géocentrique vrai se calculent à laide des formules suivantes :

          (4)

R étant le rayon terrestre.

Attention les unités utilisées doivent être cohérentes, si le coordonnées géocentriques sont exprimées en unités astronomiques, R doit également être en unités astronomiques, alors R est égal au rapport du rayon équatorial terrestre exprimé en kilomètres sur la valeur de lunité astronomique exprimée également en kilomètres ; cest donc aussi la valeur de la parallaxe horizontale équatoriale moyenne du Soleil exprimé en radians.

On peut raffiner le calcul en tenant compte de laplatissement terrestre f et de laltitude h de lobservateur ; les formules deviennent alors :

          (5)

avec

          

La hauteur h de lobservateur est exprimée dans la même unité que R.

Le calcul direct.

De nos jours, le calcul direct à laide dun ordinateur ne pose aucun problème.

Pour cela il suffit de reprendre les formules utilisées pour le calcul géocentrique en remplaçant les coordonnées géocentriques par les coordonnées topocentriques. Les vecteurs positions topocentriques MV et MS joignant le point M au centre V de Vénus et au centre S du Soleil sobtiennent par différences des vecteurs géocentriques correspondants et du vecteur OM joignant le centre de la Terre à lobservateur.

Le calcul approché.

La méthode approchée consiste à faire un certain nombre de simplifications.

Comme Da, Dd et D sont de lordre du demi-diamètre solaire, on peut dans les équations donnant les coordonnées différentielles planes tangentielles X et Y négliger les quantités du troisième ordre sans introduire sur D et w des erreurs supérieures au centième de seconde. Ce qui permet de mettre la formule (1) sous la forme simplifiée suivante :

           (6)

Le terme sin1" permet davoir les différentes coordonnées exprimées en secondes de degré, cest la valeur en radians dun angle de une seconde de degré. Si les valeurs Da et Dd sont en secondes de degré, on obtient les valeurs X et Y dans la même unité.

Calcul de la parallaxe en coordonnées différentielles tangentielles

Notons les paramètres et coordonnées topocentriques avec un indice p.

Les équations précédentes pour un lieu donné sécrivent :

           (7)

Dans les seconds termes de ces équations qui sont de lordre du carré du demi-diamètre solaire, on peut négliger la parallaxe, cest-à-dire remplacer les coordonnées topocentriques par des coordonnées géocentriques, sans quil en résulte derreur atteignant un dixième de seconde dans les valeurs de Xp et Yp. Ce qui donne pour ces équations à ce degré dappro-ximation :

           (8)

De même si lon prend la déclinaison géocentrique d du Soleil à la place de sa déclinaison topocentrique dans le calcul du cosinus du premier terme, lerreur sera inférieure à 0,02".

Les équations ont donc la forme finale suivante :

           (9)

et si lon pose :

nous aurons :

U et V sont donc les expressions de la parallaxe en coordonnées différentielles tangentielles.

Nous allons faire apparaître dans ces expressions les coordonnées cartésiennes de lobservateur. Si lon désigne de nouveau par , les coordonnées cartésiennes équatoriales du lieu dobservation, à lépoque t, nous avons par les formules de la parallaxe en coordonnées équatoriales :

Dans les deux dernières formules on peut remplacer les sinus et les cosinus de la déclinaison et de lascension droite de Vénus par ceux du Soleil sans commettre derreur supérieure à 0,2".

En substituant ces valeurs dans les expressions de U et V, et en posant :

nous trouvons :

           (10)

or on a également par définition :

doù lon tire

doù en négligeant le carré de la parallaxe solaire,

           (11)

Les coordonnées équatoriales cartésiennes équatoriales de lobservateur en négligeant laplatissement terrestre sont les suivantes :

p0 est la valeur de la parallaxe horizontale équatoriale moyenne du Soleil exprimée en secondes de degré.

En substituant ces valeurs dans la formule (6) donnant U et V, on trouve :

           (12)

Et en développant ces expressions et en regroupant les termes dépendant uniquement du lieu on obtient lexpression suivante :

           (13)

avec

           (14)

Alors la formule (11) sécrit sous la forme suivante :

           (15)

Les coefficients a, b et c sont indépendants de lobservateur, mais varient avec le temps, principalement à cause des changements rapides que subissent les angles TG et w durant le passage.

On remarque également que lon a la relation suivante :

Calcul des instants de contacts topocentriques.

Nous allons utiliser les formules précédentes pour calculer les instants des contacts topocentriques à partir des instants des contacts géocentriques.

Si lon affecte de lindice 1 les paramètres des premiers contacts géocentriques, cest à dire lun des cas suivants :

Le premier contact géocentrique extérieur :
ou le premier contact géocentrique intérieur :
t1 étant lépoque de ces contacts.

Pour un lieu quelconque de latitude j et de longitude l lépoque du contact t1+ Dt vérifie léquation approchée suivante :

r1 désignant la somme des demi-diamètres ou la différence des demi-diamètres en fonction du contact considéré.

En remarquant que et en négligeant les carrés et les produits du second ordre, on trouve  :

ce qui sécrit également

Mais on a également pour une époque t quelconque et pour le lieu M la relation :

On en déduit que :

           (16)

Donc pour un point quelconque les instants des premiers contacts se calculent à laide de la formule suivante :

           (17)

On peut faire un raisonnement identique pour les derniers contacts et lon aura les formules suivantes, lindice 2 correspondant aux derniers contacts.

           (18)

et

           (19)

Autres formulations.

On trouve également une formulation légèrement différente de ces deux relations.

En remarquant que :


et que si lon note la parallaxe de Vénus et


la parallaxe solaire exprimées en secondes de degré ; R, D et Dv doivent être exprimés en unités astronomiques, on a également

donc :

Lexpression Dt prend alors la forme

Et en intégrant

dans les coefficients a, b et c et en rappelant que X=D sinw et que Y=D cosw et en faisant apparaître la variable l'expression Dt devient  

           (20)

avec :

           (21)

Langle HG sappelle langle horaire, cest langle formé par la direction du méridien sud et la projection de la direction du corps (ici le Soleil) dans le plan équatorial (Oxy).

Pour calculer chaque contact du passage de Vénus en un lieu quelconque, il suffit de calculer les valeurs de ces trois coefficients pour chacun des contacts puis dappliquer la relation précédente donnant le Dt pour chaque contact en chaque lieu considéré. Cette méthode est donc très rapide car elle ne nécessite que la connaissance de douze coefficients pour chaque passage de Vénus. Par contre sa précision nest que de lordre du dixième de minute de temps, ce qui est suffisant pour prédire les instants des passages, mais est insuffisant si lon désire réduire des observations et en déduire une valeur de la parallaxe solaire.

Contact

Pour tous renseignements concernant ces pages, veuillez contacter ou le Service de renseignements de l'IMCCE.
Caution : This Website was created with the Ministère de l'Education Nationale, the CNRS and the CNES support. Any use of the data published on this website requires the IMCCE agreement..