A HHO formulation for variable density incompressible flows where the density is purely advected

Cet article propose une formulation Hybrid High-Order (HHO) pour les écoulements incompressibles à densité variable qui assure une conservation exacte du volume et une pure advection de la densité, tout en offrant des propriétés avantageuses telles que la robustesse vis-à-vis de la pression et la conservation des bornes de densité, validées par des tests de convergence et l'étude de l'instabilité de Rayleigh-Taylor.

L'essentiel

Imaginez que vous essayez de simuler le comportement de deux liquides qui ne se mélangent pas, comme de l'huile et de l'eau, ou de l'air chaud et froid, dans un ordinateur. C'est un défi énorme pour les mathématiciens et les ingénieurs. Si vous utilisez les bons outils, vous pouvez prédire comment ces fluides vont tourbillonner, se séparer ou se mélanger sous l'effet de la gravité.

C'est exactement ce que font Lorenzo Botti et Francesco Carlo Massa dans leur article. Ils ont créé une nouvelle méthode de calcul très puissante, appelée HHO (Hybrid High-Order), spécialement conçue pour gérer des fluides dont la densité change (comme un mélange d'air et d'eau).

Voici une explication simple de leur travail, avec quelques images pour aider à visualiser :

1. Le problème : Mélanger des fluides sans "fuite"

Dans la vraie vie, si vous avez un seau d'eau, le volume d'eau reste le même, peu importe comment vous le secouez. En informatique, simuler cela est difficile. Souvent, les méthodes numériques font des "fuites" invisibles : l'eau semble disparaître ou apparaître par magie, ce qui fausse tout le résultat.

De plus, quand on a deux fluides de densités différentes (l'un lourd, l'autre léger), la densité doit simplement "glisser" avec le courant, comme une feuille sur un ruisseau. Si la méthode de calcul est imparfaite, cette feuille peut se déformer bizarrement ou se transformer en quelque chose d'impossible.

2. La solution : La méthode "HHO" (L'architecte hybride)

Les auteurs proposent une méthode qui agit comme un architecte très méticuleux.

• L'approche hybride : Imaginez que pour construire une maison (simuler le fluide), vous avez besoin de deux types d'ouvriers : ceux qui travaillent à l'intérieur de chaque pièce (les cellules du maillage) et ceux qui travaillent sur les murs entre les pièces (les arêtes). La méthode HHO utilise les deux. Cela permet d'avoir une précision incroyable tout en gardant le contrôle sur chaque pièce individuellement.
• La conservation parfaite : Grâce à cette méthode, le volume d'eau est conservé à la perfection, jusqu'à la dernière décimale de l'ordinateur. C'est comme si chaque goutte d'eau avait un GPS et que l'ordinateur savait exactement où elle est, sans jamais en perdre une seule.
• La densité pure : Puisque le volume est parfaitement conservé, la densité (le poids du fluide) se déplace simplement comme un passager dans un bus. Elle ne change pas de poids, elle ne fait que suivre le mouvement. C'est ce qu'ils appellent "l'advection pure".

3. La robustesse : Le bouclier contre les erreurs

L'un des plus grands problèmes dans ces simulations est que la pression (la force qui pousse le fluide) peut parfois "polluer" le calcul de la vitesse. C'est comme si vous essayiez de mesurer la vitesse d'une voiture, mais que le compteur de vitesse était faussé par la température extérieure.

Leur méthode est "robuste à la pression". Cela signifie que peu importe à quel point le calcul de la pression est complexe ou approximatif, la vitesse du fluide et la densité restent parfaitement justes. C'est comme avoir un GPS qui continue de vous indiquer le bon chemin même si la carte routière (la pression) est un peu floue.

4. Le moteur : Le chronomètre intelligent (ESDIRK)

Pour simuler le temps qui passe, ils utilisent un type de chronomètre très sophistiqué appelé ESDIRK.
Imaginez que vous devez filmer un match de football en très haute définition. Vous ne pouvez pas juste prendre une photo toutes les secondes (trop flou) ni une photo à chaque milliseconde (trop lourd pour l'ordinateur).
Leur chronomètre est intelligent : il prend des photos à des moments stratégiques, en anticipant le mouvement, ce qui permet d'avoir un film très fluide et précis sans surcharger l'ordinateur.

5. L'économie d'énergie : La condensation statique

Simuler ces fluides demande beaucoup de mémoire (comme remplir une bibliothèque de livres). Pour éviter que l'ordinateur ne plante, ils utilisent une astuce appelée condensation statique.
C'est comme si, au lieu de ranger chaque livre individuellement sur une étagère, vous les empiliez dans des boîtes compactes. Vous gardez tout l'information, mais vous occupez beaucoup moins de place. Cela permet de faire des calculs très précis (haute ordre) même sur des ordinateurs standards.

6. Le test final : La Tour de Rayleigh-Taylor

Pour prouver que leur méthode fonctionne, ils l'ont testée sur un scénario célèbre : la Tour de Rayleigh-Taylor.
Imaginez un verre d'eau où vous versez délicatement de l'huile par-dessus. Si vous secouez le verre, l'huile (plus légère) essaie de monter et l'eau (plus lourde) essaie de descendre, créant de jolis motifs de tourbillons.

• Cas facile (faible différence de poids) : Ils ont simulé cela avec une grande précision. Même en utilisant des calculs très fins (haute précision) sur un maillage grossier, ils ont obtenu le même résultat que des calculs très lourds sur un maillage fin. C'est comme obtenir une photo HD en utilisant moins de pixels.
• Cas difficile (forte différence de poids) : Quand la différence de poids est énorme, les fluides deviennent très instables. Là, leur méthode a prouvé qu'elle ne créait pas de résultats impossibles (comme une densité négative, ce qui n'existe pas dans la réalité). Elle reste stable, comme un pilote d'avion expert qui garde l'avion en vol même dans une tempête.

En résumé

Botti et Massa ont créé un outil de simulation qui est :

1. Précis : Il ne perd aucune goutte d'eau.
2. Fiable : Il ne se trompe pas à cause de la pression.
3. Économe : Il utilise moins de mémoire que les méthodes classiques.
4. Puissant : Il peut simuler des mélanges d'huiles et d'eau complexes, du calme à la tempête.

C'est une avancée majeure pour comprendre comment les fluides se comportent, que ce soit pour la météo, la conception de moteurs, ou l'étude des océans.

Résumé technique

Voici un résumé technique détaillé de l'article proposé par Botti et Massa, rédigé en français.

Titre de l'étude

Formulation HHO pour les écoulements incompressibles à densité variable où la densité est purement advectée.

1. Problématique

L'article aborde la simulation numérique des écoulements de fluides incompressibles à densité variable, spécifiquement pour des mélanges de fluides immiscibles. Le défi principal réside dans la conservation exacte du volume (incompressibilité) et l'advection pure de la densité sans diffusion numérique excessive ni oscillations non physiques. Les méthodes existantes peinent souvent à concilier haute précision, robustesse dans les régimes dominés par la convection, et conservation stricte des propriétés physiques (comme la positivité de la densité) tout en gérant efficacement les coûts de calcul pour des ordres d'approximation élevés.

2. Méthodologie

Les auteurs proposent une extension des méthodes Hybrid High-Order (HHO) aux équations de Navier-Stokes incompressibles à densité variable. La méthodologie repose sur plusieurs piliers fondamentaux :

• Discrétisation Spatiale Hybride :

  • Toutes les variables (vitesse, pression, densité) sont discrétisées de manière hybride, c'est-à-dire avec des degrés de liberté définis à la fois sur les cellules du maillage et sur les faces (squelette du maillage).
  • L'utilisation d'un espace de pression hybride permet de reconstruire une approximation de la vitesse conforme à l'espace $H(\text{div})$ à partir d'espaces polynomiaux non conformes.
  • Cette conformité $H(\text{div})$ garantit une vitesse exactement divergence-free (incompressibilité ponctuelle), ce qui permet de formuler l'équation de conservation de la masse comme une dérivée temporelle matérielle de la densité ($\partial_t \rho + \mathbf{u} \cdot \nabla \rho = 0$), assurant ainsi une advection pure de la densité.

• Stabilisation et Robustesse :

  • Robustesse vis-à-vis de la pression (Pressure-robustness) : Grâce à la vitesse $H(\text{div})$, les forces conservatives (irrotationnelles) n'affectent que le champ de pression, rendant les erreurs de vitesse et de densité indépendantes de la précision de l'approximation de la pression.
  • Régime dominé par la convection : Des stabilisations de type "upwind" sont introduites pour le transport de masse et de quantité de mouvement. Deux variantes sont proposées :
    • HHO-$\pi^k$ : Utilise des projections $L^2$ de degré $k$ pour les vitesses advectées, garantissant une stabilité prouvée (formes skew-symétriques).
    • HHO-$\pi^{k+1}$ : Utilise des projections de degré $k+1$, offrant une meilleure précision dans les régimes à fort nombre de Reynolds, bien que la stabilité soit moins triviale.
  • Stabilisation temporelle : Un terme de stabilisation supplémentaire est ajouté pour la densité afin de garantir la robustesse lorsque la composante normale de la vitesse s'annule sur une face.

• Discrétisation Temporelle (ESDIRK) :

  • L'intégration temporelle utilise des schémas ESDIRK (Explicit Singly Diagonal Implicit Runge-Kutta) d'ordre élevé, conçus spécifiquement pour les Équations Différentielles Algébriques (DAE) d'indice 2.
  • Le schéma est entièrement implicite, permettant des pas de temps flexibles et une grande stabilité.

• Aspects Algébriques et Efficacité :

  • Condensation Statique : Les inconnues internes aux cellules sont éliminées localement, ne laissant que les inconnues sur les faces comme inconnues globales. Cela réduit considérablement la taille et la bande de la matrice jacobienne globale.
  • Préconditionnement : Une stratégie de préconditionnement p-multigrid est utilisée pour accélérer la convergence des solveurs linéaires itératifs (FGMRES), exploitant la structure hiérarchique des espaces polynomiaux.

3. Contributions Clés

1. Conservation Exacte du Volume : La formulation impose la conservation du volume cellule par cellule jusqu'à la précision machine, assurant que la densité est strictement une quantité advectée.
2. Robustesse et Stabilité : La méthode est prouvée stable (dissipative) pour des conditions aux limites de Dirichlet homogènes, avec des termes de stabilisation adaptés pour les régimes à convection dominante et les faibles viscosités.
3. Préservation des Bornes (Low-Order) : Pour les ordres polynomiaux faibles ($k=0$), la méthode préserve les bornes de la densité, évitant les densités négatives, ce qui est crucial pour les rapports de densité élevés.
4. Efficacité Computationnelle : La combinaison de la condensation statique et du préconditionnement p-multigrid permet d'obtenir des solutions de haute précision sur des maillages grossiers, réduisant le coût de calcul par rapport aux méthodes de bas ordre sur des maillages très fins.

4. Résultats Numériques

Les auteurs valident la méthode sur plusieurs cas tests :

• Cas de Kovasznay (Écoulement stationnaire) :

  • Convergence spatiale observée conforme aux théories : taux $k+1$ pour la pression et $k+2$ pour la vitesse (en norme $L^2$) dans le régime diffusif.
  • La variante HHO-$\pi^{k+1}$ montre une précision supérieure pour $k=0$ (convergence d'ordre 2 pour la vitesse).
  • La densité reste constante (égale à 1) avec une erreur négligeable, confirmant l'advection pure.

• Solution Manufacturée de Guermond et Quartapelle (Écoulement instationnaire) :

  • Validation des taux de convergence spatiale et temporelle.
  • Convergence temporelle d'ordre 2 avec un schéma ESDIRK à 3 étapes (TR-BDF2) et d'ordre 4 (pour la densité) / 3 (pour la vitesse) avec un schéma à 6 étapes.
  • Confirmation de la robustesse vis-à-vis de la pression : l'erreur de vitesse est indépendante de l'erreur de pression, même pour des champs de pression complexes.

• Instabilité de Rayleigh-Taylor (RTI) :

  • Faible nombre d'Atwood (At=0.5) : Comparaison entre une simulation haute ordre ($k=6$) sur un maillage grossier et une simulation basse ordre ($k=1$) sur un maillage fin. Les résultats sont en excellent accord pour $Re=1000$, démontrant que la haute précision spatiale permet de réduire drastiquement le nombre de degrés de liberté (réduction de 14 fois de la dimension de la matrice jacobienne).
  • Fort nombre d'Atwood (At=0.75) : Simulation avec $k=0$ sur un maillage très fin. La méthode préserve les bornes de densité (entre 1 et 7) et reste stable, bien que la précision spatiale soit limitée pour capturer les interfaces très fines sans stabilisation par viscosité d'entropie.

5. Signification et Conclusion

Ce travail présente une avancée significative dans la modélisation numérique des écoulements multiphasiques incompressibles. La formulation HHO proposée réussit à combiner :

• Une précision spectrale (haute ordre) dans l'espace et le temps.
• Des propriétés de conservation rigoureuses (volume, advection pure).
• Une efficacité computationnelle grâce à la condensation statique et aux solveurs multigrilles.

La capacité à simuler des instabilités complexes comme Rayleigh-Taylor avec des rapports de densité élevés tout en maintenant la stabilité et la précision, même sur des maillages grossiers avec des ordres élevés, ouvre la voie à des simulations de mélanges de fluides immiscibles plus réalistes et moins coûteuses. Les auteurs prévoient d'étendre cette approche au système Cahn-Hilliard-Navier-Stokes pour traiter plus spécifiquement les interfaces de phase.