Détermination des coefficients de transport turbulent et analyse des cycles magnétiques produits dans un modèle dynamo en champ moyen avec et sans rétroaction magnétique

Avec les récents développements obtenus grâce aux modèles globaux magnétohydrodynamiques en trois dimensions de la convection solaire, il est désormais possible de simuler des champs magnétiques structurés à grande échelle et présentant des inversions de polarité bien synchronisées dans chaque hémisphère. Ces modèles qui n'incluent, pour la plupart, aucune modélisation de la surface du Soleil génèrent donc leur dynamo avec l'action de la force électromotrice turbulente (FEM) et de la rotation différentielle uniquement. À partir de cette FEM, différentes techniques peuvent être utilisées pour extraire les coefficients de transport turbulent. Notamment, différents auteurs ont obtenu un tenseur alpha (coefficient du premier ordre) dont les 9 composantes présentent des amplitudes du même ordre, remettant en doute l'approximation faite dans le cas des modèles dynamo de type alphaOmega qui ne tient en compte qu’une de composante du tenseur. À partir d'un code d'analyse par décomposition en valeurs singulières pour évaluer les coefficients du tenseur alpha, nous avons généralisé la procédure pour extraire 18 des composantes du tenseur de deuxième ordre (tenseur beta). Les tenseurs alpha et beta obtenus par cette nouvelle procédure tel qu'appliquée aux sorties du modèle global EULAG-MHD, sont similaires aux tenseurs alpha et beta équivalant obtenus en utilisant l'approximation « Second Order Correlation Approximation ».  À l'aide des coefficients de transport turbulent du premier ordre introduit dans un modèle dynamo en champ moyen, nous avons ensuite étudié certaines solutions magnétiques présentant des doubles dynamos. Cette analyse avait pour but de comparer les résultats obtenus par ce modèle simplifié caractérisé par la FEM provenant de EULAG-MHD aux résultats de EULAG-MHD directement. Cette preuve de concept nous a permis de démontrer que l'oscillation observée dans le champ magnétique en surface de EULAG-MHD pouvait provenir de l'action inductive d'une seconde dynamo. Une oscillation biennale est également observée dans plusieurs indices d'activité solaire dont l'origine n'est toujours pas établie. Il est évident que les deux modèles décrits ci-haut et le Soleil opèrent dans des régimes physiques différents. Toutefois, malgré leurs différences, le fait qu'il soit relativement facile de produire une seconde dynamo dans EULAG-MHD et dans le modèle en champ moyen indique que l'action inductive de la FEM peut facilement générer deux dynamos.  Finalement, dans le but d'étudier les périodes de grands minima, phénomène encore non reproduit par les modèles globaux, nous avons ajouté une rétroaction magnétique sur l'écoulement azimutal au modèle dynamo cinématique en champ moyen décrit ci-haut. En analysant les solutions de ce modèle dynamo de type alpha2Omega non cinématique, nous avons pu reproduire la tendance observée jusqu'ici uniquement dans les modèles de type alphaOmega selon laquelle le nombre de nombre de Prandtl magnétique contrôle le rapport des périodes générées. De plus, en analysant une solution sur 50 000 ans présentant des périodes de grands minima et maxima non périodiques, nous avons obtenu une distribution de temps de séparation des grands minima presque exponentielle, caractéristique observée dans les reconstructions de l'activité solaire. La rotation différentielle associée à ces périodes de grands minima présente un niveau de fluctuation de 1% par rapport au profil moyen. Ce niveau de fluctuation est d'ailleurs comparable avec les reconstructions historiques de la rotation différentielle en surface obtenues lors du grand minimum de Maunder. The recent developments achieved by tri-dimensionals magnetohydrodynamic (3D-MHD) global simulations of solar convection allow us to generate an organized large-scale magnetic fields with well-synchronized hemispheric polarity reversal.  Because the vast majority of these simulations do not include a modelization of the Sun's surface layer, the generation of their dynamo is thus solely due to the action of the turbulent electromotive force (EMF) in conjunction with differential rotation. From this EMF, different methods can be used to extract the turbulent transport coefficients. In particular, various authors found a full 9 component alpha-tensor (first order coefficients) where all the components are of the same order of magnitude. This finding calls into question the alphaOmega approximation made by the vast majority of mean field dynamo models. We generalized a first order (alpha-tensor) singular value decomposition (SVD) analysis procedure to extract the 18 additional components of the  second order tensor (beta-tensor). The alpha and beta tensors obtained by this new procedure as applied to the EULAG-MHD outputs, are similar to the equivalent alpha and beta tensors obtained using the second order correlation approximation (SOCA).  By introducing the first order turbulent transport coefficients in a mean field dynamo model, we study the magnetic solutions where double dynamo modes were observed. This analysis allows us to compare the mean field dynamo solutions produced with the EMF, as extracted from EULAG-MHD, with the real magnetic output of EULAG-MHD. This proof of concept demonstrated that the quasi-biennal oscillation observed in the surface toroidal magnetic field in EULAG-MHD can be produced by the inductive action of a secondary dynamo. A similar quasi-biennal oscillation signal is also observed in multiple proxies of the solar activity whose origin is still not confirmed.  Although the physical set of properties under which the two numerical models described above operate are different from the Sun, the fact that both models can reproduce a secondary dynamo shows us that the inductive action of the EMF can easily produce two dynamos.  Finally, in order to study epochs of grand minima that still cannot be reproduced in global 3D-MHD simulations of convection, we added a magnetic feedback on the mean azimutal flow in our kinematic mean field model. This non-kinematic alpha2Omega model was able to reproduce the tendency of the Prandtl number (Pm) to control the ratio of the modulation period. More specifically, we found an inverse relation between Pm and the ratio of the main magnetic cycle period to the grand minima occurrence period. Moreover, by analyzing a simulation of a length of 50,000 years, where aperiodic periods of grand minima and maxima are observed, we found a waiting time distribution (WTD) of the grand minima close to an exponential, a characteristic also observed in the reconstruction of the solar activity. Finally, the level of fluctuation in the surface differential rotation associated with epochs of grand minima is ~1%. This level of fluctuation was also observed in historical reconstructions of the surface differential rotation during the Maunder minimum.

Keyphrases: Effet Alpha, Diffusivité turbulente, Dynamo solaire, Modèle dynamo en champ moyen, Coefficient de transport turbulent, Grand minimum, Alpha effect, Turbulent diffusivity, Solar dynamo, Mean field dynamo model, Turbulent transport coefficients, Grand minima