Dynamique des grandes échelles de la convection dans la photosphère solaire

Dynamique des grandes échelles de la convection dans la photosphère solaire

Granulation solaire

 Comment caractériser quantitativement les propriétés de la convection photographique en termes de dynamique des fluides, en se basant sur les observations et les simulations ? Diverses estimations s’avèrent utiles pour la compréhension de ces processus physiques. La valeur de la diffusivité thermique 7 à proximité de la surface est κ ‘ 2 108 m2 s −1 . Cette donnée permet de calculer : – le temps de thermalisation d’un granule, τth = L 2 g κ ‘ 103 s ; (1.3) – le temps dynamique d’un granule, qui correspond au temps mis par une particule fluide pour parcourir Lg, τdyn = Lg Vg ‘ 103 s .On suppose que l’approximation de diffusion est valable dans le domaine considéré. Cette hypothèse est de toute évidence mise en défaut dans les zones supérieures de la photosphère solaire où le libre parcours moyen des photons est important. Une tentative de synthèse 41 On voit que ces deux temps sont du même ordre de grandeur, et proches du temps de vie caractéristique d’une structure granulaire. Cela se traduit par Pe ≡ τth τdyn = LgV κ = O(1), (1.5) en introduisant le nombre de Peclet. Lorsque Pe < 1, la turbulence n’est pas assez vigoureuse pour « apporter » des fluctuations de température en un point, celles-ci étant rapidement atténuées par la diffusion thermique. Au contraire, pour Pe > 1, la diffusion thermique agit sur des échelles de temps trop importantes pour gommer les fluctuations turbulentes de température (Pe est l’analogue thermique du nombre de Reynolds). Autrement dit, le granule est avant tout une structure issue de l’équilibre entre les processus de diffusion et la force d’Archimède, qui assure l’interface entre le transport convectif et le transport radiatif. La plus petite échelle de longueur à laquelle il est possible d’observer des fluctuations de température significatives est donc de l’ordre de Lg = κ/Vg. De manière analogue, on peut former le nombre de Reynolds d’un granule, en utilisant la viscosité du gaz, estimée à ν ∼ 1 m2 s −1 en surface. On a Re = LgV ν ‘ 109 (1.6) Autrement dit, la granulation est laminaire du point de vue radiatif mais extrêmement turbulente du point de vue purement hydrodynamique. Le nombre de Prandtl, défini par Pr \ ν κ (1.7) est de l’ordre de 10−9 en surface 8 . Remarque Le nombre de Peclet est proche de 1 à la surface, mais augmente dès qu’on pénètre à τ > 1, pour atteindre des valeurs de l’ordre de 107 dans les profondeurs de la zone convective où le transfert radiatif redevient prédominant et les échelles de vitesse et de longueur caractéristique sont très grandes, contrairement à la surface (Brandenburg et al. 2000). En profondeur on a simultanément Pe  1 et Re  1, tandis que le nombre de Prandtl reste très petit (environ 10−6 ). 

Le modèle standard de la convection à la surface du Soleil

 Les observations des années 60 ont conduit Simon et Leighton à proposer le modèle physique « historique » de la supergranulation (Simon et Leighton 1964). Dans ce modèle, la granulation est identifiée à la profondeur de recombinaison de H+, tandis que la supergranulation trouve son origine dans la recombinaison en 8 On peut noter que ce paramètre différencie complètement les comportements de la convection solaire et de la convection dans le manteau terrestre, caractérisée par Pr → ∞, puisque le fluide qui compose celui-ci est hypervisqueux. profondeur de He++. Dans le cadre de l’étude linéaire de l’instabilité convective, des cellules peuvent en effet se développer sur des échelles horizontales du même ordre de grandeur que la hauteur de la couche de fluide. Demeure alors la question de la profondeur intermédiaire de recombinaison de He+, qui devrait elle aussi faire émerger une échelle distincte à la surface (figure 1.6). La découverte de la méso granulation par November et al. (1981), avec une échelle horizontale de 8 000 km, a donné un argument aux défenseurs de cette théorie, qui y ont vu la manifestation de la recombinaison de l’hélium ionisé une fois.  FIG. 1.6 – Modèle standard de la convection photographique. Ce scénario est principalement basé sur une relation linéaire (i. e. une superposition) entre les différentes échelles observées. 1.5.c Le modèle en question : et la turbulence ? Plusieurs ombres viennent cependant noircir ce tableau. Les travaux de Beckers (1968), Worden et Simon (1976), Lin et Kuhn (1992), pour ne citer que ceux ci, montrent que les faibles fluctuations de température observées s’identifient surtout au réseau magnétique et peuvent difficilement être attribuées à une structure convective transportant un flux d’énergie significatif. Rast (2003b) soutient également que les profondeurs d’injection d’énergie ne sont pas aussi nettes que dans la description du modèle standard et que dans tous les cas cet apport d’énergie tendrait plutôt à déstabiliser des structures cohérentes et à faire apparaître des petites échelles. En réalité, la description classique fait totalement abstraction de la nature turbulente de la convection photosphérique et des interactions entre les diverses échelles qui font partie d’un continuum, comme le montre le spectre de puissance de la figure 1.5. Nous avons vu au paragraphe 1.5.a que l’échelle de la granulation Une tentative de synthèse 43 est intrinsèquement liée au processus de convection et au refroidissement radiatif. En utilisant des arguments généraux de turbulence, il est alors tout à fait possible de faire apparaître les bons ordres de grandeurs pour les grandes échelles (Rieutord et al. 2000) en utilisant la granulation pour normaliser le spectre turbulent. Par exemple, la durée de vie caractéristique d’une échelle k varie comme τ ∼ k α−3 2 , (1.8) en faisant l’hypothèse de similarité pour le spectre de puissance de la turbulence E(k) ∼ k −α . Sachant que pour une turbulence 3D on doit avoir α < 3, les grandes structures doivent vivre plus longtemps. On peut ainsi estimer la valeur de l’échelle super granulaire pour α = 5/3 (turbulence universelle tridimensionnelle homogène et isotrope), connaissant les temps caractéristiques de la granulation et de la supergranulation. On trouve Lsg ‘ Lg  τsg τg 3/2 ‘ 108 m , (1.9) Cet exemple n’a bien sûr pas valeur de référence en lui-même, il montre simplement qu’une alternative à l’explication historique de Simon et Leighton (1964) et November et al. (1981), basée sur des hypothèses de dynamique des fluides plus réalistes, est possible et doit être considérée attentivement, au regard des observations et simulations récentes : les études de Rieutord et al. (2000), Ploner et al. (2000), Cattaneo et al. (2001), Roudier et al. (2003), Rast (2003b) suggèrent majoritairement des mécanismes d’interactions non-linéaires entre granules pour expliquer l’émergence de structures de taille supérieure. Le scénario turbulent possède malgré tout ses propres inconvénients, le premier d’entre eux étant qu’il n’explique pas actuellement le maximum de puissance observé à l’échelle super granulaire, et dans une moindre mesure à l’échelle méso granulaire. Pour éviter de pécher par naïveté dans ce sens et illustrer la complexité du problème, rappelons simplement quelques arguments qui suggèrent qu’un spectre en k −5/3 , par exemple, n’est pas approprié à la convection photosphérique, et que l’image de la turbulence universelle doit être utilisée avec précaution pour un tel écoulement. Au vu de la figure 1.5, il est difficile d’identifier le spectre de la convection solaire avec le spectre universel de Kolmogorov de la turbulence homogène isotrope avec une pente en k −5/3 . Nordlund et al. (1997) soulignent que la granulation apparaît à la fois observationnelle ment et numériquement comme un écoulement laminaire, en particulier dans sa partie montante, et relèvent que les hypothèses d’isotropie et d’homogénéité utilisées en turbulence classique sont complètement isolées dans le cadre de la photosphère, puisque l’échelle de hauteur de densité au niveau de la surface ne constitue qu’un dixième de la taille caractéristique d’un granule. Finalement nous conclurons ce paragraphe par une question : comment définir la zone inertielle de la convection photographique ? Si la granulation s’identifie à l’échelle d’injection d’énergie dans la turbulence, la méso granulation et la supergranulation appartiendraient pas à la zone inertielle, dont les propriétés constituent, avec l’échelle de dissipation de Kolmogorov, les principales prédictions de la théorie K41. Constatant ainsi que la dynamique au-delà de l’échelle d’injection est largement inconnue, la quête d’une explication simple à l’existence des méso granules et des super granules basée uniquement sur des arguments académiques de stabilité linéaire ou de turbulence semble malheureusement partiellement compromise…

Table des matières

Remerciements
Notations utilisées
Problématique générale
1 Observations et modèles de la photosphère solaire
1.1 Introduction
1.2 Phénoménologie de la zone convective du Soleil
1.2.a Convection en profondeur
1.2.b Convection photosphérique
1.3 Histoire de l’observation de la dynamique photosphérique
1.3.a Méthodes d’observation et d’analyse
1.3.b Données relatives à la granulation
1.3.c La dynamique de la supergranulation
1.3.d Le problème de la mésogranulation
1.3.e Magnétisme et dynamique dans la photosphère
1.4 Expériences numériques sur la convection turbulente
1.4.a Les pionniers
1.4.b Propriétés des modèles détaillés de la granulation
1.4.c Propriétés des simulations d’atmosphères polytropiques
1.4.d Simulations à grands rapports d’aspects
1.5 Une tentative de synthèse
1.5.a Granulation solaire
1.5.b Le modèle standard de la convection à la surface du Soleil
1.5.c Le modèle en question : et la turbulence ?
1.5.d Que nous apportent les simulations numériques ?
1.6 Conclusion
2 Équations MHD dans un milieu fortement stratifié
2.1 Introduction
2.2 Système physique étudié
2.2.a Cadre général
2.2.b Paramètres caractéristiques
2.3 Équations sans dimension
2.4 Différents équilibres hydrostatiques
2.4.a Les polytropes
2.4.b L’atmosphère isotherme
2.4.c Échelles de hauteur, temps caractéristiques
2.5 Conditions aux limites
2.6 Définitions diverses et variées
2.6.a Nombre de Rayleigh
2.6.b Équation pour l’énergie
2.6.c Nombre de Nusselt
3 Approche linéaire de la supergranulation solaire
3.1 Introduction
3.2 Justification des hypothèses de l’étude
3.2.a Utilisation de diffusivités turbulentes
3.2.b Conditions aux limites sur la température
3.3 Détails du modèle linéaire
3.3.a Linéarisation pour un fluide complètement compressible
3.3.b Forme anélastique
3.3.c Équation cubique dans l’approximation de Boussinesq
3.4 Résultats
3.4.a Méthode numérique
3.4.b Diverses configurations
3.5 Développements analytiques dans l’approximation de Boussinesq
3.5.a Résolution de l’équation cubique
3.5.b Dépendance de Rc vis-à-vis de Q pour a = 0
3.5.c Transition vers une longueur d’onde finie
3.6 Discussion
4 Simulations numériques à grand rapport d’aspect
4.1 Introduction
4.2 La simulation
4.2.a Paramètres physiques et numériques
4.2.b Commentaires
4.2.c Définition des spectres
4.3 Description de l’écoulement
4.3.a Évolution temporelle
4.3.b Régime quasi-stationnaire
4.4 Analyse dans l’espace de Fourier
4.4.a Quelle loi de puissance pour les spectres ?
4.4.b Dynamique et transfert non-linéaire dans l’espace spectral
4.5 Les simulations et la photosphère solaire
4.5.a Comparaison avec la convection radiative
4.5.b Interprétation dans le cadre photosphérique
4.5.c Et la supergranulation ?
4.6 Récapitulatif des résultats
5 Théories de champ moyen
5.1 Introduction
5.2 Effet AKA et viscosité turbulente
5.2.a Phénoménologie
5.2.b Expression des coefficients
5.2.c Écoulements non périodiques et instabilités de phase
5.3 Méthode numérique de résolution
5.3.a Principe général du code
5.3.b Disposition des données
5.3.c Blocs matriciels
5.3.d Conditions aux limites
5.3.e Algorithmes utilisés
5.4 Calculs de viscosités turbulentes dans des écoulements simples
5.4.a Écoulement de Kolmogorov
5.4.b Écoulement hexagonal simple
5.4.c Écoulement hexagonal décoré
5.5 Conclusions
Conclusions et perspectives
A Méthodes numériques et tests
A.1 Historique et objectifs
A.2 Fonctionnement général du code
A.2.a Organisation des données
A.2.b Schémas numériques de différentiation spatiale
A.2.c Schéma temporel
A.2.d Le dilemme compressible/incompressible
A.2.e Interface de traitement
A.3 Tests numériques
A.3.a Pas de temps et résolution spatiale
A.3.b Tests de performances et de parallélisation
A.4 Tests du code sur quelques problèmes de physique
A.4.a Ondes sonores
A.4.b Convection dans un polytrope sans champ magnétique .
A.4.c Instabilité de magnétoconvection – Overstabilité
A.5 Conclusions
B Oscillations alfvéniques d’une coquille sphérique
B.1 Introduction
B.1.a Avant-propos
B.1.b Oscillations magnétiques stellaires
B.2 Oscillations magnétiques en géométrie sphérique
B.2.a Modèle et équations
B.2.b Décomposition harmonique des équations MHD
B.2.c Symétries et classification des modes propres
B.2.d Les ondes d’Alfvén en résumé
B.3 Modes axisymétriques poloïdaux
B.3.a Modes stationnaires
B.3.b Modes oscillants
B.3.c Couches de cisaillement internes
B.4 Modes axisymétriques toroïdaux
B.4.a Spectre
B.4.b Singularité et couches limites
B.5 Solutions asymptotiques à faible diffusivité
B.6 Conclusions et extensions possibles
B.6.a Principaux résultats
B.6.b Perspectives
Bibliographie
Liste des publications
Liste des communications

projet fin d'etudeTélécharger le document complet

Télécharger aussi :

Laisser un commentaire

Votre adresse e-mail ne sera pas publiée. Les champs obligatoires sont indiqués avec *