No evidence of vorticity production from initially irrotational turbulent gravitational collapse

Résumé technique : Absence de preuve de production de vorticité à partir d'un effondrement gravitationnel turbulent initialement irrotationnel

Énoncé du problème
L'effondrement gravitationnel est largement considéré comme un générateur de turbulence, laquelle pourrait potentiellement amorcer une dynamo à petite échelle pour amplifier les champs magnétiques sous l'échelle de Jeans. Cependant, les dynamos à petite échelle nécessitent fondamentalement une turbulence vorticale (rotationnelle), alors que l'effondrement gravitationnel produit principalement des écoulements irrotationnels (potentiels). En l'absence de champs magnétiques ou de forces d'entraînement externes, et en employant une équation d'état barotrope (où les gradients de pression et de densité sont parallèles), le terme barocline ne peut générer de vorticité. Par conséquent, le seul mécanisme disponible pour produire de la vorticité à partir d'un écoulement initialement irrotationnel passe par l'interaction avec la viscosité. Des études antérieures ont suggéré que la compression induite par l'effondrement pourrait stimuler l'action de la dynamo, mais il reste incertain si l'effondrement lui-même génère la vorticité nécessaire ou si les effets de dynamo observés ne sont que la décroissance d'une turbulence initiale préexistante.

Méthodologie
Les auteurs réalisent des simulations numériques directes (DNS) d'effondrement gravitationnel à l'aide du code Pencil, employant des différences centrées du sixième ordre et un schéma de discrétisation temporelle du troisième ordre. L'étude utilise un domaine cubique périodique avec une équation d'état isentropique ($\gamma = 5/3$) et un profil de densité initial basé sur une sphère de Bonnor-Ebert isotherme modifiée.

Les choix méthodologiques clés incluent :

• Conditions initiales : Le champ de vitesse est construit pour être strictement irrotationnel dans l'espace de Fourier, avec des phases aléatoires et une pente spectrale spécifique. Aucun champ magnétique n'est inclus afin d'isoler les mécanismes de génération de vorticité.
• Viscosité : Au lieu de la viscosité de choc (qui peut masquer les résultats), les auteurs utilisent une viscosité cinématique constante ($\nu$).
• Résolution et paramètres : Les simulations sont exécutées à des résolutions allant de $512^3$ à $2048^3$ points de grille. Les auteurs font varier l'amplitude de la vitesse initiale ($u_{ini}$) et la viscosité ($\nu$) pour déterminer les limites numériquement fiables.
• Analyse : L'étude suit l'évolution de l'enstrophie ($\langle \omega^2 \rangle/2$) et décompose sa dérivée temporelle en trois termes : un terme de type dynamo ($T_{dyn}$), un terme de génération piloté par la viscosité et les gradients de densité ($T_{gen}$), et un terme de dissipation ($T_{dis}$). Les auteurs analysent également les spectres de vitesse, de densité et de vorticité, et examinent la mise à l'échelle de la vorticité avec l'amplitude du flux initial.

Résultats clés

1. Dépendance à la viscosité : La production de vorticité est considérée comme un effet secondaire piloté par l'amplification des divergences du flux initial via la viscosité. L'étude montre que le rapport entre la génération de vorticité et la dissipation ($R_{gen}$) n'est pas proportionnel à la viscosité ; au contraire, la génération de vorticité est plus significative durant la phase initiale précédant l'effondrement principal, et non comme une conséquence de l'écoulement de l'effondrement lui-même.
2. Mise à l'échelle avec les conditions initiales : La vorticité quadratique moyenne ($\omega_{rms}$) suit une mise à l'échelle en carré de l'amplitude de la vitesse initiale ($u_{ini}^2$), tandis que la divergence de vitesse quadratique moyenne suit une loi linéaire en $u_{ini}$. Cette mise à l'échelle persiste tout au long de la simulation, suggérant que l'effondrement gravitationnel ne pilote pas directement la production de vorticité. Au lieu de cela, l'effondrement amplifie les divergences du flux initial, que la viscosité convertit ensuite en vorticité.
3. Fiabilité numérique : Les auteurs établissent que les résultats deviennent peu fiables lorsque le nombre de Reynolds de maille ($Re_{mesh}$) dépasse des valeurs comprises entre 0,7 et 2,0, selon la simulation spécifique. Dans le régime fiable, le terme de type dynamo ($T_{dyn}$) reste subdominant, contribuant à moins de 20 % à l'évolution de l'enstrophie.
4. Analyse spectrale : Le spectre de vorticité présente une plage subinertielle en $k^4$, analogue au comportement des champs magnétiques en MHD, mais les auteurs attribuent cela à l'instabilité de Jeans pilotant des fluctuations de vitesse aléatoires plutôt qu'à un processus de dynamo. Le spectre de vitesse présente une plage inertielle abrupte sans effet de « goulot d'étranglement » (bottleneck).
5. Vorticité critique pour la dynamo : En réexaminant le nombre de Reynolds de vorticité magnétique critique ($Pr_M Re_\omega$) pour l'amorçage de la dynamo, les auteurs confirment une bifurcation autour d'une valeur de 300. Cela implique que les écoulements purement irrotationnels, qui ne peuvent soutenir une vorticité suffisante sans pilotage externe ou conditions initiales spécifiques, ne fonctionneront pas comme des dynamos.

Signification et affirmations
L'article conclut que l'effondrement gravitationnel, sous les conditions d'une équation d'état barotrope et d'une turbulence initialement irrotationnelle, ne produit pas efficacement de vorticité. La vorticité observée est un effet secondaire résultant de la conversion visqueuse des divergences de flux initiales, et non une conséquence directe de la dynamique de l'effondrement.

Les auteurs affirment que leurs conclusions remettent en question les affirmations précédentes concernant l'action de la dynamo à petite échelle induite par l'effondrement. Ils suggèrent que les observations antérieures d'amplification de champ magnétique dans les écoulements en effondrement étaient probablement alimentées par la décroissance de la turbulence initiale plutôt que par l'effondrement lui-même. L'étude souligne la nécessité de simulations numériques directes sans modélisation de sous-maille pour éviter la génération artificielle de vorticité. Ces travaux impliquent qu'atteindre une dynamo auto-entretenue dans de tels scénarios d'effondrement est difficile et que l'écoulement d'effondrement seul est insuffisant pour générer la turbulence vorticale requise pour l'action de la dynamo.