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

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: Alpha effect, Turbulent diffusivity, Solar dynamo, Mean field dynamo model, Turbulent transport coefficients, Grand minima