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

Sintesi Tecnica: Una famiglia di metodi del secondo ordine, lineari e incondizionatamente stabili per le equazioni di Cahn-Hilliard-Navier-Stokes

Problematica
La simulazione numerica accurata di dinamiche dei fluidi multifase complesse, specificamente quelle modellate dalle equazioni di Cahn-Hilliard-Navier-Stokes (CHNS), presenta sfide significative dovute alla rigidità matematica dell'operatore di Cahn-Hilliard del quarto ordine e al suo accoppiamento non lineare con le equazioni di Navier-Stokes incomprimibili. Gli standard schemi di avanzamento temporale espliciti affrontano severi vincoli sul passo temporale (ad esempio, $\Delta t \sim O(\Delta x^4)$) per mantenere la stabilità. Sebbene i metodi completamente impliciti o di scissione convessa offrano una stabilità energetica incondizionata, essi richiedono spesso la risoluzione di grandi sistemi non lineari accoppiati ad ogni passo temporale, aumentando il costo computazionale. È necessario un algoritmo che combini l'accuratezza temporale del secondo ordine, la stabilità incondizionata e la risolvibilità lineare senza sacrificare le proprietà di conservazione fisica del sistema.

Metodologia
Gli autori propongono una famiglia di metodi IMEX (implicito-esplicito) del secondo ordine, lineari e incondizionatamente stabili per flussi bifase a densità accoppiata. Lo schema semi-discreto integra tre componenti primarie:

1. Estrapolazione dei termini non lineari: I termini convettivi e di accoppiamento non lineari sono trattati esplicitamente utilizzando quantità estrapolate (indicate con $H_{n+\theta}$), evitando la necessità di iterazioni non lineari.
2. Formulazione con variabile ausiliaria: Il termine di energia libera non lineare è gestito tramite una variabile ausiliaria $q$ (derivata dal potenziale a doppia buca $F(\phi)$), riformulando il problema per mantenere la linearità.
3. Regolarizzazione della curvatura temporale: Viene introdotto un meccanismo di stabilizzazione controllato da un parametro $\epsilon \geq 0$. Questo comporta un operatore di interpolazione specifico $J^\epsilon_{n+\theta}$ che incorpora la curvatura temporale discreta della soluzione, un concetto adattato da precedenti lavori sulla regolarizzazione di Navier-Stokes.

L'algoritmo risultante (Algoritmo 2.1) richiede la risoluzione solo di sistemi lineari a coefficienti costanti ad ogni passo temporale. Lo schema è parametrizzato da $\theta \in (1/2, 1]$ e $\epsilon \geq 0$.

Contributi Chiave e Risultati Teorici

• Stabilità Incondizionata: Il documento stabilisce una rigorosa stima dell'energia discreta provando che lo schema proposto è incondizionatamente stabile nel lungo periodo per $\theta \in (1/2, 1]$ e $\epsilon \geq 0$. La prova utilizza matrici definite positive simmetriche per definire le norme discrete e dimostra che lo schema soddisfa una legge di dissipazione dell'energia discreta senza restrizioni sul passo temporale.
• Risolvibilità Lineare: A differenza degli approcci completamente impliciti o di scissione convessa che possono richiedere solver non lineari, questo metodo produce sistemi lineari ad ogni passo temporale, riducendo significativamente la complessità computazionale.
• Accuratezza del Secondo Ordine: L'analisi teorica e gli esperimenti numerici confermano che il metodo mantiene un'accuratezza temporale del secondo ordine.

Risultati Numerici
Gli autori validano il metodo attraverso una serie di computazioni di benchmark:

• Analisi di Convergenza: Utilizzando il Metodo delle Soluzioni Manifatturate (MMS), lo schema dimostra una convergenza temporale approssimativamente del secondo ordine per velocità, pressione e variabili di campo di fase attraverso diverse configurazioni di elementi finiti ($P2-P1-P2-P2$ e $P3-P2-P3-P3$).
• Robustezza e Regolarizzazione: Nel problema della cavità guidata da una parete mobile (lid-driven cavity) a due fasi, lo studio esamina l'interazione tra $\theta$ ed $\epsilon$. Si osserva che per valori di $\theta$ vicini al limite di Crank-Nicolson (ad esempio, $\theta=0.51$) con $\epsilon=0$, lo schema può diventare instabile a causa di oscillazioni temporali debolmente smorzate. L'introduzione di una regolarizzazione della curvatura positiva ($\epsilon > 0$) smorza efficacemente queste oscillazioni, ripristinando l'limitatezza e la robustezza anche per passi temporali più ampi.
• Benchmark Fisici:
  • Decomposizione Spinodale: Il metodo cattura accuratamente la dissipazione dell'energia e la conservazione della massa in simulazioni a lungo termine ($T=100$), mostrando una separazione di fase e un ingrossamento (coarsening) efficaci.
  • Rilassamento di una Goccia: L'algoritmo risolve con successo il moto dell'interfaccia guidato dalla tensione superficiale, passando da una goccia quadrata a un equilibrio circolare preservando la massa.
  • Cavità Guidata e Instabilità di Rayleigh-Taylor: Lo schema gestisce forti sforzi di taglio, deformazioni interfacciali e instabilità guidate dalla densità (incluse le rotazioni di Kelvin-Helmholtz in regimi a bassa viscosità) senza oscillazioni spurie visibili.

Significatività
Il documento sostiene che la famiglia di metodi proposta offre un'alternativa pratica ed efficiente alle discretizzazioni completamente non lineari per i sistemi CHNS. Decoppiando le non linearità attraverso l'estrapolazione e le variabili ausiliarie, e stabilizzando la discretizzazione temporale tramite la regolarizzazione della curvatura, il metodo raggiunge un equilibrio tra accuratezza del secondo ordine, stabilità incondizionata e costo computazionale lineare. Gli autori sottolineano che lo schema preserva le principali strutture fisiche del modello, inclusa la dissipazione dell'energia e la conservazione della massa, rendendolo adatto alla simulazione di flussi interfacciali complessi in vari regimi.