Sub-Kolmogorov Intermittency and Multifractal Dissipation in Multiphase Turbulence

Par des simulations numériques directes, cette étude révèle que dans la turbulence multiphasique, la rupture et la coalescence des interfaces induisent une organisation multifractale distincte de la dissipation, provoquant une extension des événements de dissipation d'énergie intense profondément dans la gamme sous-Kolmogorov et élargissant considérablement la coupure dissipative locale par rapport à la turbulence monophasique.

L'essentiel

Imaginez une casserole d'eau bouillante sur une cuisinière. Si vous n'avez que de l'eau (monophasique), les bulles et les tourbillons sont chaotiques, mais ils suivent un schéma relativement prévisible de taille et d'énergie. Maintenant, imaginez que vous ajoutez de l'huile à cette eau et que vous remuez vigoureusement. Vous obtenez un mélange désordonné de gouttelettes, de courants et de bulles qui se forment, fusionnent et se divisent constamment. C'est la turbulence multiphasique.

Cet article étudie ce qui se passe aux niveaux les plus infimes et les plus invisibles de ce mélange chaotique. Les chercheurs voulaient comprendre pourquoi les « plus petits tourbillons » d'un mélange de liquides se comportent si différemment — et de manière plus violente — que dans un liquide unique.

Voici l'histoire de leur découverte, décomposée en concepts simples :

1. Le « filet de sécurité » qui n'existe pas

En physique des fluides classique, il existe un « filet de sécurité » théorique appelé échelle de Kolmogorov. Voyez cela comme la taille minimale qu'un tourbillon peut atteindre avant que la viscosité naturelle du fluide ne l'atténue et n'en tue l'énergie. Dans un liquide unique, l'énergie s'arrête là.

Cependant, les chercheurs ont découvert que dans un mélange de liquides (comme l'huile et l'eau), ce filet de sécurité est brisé.

• L'analogie : Imaginez un acrobate de trapèze (l'énergie) qui se balance. Dans un liquide unique, il s'arrête de se balancer à une certaine hauteur. Dans un mélange de liquides, l'acrobate continue de se balancer bien plus bas, profondément dans une zone où la physique dit qu'il aurait dû s'arrêter.
• La découverte : La « coupure dissipative » (le point où l'énergie meurt) ne s'arrête pas simplement à la limite habituelle ; elle s'étend profondément dans une plage « sub-Kolmogorov ». Les fluctuations d'énergie deviennent beaucoup plus intenses et extrêmes que ce que l'on attendait.

2. Les coupables : Rupture et fusion

Pourquoi cela se produit-il ? L'article identifie les « scènes de crime » spécifiques où cette énergie extrême est générée.

• L'analogie : Pensez à une foule de personnes se déplaçant de manière aléatoire. Si deux personnes se cognent et fusionnent, ou si un groupe se divise, cela provoque un choc soudain et chaotique.
• La découverte : Les bouffées d'énergie les plus intenses à petite échelle se produisent spécifiquement aux interfaces où les liquides se rencontrent. Plus précisément, elles se produisent lorsque les gouttelettes se brisent (rupture) ou s'entrechoquent et fusionnent (coalescence).
• Ces événements créent des courbes abruptes et des changements de vitesse soudains que le fluide ne peut pas lisser facilement, forçant l'énergie à descendre de plus en plus profondément dans le domaine microscopique.

3. La géométrie « fractale » du chaos

Les chercheurs ont utilisé un outil mathématique appelé analyse multifractale.

• L'analogie : Imaginez regarder un littoral. De loin, il ressemble à une ligne. De plus près, il est dentelé. De plus près encore, il est rempli de baies et de rochers. Un « fractal » est une forme qui semble complexe à chaque niveau de zoom.
• La découverte : Dans un liquide unique, la « rugosité » de la distribution de l'énergie est assez constante. Mais dans un mélange de liquides, la géométrie du chaos change complètement aux plus petites échelles.
  • La « rugosité » devient beaucoup plus extrême.
  • Les événements énergétiques les plus violents ne sont pas répartis uniformément ; ils sont concentrés sur des structures filamenteuses très fines (comme le col d'une gouttelette juste avant qu'elle ne se rompe).
  • L'article décrit ces événements intenses comme étant supportés par des « structures éparses », ce qui signifie qu'ils sont rares, isolés et incroyablement tranchants, plutôt qu'un brouillard général de turbulence.

4. La machine de prédiction

Les chercheurs n'ont pas seulement observé cela ; ils ont prouvé qu'ils pouvaient le prédire.

• Ils ont utilisé la « forme » mathématique du chaos (le spectre de singularité) pour prédire exactement la fréquence de ces événements minuscules et extrêmes.
• Le résultat : Lorsqu'ils ont examiné les zones « proches » et « sub-Kolmogorov » (les échelles profondes et minuscules), leurs prédictions correspondaient parfaitement aux simulations informatiques. Cela confirme que ce comportement étrange et extrême est le résultat direct de la rupture et de la fusion des interfaces liquides.

L'essentiel

L'article conclut que lorsque vous mélangez deux liquides, la turbulence ne devient pas seulement « un peu plus désordonnée ». L'acte de rupture et de fusion des gouttelettes réécrit fondamentalement les règles des plus petites échelles. Cela crée un nouveau type de géométrie chaotique distinct, où les événements énergétiques les plus violents sont verrouillés dans des régions filamenteuses et fines à l'interface.

En bref : La rupture et la fusion des gouttelettes ne font pas que perturber le flux ; elles impriment un motif de chaos unique, extrême et hautement organisé sur les échelles les plus infimes du fluide.

Résumé technique

Résumé Technique : Intermittence Sub-Kolmogorov et Dissipation Multifractale dans la Turbulence Multiphasique

Énoncé du Problème
Les écoulements turbulents multiphasiques présentent une intermittence plus forte que leurs homologues monophasiques, pourtant l'origine physique et l'organisation géométrique de leurs fluctuations les plus intenses à petite échelle restent mal comprises. Bien qu'il soit établi que la tension superficielle assure un échange dépendant de l'échelle entre l'énergie cinétique turbulente et l'énergie interfaciale — modifiant ainsi le cascade d'énergie classique — l'impact spécifique des événements de changement de topologie interfaciale (rupture et coalescence) sur les statistiques du champ dissipatif et la coupure dissipative locale reste incertain. Une question cruciale demeure : comment ces événements influencent-ils les plus petites échelles dynamiquement actives, là où une fraction significative du mélange et de la dissipation se produit, et permettent-ils de modifier l'organisation multifractale de l'écoulement ?

Méthodologie
Les auteurs utilisent des simulations numériques directes (DNS) des équations de Navier-Stokes incompressibles couplées à la tension superficielle. La dynamique de l'interface est capturée via une méthode de Volume de Fluide (VOF) implémentée dans le code FLUTAS.

• Configuration de la Simulation : Les simulations sont conduites dans une boîte périodique triple ($L=2\pi$) avec une résolution de grille de $512^3$ points. La turbulence est maintenue par un forçage à grande échelle.
• Configurations : Deux cas sont comparés : un écoulement de référence monophasique ($Re_\lambda \approx 50$, $\eta_K \approx 6.8 \Delta x$) et un écoulement multiphasique présentant une densité et une viscosité identiques pour les phases continue et dispersée ($\psi = 0.1$, $We(d_v) = 690$). L'appariement de la densité et de la viscosité isole les effets de la dynamique interfaciale et de la tension superficielle.
• Cadre Analytique : L'étude combine une analyse des échelles dissipatives fluctuantes avec un formalisme multifractal.
  • Échelle Dissipative Locale : Définie via un nombre de Reynolds local dépendant de l'échelle, $Re(x, \eta) \equiv \eta |\delta u_\eta(x)| / \nu \sim 1$, pour caractériser la fonction de densité de probabilité (PDF), $Q(\eta)$, de la coupure dissipative locale.
  • Analyse Multifractale : Le champ de dissipation est traité comme une mesure singulière. L'étude calcule les dimensions de Rényi généralisées ($D_q$) et les spectres de singularité ($f(\alpha)$) sur deux intervalles d'échelle distincts : Au-dessus de Kolmogorov (AK, $2.4 < r/\eta_K < 20$) et En-dessous de Kolmogorov (BK, $0.3 < r/\eta_K < 2.4$).
  • Validation : Le spectre multifractal est utilisé pour reconstruire théoriquement la distribution des longueurs de coupure, $Q(\eta)$, afin de vérifier la capacité prédictive de l'approche.

Résultats Clés

1. Élargissement de la Coupure Dissipative : Dans la turbulence multiphasique, la PDF de l'échelle dissipative locale, $Q(\eta)$, est nettement plus large que dans la turbulence monophasique. Plus précisément, la queue gauche est substantiellement renforcée, indiquant une haute probabilité d'événements où l'échelle dissipative locale $\eta$ est significativement plus petite que l'échelle de Kolmogorov globale ($\eta \ll \eta_K$). Cela confirme que les interfaces amplifient l'intermittence de la dynamique dissipative, permettant aux fluctuations intenses de pénétrer profondément dans la plage sub-kolmogorovienne.
2. Concentration Spatiale des Événements : Les événements sub-kolmogoroviens intenses ne sont pas distribués de manière aléatoire mais sont spatialement concentrés autour des régions interfaciales changeant de topologie, spécifiquement lors de la rupture et de la coalescence. Ces régions sont caractérisées par une grande courbure interfaciale et des fluctuations de vitesse extrêmes à petite échelle.
3. Divergence du Spectre Multifractal :
   • Au-dessus de Kolmogorov (AK) : Le spectre de singularité pour le cas multiphasique reste proche du cas monophasique pour les valeurs modérées et élevées de $\alpha$, suggérant que les régions faiblement dissipatives sont seulement légèrement affectées par l'interface.
   • En-dessous de Kolmogorov (BK) : Le cas multiphasique présente un comportement qualitativement différent. Le spectre développe un spectre de singularité nettement plus large avec des queues gauches fortes, supportant des valeurs très petites de $\alpha$ sur des ensembles possédant une dimension fractale finie. La valeur $f(\alpha \approx 0) \approx 1$ suggère que les structures dissipatives les plus singulières sont supportées par des ensembles approximativement unidimensionnels (par exemple, des ligaments amincis et des filaments en reconnexion).
4. Capacité Prédictive : En utilisant le spectre de singularité dérivé de la zone BK, la reconstruction théorique de $Q(\eta)$ concorde avec les données DNS avec un excellent accord. Inversement, l'utilisation du spectre AK échoue à capturer les événements sub-kolmogoroviens, démontrant que l'approche multifractale n'est prédictive que lorsque les intervalles d'échelle appropriés (spécifiquement ceux capturant les effets interfaciaux) sont sélectionnés.

Signification et Revendications
Cet article établit que la rupture et la coalescence dans la turbulence multiphasique ne font pas que perturber l'écoulement localement ; elles impriment une géométrie multifractale distincte sur le champ de dissipation. Les résultats démontent que la dynamique interfaciale soutient un régime dissipatif fortement multifractal au voisinage et en-dessous de l'échelle de Kolmogorov, ce qui est fondamentalement différent de l'intermittence de la gamme inertielle de la turbulence monophasique. En liant les fluctuations du nombre de Reynolds local au spectre de singularité, l'étude fournit un lien direct entre les changements de topologie interfaciale et l'activité extrême à petite échelle de la turbulence multiphasique. Ce travail étend le formalisme multifractal aux écoulements multiphasiques, offrant un cadre pour décrire comment le transfert d'énergie médié par l'interface remodèle la géométrie statistique de la dissipation.