Sharp-interface Simulations of Energetic Multiphase Flows with Large Density and Viscosity Ratios

Questo articolo propone un framework di flussi di Momento Sincronizzati tra Regione Donatrice e Regione Ricevente (SynDRoM) combinato con un limitatore di viscosità per migliorare la robustezza numerica e la fedeltà fisica delle simulazioni a interfaccia netta per flussi multifase energetici con elevati rapporti di densità e viscosità.

L'essenziale

Il Quadro Generale: Simulare un Mare in Tempesta

Immagina di voler creare una simulazione al computer di una violenta tempesta oceanica. Vuoi vedere come si infrangono le onde, come l'aria viene risucchiata nell'acqua e come si formano le bolle. Questo è complicato perché l'acqua è pesante e densa, mentre l'aria è leggera e sottile. In termini fisici, hanno una massiccia differenza di "densità".

Quando i computer cercano di simulare questo fenomeno, spesso vanno in crash o producono risultati strani e impossibili (come l'acqua che improvvisamente diventa un fantasma o l'aria che attraversa l'acqua come un proiettile). Questo articolo introduce un nuovo insieme di regole (algoritmi) per rendere queste simulazioni stabili, accurate e fisicamente realistiche, anche quando le onde si infrangono violentemente.

Il Problema: Il "Fantasma" e lo "Shock"

Gli autori spiegano che i vecchi metodi per simulare questi flussi presentano due difetti principali:

1. Il Problema del "Fantasma" (Penetrazione della Velocità):
   Immagina un camion pesante (acqua) e una piuma (aria) che si muovono l'uno accanto all'altro. Nelle vecchie simulazioni, il "vento" della piuma poteva a volte spingere il camion all'indietro, o il camion poteva spingere la piuma attraverso il proprio corpo. Questo è chiamato "penetrazione della velocità". Crea forme non fisiche e finte nell'acqua, come un "corno del diavolo" che spunta da un'onda.

2. Il Problema dello "Shock" (Picchi di Momento):
   Per risolvere il problema del fantasma, gli scienziati hanno provato un nuovo metodo chiamato CMOM (Consistenza Massa-Momento). È come tenere un registro rigoroso di quanto "slancio" (momento) ogni goccia d'acqua possiede. Tuttavia, questo metodo ha un effetto collaterale. Quando un piccolo frammento di acqua pesante si muove in una cella piena d'aria, la matematica si confonde. È come dividere un numero enorme per un numero minuscolo, risultando in un picco di velocità enorme e impossibile. Questo crea "blob di velocità": sacche d'aria finte che si muovono a velocità supersoniche e che non dovrebbero esistere.

La Soluzione: Il Metodo "SynDRoM"

Gli autori propongono una nuova soluzione chiamata SynDRoM (Flusso di Momento della Regione Donatrice Sincronizzata). Ecco come funziona, usando un'analogia:

L'Analogia: Il Nastro Trasportatore in Movimento
Immagina un nastro trasportatore che trasporta scatole.

• Il Vecchio Modo: Conti le scatole (massa) e il peso delle scatole (momento) separatamente. Se una scatola si muove, potresti accidentalmente contare il suo peso in un punto dove la scatola non è ancora effettivamente arrivata. Questo causa lo "shock" o il "picco" di velocità.
• Il Modo SynDRoM: Questo metodo agisce come una squadra sincronizzata. Prima di spostare il peso, guardi esattamente quale parte del nastro trasportatore il peso sta provenendo.
  • Chiede: "Se sto muovendo questo specifico pezzo di aria, esattamente a quale pezzo di momento è attaccato?"
  • Assicura che il momento venga spostato solo se la massa è effettmente lì per trasportarlo.
  • Il Risultato: Niente più picchi di velocità falsi. L'aria rimane lenta e l'acqua rimane pesante, proprio come nella realtà. La simulazione rimane fluida e non "esplode".

Il Secondo Problema: La Viscosità "Scivolosa"

L'articolo affronta anche un secondo problema: la Viscosità (quanto un fluido è denso o appiccicoso).

• Il Problema: L'acqua è appiccicosa; l'aria è scivolosa. Quando si mescolano in un confine netto (come un'onda che si infrange), il computer cerca di indovinare la "viscosità" nel mezzo. Se indovina male, la matematica diventa instabile, come cercare di bilanciare una matita sulla sua punta.
• La Soluzione: Gli autori introducono un Limitatore di Viscosità.
  • L'Analogia: Immagina un segnale stradale di limite di velocità. Anche se la matematica prova a calcolare una "appiccicosità" che farebbe muovere il fluido in modo impossibile (instabile), il limitatore dice: "No, non puoi andare più veloce della velocità del fluido più sottile qui presente". Esso blocca il calcolo per evitare che la simulazione vada in crash, senza cambiare la fisica reale dell'acqua o dell'aria.

La Prova: Funziona?

Gli autori hanno testato le loro nuove regole in tre modi:

1. La Rottura della Diga (Dam Break): Hanno simulato un muro d'acqua che crolla.
   • Vecchi metodi: L'acqua appariva distorta con picchi falsi.
   • SynDRoM: L'acqua si è infranta naturalmente e l'aria non è stata risucchiata nell'acqua in modi strani.
2. L'Instabilità di Kelvin-Helmholtz: Questa è quando il vento soffia sull'acqua, creando onde rotanti (come le nuvole).
   • Risultato: La simulazione ha mostrato correttamente le onde che si arrotolavano e crescevano, senza che il computer aggiungesse energia falsa o smorzasse le onde. Ha dimostrato che il metodo rispetta le leggi della fisica.
3. L'Onda che si Infrange: Hanno simulato un'enorme onda diagonale che si infrange.
   • Risultato: L'onda si è infranta, ha schizzato e ha creato schiuma proprio come un vero oceano. L'energia totale del sistema è rimasta bilanciata (non è sparita magicamente né è esplosa). Anche aggiungendo "appiccicosità" (viscosità), la simulazione è rimasta stabile.

In Sintesi

Questo articolo presenta un nuovo "poliziotto del traffico" per le simulazioni al computer di acqua e aria.

• Impedisce all'aria di passare come un fantasma attraverso l'acqua.
• Impedisce all'acqua di creare impossibili picchi di velocità.
• Mantiene i calcoli della "appiccicosità" impedendo che la matematica si rompa.

Sincronizzando esattamente cosa si muove con dove si sta muovendo, gli autori hanno creato uno strumento di simulazione molto più robusto e affidabile per studiare eventi oceanici violenti, come quelli di cui gli ingegneri navali hanno bisogno per la progettazione delle navi.

Sommario tecnico

Sintesi Tecnica: Simulazioni a interfaccia netta di flussi multifase energetici con elevati rapporti di densità e viscosità

Definizione del Problema
La simulazione di flussi multifase energetici con rapporti di densità estremi (ad esempio, sistemi aria-acqua con rapporti $\sim 10^3$ o $10^4$) rimane una sfida significativa nell'idrodinamica navale. Le tradizionali formulazioni basate sulla velocità soffrono spesso di divergenza numerica a causa del disallineamento tra i gradienti di pressione e le discontinuità di densità, portando a trasferimenti di quantità di moto spuri e a esplosioni di energia. Sebbene i framework di trasporto a Consistenza Massa-Quantità di Moto (CMOM) migliorino la robustezza garantendo la conservazione della quantità di moto e la conservazione semi-discreta dell'energia, essi introducono nuove difficoltà quando combinati con metodi a interfaccia netta. Nello specifico, il CMOM può generare shock di quantità di moto severi, provocando oscillazioni non fisiche della velocità e superamenti (overshoot) che violano la condizione di Variazione Totale Limitata (TVD), quando le interfacce nette attraversano i confini delle celle. Inoltre, una stima impropria della viscosità vicino alle interfacce, in particolare dividendo la viscosità dinamica per la densità su griglie sfasate (staggered), può introdurre instabilità numeriche a passi temporali finiti, compromettendo la robustezza del solutore.

Metodologia
Il documento propone una serie di schemi numerici preservanti la fisica, progettati per conciliare la fedeltà fisica con la robustezza numerica per flussi con elevati rapporti di densità.

1. Flusso di Quantità di Moto con Regione di Donatore Sincronizzata (SynDRoM):
   Per affrontare le oscillazioni di velocità inerenti al CMOM standard, gli autori introducono il metodo SynDRoM. Questo metodo impone la monotonicità del campo di velocità trasportato sincronizzando il flusso di quantità di moto con il flusso di massa.

• Meccanismo: A differenza degli esistenti limitatori di pendenza sensibili all'interfaccia che si concentrano esclusivamente sulla ricostruzione della velocità, SynDRoM ricostruisce una coppia distribuzione densità-velocità all'interno di una cella per soddisfare i vincoli medi di cella e i limiti TVD. Esso evolve poi la quantità di moto utilizzando un concetto di "regione donatrice" determinato strettamente dal flusso di massa. La velocità addetta viene calcolata come una media ponderata su questa specifica regione donatrice, garantendo che la velocità aggiornata rimanga limitata.
• Implementazione: Il metodo opera all'interno di un framework di splitting direzionale coerente con l'advezione cVOF (Volume of Fluid conservativo). Utilizza uno schema di integrazione temporale Runge-Kutta di secondo ordine (RK2) pesato sulla densità per gestire l'evoluzione non lineare della velocità rispetto alla densità e alla quantità di moto.

2. Limitatore di Viscosità Cinematica Limitata:
   Per affrontare le instabilità derivanti dai salti di viscosità alle interfacce, il documento introduce un limitatore di viscosità.

• Meccanismo: Gli autori identificano che l'interpolazione standard (aritmetica o armonica) della viscosità dinamica ai vertici della griglia sfasata, seguita dalla divisione per la densità, può produrre una viscosità cinematica effettiva illimitata, violando i vincoli CFL viscosi. Il limitatore proposto vincola la viscosità dinamica effettiva utilizzata nei termini di differenza incrociata tramite la densità minima delle celle di quantità di moto circostanti e un limite superiore basato sull'occupazione di fase della viscosità cinematica. Ciò assicura che la viscosità cinematica effettiva rimanga entro un intervallo di stabilità.

Contributi Chiave

• Algoritmo SynDRoM: Un metodo generalizzato per il CMOM che elimina le oscillazioni spurie della velocità senza sacrificare la conservazione della quantità di moto o la nitidezza dell'interfaccia. Tiene esplicitamente conto della distribuzione della quantità di moto e della corrispondente regione donatrice basata sul flusso di massa, evitando la necessità di una completa ricostruzione dell'interfaccia su griglie sfasate.
• Stabilizzazione della Viscosità: Un limitatore di viscosità semplice ma efficace che previene la divergenza numerica nei flussi con elevi rapporti di viscosità, limitando la viscosità cinematica effettiva, preservando così la stabilità senza compromettere la fedeltà fisica.
• Coerenza di Splitting Direzionale: Una formulazione che assicura che gli aggiornamenti della quantità di moto rimangano coerenti con l'advezione di massa a splitting direzionale, prevenendo flussi fittizi che potrebbero violare il limite della velocità.

Risultati
I metodi proposti sono stati validati attraverso una serie di casi di test:

• Dam Break 2D: La formulazione SynDRoM ha risolto con successo la deformazione dell'interfaccia ed eliminato le strutture di velocità spurie (come le distorsioni "devil-horn" viste nelle formulazioni di velocità e i "blob di velocità" nel CMOM originale) mantenendo la limitatezza dell'energia.
• Advezione di Goccia 1D: In un test quasi monodimensionale con un elevato rapporto di densità ($10^{-3}$), SynDRoM ha mantenuto la monotonicità della velocità e le proprietà TVD, mentre il CMOM originale e le varianti Favre-averaged hanno prodotto overshoot non fisici.
• Instabilità di Kelvin-Helmholtz: Lo schema ha catturato correttamente la crescita lineare delle perturbazioni e la successiva evoluzione non lineare (accoppiamento di vortici, roll-up) senza smorzamento artificiale o esplosione di energia. La quantità di moto orizzontale è stata conservata entro la tolleranza del solutore e le regioni ad alta velocità sono state correttamente confinate alla fase del fluido leggero.
• Flusso di Poiseuille Stratificato: Il limitatore di viscosità ha impedito con successo la divergenza numerica osservata con la media aritmetica standard e ha corretto i gradienti di velocità eccessivi visti con la media armonica, fornendo soluzioni stabili e fisicamente accurate.
• Onda Stazionaria Rompente 3D: Le simulazioni di onde rompenti stazionarie diagonali hanno dimostrato la capacità del solutore di gestire l'estrema non linearità, la formazione di getti e l'intrappolamento dell'aria. Il metodo ha mantenuto l'energia meccanica limitata anche in assenza di viscosità esplicita, con la dissipazione dell'energia che avviene naturalmente attraverso il meccanismo iLES (Large Eddy Simulation implicita) man mano che il flusso diventa irregolare. Gli studi di convergenza della griglia hanno confermato tendenze coerenti nell'evoluzione dell'energia e nei tassi di dissipazione.

Significatività
Il documento sostiene che le modifiche proposte agli schemi di advezione e viscosità migliorano significativamente la fedeltà e la robustezza delle simulazioni di flussi multifase, in particolare per i flussi a superficie libera energetici incontrati nell'idrodinamica navale. Garantendo che i metodi CMOM soddisfino la condizione TVD vicino alle celle di interfaccia in svuotamento e limitando la viscosità cinematica, l'approccio rimuove gli artefatti numerici che potrebbero innescare falsi eventi di rottura. Gli autori sottolineano che SynDRoM e il limitatore di viscosità sono ortogonali al metodo specifico di cattura dell'interfaccia utilizzato (dimostrato qui con il VOF geometrico), rendendo l'approccio ampiamente applicabile. Il lavoro fornisce una via per simulare flussi complessi con elevati rapporti di densità con coerenza fisica, conservazione della quantità di moto e limitatezza empirica dell'energia, affrontando sfide di lunga data nel settore.