L'anomalie dissipative en turbulence compressible dépend du nombre de Mach ✛ ✎

Original: "Mach-number-dependent dissipative anomaly in isothermal compressible turbulence"
arXiv:2609.03945v1 · 2026-09-03 · arXiv non-exclusive · analyse indépendante · 4 min
En turbulence compressible isotherme, le taux de dissipation d'énergie atteint deux constantes différentes selon le nombre de Mach.

Contexte

La question de l'anomalie dissipative, c'est-à-dire l'indépendance du taux moyen de dissipation d'énergie cinétique par rapport à la viscosité à grands nombres de Reynolds, a longtemps été étudiée pour la turbulence incompressible. Pour les écoulements compressibles, les données sont rares et limitées au régime subsonique. Notre travail se démarque : les travaux les plus proches dans les archives concernent soit la turbulence incompressible, soit les champs magnétiques et la superfluidité, mais n'abordent pas le comportement de la dissipation lors de la transition à travers la barrière du son. Nous avons, pour la première fois, suivi systématiquement comment l'anomalie change lorsque le nombre de Mach varie de 0.1 à 10.

Méthodes

Nous avons réalisé une série de simulations numériques directes de turbulence compressible isotherme tridimensionnelle dans un domaine cubique avec des conditions aux limites périodiques. La turbulence était maintenue par une force stochastique purement solénoïdale à grande échelle. Le nombre de Mach variait de 0.1 à 10, et le nombre de Reynolds de 100 à 2400. Pour l'analyse, nous avons utilisé trois approches : les équations de bilan pour le taux de dissipation avec séparation des contributions des mouvements solénoïdaux et dilatationnels, l'extension du formalisme de Duchon–Robert aux écoulements compressibles pour les précurseurs locaux de l'anomalie dissipative, et l'analyse géométrique des régions de dissipation intense à l'aide de la méthode de comptage de boîtes et du formalisme multifractal.

Résultats

Le taux moyen de dissipation d'énergie cinétique normalisé, ainsi que ses composantes solénoïdale et dilatationnelle, atteignent des valeurs asymptotiques finies lorsque le nombre de Reynolds augmente. Cependant, ces valeurs ne sont pas universelles : les données se répartissent en deux branches. Une branche correspond aux régimes subsonique et transsonique (nombre de Mach quadratique moyen inférieur ou environ égal à 1), l'autre au régime fortement supersonique (nombre de Mach supérieur ou égal à 3). Pour la dissipation totale, la différence est la plus marquée. Pour les composantes solénoïdale et dilatationnelle, la dépendance au nombre de Reynolds est en outre sensible au choix de l'échelle intégrale et de la vitesse quadratique moyenne. Au nombre de Mach 10, les régions de dissipation maximale deviennent des feuilles bidimensionnelles, ce qui est confirmé par une dimension de comptage de boîtes d'environ 2 et par l'analyse multifractale. Au nombre de Mach 0.1, la dissipation intense est répartie sur un réseau de couches de cisaillement avec une dimension effective d'environ 2.5. Les mécanismes de génération changent également : aux faibles nombres de Mach, l'étirement tourbillonnaire domine ; aux nombres de Mach élevés, c'est l'interaction entre les mouvements tourbillonnaires et dilatationnels au sein des couches de choc.

Implications

Ce travail montre que l'anomalie dissipative en turbulence compressible n'est pas unique : elle dépend du nombre de Mach. Cela la distingue du cadre classique des écoulements incompressibles, où l'anomalie est unique. Contrairement aux travaux voisins, qui étudiaient soit la turbulence incompressible, soit des aspects particuliers de la compressibilité, nous avons, pour la première fois, relié les deux branches de l'anomalie à la réorganisation spatiale des structures dissipatives. Ce résultat est important pour construire des modèles de turbulence corrects en astrophysique et dans les applications d'ingénierie où la compressibilité est essentielle.

Développements futurs

Une suite naturelle consiste à inclure un forçage dilatationnel pour vérifier comment il affecte les branches de l'anomalie. On peut également sortir du cadre isotherme, par exemple en utilisant une équation d'état polytropique avec un exposant supérieur à un, ou en résolvant l'équation d'énergie complète. Cela permettra d'évaluer la robustesse des deux branches face aux changements de la thermodynamique. L'augmentation du nombre de Reynolds dans de futures simulations aidera à distinguer l'anomalie faible de l'anomalie forte.

Impact

Les résultats influenceront les modèles de turbulence en astrophysique, où les écoulements compressibles supersoniques sont courants, et les calculs d'ingénierie des flux à haute vitesse.

Prochaines étapes

Effectuer des simulations avec forçage dilatationnel et avec des nombres de Reynolds plus élevés afin de préciser la nature de l'anomalie.

Problèmes ouverts clés

En hydrodynamique, la question de la preuve rigoureuse de l'existence de l'anomalie dissipative pour les solutions des équations d'Euler et de Navier–Stokes reste ouverte. Nos données fournissent une base empirique pour les constructions mathématiques, en particulier pour les méthodes d'intégration convexe, et montrent que dans le cas compressible, l'anomalie peut ne pas être unique.

Le regard de la machine du savoir

Cette section s'adresse à l'auteur du travail et à ses collègues. Les travaux voisins ont été trouvés par une recherche sémantique dans notre archive, et non par correspondance de mots. Nous n'évaluons pas : nous montrons le champ autour.

Ce que fait ce travail

L'article examine l'anomalie dissipative dans la turbulence compressible isotherme forcée de manière solénoïdale, en s'appuyant sur des simulations numériques directes avec une résolution allant jusqu'à 2048³. Il est montré que le taux de dissipation d'énergie cinétique normalisé atteint des valeurs finies lorsque le nombre de Reynolds augmente, mais ces valeurs diffèrent entre les régimes subsonique/transsonique et fortement supersonique. L'analyse comprend les équations de bilan, une extension du formalisme de Duchon–Robert et une analyse géométrique des régions de dissipation intense.

Où est sa force

La force de ce travail réside dans le suivi systématique de la dépendance de l'anomalie dissipative en fonction du nombre de Mach, jusqu'au régime fortement supersonique, ce qui révèle deux branches de l'anomalie et les relie à la réorganisation des structures de dissipation, passant de tubes de vorticité à des couches de choc.

Vers où avancer

À essayer

  • Tester l'analyse des précurseurs locaux de dissipation sur les données de l'article sur la goutte dans la turbulence active, afin de vérifier si le lien entre la géométrie des structures dissipatives et l'anomalie persiste en présence de particules.
  • Comparer les spectres multifractals des régions de dissipation intense à différents nombres de Mach avec les spectres correspondants en turbulence incompressible, en utilisant les méthodes de l'article sur les mesures d'énergie, pour préciser l'universalité des propriétés géométriques.

Pourquoi cela peut compter

Comprendre la dépendance de l'anomalie dissipative au nombre de Mach est important pour les applications astrophysiques et d'ingénierie, où la turbulence est souvent supersonique, et cela peut aider à construire des modèles de dissipation sous-maille pour les simulations numériques.

Au nombre de Mach 10, les régions de dissipation maximale ont une dimension d'environ 2, c'est-à-dire qu'elles sont presque parfaitement plates, comme des feuilles de papier, alors que l'écoulement lui-même est tridimensionnel.

\frac{\langle \epsilon \rangle}{u_{\mathrm{rms}}^3/L} = \mathrm{const}
Le rapport de la dissipation moyenne au cube de la vitesse quadratique moyenne divisé par l'échelle intégrale tend vers une constante aux grands nombres de Reynolds.
M_{\mathrm{rms}} = \frac{u_{\mathrm{rms}}}{c_s}
Rapport de la vitesse quadratique moyenne à la vitesse du son.
Re = \frac{u_{\mathrm{rms}} L_f}{\nu}
Rapport des forces d'inertie aux forces visqueuses, où L_f est l'échelle caractéristique de forçage et ν la viscosité cinématique.

Chiffres clés

  • nombre de Mach maximal: 10
  • nombre de Reynolds maximal: 2400
  • résolution de la grille: jusqu'à 2048³ points
  • plage de nombres de Mach: 0.1–10
  • dimension effective des régions de dissipation à M=10: environ 2
Notions
nombre de Mach nombre de Reynolds