Trois traducteurs et la danse du tourbillon ✛

Original: "Effects of Turbulence Modeling and Parcel Approach on Dispersed Two-Phase Swirling Flow"
· Osama A. Marzouk, E. David Huckaby
arXiv:2501.00037 · 2024-12-24 · CC BY 4.0 · 4 min
La comparaison de trois modèles k-ε et de deux approches pour les particules dans un écoulement tourbillonnaire a montré que le modèle standard est plus précis, et que le regroupement des particules n'accélère presque pas le calcul.

Un tube vertical avec un élargissement brusque est un terrain d'essai classique pour la turbulence. L'air, tourbillonnant comme une tornade, s'engouffre dans une chambre large et engendre des tourbillons, des écoulements inverses et des zones où des particules de verre de 45 µm de diamètre, soit la moitié de l'épaisseur d'un cheveu, décrivent des trajectoires complexes. Martin Sommerfeld et son collègue Hui Qiu ont mesuré les vitesses dans un tel écoulement dès les années 1990, et depuis, leur expérience est devenue une référence pour valider les modèles d'écoulements diphasiques. Trois versions populaires du modèle k-ε (une façon d'estimer les fluctuations turbulentes via l'énergie et la dissipation) prétendent décrire avec précision cette danse des tourbillons. Chacune est comme un traducteur qui entend la même musique mais la restitue à sa manière.

Les traducteurs divergent dans les détails. Le modèle k-ε standard, le plus ancien et le plus simple, s'avère le plus fidèle : ses profils de vitesse coïncident avec l'expérience dans les quatre sections (x = 3, 52, 112 et 155 mm). Le modèle réalisable, créé plus tard pour corriger les défauts du standard, se trompe sur la vitesse radiale : loin des jets, il prédit un flux vers l'intérieur trop fort, et près de la paroi, un flux vers l'extérieur, alors que dans l'expérience le flux était dirigé vers l'intérieur. Le modèle RNG, dérivé de la mécanique statistique, génère des zones de recirculation supplémentaires au début du bulbe central, des tourbillons factices absents de l'expérience. En termes de temps de calcul pour un intervalle de 0,5 s, le modèle standard s'avère aussi le plus rapide : 17,13 heures contre 19,44 pour le RNG et 24,81 pour le réalisable (soit 45% de plus).

Les particules de verre de 45 µm de diamètre représentent environ la moitié de l'épaisseur d'un cheveu humain. Avec un chargement de 0,034 kg/kg, le flux contenait environ 12 000 particules individuelles, mais leur masse totale représentait moins de 4% de la masse de l'air.

Les particules peuvent aussi être décrites de deux manières : suivre chacune individuellement (approche lagrangienne, qui consiste à suivre chaque particule le long de sa trajectoire) ou les regrouper en paquets, appelés parcels, où une centaine de particules identiques se déplacent comme une seule. La logique suggère que le regroupement fera gagner du temps, car il y aura moins d'équations du mouvement. Mais dans ce test, l'économie n'est que de 4,2% (31,15 heures contre 29,84), soit presque rien. La raison est simple : la charge de calcul principale provient de la résolution des équations de la phase gazeuse, et les particules sont peu nombreuses. Les profils moyens de vitesse du gaz avec un regroupement jusqu'à 100 particules par paquet ne changent pas du tout, la différence est imperceptible à l'œil.

Ce travail donne aux ingénieurs un repère clair : pour les écoulements diphasiques tourbillonnaires avec une charge modérée de particules, le modèle k-ε standard reste un compromis raisonnable entre précision et coût de calcul. Le modèle réalisable, malgré des améliorations théoriques, peut donner une vitesse radiale qualitativement incorrecte, et le RNG peut produire de fausses structures tourbillonnaires. L'approche par paquets pour les particules n'est justifiée que lorsque la concentration de la phase dispersée est élevée, car l'économie sur le nombre d'équations du mouvement des particules devient alors significative. Aujourd'hui, cette connaissance s'applique à la conception de foyers à charbon pulvérisé, de séparateurs cycloniques et de systèmes de gazéification avec capture du CO2, où il est important de prédire la distribution des particules et les zones de recirculation.

La prochaine étape naturelle consiste à répéter la comparaison avec une charge de particules plus élevée (supérieure à 0,034 kg/kg) et un taux de swirl plus fort (nombre de swirl supérieur à 0,6), caractéristique des brûleurs réels. Il faudra peut-être adapter les fonctions de paroi (description approchée de l'écoulement près de la paroi) ou passer à des modèles résolvant la sous-couche visqueuse pour mieux décrire les écoulements secondaires. Mais une chose est déjà claire : dans le monde des écoulements diphasiques turbulents, un traducteur fiable vaut mieux qu'un traducteur à la mode.

Les particules de verre de 45 μm de diamètre représentent environ la moitié de l'épaisseur d'un cheveu humain. Avec une charge de 0,034 kg/kg, il y avait environ 12 000 particules individuelles dans l'écoulement, mais leur masse totale représentait moins de 4 % de la masse de l'air.

Scientifiques
Daniel BernoulliLeonhard EulerGeorge Gabriel StokesLudwig Boltzmann
Notions
Équations de Navier-Stokes mélange turbulent turbulence équation de Poiseuille coupled general circulation model dynamique des fluides à motifs coagulation turbulente action simulation hydrodynamique marche aléatoire équation de Fokker-Planck