Dans ce travail, nous avons abordé le phénomène de la Convection Thermique Turbulente sous des angles nouveaux et innovants. <br /><br />Le premier système que nous abordons expérimentalement est un Canal Vertical dit Infini au sein duquel règne un gradient de température moyen constant. Le flux de masse dans ce canal est nul. L'image qui se dégage de nos mesures est celle d'un écoulement majoritairement inertiel, où les coefficients dissipatifs (la viscosité en l'occurrence) n'interviennent que pour fixer une longueur de cohérence L. Cette longueur est celle sur laquelle les panaches thermiques peuvent êtres considérés comme en « chute libre ». Le transport horizontal (d'impulsion et de chaleur) est entièrement dû aux fluctuations. La « longueur de mélange » associée est petite devant la largeur du canal. Par contre, le transport de chaleur vertical est dû à des structures cohérentes, les panaches. <br />Ces panaches, nous les retrouvons dans une étude Lagrangienne de l'écoulement au centre d'une cellule de Rayleigh Bénard. La sonde est une sphère de 2 cm de diamètre qui a la même densité que le fluide que nous utilisons, équipée de thermomètres et d'un émetteur radio. Elle est transportée par les panaches, ce qui nous permet une étude statistique de ceux-ci.

Le confinement magnétique dans les tokamaks est à l'heure actuelle la voie la plus avancée pour produire de l'énergie par fusion thermonucléaire. Des études théoriques et expérimentales ont montré que la génération de rotation permet d'en augmenter les performances par la réduction du transport turbulent à l'œuvre dans les plasmas de tokamaks. L'influence de la rotation sur les flux turbulents de chaleur et de particules ainsi que le transport du moment angulaire sont étudiés par simulation numérique dans le cadre du code gyro-cinétique, quasi-linéaire QuaLiKiz. A cette occasion, le code QuaLiKiz est modifié pour prendre en compte la rotation du plasma et calculer le flux de moment angulaire. Il est montré que le cadre de travail de QuaLiKiz permet de calculer le flux de moment angulaire y compris le stress résiduel induit par le cisaillement du champ électrique radial ainsi que l'effet de la rotation sur les flux de chaleur et de particules. Les approximations majeures du formalisme utilisé, en particulier la représentation de ballonnement à son ordre le plus bas et l'utilisation de fonctions propres analytiques calculées dans la limite hydrodynamiques, sont analysées en détail et leur validité vérifiée. La construction des flux quasi-linéaires est ensuite détaillée et le flux quasi-linéaire de moment angulaire dérivé. Les différentes contributions au flux turbulent de moment angulaire sont étudiées et comparées avec succès à la fois aux données de simulations gyro-cinétiques non-linéaires ainsi qu'aux données expérimentales.

L'approche de simulation numérique directe est utilisée pour la simulation d'un écoulement en canal pleinement turbulent afin d'étudier l'influence des grandes échelles de l'écoulement ainsi que la dynamique du transport des contraintes de Reynolds et de la vorticité. Les simulations sont réalisées sur un domaine de calcul de grande taille afin de pouvoir capturer l'intégralité des grandes structures de l'écoulement, et portent sur une gamme relativement étendue de nombres de Reynolds (Reτ =180, 395, 590 et 1100) allant des écoulements faiblement turbulents à des écoulements modérément turbulents. L'invariance remarquable des fluctuations de vorticité normale est expliquée à travers une analyse spectrale de la vorticité. L'étude des différents termes du transport de l'intensité turbulente de la vorticité révèle par ailleurs que le pic de production de la vorticité transverse est situé à proximité immédiate de la paroi et pourrait ouvrir la voie à des stratégies de réduction de la traînée basées sur la réduction de la vorticité transverse. Le transport des contraintes de Reynolds dans la couche interne et dans la couche de recouvrement est également étudié. A proximité des parois, la dépendance des termes de transport avec le nombre de Reynolds dans les échelles internes montre que ces dernières ne suffisent pas à caractériser la dynamique des contraintes de Reynolds dans cette zone. Cette insuffisance des échelles internes nous a amenés à nous intéresser plus particulièrement au processus de production à travers les statistiques de la production conditionnées par le passage par niveau des fluctuations de la vitesse normale ou longitudinale. Cette étude nous a permis d'identifier les fluctuations qui contribuent le plus à la production et celles qui sont à l'origine de la dépendance avec le nombre de Reynolds. / The direct numerical simulations of a fully turbulent channel flow are investigated to study the large scales effects on the flow quantities such as the Reynolds stresses and vorticity transport processes. Large computational domains are used so as to cover the largest scales of the flow. The simulations are performed in a wide range of Reynolds numbers (Reτ=180, 395, 590 and 1100) going from weakly to moderately high Reynolds number turbulent flows. The invariance of the wall-normal vorticity fluctuations scaled in wall variables in the inner layer versus the Reynolds number is analyzed using a spectral analysis. The vorticity transport equations are investigated in detail, presumably for the first time. The transport mechanism of the Reynolds shear stresses are subsequently analyzed in the inner layer and the overlapping zone. In the wall layer, different terms of the Reynolds stresses transport expressed in inner scales depend on the Reynolds number. This scaling failure lead us to focus on the statistics of the production when the streamwise or normal velocity fluctuations cross a given level, through the conditional Palm statistics. The main aim is to identify those amplitudes of the fluctuations that contribute more to the production and those which are responsible for the production Reynolds dependence.

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.