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

Technische Samenvatting: Een familie van tweede-orde, lineaire, onvoorwaardelijk stabiele methoden voor de Cahn-Hilliard-Navier-Stokes-vergelijkingen

Probleemstelling
De nauwkeurige numerieke simulatie van complexe meerfasige fluïdumdynamica, specifiek die gemodelleerd door de Cahn-Hilliard-Navier-Stokes (CHNS) vergelijkingen, brengt aanzienlijke uitdagingen met zich mee vanwege de wiskundige stijfheid van de vierde-orde Cahn-Hilliard-operator en de nietlineaire koppeling met de incompressibele Navier-Stokes-vergelijkingen. Standaard expliciete tijdstap-schema's kampen met strikte beperkingen voor de tijdstap (bijv. $\Delta t \sim O(\Delta x^4)$) om stabiliteit te handhaven. Hoewel volledig impliciete of convex-splitting methoden onvoorwaardelijke energie-stabiliteit bieden, vereisen ze vaak het oplossen van grote, gekoppelde nietlineaire systemen bij elke tijdstap, wat de computationele kosten verhoogt. Er is behoefte aan algoritmen die tweede-orde temporele nauwkeurigheid, onvoorwaardelijke stabiliteit en lineaire oplosbaarheid combineren zonder de fysische conserverende eigenschappen van het systeem op te offeren.

Methodologie
De auteurs stellen een familie van tweede-orde, lineaire, onvoorwaardelijk stabiele IMEX (implicit-explicit) methoden voor twee-fasige stromingen met gelijke dichtheid voor. Het semi-discrete schema integreert drie primaire componenten:

1. Extrapolatie van nietlineaire termen: Nietlineaire convectieve en koppelingstermen worden expliciet behandeld met behulp van geëxtrapoleerde grootheden (aangeduid als $H_{n+\theta}$), waardoor nietlineaire iteraties worden vermeden.
2. Auxiliaire variabele formulering: De nietlineaire vrije-energie term wordt afgehandeld via een hulpvariabele $q$ (afgeleid van de double-well potentiaal $F(\phi)$), wat het probleem herformuleert om lineariteit te behouden.
3. Temporele-kromming regularisatie: Een stabilisatiemechanisme gecontroleerd door een parameter $\epsilon \geq 0$ wordt geïntroduceerd. Dit omvat een specifieke interpolatie-operator $J^\epsilon_{n+\theta}$ die de discrete temporele kromming van de oplossing incorporeert, een concept dat is aangepast van eerder werk op het gebied van Navier-Stokes regularisatie.

Het resulterende algoritme (Algoritme 2.1) vereist bij elke tijdstap enkel het oplossen van lineaire systemen met constante coëfficiënten. Het schema is geparametriseerd door $\theta \in (1/2, 1]$ en $\epsilon \geq 0$.

Belangrijkste bijdragen en theoretische resultaten

• Onvoorwaardelijke stabiliteit: Het artikel stelt een rigoureuze discrete energie-inschatting vast die bewijst dat het voorgestelde schema onvoorwaardelijk langetermijnstabiel is voor $\theta \in (1/2, 1]$ en $\epsilon \geq 0$. Het bewijs maakt gebruik van symmetrische positief-definiete matrices om discrete normen te definiëren en toont aan dat het schema voldoet aan een discrete energie-dissipatiewet zonder beperkingen voor de tijdstap.
• Lineaire oplosbaarheid: In tegenstelling tot volledig impliciete of convex-splitting benaderingen die mogelijk nietlineaire solvers vereisen, levert deze methode bij elke tijdstap lineaire systemen op, wat de computationele complexiteit aanzienlijk vermindert.
• Tweede-orde nauwkeurigheid: Theoretische analyse en numerieke experimenten bevestigen dat de methode een tweede-orde temporele nauwkeurigheid behoudt.

Numerieke resultaten
De auteurs valideren de methode via een reeks benchmark-berekeningen:

• Convergentieanalyse: Met behulp van de Method of Manufactured Solutions (MMS) vertoont het schema een benadering van tweede-orde temporele convergentie voor snelheid, druk en faseveld-variabelen over verschillende eindige elementen-configuraties ($P2-P1-P2-P2$ en $P3-P2-P3-P3$).
• Robuustheid en regularisatie: In het twee-fasige lid-driven cavity probleem onderzoekt de studie de interactie tussen $\theta$ en $\epsilon$. Er is waargenomen dat voor $\theta$-waarden nabij de Crank-Nicolson limiet (bijv. $\theta=0.51$) met $\epsilon=0$, het schema instabiel kan worden door zwak gedempte temporele oscillaties. De introductie van positieve kromming-regularisatie ($\epsilon > 0$) dempt deze oscillaties effectief, wat de begrensdheid en robuustheid herstelt, zelfs voor grotere tijdstappen.
• Fysische benchmarks:
  • Spinodale decompositie: De methode legt de energie-dissipatie en massaconservering accuraat vast tijdens langdurige simulaties ($T=100$), waarbij effectieve fase-scheiding en coarsening zichtbaar zijn.
  • Druppelrelaxatie: Het algoritme lost de oppervlaktespanning-gedreven interfacebeweging succesvol op, waarbij een vierkante druppel transformeert naar een circulaire evenwichtstoestand met behoud van massa.
  • Lid-Driven Cavity & Rayleigh-Taylor Instabiliteit: Het schema handelt sterke afschuiving, interface-deformatie en dichtheidsgedreven instabiliteiten (inclusief Kelvin-Helmholtz roll-up in regimes met lage viscositeit) af zonder zichtbare spuria true oscillaties.

Betekenis
Het artikel beweert dat de voorgestelde familie van methoden een praktische en efficiënte alternatief biedt voor volledig nietlineaire discretisaties voor CHNS-systemen. Door de nietlineariteiten te ontkoppelen via extrapolatie en hulpvariabelen, terwijl de temporele discretisatie wordt gestabiliseerd via kromming-regularisatie, bereikt de methode een balans tussen tweede-orde nauwkeurigheid, onvoorwaardelijke stabiliteit en lineaire computationele kosten. De auteurs benadrukken dat het schema de belangrijkste fysische structuren van het model behoudt, inclus_ief energie-dissipatie en massaconservering, wat het geschikt maakt voor het simuleren van complexe interface-stromingen in diverse regimes.