A family of second order, linear, unconditionally stable methods for the Cahn-Hilliard-Navier-Stokes equations

Résumé Technique : Une famille de méthodes de second ordre, linéaires et inconditionnellement stables pour les équations de Cahn-Hilliard-Navier-Stokes

Énoncé du Problème
La simulation numérique précise de la dynamique des fluides multiphasiques complexes, spécifiquement celles modélisées par les équations de Cahn-Hilliard-Navier-Stokes (CHNS), présente des défis importants en raison de la rigidité mathématique de l'opérateur de Cahn-Hilliard d'ordre quatrième et de son couplage non linéaire avec les équations de Navier-Stokes incompressibles. Les schémas classiques de pas de temps explicites font face à de sévères restrictions sur le pas de temps (par exemple, $\Delta t \sim O(\Delta x^4)$) pour maintenir la stabilité. Bien que les méthodes entièrement implicites ou de division convexe (convex-splitting) offrent une stabilité énergétique inconditionnelle, elles nécessitent souvent la résolution de grands systèmes non linéaires couplés à chaque pas de temps, ce qui augmente le coût computationnel. Il existe un besoin pour des algorithmes combinant une précision temporelle de second ordre, une stabilité inconditionnelle et une résolubilité linéaire sans sacrifier les propriétés de conservation physique du système.

Méthodologie
Les auteurs proposent une famille de méthodes IMEX (implicite-explicite) de second ordre, linéaires et inconditionnellement stables pour les écoulements biphasiques à densité égale. Le schéma semi-discret intègre trois composantes primaires :

1. Extrapolation des termes non linéaires : Les termes de convection et de couplage non linéaires sont traités explicitement à l'aide de quantités extrapolées (notées $H_{n+\theta}$), évitant ainsi le recours à des itérations non linéaires.
2. Formulation par variable auxiliaire : Le terme d'énergie libre non linéaire est géré via une variable auxiliaire $q$ (dérivée du potentiel de double puits $F(\phi)$), reformulant le problème pour maintenir la linéarité.
3. Régularisation de la courbure temporelle : Un mécanisme de stabilisation contrôlé par un paramètre $\epsilon \geq 0$ est introduit. Cela implique un opérateur d'interpolation spécifique $J^\epsilon_{n+\theta}$ qui incorpore la courbure temporelle discrète de la solution, un concept adapté de travaux antérieurs sur la régularisation de Navier-Stokes.

L'algorithme résultant (Algorithme 2.1) nécessite la résolution de seuls systèmes linéaires à coefficients constants à chaque pas de temps. Le schéma est paramétré par $\theta \in (1/2, 1]$ et $\epsilon \geq 0$.

Contributions Clés et Résultats Théoriques

• Stabilité Inconditionnelle : L'article établit une estimation d'énergie discrète rigoureuse prouvant que le schéma proposé est inconditionnellement stable sur le long terme pour $\theta \in (1/2, 1]$ et $\epsilon \geq 0$. La preuve utilise des matrices définies positives symétriques pour définir les normes discrètes et démontre que le schéma satisfait une loi de dissipation d'énergie discrète sans restrictions sur le pas de temps.
• Résolubilité Linéaire : Contrairement aux approches entièrement implicites ou de division convexe qui peuvent nécessiter des solveurs non linéaires, cette méthode produit des systèmes linéaires à chaque pas de temps, réduisant considérablement la complexité computationnelle.
• Précision de Second Ordre : L'analyse théorique et les expériences numériques confirment que la méthode conserve une précision temporelle de second ordre.

Résultats Numériques
Les auteurs valident la méthode à travers une série de calculs de référence :

• Analyse de Convergence : En utilisant la Méthode des Solutions Manufacturées (MMS), le schéma démontre une convergence temporelle d'environ second ordre pour la vitesse, la pression et les variables de champ de phase à travers différentes configurations d'éléments finis ($P2-P1-P2-P2$ et $P3-P2-P3-P3$).
• Robustesse et Régularisation : Dans le problème de la cavité entraînée par un couvercle (lid-driven cavity) à deux phases, l'étude examine l'interaction entre $\theta$ et $\epsilon$. Il est observé que pour des valeurs de $\theta$ proches de la limite de Crank-Nicolson (par exemple, $\theta=0,51$) avec $\epsilon=0$, le schéma peut devenir instable en raison d'oscillations temporelles faiblement amorties. L'introduction d'une régularisation de courbure positive ($\epsilon > 0$) amortit efficacement ces oscillations, restaurant la bornitude et la robustesse, même pour des pas de temps plus importants.
• Benchmarks Physiques :
  • Décomposition Spinodale : La méthode capture avec précision la dissipation d'énergie et la conservation de la masse sur des simulations de longue durée ($T=100$), montrant une séparation de phase et un mûrissement efficaces.
  • Relaxation de Gouttelette : L'algorithme résout avec succès le mouvement d'interface piloté par la tension superficielle, faisant passer une gouttelette carrée à un équilibre circulaire tout en préservant la masse.
  • Cavité Entraînée par un Couvercle & Instabilité de Rayleigh-Taylor : Le schéma gère les cisaillements forts, la déformation d'interface et les instabilités pilotées par la densité (incluant le roulement de Kelvin-Helmholtz dans les régimes de faible viscosité) sans oscillations parasites visibles.

Signification
L'article affirme que la famille de méthodes proposée offre une alternative pratique et efficace aux discrétisations entièrement non linéaires pour les systèmes CHNS. En découplant les non-linéarités par extrapolation et variables auxiliaires, tout en stabilisant la discrétisation temporelle via la régularisation de la courbure, la méthode atteint un équilibre entre précision de second ordre, stabilité inconditionnelle et coût de calcul linéaire. Les auteurs soulignent que le schéma préserve les structures physiques principales du modèle, notamment la dissipation d'énergie et la conservation de la masse, ce qui le rend adapté à la simulation d'écoulements interfaciaux complexes dans divers régimes.