Asymptotics of spherical dynamos exhibiting a small-scale MAC balance

Résumé Technique : Asymptotique des dynamos sphériques présentant un équilibre MAC à petite échelle

Énoncé du problème
Comprendre la dynamique des noyaux planétaires terrestres, tels que le noyau externe de fer liquide de la Terre, nécessite l'extrapolation des modèles de dynamo numériques vers des régimes de paramètres physiques qui sont actuellement inaccessibles par calcul direct. Plus précisément, le nombre d'Ekman ($Ek$) dans le noyau terrestre est d'environ $10^{-15}$, alors que les simulations numériques directes (DNS) de dynamos sphériques atteignent typiquement seulement $Ek \sim 10^{-7}$. Un défi critique consiste à combler cet écart en identifiant les lois d'échelle asymptotiques correctes pour les quantités physiques lorsque $Ek \to 0$. Bien que la convection en rotation rapide non magnétique soit bien décrite par la théorie quasi-géostrophique (QG), la présence de champs magnétiques forts introduit la force de Lorentz, ce qui peut altérer l'équilibre des forces dominantes. Il reste incertain de savoir comment le comportement asymptotique des dynamos sphériques change lorsqu'un système passe d'un équilibre QG (Coriolis-pression) à un équilibre Magnétique-Archimède-Coriolis (MAC), où la force de Lorentz est comparable aux forces de Coriolis et de flottabilité.

Méthodologie
Les auteurs utilisent le code pseudo-spectral Rayleigh pour résoudre les équations de Navier-Stokes non dimensionnelles pour un fluide de Boussinesq auto-gravitant dans une géométrie de couche sphérique en rotation. L'étude emploie un balayage systématique des paramètres pour étudier la dépendance asymptotique vis-à-vis du nombre d'Ekman ($Ek$).

• Paramètres : Les simulations sont menées avec un nombre de Prandtl fixé ($Pr = 1$) et un nombre de Prandtl magnétique primaire ($Pm = 2$), avec un sous-ensemble variant $Pm$ de 0,2 à 5. Le nombre d'Ekman varie de $3 \times 10^{-6}$ à $10^{-4}$. Le nombre de Rayleigh ($Ra$) est fait varier jusqu'à 50 fois la valeur critique pour l'apparition de la convection non magnétique.
• Approche comparative : L'étude compare explicitement les cas de dynamo (où $B \neq 0$) à des cas équivalents de convection rotative non magnétique ($B = 0$) afin d'isoler les effets du champ magnétique sur les échelles asymptotiques.
• Analyse : Les auteurs analysent les forces fluctuantes intégrées par volume (Coriolis, pression, flottabilité, visqueuse, advection et Lorentz), les nombres de Reynolds, les spectres d'énergie et diverses échelles de longueur (échelles de Taylor, échelles spectrales de pic et échelles de dissipation). Ils se concentrent sur la composante radiale de l'équation de quantité de mouvement fluctuante pour déterminer les équilibres de force de premier ordre.

Contributions clés et résultats

1. Équilibre des forces et régimes d'échelle :

   • Cas non magnétiques : Confirment l'échelle QG standard où les forces de Coriolis et de pression dominent au premier ordre. La force de flottabilité fluctuante suit une échelle de $O(Ek^{-1})$, ce qui est asymptotiquement plus petit que la force de Coriolis ($O(Ek^{-4/3})$).
   • Cas de dynamo (Équilibre MAC) : Pour $Pm=2$ et un $Ek$ suffisamment bas, les simulations présentent un équilibre MAC où les forces de Coriolis, de pression, de flottabilité et de Lorentz entrent toutes au premier ordre.
   • Échelle de la flottabilité : Une découverte majeure est que dans le régime MAC, la température fluctuante devient de l'ordre de l'unité ($O(1)$) et indépendante de $Ek$. Par conséquent, la force de flottabilité fluctuante suit une échelle de $O(Ek^{-4/3})$, correspondant à la force de Coriolis. Cela contraste nettement avec les cas non magnétiques où la flottabilité est asymptotiquement plus faible. Ce comportement est attribué à la dominance de la dissipation ohmique sur la dissipation visqueuse dans le bilan énergétique, permettant des fluctuations de température plus importantes.

2. Échelle de la vitesse et du nombre de Reynolds :

   • Malgré le passage à un équilibre de force MAC, le nombre de Reynolds convectif ($Re_c$) suit la prédiction QG de $O(Ek^{-1/3})$.
   • L'échelle de vitesse du flux est cohérente avec un équilibre entre les forces de Coriolis et de flottabilité (échelle CIA), produisant $Re_c \sim Ra Ek / Pr$. Cela suggère que même dans un régime MAC, l'échelle de vitesse reste sans diffusion et identique à la théorie QG.

3. Échelles de longueur et spectres :

   • Champ de vitesse : Les spectres d'énergie cinétique révèlent un double comportement d'échelle. Les structures à petite échelle (degrés harmoniques sphériques élevés) suivent l'échelle de longueur prédite par la QG de $O(Ek^{1/3})$. Cependant, le pic du spectre d'énergie cinétique et les structures à grande échelle présentent une échelle de $O(Ek^{1/6})$.
   • Champ magnétique : Le spectre d'énergie magnétique et les échelles de longueur magnétique (microéchelle de Taylor et échelle de la force de Lorentz) sont uniformément caractérisés par une échelle de $O(Ek^{1/6})$. Ceci est cohérent avec une échelle du nombre de Reynolds magnétique de $Rm^{-1/2}$.
   • Rôle de la viscosité : Il est démontré que la présence de l'échelle de longueur $O(Ek^{1/3})$ dans le champ de vitesse provient de la nature asymptotique des termes d'advection et d'inertie, et non nécessairement de la dominance de la force visqueuse.

4. Dissipation :

   • Dans le régime MAC, la dissipation ohmique suit une échelle de $O(Ek^{-5/3})$, ce qui est nettement supérieur à l'échelle de dissipation visqueuse de $O(Ek^{-4/3})$ observée dans les cas non magnétiques. Cela confirme que les simulations de dynamo possèdent une dissipation asymptotique totale plus élevée que leurs homologues non magnétiques à bas $Ek$.

5. Rôle de l'advection et de l'inertie :

   • L'étude montre que les termes d'advection et d'inertie restent comparables à la force visqueuse et suivent une échelle de $O(Ek^{-1})$. Cela implique que les ondes de Rossby convectives non linéaires, caractérisées par une échelle de temps de $O(Ek^{2/3})$, jouent un rôle significatif dans la dynamique, même au sein d'un régime MAC, remettant en question l'idée que l'advection est négligeable dans les modèles de géodynamo.

Signification et affirmations
L'article affirme fournir une mesure quantitative du comportement asymptotique des dynamos sphériques, offrant un cadre pour extrapoler les résultats numériques aux conditions des noyaux planétaires. Les affirmations clés sont :

• Ordres asymptotiques distincts : La différence la plus notable entre les dynamos MAC et la convection non magnétique est l'ordre asymptotique de la force de flottabilité. Dans les dynamos MAC, la flottabilité est de premier ordre ($O(Ek^{-4/3})$), tandis que dans les cas non magnétiques, elle est de second ordre ($O(Ek^{-1})$).
• Robustesse de l'échelle de vitesse : L'échelle de vitesse ($O(Ek^{-1/3})$) est robuste à travers les régimes magnétiques et non magnétiques, suggérant que les propriétés de transport fondamentales du flux sont moins sensibles à l'intensité du champ magnétique que l'équilibre des forces lui-même.
• Clarification des échelles de longueur : Les auteurs soutiennent que l'observation d'échelles de $O(Ek^{1/3})$ dans les simulations n'implique pas intrinsèquement que la viscosité est la force dominante ; c'est une conséquence de la structure asymptotique des termes d'advection et d'inertie.
• Limites des simulations actuelles : L'étude note que bien que la force de Lorentz suive une échelle plus forte ($O(Ek^{-3/2})$) que la force de Coriolis dans la plage de paramètres actuelle, il n'est pas certain si cette tendance persiste vers l'extrême $Ek$ du noyau terrestre ou si l'échelle doit changer. Les résultats suggèrent que les investigations futures doivent explorer la limite double de $Ek \to 0$ et $Pm \to 0$ pour résoudre pleinement l'équilibre des forces dans les noyaux planétaires.

Les auteurs concluent que leurs découvertes aident à expliquer le comportement des forces dans les études antérieures et fournissent une base rigoureuse pour distinguer les régimes QG et MAC dans les modèles de dynamo numériques.