Generalization of the viscous stress tensor to the case of non-small gradients of hydrodynamic velocity: a path to numerical modeling of turbulence non-locality

Ce papier généralise la méthode de Chapman-Enskog pour dériver une représentation intégrale du tenseur des contraintes visqueuses pour de grands gradients de vitesse, permettant la modélisation numérique de la non-localité de la turbulence et de phénomènes tels que les discontinuités tangentielles que les formulations standard de Navier-Stokes peinent à capturer.

L'essentiel

Le Grand Problème : Le « Smoothie » contre la « Tempête »

Imaginez que vous essayez de prédire comment un fluide (comme l'eau ou l'air) se déplace. Depuis plus d'un siècle, les scientifiques utilisent un célèbre ensemble de règles appelé les équations de Navier-Stokes. Ces règles reposent sur un ingrédient spécifique appelé le tenseur des contraintes visqueuses.

Considérez ce tenseur comme un « calculateur de frottement ». Dans la version standard, il suppose que si vous poussez un fluide, la résistance (frottement) qu'il ressent dépend uniquement de la vitesse du fluide immédiatement à côté du point que vous observez. C'est comme supposer que si vous remuez une tasse de café, la résistance que vous ressentez est déterminée uniquement par les molécules de café touchant votre cuillère.

Le Défaut :
L'auteur, A.B. Kukushkin, souligne que ce « calculateur de frottement » standard s'effondre lorsque les choses deviennent chaotiques, comme dans une tempête turbulente ou lorsque deux courants de fluide glissent l'un à côté de l'autre à grande vitesse (discontinuités tangentielles).

• L'Analogie : Imaginez une foule de personnes marchant dans un couloir. Le modèle standard suppose que chacun ne heurte que la personne immédiatement à côté de lui. Mais dans une vraie foule (turbulence), une personne peut être poussée par quelqu'un situé trois rangées en arrière, ou une vague de mouvement peut traverser toute la pièce. Le modèle standard ignore ces interactions « à longue distance ».
• Le Paradoxe : Les mathématiques standard conduisent également à un résultat étrange : elles suggèrent que si les particules entrent en collision plus souvent (comme dans un brouillard épais), le fluide devrait en réalité s'écouler plus facilement (viscosité plus faible). Cela va à l'encontre de notre intuition.

La Solution : Examiner l'Image Globale

Kukushkin propose une nouvelle façon de calculer ce frottement. Au lieu de ne regarder que le voisinage immédiat, sa nouvelle formule examine l'ensemble de l'histoire et de la localisation du mouvement du fluide.

La Nouvelle Approche :

1. Abandonner la Règle des « Petits Pas » : L'ancienne mathématique (méthode de Chapman-Enskog) ne fonctionne que si le fluide change très lentement et doucement. Kukushkin élimine cette règle. Il autorise des changements soudains et brutaux de vitesse, qui se produisent dans la vraie turbulence.
2. L'Analogie du « Messager » : Au lieu de ne regarder le fluide qu'à un seul endroit, imaginez que le fluide est rempli de petits « messagers » (tourbillons ou perturbations).
   • Dans l'ancien modèle, un messager ne parle qu'à son voisin.
   • Dans le nouveau modèle de Kukushkin, un messager naît à un endroit, traverse la pièce et livre son message à un endroit éloigné avant de s'arrêter.
3. La Formule Intégrale : Les nouvelles mathématiques sont une intégrale (une somme sur une grande zone). Elle calcule la contrainte (frottement) à un point spécifique en additionnant les effets de tous ces messagers voyageant de partout ailleurs dans le fluide vers ce point.

Pourquoi Cela Importe

1. Réparer le « Paradoxe » :
En permettant à ces messagers de parcourir de longues distances, la nouvelle formule corrige le paradoxe étrange concernant la viscosité. Elle explique pourquoi les fluides se comportent comme ils le font même lorsque les particules entrent fréquemment en collision. Les « longs voyages » des messagers rendent compte de la résistance d'une manière que l'ancien modèle de « petits pas » ne pouvait pas.

2. Relier au Chaos Réel (Loi de Richardson) :
L'article mentionne une observation célèbre appelée la loi en $t^3$ de Richardson.

• L'Analogie : Si vous lâchez deux feuilles dans une rivière turbulente, le modèle standard prédit qu'elles dériveront lentement l'une de l'autre (comme en $t^2$). Mais en réalité, elles s'éloignent beaucoup plus vite (comme en $t^3$).
• Le Lien : Ce nouveau modèle « à longue distance » explique naturellement pourquoi les particules se séparent si rapidement. Les messagers voyagent loin et vite, transportant la perturbation à travers le fluide, ce qui correspond à l'observation réelle de la propagation de la turbulence.

3. Un Pont vers de Meilleures Simulations Informatiques :
Actuellement, les simulations informatiques de la turbulence doivent souvent utiliser des « trucs » ou des nombres inventés parce que les mathématiques standard échouent aux bords nets (comme là où une aile se sépare de l'air).

• La nouvelle formule de Kukushkin fournit un pont mathématique. Elle transforme les « trucs » en un calcul rigoureux basé sur les premiers principes. Elle permet aux ordinateurs de modéliser la turbulence en additionnant ces interactions à longue distance, plutôt que de simplement deviner.

Résumé en Bref

L'article soutient que l'ancienne façon de calculer le frottement des fluides est comme essayer de comprendre une conversation en n'écoutant que la personne debout à côté de vous. Elle manque la vue d'ensemble.

Kukushkin a écrit un nouveau code de règles qui écoute toute la pièce. En tenant compte de la façon dont les perturbations voyagent à travers tout le fluide (non-localité), ces nouvelles mathématiques :

• Réparent un paradoxe logique concernant l'épaisseur ou la fluidité d'un fluide.
• Expliquent pourquoi les particules dans une tempête s'éloignent si rapidement.
• Offrent une voie pour que les ordinateurs simulent des écoulements complexes et chaotiques (comme le vent autour d'un avion ou l'eau dans un tuyau) beaucoup plus précisément, sans avoir besoin de se fier à des suppositions.

L'auteur note que cette même logique pourrait éventuellement être appliquée au transfert de chaleur et même à la physique des plasmas, mais le principal accomplissement ici est la réécriture des règles du frottement des fluides pour gérer la réalité « désordonnée » de la turbulence.

Résumé technique

Résumé technique : Généralisation du tenseur des contraintes visqueuses pour des gradients non faibles

Énoncé du problème
L'article aborde les limitations fondamentales de la description standard de la turbulence hydrodynamique par les équations de Navier-Stokes, concernant spécifiquement le tenseur des contraintes visqueuses ($\sigma_{ik}$). Le tenseur conventionnel, dérivé par la méthode de Chapman-Enskog à partir de l'équation de Boltzmann, repose sur l'hypothèse de gradients faibles de la vitesse hydrodynamique. L'auteur soutient que cette hypothèse rend le tenseur standard inadéquat pour décrire des phénomènes d'écoulement complexes tels que les discontinuités tangentielles, les écoulements décollés et la non-localité inhérente à la turbulence.

De plus, l'article met en évidence un paradoxe théorique dans la dérivation standard : le coefficient de viscosité dynamique ($\eta$) s'avère inversement proportionnel à la section efficace de collision des particules ($\sigma_{coll}$), contrairement aux attentes intuitives. L'auteur attribue ce paradoxe à la condition restrictive de « gradient faible » de la méthode de Chapman-Enskog, qui ne devient de plus en plus réalisable qu'à mesure que la section efficace de collision augmente, masquant ainsi le comportement des tourbillons localisés de longue durée.

Les modèles de turbulence existants (par exemple, Spalart-Allmaras, RANS) ont souvent recours à des corrections phénoménologiques ou à des équations différentielles locales qui échouent à capturer la nature non locale de la turbulence, comme en témoigne empiriquement la loi en $t^3$ de Richardson pour les corrélations de paires dans les milieux turbulents. Bien que des modèles de transport non local (tels que les vols/promenades de Lévy) aient été proposés, ils manquent souvent d'une dérivation rigoureuse à partir de premiers principes au sein de l'hydrodynamique des gaz.

Méthodologie
L'auteur propose une généralisation du tenseur des contraintes visqueuses en assouplissant la condition de gradients de vitesse faibles dans le cadre de Chapman-Enskog. La méthodologie procède comme suit :

1. Fondement cinétique : La dérivation commence par l'équation cinétique de Boltzmann pour un gaz à une seule espèce. La fonction de distribution $f(\mathbf{p}, \mathbf{r}, t)$ est décomposée en une partie d'équilibre thermodynamique local (distribution de Maxwell $f^{(0)}$) et un terme de correction $\varphi$.
2. Assouplissement des contraintes de gradient : Contrairement à l'approche standard où les dérivées spatiales du terme de correction sont négligées, l'auteur conserve les dérivées de $\varphi$ du côté gauche de l'équation de transport. Cela reconnaît que dans les écoulements cisaillés et les tourbillons localisés, les gradients du terme de correction peuvent être significativement plus grands que ceux de la distribution d'équilibre lisse.
3. Formulation de l'équation de transport : Une intégrale de collision simplifiée (de type BGK) est utilisée, conduisant à une équation de transfert cinétique pour $\varphi$. Cette équation présente un terme source dépendant du gradient de la vitesse massique moyenne et un terme puits représentant l'absorption des porteurs (libre parcours inverse).
4. Solution intégrale : L'équation de transport est résolue en supposant que le mouvement des porteurs est significativement plus rapide que le mouvement du milieu global. La solution pour $\varphi$ est exprimée comme une intégrale sur les trajectoires des porteurs, depuis la limite du milieu jusqu'au point d'observation, incorporant un facteur de décroissance exponentielle basé sur la longueur du libre parcours.
5. Dérivation du tenseur des contraintes : Le tenseur des contraintes visqueuses est calculé en intégrant le produit du terme de correction, de la distribution d'équilibre et du flux de quantité de mouvement sur l'espace des phases. En transformant l'intégration sur les trajectoires en une intégration volumique sur le milieu, l'auteur dérive une représentation intégrale non locale du tenseur des contraintes.

Contributions et résultats clés
Le résultat principal est la dérivation d'un tenseur des contraintes visqueuses généralisé $\sigma_{ik}$ qui est une intégrale sur les coordonnées spatiales plutôt qu'un opérateur différentiel local.

• Représentation non locale : La nouvelle expression de $\sigma_{ik}$ (Équation 2.13) dépend explicitement des gradients de vitesse en tous les points du milieu, pondérés par un noyau impliquant la distance entre les points et la longueur du libre parcours. Cette structure reflète les théories de transport non local, telles que le modèle de Biberman-Holstein pour le transfert de radiation résonnante.
• Récupération de la théorie standard : Dans la limite des faibles longueurs de parcours (faibles gradients), l'expression intégrale se réduit au tenseur des contraintes visqueuses local standard (Équation 1.2), avec le coefficient de viscosité $\eta = P\tau$. Cela assure la cohérence avec la théorie établie dans le régime approprié.
• Résolution du paradoxe : La forme généralisée résout la dépendance inverse de la viscosité par rapport à la section efficace de collision. L'auteur postule que le paradoxe découle du fait que la dérivation standard suppose des gradients faibles, une condition incompatible avec les tourbillons localisés de longue durée qui agissent comme porteurs de perturbations dans les écoulements turbulents.
• Lois de conservation : La dérivation démontre que la correction à la fonction de distribution préserve le nombre total de particules dans le milieu, représentant un échange dynamique entre le milieu principal et l'ensemble des porteurs de perturbations.

Signification et affirmations
L'article affirme que ce travail constitue une étape vers une théorie non locale de la turbulence hydrodynamique dérivée de premiers principes.

• Modélisation des discontinuités : Le tenseur généralisé fournit un opérateur mathématique capable de décrire les discontinuités tangentielles et les écoulements décollés où les dérivées de vitesse sont non bornées, comblant une lacune de longue date dans la modélisation numérique.
• Lien avec la loi de Richardson : La structure non locale du tenseur dérivé est présentée comme une fondation nécessaire pour la modélisation numérique de la non-localité de la turbulence, spécifiquement pour reproduire la loi empirique en $t^3$ de Richardson pour les corrélations de paires.
• Validation de l'approche de Kolmogorov : L'auteur suggère que cette généralisation peut expliquer le succès du raisonnement dimensionnel de Kolmogorov dans la prédiction des phénomènes de turbulence, malgré le fait que ces résultats aient historiquement été difficiles à dériver directement de l'équation standard de Navier-Stokes.
• Applicabilité plus large : La méthodologie est notée comme potentiellement applicable au transfert de chaleur dans les gaz et les liquides, ainsi qu'à des milieux plus complexes comme le plasma, offrant une approche unifiée des phénomènes de transport non local.

L'article conclut que, bien que l'expression dérivée constitue une avancée théorique significative, la démonstration complète de sa capacité à reproduire la loi de Richardson nécessite une modélisation numérique adéquate ultérieure.