Voir la source de Effet Yarkovsky

{{À déjargoniser|date=novembre 2023}}[[Image:YarkovskyEffect.svg|vignette|300px|Effet Yarkovsky<br />1. Radiation de la surface de l'astéroïde <br />2. Rotation de l'astéroïde<br />3. Orbite de l'astéroïde <br />4. Radiation du soleil]]
L’'''effet Yarkovsky''' est la conséquence du [[rayonnement]] [[Thermodynamique|thermique]] reçu par un [[astéroïde]] sur son déplacement.

Les astéroïdes concernés ont un diamètre inférieur à environ {{nombre|20|[[Kilomètre|km]]}}, et ceux de la [[Ceinture d'astéroïdes|ceinture principale]] peuvent être mis en [[résonance]], {{quoi|et regroupés en [[Objet géocroiseur|objets géocroiseurs]]}}<ref>Bottke {{Jr.}} {{et al.}}, 2002.</ref>.

== Notations ==
Soient les notations suivantes concernant l'astéroïde :
* ''R'' : rayon moyen ;
* ''D'' : diamètre moyen ;
* ''E'' : [[obliquité]] de l'axe de rotation ;
* ''ω'' : vitesse angulaire, ou « pulsation », de la rotation ;
* ''P'' : période de rotation ;
** ''ω'' = 2π/''P''
* ''κ'' : [[conductivité thermique]] ;
* ''ρ'' : [[masse volumique]] ;
* ''C''{{ind|p}} : [[chaleur spécifique]] à [[pression]] constante ;
* ''T'' : température de surface ;
* ''A'' : [[albédo]] de [[Albédo de Bond|Bond]] ;
* ''m'' : masse ;
* ''S'' : surface ;
* ''ε'' : [[Émissivité|émissivité thermique]] de la surface ;
ainsi que les constantes physiques suivantes :
* ''c'' : [[vitesse de la lumière]] ;
* ''σ'' : [[constante de Stefan-Boltzmann]].

== Généralités ==
[[Image:yarkov-force.png|250px|vignette|droite|{{abréviation discrète|Fig.|figure}} 1 – Schéma des composantes de la force résultante de l'effet Yarkovsky. C'est la composante le long de la trajectoire qui est responsable du changement du demi-grand axe de l'orbite de l'astéroïde.]]

Lors de la révolution d'un astéroïde autour du [[Soleil]], un côté va être plus illuminé que l'autre (et donc recevra plus de photons). Son insolation est donc asymétrique. Il faudra alors un certain temps τ pour réémettre ces [[photon]]s, pendant lequel l'astéroïde va tourner sur lui-même. Il en résulte donc une force opposée à sa direction de réémission (voir [[#Figure_3|{{nobr|{{abréviation discrète|Fig.|figure}} 3}}]], et la [[pression de radiation]]). L'effet Yarkovsky est donc la force qui résulte de cette réémission du photon : c'est un effet différé dans le temps, contrairement à la [[pression de radiation]] classique.

Cette force peut se décomposer, selon l'axe porté par la droite passant par le Soleil et l'astéroïde, et la droite qui lui est perpendiculaire (voir {{nobr|{{abréviation discrète|Fig.|figure}} 1}}). C'est la composante le long de la trajectoire<ref>Dans le cas d'une orbite circulaire, elle est confondue avec la composante tangentielle, mais pour une [[orbite]] excentrique, ce n'est pas le cas.</ref> qui a pour effet une augmentation séculaire du demi-grand axe<ref>Car on raisonne sur une révolution de l'astéroïde : pendant une révolution la composante radiale va être en moyenne nulle, tandis que la composante le long de la trajectoire va avoir un effet moyenné non nul.</ref> pour un sens de rotation prograde de l'astéroïde, et ''une diminution'' pour un sens rétrograde. Cette force modifie aussi l'[[Excentricité orbitale|excentricité]] et l'inclinaison  de l'orbite, et la vitesse de rotation de l'astéroïde.

L'effet Yarkovsky dépend des variables et paramètres suivants :
* distance au Soleil,
* inclinaison de l'axe de rotation,
* taille de l'astéroïde (diamètre moyen ''D''),
* forme de l'astéroïde,
* propriétés thermiques, et en particulier l'[[inertie thermique]],
* vitesse angulaire de rotation ''ω'' <ref>Si ''ω'' est nul ou tend vers l'infini, il n'y a pas d'effet Yarkovsky.</ref>.

On distingue deux effets Yarkovsky différents (voir {{nobr|{{abréviation discrète|Fig.|figure}} 2}}<!-- légende de la figure 2 disponible sur https://fr.wikipedia.org/w/index.php?title=Effet_Yarkovsky&oldid=5853138 mais je n’ai pas trouvé trace de la figure… -->) :
# '''l'effet diurne''', qui se caractérise par :
#* l'axe de rotation perpendiculaire au plan orbital,
#* une échelle de temps de l'inertie thermique de l'ordre de la période de rotation,
#* une dépendance du sens de rotation ;
# '''l'effet saisonnier''', qui se caractérise par :
#* l'axe de rotation dans le plan orbital,
#* une échelle de temps de l'[[inertie thermique]] de l'ordre de la période de révolution,
#* une indépendance du sens de rotation, ce qui a pour effet de circulariser les orbites.

== Brève théorie ==
Il y a deux étapes dans le calcul de la force Yarkovsky [Bottke Jr. {{et al.}}, 2002] :
# détermination de la distribution de température à la surface,
# évaluation de la force de recul due à la radiation thermique.

=== Température de surface ===
Pour déterminer la distribution de température à la surface de l'astéroïde, on doit prendre en compte l'inertie thermique de celui-ci, la vitesse de rotation, et la [[Conduction thermique|diffusion]] de chaleur. Cette dernière est trop complexe pour pouvoir la traiter analytiquement sans simplification. On la linéarise donc, ce qui revient à considérer des petites variations de chaleur [Vokrouhlicky, 1999].

{{Ancre|Figure_3}}
[[Image:yarkov-princ-diu.png|300px|vignette|centre|{{abréviation discrète|Fig.|figure}} 3 - L'effet Yarkovsky diurne. Un photon incident frappe la surface de l'astéroïde (en rouge), puis est réémis après un certain temps τ (caractéristique de l'inertie thermique) perpendiculairement à la surface (voir [[pression de radiation]]) de l'astéroïde (en vert). En moyenne (sur la surface ensoleillé) la force sera donc orientée d'un angle dα(τ) par rapport à la direction d'incidence des photons (en magenta).]]

Pour des astéroïdes de forme complexe et/ou de paramètres thermiques hétérogènes, le calcul nécessite un traitement plus sophistiqué pour éviter des calculs trop gourmands en temps. Ainsi, on déterminera la température de surface en utilisant le flux d'énergie à l'intérieur de l'astéroïde :
:<math>\vec{\nabla} \cdot (\kappa\ \vec{\nabla}T) = \rho\ C_{p} \frac{\partial T}{\partial t}</math>, 
et la condition limite à la surface de l'astéroïde :
:<math>(\kappa\ \vec{\nabla}T \cdot \vec{n}_{\perp}) + \varepsilon\ \sigma\ T^{4} = \alpha\ E</math> .

La dernière équation se réfère à un [[élément de surface]] orientée par le vecteur unitaire <math>\vec{n}_{\perp}</math> ; ''E'' est le flux solaire à travers cette surface, et <math>\alpha\ = 1 - A</math>. Cette équation est une condition limite pour la détermination de la température ''T'' mais requiert la connaissance de la forme et l'état de rotation de l'astéroïde par le biais de ''E'', ce qui implique de longs temps de calcul.

Dans le cas de petites variations thermiques, on peut dire que la température ''T'' autour du corps est proche de ''T''{{ind|0}} (c'est-à-dire <math>T = T_{0} + \Delta T</math>, la variation de température étant petite devant ''T''{{ind|0}} :  <math>\delta = \frac{\Delta T}{T_{0}} \ll 1</math>). Si ''T''{{ind|0}} est constant, ces deux équations peuvent être réécrites pour la variable δ, tandis qu'on peut écrire le terme de température de la dernière équation comme suit [Vokrouhlicky, 1998] :
:<math>T^{4} \sim T_{0}^{4}\ (1 + 4\delta + \ldots)</math>.

=== Évaluation de la force ===
L'astéroïde conserve presque toujours une asymétrie dans son [[irradiation]] (à l'exception des très petits corps de moins de {{unité|10|centimètres}} de diamètre, qui sont gardés isothermes par une conduction de chaleur interne efficace). Il en résulte donc une direction privilégiée vers laquelle la force va agir [Vokrouhlicky {{et al.}}, 2005].

À ce stade, il parait utile de mettre en évidence quelques grandeurs caractéristiques. Pour un terme de Fourier donné à la fréquence ν de la décomposition de ''E'', deux paramètres fondamentaux apparaissent [Vokrouhlicky {{et al.}}, 2000 ; Bottke {{Jr.}} {{et al.}}, 2002] :
* la profondeur de pénétration de l'onde thermique dans l'astéroïde : <math>l_{\nu} = \sqrt{\frac{\kappa}{\rho\ C_{p}\ \nu}}</math>,
* le paramètre thermique <math>\Theta_{\nu} = \frac{\sqrt{\kappa\ \rho\ C_{p}\ \nu}}{\varepsilon \sigma T_{\star}^{3}}</math>, qui permet de mesurer la relaxation entre l'absorption et la réémission à la fréquence ν d'un photon, ''T''{{ind|*}} étant la température subsolaire définie comme <math>\varepsilon \sigma T_{\star}^{4} = \alpha E_{\star}</math>, où <math>E_{\star}</math> est le flux de radiation solaire à la distance de l'astéroïde.

Pour les astéroïdes qui viennent de la ceinture principale, on pourra considérer une excentricité ''e'' petite, ce qui permet de linéariser le terme δ. Enfin, on supposera une émission isotropique des photons, ce qui correspond à la [[loi de Lambert]].

Avec toutes les simplifications précédentes, on obtient l'expression de la force due à l'effet Yarkovsky suivante : 
:<math>{\mathbf f} = -\ \frac{2}{3}\ \frac{\varepsilon \sigma}{m c} \int_{S} dS(u,v)\ T^{4}\ {\mathbf n}_{\perp}</math>,
où '''f''' est la force par unité de masse, et (''u'',&nbsp;''v'') les coordonnées de l'[[élément de surface]].

Si on adopte un système local de [[Système de coordonnées|coordonnées]] avec ''z'' selon l'axe de rotation de l'astéroïde et ''x'',&nbsp;''y'' dans le plan équatorial, on différencie alors deux variantes de <math>\vec{f}</math> :
* l'effet diurne : ''f{{ind|x}}'', ''f{{ind|y}}'' dépendent de la fréquence de rotation ω <ref>avec une coupure sans grande importance à <math>\omega = \pm\ \nu_{mm}</math>, ν{{ind|mm}} étant la fréquence de [[mouvement moyen]], ν{{ind|mm}}&nbsp;=&nbsp;2π/''P''.</ref>,
* l'effet saisonnier : ''f{{ind|z}}'' dépend uniquement de ν{{ind|mm}}&nbsp;=&nbsp;2π/''P'', la fréquence de mouvement moyen.
Dans un formalisme plus complet (c'est-à-dire sans linéarisation de la température), les deux effets sont couplés.

Ces forces agissent sur le demi-grand axe de l'orbite de l'astéroïde. En moyennant ces petites perturbations sur une révolution, on a :
:<math>\left(\frac{{\rm d}a}{{\rm d}t}\right)_{\rm diurne} = -\ \frac{8\alpha}{9}\ \frac{\phi}{\nu_{mm}}\ F_{\omega}(R',\Theta_{\omega}) \cos\epsilon\ +\ \mathcal{O}(e)</math>,
:<math>\left(\frac{{\rm d}a}{{\rm d}t}\right)_{\rm saisonnier} = \frac{4\alpha}{9}\ \frac{\phi}{\nu_{mm}}\ F_{\nu_{mm}}(R',\Theta_{\nu_{mm}}) \sin^{2}\epsilon\ +\ \mathcal{O}(e)</math>,
où <math>\phi = \pi R^{2}\mathcal{E}_{0}/(mc)</math> est le coefficient de pression de radiation, <math>R' = R/l_{\nu}</math> une longueur normalisée, et <math>F_{\nu}</math> la fonction définie comme suit<ref>Cette fonction est donnée pour mieux comprendre les équations précédentes. De plus amples détails peuvent être trouvés dans la littérature.</ref>:
:<math>F_{\nu}(R',\Theta) = -\ \frac{K_{1}(R')\ \Theta_{\nu}}{1 + 2K_{2}(R')\ \Theta_{\nu} + K_{3}(R')\ \Theta_{\nu}^{2}}</math>.

Les paramètres <math>K_{i}(R')</math> sont des fonctions analytiques de <math>R'</math>. La dernière équation est la même pour les deux effets Yarkovsky, en remplaçant <math>\nu = \omega</math> pour l'effet diurne, et <math>\nu = \nu_{\rm mm}</math> pour l'effet saisonnier. On peut remarquer que le paramètre <math>\Theta_{\nu}</math> est le paramètre qui fait la plus grande différence entre les deux effets.

Les paramètres importants dans le changement du demi-grand axe par l'effet Yarkovsky sont donc :
* '''l'obliquité et la rotation''' : l'effet saisonnier produit toujours une diminution de ''a'' et est maximum à <math>\epsilon = 90^{\circ}</math>, alors que l'effet diurne peut produire une augmentation (<math>\epsilon < 90^{\circ}</math>) et une diminution (<math>\epsilon > 90^{\circ}</math>) de ''a'', et est maximum pour <math>\epsilon = 0^{\circ}</math>. De plus, l'effet Yarkovsky est négligeable dans la limite des rotations à vitesse infinie ou nulle car les variations de température sont nulles le long d'une même longitude.
* '''la taille''' : l'effet Yarkovsky est négligeable pour les très petits astéroïdes car l'onde thermique chauffe de la même manière toute sa surface (donc typiquement quand ''D'' est de l'ordre de <math>l_{\nu}</math>), et pour les grands objets car la gravitation l'emporte (à partir de ''D'' égal à {{unité|20|kilomètres}}).
* '''conductivité thermique''' : pour les corps très poreux ou ayant une surface semblable au régolithe, κ est très faible (de l'ordre de 10{{exp|-3}}&nbsp;W·m{{exp|-1}}·lK{{exp|-1}}). Pour les corps rocheux (chondrites ou glace), κ est intermédiaire ({{unité|1 W m-1 K-1}}), alors que pour les corps ferreux κ est très élevé ({{unité|40 W m-1 K-1}}). Les variations de κ modifient les paramètres <math>l_{\nu}</math> et <math>\Theta_{\nu}</math> : un κ petit implique que <math>l_{\nu} \approx 0</math>, et donc <math>\dot{a} \approx \Theta</math> et l'effet Yarkovsky disparaît ; un κ grand implique que <math>R'\rightarrow 0</math> et <math>l_{\nu}\rightarrow \infty</math>, et donc <math>\dot{a} \approx \frac{R'^{2}}{\Theta}</math>, ce qui correspond à un équilibre thermique de l'astéroïde.
* '''distance au Soleil''' : l'effet Yarkovsky diminue quand la distance de l'astéroïde au Soleil augmente.

Les valeurs caractéristiques de changement du [[Grand axe|demi-grand axe]] sont de l'ordre de {{nombre|0,1|[[unité astronomique]]}} pour les objets sub-kilométriques, et d'environ {{nombre|0,01|ua}} pour les objets de diamètre supérieur au kilomètre{{quoi}}. Par ailleurs, le paramètre κ a de moins en moins d'importance au fur et à mesure que la taille de l'astéroïde augmente.

Ainsi, l'effet Yarkovsky est le premier mécanisme pour les astéroïdes ayant un diamètre inférieur à {{unité|20|kilomètres}} qui leur permette d'échapper à la ceinture principale et est susceptible de les amener sur des orbites de type [[Objet géocroiseur|géocroiseur]].

Outre le changement du demi-grand axe, qui est l'effet le plus important et le plus prononcé de l'effet Yarkovsky, celui-ci se traduit aussi par des changements de l'excentricité et l'[[Inclinaison orbitale|inclinaison]] de l'orbite de l'astéroïde.

== Étude de l'effet Yarkovsky ==
{{...}}

L'étude de l'effet Yarkovsky est l'un des objectifs de la mission [[OSIRIS-REx]] de la [[National Aeronautics and Space Administration|NASA]]. À partir de la fin de l'année 2018, la sonde OSIRIS-REx orbitera autour de l'astéroïde [[(101955) Bénou|Bénou]], un astéroïde [[géocroiseur]], pour en étudier avec plus de précision l'effet Yarkovsky afin de mieux déterminer la trajectoire de ce corps céleste. Cet effet avait déjà, au préalable, été étudié sur cet astéroïde à l'aide de télescopes terrestres et spatiaux montrant que son orbite s'était décalée de {{nombre|160|km}} depuis sa découverte en 1999<ref>{{lien|lang=en|trad=Dante Lauretta}}, ''[[Pour la science|Pour la Science]]'', {{date-|octobre 2016}}, {{p.|38-45}}.</ref>. 

== Notes et références ==
{{Références}}

== Voir aussi ==
=== Bibliographie ===
* {{en}} [Bottke {{Jr.}} {{et al.}}, 2002] [[William F. Bottke|W.F. Bottke {{Jr.}}]], {{lien|lang=de|trad=David Vokrouhlický|texte=D. Vokrouhlicky}}, D.P. Rubincam et M. Broz (2002). ''{{lang|en|The effect of Yarkovsky thermal forces on the dynamical evolution of asteroids and meteoroids}}'', in {{lang|en|{{nobr romains|Asteroids III}}}} (W.F. Bottke {{Jr.}} {{et al.}}, eds.), {{p.|395–408}}. [[Université de l'Arizona|{{lang|en|Univ. of Arizona}}]], ([[Tucson]]). {{ads|2002aste.conf..395B}}.
* {{en}} [Vokrouhlicky, 1998] D. Vokrouhlicky (1998). ''{{lang|en|Diurnal Yarkovsky effect as asource of mobility of meter-sized asteroidal fragments}}''. {{lang|en|[[Astronomy & Astrophysics]]}}, {{vol.|335}}, {{p.|1093–1100}}. {{ads|1998A&A...335.1093V}}.
* {{en}} [Vokrouhlicky, 1999] D. Vokrouhlicky (1999). ''{{lang|en|A complete linear model for the Yarkovsky thermal force on spherical asteroid fragments}}''. {{lang|en|[[Astronomy & Astrophysics]]}}, {{vol.|344}}, {{p.|362–366}}. {{ads|1999A&A...344..362V}}.
* {{en}} [Vokrouhlicky {{et al.}}, 2000] D. Vokrouhlicky, A. Milani et S.R. Chesley (2000). ''{{lang|en|Yarkovsky effect on small near-Earth asteroids : mathematical formulation and examples}}''. [[Icarus (journal)|Icarus]], {{vol.|148}}, {{p.|118–138}}. {{ads|2000Icar..148..118V}} {{DOI|10.1006/icar.2000.6469}}.
* {{en}} [Vokrouhlicky {{et al.}}, 2005] D. Vokrouhlicky, M. Broz, W.F. Bottke, D. Nesvorny et A. Morbidelli (2005). ''{{lang|en|Non-gravitational perturbations and evolution of the asteroid main belt}}'', in {{lang|en|Dynamics of Populations of Planetary Systems}} (Z. Knezevic et A. Milani, eds.), {{numéro|197}} in {{lang|en|Proceedings IAU Colloquium}}, {{p.|145–156}}. [[Union astronomique internationale|{{lang|en|International Astronomical Union}}]], {{lang|en|Cambridge Academic Press}}, ([[Cambridge]]). {{ads|2005dpps.conf..145V}} {{DOI|10.1017/S1743921304008609}}.

=== Articles connexes ===
* [[Effet YORP]]
* [[Effet Magnus]] (bien que concerné par des forces de nature très différente, il concerne aussi un changement de trajectoire liée à une rotation de l'objet en déplacement dans un fluide, donc hors vide spatial)
* [[Ivan Yarkovsky]]
* [[Paramètres d'accélération non gravitationnelle]]
* [[Trou de serrure gravitationnel]]
* [[:Catégorie:Planète mineure où l'effet Yarkovsky a été détecté|Planètes mineures où l'effet Yarkovsky a été détecté]]

=== Liens externes ===
* {{lien brisé|url=http://www.hypparque.fr/spip.php?article14 |titre=Rapport de stage sur les effets des perturbations non gravitationnelles sur l'évolution du spin des NEAs}}
* {{en}} Article arXiv sur l'influence de la dérive due à l'effet Yarkovsky  et prédiction sur les astéroïdes géocroiseurs connus (C. R. Nugent, A. Mainzer, J. Masiero, T. Grav et J. Bauer) [https://arxiv.org/pdf/1206.4299.pdf {{pdf}}]
* {{site web|langue=en |auteur1={{lien|lang=de|trad=David Vokrouhlický}} |auteur2=[[William F. Bottke]] |titre=Yarkovsky and YORP effects |url=http://www.scholarpedia.org/article/Yarkovsky_and_YORP_effects|site=[[Scholarpedia]]}} 

{{Palette|Astéroïde}}
{{Portail|astronomie}}

[[Catégorie:Effet Yarkovsky|*]]