A linear, fully decoupled, and unconditionally energy-stable SAV-FEM for the Cahn--Hilliard--Navier--Stokes model

Probleemstelling
Het artikel behandelt de numerieke simulatie van het Cahn–Hilliard–Navier–Stokes (CHNS) systeem, een gekoppeld nietlineair partiëel differentiaalvergelingsmodel dat twee-fasen incompressibele stromingen beschrijft. Het systeem brengt aanzienlijke uitdagingen met zich mee voor numerieke discretisatie, waaronder:

1. Nietlineaire Koppeling: Het snelheidsveld en de fasevariabele zijn gekoppeld via convectieve termen in de fasevergelijking en stress-termen in de momentumvergelijking.
2. Nietlineariteit en Stijfheid: Het systeem bevat een nietlineair energiepoteential (Ginzburg–Landau vrije energie) en stijfheid veroorzaakt door kleine interface-breedteparameters.
3. Energie-stabiliteit: Het ontwerpen van schema's die de fysieke energie-dissipatiewet op discreet niveau behouden, is moeilijk, vooral wanneer men streeft naar lineaire en volledig ontkoppelde formuleringen.

Bestaande methoden vertrouwen vaak op multiple scalar auxiliary variables (MSAV) of combinaties van hulpvariabelen met zero-energy-contribution (ZEC) technieken om deze nietlineariteiten te behandelen. Hoewel effectief, verhogen deze benaderingen de computationele complexiteit doordat het nietlineaire systeem meerdere keren moet worden gedecomposeerd en meerdere hulp-gewone differentiaalvergelijkingen moeten worden opgelost.

Methodologie
De auteurs stellen een lineair, volledig ontkoppeld en onvoorwaardelijk energie-stabiel volledig discreet eindige elementen-schema voor, gebaseerd op de Scalar Auxiliary Variable (SAV) benadering in combinatie met een Implicit–Explicit (IMEX) Euler tijd-discretisatie.

• Herformulering: Een enkele scalaire hulpvariabele, $r(t) = \sqrt{E_1[\phi(t)]}$, wordt geïntroduceerd om alle nietlineaire termen die voortvloeien uit zowel de Cahn–Hilliard als de Navier–Stokes vergelijkingen simultaan te herformuleren. In tegenstelling tot standaard SAV-benaderingen, introduceren de auteurs een nieuwe update-vergelijking voor deze enkele variabele die specifieke correctietermen bevat. Dit maakt de behandeling van meerdere nietlineariteiten mogelijk zonder ze in afzonderlijke hulpvariabelen te splitsen.
• Discretisatie:
  • Temporeel: Het IMEX Euler-schema behandelt lineaire termen impliciet en nietlineaire termen expliciet.
  • Ruimtelijk: Een standaard eindige elementenmethode (FEM) wordt gebruikt. De snelheid- en drukruimtes voldoen aan de inf-sup conditie (bijv. MINI-elementen of Taylor-Hood elementen), en het faseveld gebruikt continue stuksgewijze lineaire elementen.
• Algoritmische Structuur: Het resulterende volledig discrete schema deelt zich op in drie onafhankelijke subproblemen per tijdstap:
  1. Een lineair Cahn–Hilliard-achtig systeem.
  2. Een lineair Stokes-achtig systeem.
  3. Een scalaire kwadratische algebraïsche vergelijking voor de hulpvariabele.
     Cruciaal is dat de snelheid en druk worden opgelost als een gekoppeld Stokes-systeem op elke stap (in plaats van een druk-correctie projectie), wat de foutanalyse vereenvoudigt.

Belangrijkste Bijdragen en Resultaten

1. Efficiëntie van Enkele Hulpvariabele: De primaire bijdrage is de demonstratie dat een enkele scalaire hulpvariabele, gekoppeld aan een specifieke update-strategie, voldoende is om alle nietlineariteiten in het CHNS-systeem te behandelen. Dit vermindert de computationele overhead vergeleken met MSAV of SAV-ZEC methoden die meerdere hulpvariabelen vereisen.
2. Onvoorwaardelijke Energie-stabiliteit: De auteurs bewijzen dat het voorgestelde schema voldoet aan een onvoorwaardelijke discrete energie-dissipatiewet. De gemodificeerde discrete energie neemt bij elke tijdstap monotoon af, ongeacht de grootte van de tijdstap.
3. Optimale Foutschattingen: Onder standaard regulariteitsveronderstellingen ($u, \phi \in H^2$), stelt het artikel optimale $L^2$ foutschattingen vast voor de volledig discrete benadering.
   • Het faseveld, de chemische potentiaal en de druk bereiken optimale convergentiegraden.
   • De snelheid bereikt een optimale $L^2$ convergentie.
   • Een belangrijke technische prestatie is het behandelen van de niet-globaal Lipschitz-aard van de nietlineaire potentiaal $g(\phi) = F'(\phi)$. De auteurs bewijzen de uniforme $L^\infty$ begrensdheid van de numerieke oplossing, waardoor de foutanalyse kan voortgaan zonder restrictieve aannames over de vrije energie.
4. Numerieke Verificatie: Numerieke experimenten bevestigen de theoretische convergentiegraden (waarbij $O(h^2)$ in $L^2$ en $O(h)$ in $H^1$ wordt waargenomen voor ruimtelijke discretisatie) en demonstreren de monotone afname van de gemodificeerde discrete energie in diverse testgevallen, inclusief fase-scheiding en coarsening-dynamica.

Betekenis
Het artikel beweert dat de benadering een efficiënter en eenvoudiger alternatief biedt voor bestaande SAV-gebaseerde methoden voor het CHNS-systeem. Door gebruik te maken van een enkele hulpvariabele en een verenigde update-strategie, vermijdt de methode de complexiteit van het meervoudig decomponeren van het systeem. Bovendien vereenvoudigt de directe oplossing van het Stokes-subsysteem (zonder projectiestappen) de theoretische foutanalyse, wat optimale foutschattingen mogelijk maakt onder zwakkere regulariteitsveronderstellingen ($\phi \in H^2$) vergeleken met sommige projectie-gebaseerde methoden die hogere regulariteit vereisen. Het werk biedt een robuust kader voor het simuleren van twee-fasen stromingen dat zowel computationeel efficiënt (lineair en ontkoppeld) als theoretisch rigoureus is (onvoorwaardelijk stabiel met optimale foutgrenzen).