Finite-Time Relaxation of Inertial Particle Clustering in Non-Equilibrium Turbulence

Cette étude démontre que l'approximation d'équilibre instantané échoue à prédire avec précision le regroupement de particules inertielles dans une turbulence hors équilibre en raison d'effets de relaxation à temps fini, et propose un nouveau modèle de relaxation linéaire intégrant une loi d'échelle spécifique qui réduit considérablement les erreurs de prédiction par rapport aux méthodes traditionnelles.

L'essentiel

Imaginez que vous êtes sur une piste de danse bondée (la turbulence). Au milieu de cette foule, il y a des milliers de danseurs minuscules et lourds (les particules inertielles, comme des gouttelettes d'eau dans un nuage ou de la poussière dans l'air).

Parce que ces danseurs lourds possèdent une quantité de mouvement, ils ne peuvent pas tourner instantanément comme les gens légers et agiles qui les entourent. Au lieu de cela, ils sont projetés hors des tourbillons en rotation et poussés dans les couloirs droits et étirés. Cela les amène à se regrouper en des endroits spécifiques, formant de petits groupes serrés ou des amas.

Les scientifiques savent depuis longtemps que ces amas sont importants car ils font que les danseurs lourds se heurtent plus souvent les uns aux autres. Cependant, la plupart des anciennes règles pour prédire ces amas ont été écrites pour une piste de danse où la musique et l'énergie de la foule ne changent jamais (un état statistiquement stationnaire).

Le Problème : L'Erreur de l'« Instantané »

La grande question que pose cet article est : Que se passe-t-il lorsque la musique change soudainement ?

Imaginez que le DJ passe brusquement d'un rythme lent et doux à un morceau énergique et rapide, puis revient en arrière.

• L'Ancienne Hypothèse : Les scientifiques supposaient autrefois que les danseurs lourds se réorganiseraient instantanément en de nouveaux groupes dès que le rythme changerait. Ils pensaient que l'agrégation était un « équilibre instantané ».
• La Réalité : Les auteurs de cet article ont découvert que cette hypothèse est fausse. Tout comme un danseur lourd met quelques secondes à arrêter de tourner et à courir vers un nouvel endroit, les amas de particules prennent du temps pour réagir à l'énergie changeante de la turbulence. Ils ne s'adaptent pas immédiatement à la nouvelle forme ; ils s'y « relaxent » sur une période finie.

L'Expérience : Une Piste de Danse avec un Rythme

Pour le prouver, les chercheurs ont utilisé un superordinateur pour simuler une piste de danse en 3D. Ils n'ont pas simplement laissé la musique jouer au hasard ; ils ont programmé l'injection d'énergie pour qu'elle pulse de haut en bas selon un rythme parfait (comme un battement de cœur).

Ils ont testé différentes vitesses pour ce rythme :

1. Rythme Rapide : Le rythme changeait si vite que les danseurs lourds ne pouvaient pas du tout suivre.
2. Rythme Lent : Le rythme changeait assez lentement pour que les danseurs aient le temps de réagir, mais pas assez lentement pour qu'ils soient parfaitement synchronisés.

Ce qu'ils ont découvert :
Lorsque le rythme était suffisamment lent (spécifiquement, lorsque le temps entre les battements était plus long que le temps qu'il faut à un grand tourbillon dans la foule pour faire un tour complet), les amas ont montré un phénomène appelé hystérésis.

Pensez à l'hystérésis comme à une porte collante.

• Si vous poussez la porte pour l'ouvrir (augmentation de l'énergie), elle s'ouvre à un certain point.
• Si vous tirez pour la fermer (diminution de l'énergie), elle ne se ferme pas exactement au même point ; elle reste ouverte un peu plus longtemps à cause de la « colle » (l'inertie).
• Dans la simulation, pour la même quantité d'énergie dans la pièce, les amas étaient totalement différents selon que l'énergie était en train de monter ou en train de descendre.
  • Lorsque l'énergie montait, les amas étaient très faibles (seulement 80 % de la taille attendue).
  • Lorsque l'énergie descendait, les amas étaient très forts (156 % de la taille attendue).

Cela a prouvé que vous ne pouvez pas simplement regarder le niveau d'énergie actuel pour savoir comment les particules sont agglomérées ; vous devez connaître l'histoire de la manière dont l'énergie y est arrivée.

La Solution : Un Nouveau Code de Règles

Les chercheurs ont réalisé que l'ancien code de règles « instantané » échouait. Alors, ils ont construit un nouveau modèle, plus simple, pour le corriger.

Ils ont traité le processus d'agrégation comme un ressort ou un amortisseur sur une voiture.

• Lorsque la route (la turbulence) change, la voiture ne s'adapte pas instantanément à la nouvelle hauteur. Elle rebondit et se stabilise sur une période de temps spécifique.
• Ils ont calculé exactement combien de temps prend ce « temps de stabilisation » (temps de relaxation). Ils ont découvert qu'il dépend de deux choses :
  1. La taille des tourbillons dans la foule (temps de rotation des grands tourbillons).
  2. Le poids des danseurs par rapport à la foule (nombre de Stokes).

Leur nouvelle formule est : Temps de Relaxation = (Taille du Tourbillon) × (Poids)^0,40.

Le Résultat : Des Prédictions Bien Meilleures

Ils ont testé ce nouveau modèle « ressort » contre leurs simulations informatiques.

• L'Ancien Modèle (Instantané) : Faisait d'énormes erreurs, parfois en étant à près de 50 % pour les particules les plus lourdes. C'était comme deviner la hauteur de la voiture sans tenir compte du rebond.
• Le Nouveau Modèle (Temps Fini) : A considérablement amélioré la précision, réduisant l'erreur à seulement 10 %. Même lorsqu'ils l'ont testé sur un ensemble de conditions complètement différent (une autre « piste de danse »), il a toujours réduit l'erreur de 76 % à 22 %.

L'Essentiel

Cet article nous dit que dans le monde chaotique de la turbulence hors équilibre (où l'énergie change constamment), les particules ne réagissent pas instantanément. Elles ont une « mémoire » et un temps de réaction. En reconnaissant ce délai et en ajoutant un simple « temps de stabilisation » à nos calculs, nous pouvons prédire comment les particules s'agglutinent avec beaucoup plus de précision. Cela est crucial pour comprendre des phénomènes comme la formation de la pluie dans les nuages, où le moment des collisions entre gouttelettes compte énormément.

Résumé technique

Résumé technique : Relaxation en temps fini du regroupement de particules inertielles dans la turbulence hors équilibre

Énoncé du problème
Les particules inertielles dans les écoulements turbulents présentent une concentration préférentielle, formant des amas qui influencent considérablement les taux de collision et les propriétés de transport. Les modèles actuels de regroupement de particules, en particulier en microphysique des nuages, sont largement dérivés d'une turbulence statistiquement stationnaire. Lorsqu'ils sont appliqués à une turbulence variant dans le temps et hors équilibre, ces modèles supposent implicitement une « approximation d'équilibre instantané », postulant que la fonction de distribution radiale (RDF), $g(r)$, s'ajuste immédiatement à l'état turbulent instantané (par exemple, le taux de dissipation d'énergie). Cependant, la validité de cette hypothèse dans des conditions hors équilibre reste incertaine. Des études antérieures ont suggéré que le regroupement de particules peut répondre aux changements d'état turbulent avec un délai temporel, présentant une hystérésis, mais l'intensité de cette hystérésis et les temps de relaxation associés n'ont pas été caractérisés quantitativement.

Méthodologie
Les auteurs ont employé la simulation numérique directe (DNS) d'une turbulence isotrope homogène en utilisant le simulateur de nuages lagrangien en Julia (LCS.jl). L'étude a utilisé deux stratégies de forçage principales pour générer une turbulence hors équilibre :

1. Forçage périodique instationnaire : Le taux d'injection d'énergie, $P(t)$, a été commuté alternativement entre des phases fortes ($P_s$) et faibles ($P_w$) avec des périodes variables ($T$). Cela a permis d'évaluer les réponses périodiques et d'extraire des boucles d'hystérésis par moyennage de phase sur 100 cycles.
2. Forçage instationnaire par impulsion unique : Le taux d'injection d'énergie a été commuté une fois d'un état stationnaire fort vers un état stationnaire faible. Cette réponse transitoire a été utilisée pour construire et valider un modèle de relaxation linéaire.

Les simulations couvrent une gamme de nombres de Reynolds de Taylor ($Re_\lambda \approx 84\text{--}216$) et de nombres de Stokes ($St_s = 0.25\text{--}4.0$). L'intensité du regroupement a été quantifiée à l'aide de la RDF à la distance de contact, $g \equiv g(r=2r_p)$.

Résultats clés

• Échelle hors équilibre : L'écoulement turbulent a systématiquement présenté une échelle de dissipation hors équilibre ($C_\epsilon \propto Re_\lambda^{-1}$) pour toutes les périodes de forçage, confirmant la génération d'un état turbulent hors équilibre.
• Hystérésis dans le regroupement : Une boucle d'hystérésis claire a été observée dans la relation entre le taux de dissipation d'énergie instantané ($\epsilon$) et l'intensité du regroupement ($g^*$) lorsque la période de forçage ($T^*$) dépassait plusieurs temps de rotation des grands tourbillons ($T_e$).
  • Pour $T^* = 12.0$ et $St_s = 4.0$, l'intensité du regroupement a pris deux valeurs distinctes (0,80 et 1,56 fois la référence) pour le même $\epsilon$ instantané, selon l'historique de l'écoulement.
  • Cette hystérésis dépassait les fluctuations naturelles observées dans la turbulence statistiquement stationnaire, démontrant que l'approximation d'équilibre instantané échoue dans ces conditions.
  • La direction et la forme des boucles d'hystérésis dépendaient du nombre de Stokes, reflétant la relation non monotone entre l'intensité du regroupement à l'équilibre et $\epsilon$.
• Relaxation en temps fini : L'étude a identifié que l'intensité du regroupement se relâche vers la valeur d'équilibre instantané avec une échelle de temps finie. Un modèle de relaxation linéaire a été construit :
  $$\frac{d g_{noneq}}{dt} = \frac{g_{eq}(t) - g_{noneq}(t)}{\tau_g}$$
  où le temps de relaxation $\tau_g$ évolue selon $\tau_g = 1.0 T_e(t) St(t)^{0.40}$.

Validation du modèle et performance
Le modèle de relaxation linéaire proposé (modèle noneq) a été validé par rapport à des données DNS indépendantes (Cas S2) et comparé au modèle d'équilibre instantané (modèle eq).

• Pour les particules ayant la plus grande inertie ($St_s = 4.0$), l'erreur relative maximale du modèle noneq a été réduite de 49 % (modèle eq) à 10 % dans le cas d'entraînement (S1) et de 76 % à 22 % dans le cas de validation indépendant (S2).
• L'amélioration de la précision de la prédiction était la plus significative pour les nombres de Stokes élevés, où le temps de relaxation est le plus long. Pour les nombres de Stokes plus faibles, l'approximation d'équilibre instantané restait relativement précise en raison de temps de relaxation plus courts.

Signification
Cette étude démontre que l'approximation d'équilibre instantané est inappropriée pour décrire le regroupement de particules inertielles dans une turbulence hors équilibre lorsque la période de forçage dépasse plusieurs temps de rotation des grands tourbillons. En caractérisant la réponse du regroupement comme un processus de relaxation en temps fini, les auteurs fournissent une base quantitative pour améliorer les modèles de regroupement dans les écoulements turbulents variant dans le temps. La loi d'échelle identifiée pour le temps de relaxation permet de construire des modèles prédictifs plus précis pour les écoulements chargés de particules dans des conditions hors équilibre, telles que celles rencontrées en microphysique des nuages atmosphériques et dans les applications d'ingénierie.