A bound-preserving oscillation-eliminating discontinuous Galerkin scheme for compressible two-phase flow

Cet article propose un schéma de Galerkin discontinu à haute précision, préservant les bornes et éliminant les oscillations pour les écoulements diphasiques compressibles, qui surmonte les restrictions sévères de CFL induites par la rigidité grâce à une stratégie novatrice de décomposition d'opérateurs associée à un solveur implicite adaptatif, tout en garantissant rigoureusement la stabilité, la précision et le respect de la condition d'Abgrall.

L'essentiel

Imaginez que vous essayez de simuler une collision à haute vitesse entre deux fluides très différents, comme une onde de choc dans l'eau percutant une bulle d'air. Dans le monde des simulations informatiques, c'est un cauchemar. Les fluides se comportent différemment, ils s'écrasent et s'étirent à des rythmes différents, et les mathématiques régissant leur interaction sont incroyablement « rigides ».

Pensez à cette « rigidité » comme à l'effort de conduire une voiture dont les freins sont bloqués au plancher. Si vous essayez d'avancer ne serait-ce qu'un tout petit peu (un petit pas de temps dans la simulation), les freins résistent si violemment que la voiture pourrait se retourner ou le moteur pourrait exploser. En termes informatiques, cela force la simulation à prendre des pas si incroyablement petits qu'il faudrait des années pour simuler une fraction de seconde de temps réel.

Cet article introduit une nouvelle méthode, plus intelligente, pour conduire cette voiture. Voici la décomposition de leur solution en utilisant de simples analogies :

1. Le Problème : Le Frein « Rigide »

Les auteurs travaillent avec un ensemble spécifique de règles (le modèle à cinq équations de Kapila) qui décrit comment deux fluides se mélangent et se déplacent. L'ennui vient d'une règle spécifique (le terme source $\kappa$) qui gère la compression des fluides. Lorsqu'une onde de choc frappe la frontière entre l'eau et l'air, cette règle passe en surrégime.

Si l'ordinateur tente de tout résoudre d'un coup (la méthode traditionnelle), il reste bloqué. Pour éviter que les mathématiques ne s'effondrent, il doit ralentir le temps de simulation de manière si drastique que le calcul devient impossible.

2. La Solution : La Stratégie « Fractionnée »

Les auteurs proposent un astucieux tour de passe-passe appelé Décomposition d'Opérateurs. Imaginez que vous essayez de cuire un gâteau tout en réparant simultanément une fuite de tuyau. Faire les deux exactement au même moment est chaotique et susceptible d'échouer. À la place, vous les faites en étapes séparées et concentrées :

• Étape A : Réparez le tuyau (résolvez la partie de compression « rigide »).
• Étape B : Cuisez le gâteau (résolvez la partie du mouvement et de l'écoulement).

En séparant ces deux tâches, l'ordinateur peut gérer le « tuyau qui fuit » (les mathématiques rigides) en utilisant une méthode implicite spéciale, lente mais sûre, qui ne casse jamais, puis gérer la « cuisson » (l'écoulement) en utilisant une méthode rapide et de haute précision.

3. Le Filet de Sécurité « Préservant les Bornes »

Dans ces simulations, les nombres représentent des choses physiques comme la densité et la pression. Si les mathématiques se trompent, l'ordinateur pourrait calculer que l'air a une densité négative ou qu'une bulle a 150 % de son volume (ce qui est impossible). Cela provoque l'arrêt de la simulation.

Les auteurs ont construit un limiteur Préservant les Bornes (BP). Imaginez cela comme un videur dans un club. Si un nombre tente de quitter la « zone sûre » (par exemple, une fraction de volume tentant de dépasser 100 % ou de descendre en dessous de 0 %), le videur le renvoie immédiatement à l'intérieur de la zone sûre. Cela garantit que la simulation ne produit jamais de physique « absurde », même lorsque les choses deviennent chaotiques.

4. L'Amortisseur « Éliminant les Oscillations »

Lorsqu'une onde de choc frappe une bulle, elle crée des bords nets et des ondulations. Les mathématiques standards créent souvent de fausses, de saccadées « ondes fantômes » (oscillations) autour de ces bords nets, rendant l'image bruitée et erronée.

Les auteurs utilisent une technique Éliminant les Oscillations (OE). Imaginez conduire sur une route cahoteuse. Une voiture standard pourrait rebondir sauvagement. Cette nouvelle méthode agit comme un système de suspension haute technologie qui lisse la conduite sans perdre le détail des bosses. Elle élimine le faux bruit tout en gardant la physique réelle nette, et elle le fait sans avoir besoin d'effectuer des calculs complexes et lents pour déterminer la direction des ondes.

5. Le Résultat : Une Trajectoire Fluide et Rapide

Les auteurs ont testé leur nouvelle méthode sur des scénarios très difficiles :

• Choc frappant une bulle d'hélium : Comme un bang supersonique frappant une bulle de savon.
• Choc d'eau frappant une bulle d'air : Une explosion sous-marine massive frappant une poche d'air.

Dans ces tests, leur méthode a pu fonctionner rapidement (en utilisant des pas de temps standards) sans planter, tout en maintenant tous les nombres physiquement réalistes. Elle a capturé les formes complexes des bulles et des ondes de choc avec une grande précision, prouvant que l'on peut simuler ces événements extrêmes sans que l'ordinateur ne reste bloqué en « ralenti ».

En résumé : L'article présente un nouveau moteur mathématique qui divise un problème difficile en morceaux gérables, utilise un filet de sécurité pour maintenir les nombres réalistes, et lisse le bruit. Cela permet aux ordinateurs de simuler rapidement et avec précision les collisions violentes entre différents fluides, résolvant un problème qui nécessitait auparavant des quantités de puissance de calcul impossibles.

Résumé technique

1. Énoncé du problème

L'article traite de la simulation numérique des écoulements diphasiques gazeux-gazeux et gazeux-liquide compressibles régis par le modèle de Kapila à cinq équations. Ce modèle est une forme réduite du modèle de Baer-Nunziato, supposant l'équilibre mécanique et thermique (vitesses et pressions égales) entre les phases.

Défis clés :

• Rigidité du terme source $\kappa$ : L'équation de la fraction volumique contient un terme source non conservatif ($\kappa \nabla \cdot \mathbf{u}$) qui rend compte des différences de compressibilité entre les fluides. Lorsque des ondes de choc ou de détente interagissent avec des interfaces matérielles, ce terme devient extrêmement rigide.
• Restriction CFL : Les schémas d'intégration temporelle explicites traditionnels imposent des pas de temps restreints par la rigidité du terme $\kappa$, réduisant souvent le pas de temps admissible de plusieurs ordres de grandeur par rapport à la condition CFL standard. Cela rend la simulation de cas extrêmes (par exemple, les interactions choc-bulle eau-air) prohibitivement coûteuse en calcul.
• Bornes physiques et oscillations : Les méthodes de Galerkin discontinu (DG) d'ordre élevé sont sujettes à générer des oscillations parasites près des discontinuités et peuvent violer les bornes physiques (par exemple, pression négative, fractions volumiques en dehors de $[0,1]$), entraînant une rupture numérique.
• Condition d'Abgrall : Le maintien de l'équilibre de pression et de vitesse aux interfaces matérielles isolées est crucial pour éviter des artefacts non physiques, une propriété connue sous le nom de condition d'Abgrall.

2. Méthodologie

Les auteurs proposent un schéma DG à élimination d'oscillations et préservant les bornes (BP-OEDG) intégré à une nouvelle stratégie de décomposition d'opérateurs.

A. Stratégie de décomposition d'opérateurs

Pour contourner les contraintes de stabilité induites par la rigidité, le modèle de Kapila à cinq équations est découplé en deux sous-problèmes :

1. Modèle de transport (Transport à cinq équations) : Contient les lois conservatives pour la masse, la quantité de mouvement et l'énergie, plus l'advection de la fraction volumique (sans le terme source rigide).
   • Discrétisation : Résolu à l'aide d'une méthode DG quasi-conservative (basée sur Zhang et al. [8]) satisfaisant la condition d'Abgrall.
2. Terme source rigide (terme $\kappa$) : Contient uniquement le terme source rigide de l'équation de la fraction volumique.
   • Discrétisation : Résolu à l'aide d'une stratégie implicite adaptative.

B. Solveur implicite pour le terme rigide

• Schéma implicite hybride : Les auteurs proposent un hybride des méthodes d'Euler implicite (premier ordre, inconditionnellement stable) et SDIRK2 (Runge-Kutta implicite diagonalement d'ordre deux).
• Mécanisme adaptatif : Le schéma initie une étape d'Euler implicite. Si la solution prédite reste dans les bornes physiques $[0, 1]$, il procède à la deuxième étape de SDIRK2 pour maintenir une précision temporelle d'ordre deux. Si la prédiction viole les bornes, il dégénère localement en une étape robuste d'Euler implicite.
• Préservation des bornes inconditionnelle (BP) : Les auteurs prouvent rigoureusement que ce traitement implicite est inconditionnellement préservateur des bornes (BP), éliminant efficacement la restriction de pas de temps induite par la rigidité.
• Reconstruction de la divergence de vitesse : Pour restaurer une précision spatiale d'ordre élevé optimale (perdue lors de la différentiation directe des polynômes DG), une reconstruction inspirée du Galerkin discontinu local (LDG) est utilisée pour calculer $\nabla \cdot \mathbf{u}$ au sein du solveur implicite.

C. Limiteurs

• Limiteur à élimination d'oscillations (OE) : Une procédure OE est appliquée après chaque étape de Runge-Kutta. Contrairement aux limiteurs traditionnels (par exemple WENO, TVD) qui nécessitent une décomposition caractéristique, le limiteur OE supprime les oscillations parasites sans décomposer le système, réduisant considérablement la complexité computationnelle pour les écoulements multiphasiques.
• Limiteur préservant les bornes (BP) : Un limiteur d'échelle est appliqué pour garantir que les densités partielles, la pression et les fractions volumiques restent dans l'ensemble admissible physique strict $G_p = \{U : z_1\rho_1 > 0, z_2\rho_2 > 0, p(U) > 0, z_1 \in [0, 1]\}$.

D. Preuves théoriques

• Convexité de l'ensemble admissible : L'article démontre que l'ensemble $G_p$ est convexe sous l'équation d'état (EOS) de Tammann, tandis que l'ensemble plus lâche $G$ (basé sur le carré de la vitesse du son) est non convexe.
• Préservation de la condition d'Abgrall : Les auteurs prouvent que l'ensemble du cadre (décomposition, solveur implicite, limiteur OE et limiteur BP) satisfait strictement la condition d'Abgrall, garantissant que la pression et la vitesse constantes sont préservées aux interfaces matérielles isolées.

3. Contributions clés

1. Nouveau cadre DG à décomposition d'opérateurs : Première application d'une méthode DG d'ordre élevé au modèle de Kapila gaz-liquide qui contourne avec succès les restrictions CFL extrêmes causées par le terme source rigide $\kappa$ via la décomposition d'opérateurs.
2. Solveur implicite inconditionnellement BP : Développement d'un solveur implicite adaptatif (hybride Euler implicite/SDIRK2) garantissant que la fraction volumique reste dans $[0,1]$ indépendamment de la taille du pas de temps.
3. Élimination des oscillations sans décomposition caractéristique : Mise en œuvre d'un limiteur OE qui supprime efficacement les oscillations dans des systèmes multiphasiques complexes sans la décomposition caractéristique coûteuse et complexe requise par les limiteurs traditionnels.
4. Garanties théoriques rigoureuses : Preuve de la convexité de l'ensemble admissible physique $G_p$ et de la préservation de la condition d'Abgrall pour le schéma entièrement discret.
5. Précision d'ordre élevé : Restauration de la précision d'ordre $(K+1)$ optimale pour le terme de divergence via une reconstruction inspirée du LDG au sein du solveur implicite.

4. Résultats numériques

La méthode a été validée par de nombreux cas de référence 1D et 2D :

• Écoulement gazeux-gazeux lisse : Démonstration de taux de convergence d'ordre deux (P1) et d'ordre trois (P2), correspondant aux solutions de référence.
• Interface isolée : Confirmation de la propriété de non-perturbation des champs de pression et de vitesse à une interface gaz-liquide.
• Double détente et problèmes de Riemann : Démonstration de la capacité du schéma à gérer des ondes de détente fortes et des rapports de pression élevés tout en maintenant les bornes physiques (pas de pression négative).
• Tube à choc : Capture réussie des chocs, des discontinuités de contact et des ondes de détente dans un tube à choc gaz-liquide à haute pression.
• Interactions choc-bulle :
  • Choc d'air frappant une bulle d'hélium : Les résultats correspondent aux données expérimentales (Hass & Sturtevant) avec une grande fidélité et sans oscillations parasites.
  • Choc d'eau frappant une bulle d'air : Simulation réussie de la rigidité extrême des interactions eau-air, un cas souvent prohibitif en calcul pour les méthodes DG explicites, démontrant une robustesse supérieure.

5. Importance

Ce travail représente une avancée significative en dynamique des fluides computationnelle pour les écoulements multiphasiques. En combinant la décomposition d'opérateurs avec des méthodes DG d'ordre élevé, il résout le compromis de longue date entre stabilité numérique (gestion de la rigidité) et efficacité computationnelle (autorisation de pas de temps plus grands). La capacité à simuler des interactions gaz-liquide extrêmes (comme les collisions choc-bulle) avec une précision d'ordre élevé et des bornes physiques strictes fait de ce schéma un outil puissant pour des applications en aéronautique, énergie et sciences environnementales où de telles dynamiques multiphasiques complexes se produisent.

L'article reconnaît que, bien que la méthode de décomposition introduise une dégradation locale de l'ordre aux interfaces, elle est suffisante pour la dynamique de compression extrême visée. Les travaux futurs visent à développer un cadre IMEX-OEDG entièrement couplé pour éliminer totalement les erreurs de décomposition.