Sharp-interface Simulations of Energetic Multiphase Flows with Large Density and Viscosity Ratios

Cet article propose un cadre de flux de momenturs synchronisés par région de donneur (SynDRoM) combiné à un limiteur de viscosité pour améliorer la robustesse numérique et la fidélité physique des simulations d'interfaces nettes pour les écoulements multiphasiques énergétiques présentant de grands rapports de densité et de viscosité.

L'essentiel

L'idée générale : Simuler une mer déchaînée

Imaginez que vous essayez de créer une simulation informatique d'une tempête océanique violente. Vous voulez voir comment les vagues s'écrasent, comment l'air est aspiré dans l'eau et comment les bulles se forment. C'est délicat car l'eau est lourde et épaisse, tandis que l'air est léger et fin. En termes physiques, ils présentent une différence énorme de « densité ».

Lorsque les ordinateurs tentent de simuler cela, ils plantent souvent ou produisent des résultats bizarres et impossibles (comme de l'eau qui se transforme soudainement en fantôme ou de l'air qui traverse l'eau comme une balle). Ce document présente un nouvel ensemble de règles (algorithmes) pour rendre ces simulations stables, précises et physiquement réalistes, même lorsque les vagues s'écrasent violemment.

Le problème : Le « Fantôme » et le « Choc »

Les auteurs expliquent que les anciennes méthodes de simulation de ces écoulements présentent deux défauts principaux :

1. Le problème du « Fantôme » (Pénétration de la vitesse) :
   Imaginez un camion lourd (l'eau) et une plume (l'air) se déplaçant l'un à côté de l'autre. Dans les anciennes simulations, le « vent » de la plume pouvait parfois faire reculer le camion, ou le camion pouvait pousser la plume à travers son propre corps. C'est ce qu'on appelle la « pénétration de la vitesse ». Cela crée des formes non physiques et fausses dans l'eau, comme une « corne du diable » dépassant de la vague.

2. Le problème du « Choc » (Pics de quantité de mouvement) :
   Pour corriger le problème du fantôme, les scientifiques ont tenté une nouvelle méthode appelée CMOM (Mass-Momentum Consistent). C'est comme tenir un registre strict de l'« élan » (quantité de mouvement) de chaque goutte d'eau. Cependant, cette méthode a un effet secondaire. Lorsqu'un tout petit peu d'eau lourde se déplace dans une cellule remplie d'air, les mathématiques s'embrouillent. C'est comme diviser un très grand nombre par un nombre minuscule, ce qui résulte en un pic de vitesse massif et impossible. Cela crée des « blobs de vitesse » : des poches d'air factices se déplaçant à des vitesses supersoniques qui ne devraient pas exister.

La solution : La méthode « SynDRoM »

Les auteurs proposent une nouvelle correction appelée SynDRoM (Synchronized Donor Region of Momentum flux). Voici comment cela fonctionne, en utilisant une analogie :

L'analogie : Le tapis roulant en mouvement
Imaginez un tapis roulant transportant des boîtes.

• L'ancienne méthode : Vous comptez les boîtes (masse) et le poids des boîtes (quantité de mouvement) séparément. Si une boîte se déplace, vous pourriez accidentellement compter son poids dans un endroit où la boîte n'est pas encore réellement arrivée. Cela provoque le « choc » ou le « pic » de vitesse.
• La méthode SynDRoM : Cette méthode agit comme une équipe synchronisée. Avant de déplacer le poids, vous regardez exactement quelle partie du tapis roulant ce poids est en train de quitter.
  • Elle demande : « Si je déplace cette portion spécifique d'air, à quelle portion de quantité de mouvement exacte est-elle attachée ? »
  • Elle garantit que la quantité de mouvement n'est déplacée que si la masse est réellement là pour la porter.
  • Le résultat : Plus de pics de vitesse factices. L'air reste lent et l'eau reste lourde, tout comme dans la réalité. La simulation reste fluide et ne « s'emballe » pas.

Le second problème : La viscosité « glissante »

Le document traite également d'un second problème : la Viscosité (la façon dont un fluide est épais ou collant).

• Le problème : L'eau est collante ; l'air est glissant. Lorsqu'ils se mélangent à une frontière nette (comme une vague qui déferle), l'ordinateur essaie de deviner la « viscosité » au milieu. S'il se trompe, les mathématiques deviennent instables, comme essayer de faire tenir un crayon en équilibre sur sa pointe.
• La correction : Les auteurs introduisent un Limiteur de Viscosité.
  • L'analogie : Imaginez un panneau de limitation de vitesse. Même si les mathématiques tentent de calculer une « viscosité » qui ferait bouger le fluide de manière impossible (instable), le limiteur dit : « Non, vous ne pouvez pas aller plus vite que la vitesse du fluide le plus fin ici. » Il plafonne le calcul pour empêcher la simulation de planter, sans changer la physique réelle de l'eau ou de l'air.

La preuve : Est-ce que cela fonctionne ?

Les auteurs ont testé leurs nouvelles règles de trois manières :

1. La rupture de barrage (Dam Break) : Ils ont simulé l'effondrement d'un mur d'eau.
   • Anciennes méthodes : L'eau paraissait déformée avec de faux pics.
   • SynDRoM : L'eau s'est écrasée naturellement, et l'air n'a pas été aspiré dans l'eau de manière étrange.
2. L'instabilité de Kelvin-Helmholtz : C'est le phénomène où le vent souffle sur l'eau, créant des vagues qui s'enroulent (comme des nuages).
   • Résultat : La simulation a correctement montré les vagues s'enrouler et croître, sans que l'ordinateur n'ajoute de l'énergie factice ou n'amortisse les vagues. Cela prouve que la méthode respecte les lois de la physique.
3. La vague déferlante : Ils ont simulé une immense vague diagonale qui s'écrase.
   • Résultat : La vague a déferlé, a éclaboussé et a créé de l'écume comme un véritable océan. L'énergie totale du système est restée équilibrée (elle n'a pas mystérieusement disparu ou explosé). Même en ajoutant de la « viscosité », la simulation est restée stable.

L'essentiel

Ce document présente un nouveau « agent de circulation » pour les simulations informatiques d'eau et d'air.

• Il empêche l'air de traverser l'eau comme un fantôme.
• Il empêche l'eau de créer des pics de vitesse impossibles.
• Il empêche les calculs de « viscosité » de briser les mathématiques.

En synchronisant exactement ce qui se déplace avec l'endroit où il se déplace, les auteurs ont créé un outil de simulation beaucoup plus robuste et fiable pour étudier les événements océaniques violents, du genre dont les ingénieurs navals ont besoin pour la conception de navires.

Résumé technique

Résumé Technique : Simulations à interface nette de flux multiphasiques énergétiques avec de grands rapports de densité et de viscosité

Énoncé du problème
La simulation de flux multiphasiques énergétiques présentant des rapports de densité extrêmes (par exemple, les systèmes air-eau avec des rapports de l'ordre de $10^3$ à $10^4$) demeure un défi majeur en hydrodynamique navale. Les formulations traditionnelles basées sur la vitesse souffrent souvent de divergences numériques dues au désalignement des gradients de pression et des discontinuités de densité, entraînant des transferts de quantité de mouvement parasites et une explosion de l'énergie. Bien que les cadres de transport de type Consistent Mass-Momentum (CMOM) améliorent la robustesse en imposant la conservation de la quantité de mouvement et une conservation semi-discrète de l'énergie, ils introduisent de nouvelles difficultés lorsqu'ils sont combinés avec des méthodes à interface nette. Plus précisément, le CMOM peut générer des chocs de quantité de mouvement sévères, entraînant des oscillations de vitesse non physiques et des dépassements (violant la condition de Variation Totale Limitée, ou TVD) lorsque les interfaces nettes traversent les frontières de mailles. De plus, une estimation inappropriée de la viscosité à proximité des interfaces, notamment lors de la division de la viscosité dynamique par la densité sur des grilles décalées (staggered), peut introduire des instabilités numériques à des pas de temps finis, compromettant la robustesse du solveur.

Méthodologie
L'article propose une suite de schémas numériques préservant la physique, conçus pour concilier fidélité physique et robustesse numérique pour les flux à rapport de densité élevé.

1. Synchronized Donor-Region of Momentum flux (SynDRoM) :
   Pour remédier aux oscillations de vitesse inhérentes au CMOM standard, les auteurs introduisent le SynDRoM. Cette méthode impose la monotonicité du champ de vitesse transporté en synchronisant le flux de quantité de mouvement avec le flux de masse.

• Mécanisme : Contrairement aux limiteurs de pente existants axés sur l'interface qui se concentrent uniquement sur la reconstruction de la vitesse, le SynDRoM reconstruit une paire distribution densité-vitesse au sein d'une cellule pour satisfaire les contraintes de moyenne cellulaire et les bornes TVD. Il fait ensuite évoluer la quantité de mouvement via un concept de « région de donneur » (donor region) déterminé strictement par le flux de masse. La vitesse advectée est calculée comme une moyenne pondérée sur cette région de donneur spécifique, garantissant que la vitesse mise à jour reste bornée.
• Implémentation : La méthode opère dans un cadre de séparation directionnelle (directionally split) cohérent avec l'advection cVOF (Volume de Fluide Conservatif). Elle utilise un schéma d'intégration temporelle de second ordre de Runge-Kutta (RK2) pondéré par la densité pour gérer l'évolution non linéaire de la vitesse par rapport à la densité et à la quantité de mouvement.

2. Limiteur de viscosité cinématique borné :
   Pour traiter les instabilités découlant des sauts de viscosité aux interfaces, l'article introduit un limiteur de viscosité.

• Mécanisme : Les auteurs identifient que l'interpolation standard (arithmétique ou harmonique) de la viscosité dynamique aux sommets des grilles décalées, suivie de la division par la densité, peut engendrer une viscosité cinématique effective non bornée, violant les contraintes de stabilité CFL visqueuse. Le limiteur proposé contraint la viscosité cinématique effective utilisée dans les termes de différence croisée par la densité minimale des cellules de quantité de mouvement environnantes et par une borne supérieure basée sur l'occupation de phase pour la viscosité cinématique. Cela garantit que la viscosité cinématique effective reste dans une plage préservant la stabilité.

Contributions clés

• Algorithme SynDRoM : Une méthode généralisée pour le CMOM qui élimine les oscillations de vitesse parasites sans sacrifier la conservation de la quantité de mouvement ou la netteté de l'interface. Elle prend explicitement en compte la distribution de la quantité de mouvement et la région de donneur correspondante basée sur le flux de masse, évitant ainsi le besoin d'une reconstruction complète de l'interface sur les grilles décalées.
• Stabilisation de la viscosité : Un limiteur de viscosité simple mais efficace qui empêche les divergences numériques dans les flux à rapport de viscosité élevé en bornant la viscosité cinématique effective, préservant ainsi la stabilité sans compromettre la fidélité physique.
• Cohérence de séparation directionnelle : Une formulation qui garantit que les mises à jour de la quantité de mouvement restent cohérentes avec l'advection de masse par séparation directionnelle, empêchant les flux fictifs qui pourraient violer le caractère borné de la vitesse.

Résultats
Les méthodes proposées ont été validées par une série de cas de test :

• Éclatement de barrage 2D (2D Dam Break) : La formulation SynDRoM a résolu avec succès la déformation de l'interface et a éliminé les structures de vitesse parasites (telles que les distorsions de type « cornes du diable » observées dans les formulations de vitesse et les « blobs de vitesse » dans l'original CMOM) tout en maintenant la conservation de l'énergie.
• Advection de gouttelette 1D : Dans un test quasi-unidimensionnel avec un rapport de densité élevé ($10^{-3}$), le SynDRoM a maintenu la monotonicité de la vitesse et les propriétés TVD, alors que le CMOM original et les variantes de type Favre produisaient des dépassements non physiques.
• Instabilité de Kelvin-Helmholtz : Le schéma a correctement capturé la croissance linéaire des perturbations et l'évolution non linéaire subséquente (appariement de vortex, roulement) sans amortissement artificiel ni explosion d'énergie. La quantité de mouvement horizontale a été conservée selon la tolérance du solveur, et les régions de haute vitesse ont été correctement confinées à la phase de fluide léger.
• Écoulement de Poiseuille stratifié : Le limiteur de viscosité a réussi à prévenir la divergence numérique observée avec l'arithmétique standard et a corrigé les gradients de vitesse excessifs vus avec l'harmonique, produisant des solutions stables et physiquement exactes.
• Briseur de vague stationnaire 3D : Les simulations de briseurs de vagues stationnaires diagonaux ont démontré la capacité du solveur à gérer une non-linéarité extrême, la formation de jets et l'entraînement d'air. La méthode a maintenu l'énergie mécanique bornée même en l'absence de viscosité explicite, la dissipation d'énergie se produisant naturellement via le mécanisme iLES (Large Eddy Simulation implicite) lorsque le flux devient irrégulier. Les études de convergence de grille ont confirmé des tendances cohérentes dans l'évolution de l'énergie et les taux de dissipation.

Signification
L'article affirme que les modifications proposées aux schémas d'advection et visqueux améliorent considérablement la fidélité et la robustesse des simulations de flux multiphasiques, particulièrement pour les flux de surface énergétiques rencontrés en hydrodynamique navale. En garantissant que les méthodes CMOM satisfont la condition TVD près des cellules d'interface se vidant, et en bornant la viscosité cinématique, cette approche élimine les artefacts numériques qui pourraient déclencher de faux événements de rupture. Les auteurs soulignent que le SynDRoM et le limiteur de viscosité sont orthogonaux à la méthode spécifique de capture d'interface utilisée (démontrée ici avec le VOF géométrique), rendant l'approche largement applicable. Ce travail offre une voie pour simuler des flux multiphasiques complexes à rapport de densité élevé avec une cohérence physique, une conservation de la quantité de mouvement et une limitation empirique de l'énergie, répondant ainsi à des défis de longue date dans le domaine.