cfdmfFTFoam: A front-tracking solver for multiphase flows on general unstructured grids in OpenFOAM

Il documento presenta cfdmfFTFoam, un nuovo solver OpenFOAM che integra il metodo di tracciamento dell'interfaccia (FTM) su griglie non strutturate generali, colmando la carenza di software open-source per questo approccio e offrendo un'ampia gamma di algoritmi per simulazioni di flussi multifase.

L'essenziale

🌊 Il "Gestore di Fiumi e Bolle" per il Computer: cfdmfFTFoam

Immagina di voler simulare al computer cosa succede quando una goccia d'olio cade nell'acqua, o quando una bolla d'aria sale in un tubo. È un problema complicato perché l'interfaccia tra i due fluidi (il confine dove si toccano) si muove, si deforma, si allunga e a volte si spezza.

Per anni, i computer hanno usato due metodi principali per disegnare queste interfacce:

1. Il metodo "Sfocato" (VOF): Come se disegnassi il confine con un pennarello a punta grossa su un foglio a quadretti. È veloce, ma il confine diventa un po' sfocato e perde precisione quando si muove velocemente.
2. Il metodo "Preciso" (Front-Tracking o FTM): Come se avessi un filo elastico invisibile che traccia esattamente il contorno della goccia. Questo filo è fatto di piccoli triangoli che si muovono indipendentemente dai quadretti del foglio. È molto più preciso, ma è anche molto difficile da programmare, specialmente se il foglio non è fatto di quadretti perfetti, ma di forme strane e irregolari (griglie non strutturate).

🚧 Il Problema: Il "Muro" tra i Vecchi e i Nuovi

Fino a poco tempo fa, c'era un grosso problema: il metodo "filo elastico" (FTM) funzionava bene solo su griglie semplici (come i quadretti di un foglio). Ma nel mondo reale (e nei computer moderni), le forme sono complesse e le griglie sono irregolari. Inoltre, il codice per farlo funzionare era così complicato e segreto che pochi lo usavano.

✨ La Soluzione: cfdmfFTFoam

L'autore, Ehsan Amani, ha creato un nuovo "motore" per il software OpenFOAM (che è come il "motore" gratuito più famoso per simulare i fluidi) chiamato cfdmfFTFoam.

Ecco come funziona, usando un'analogia:

Immagina che il computer stia simulando una festa in una stanza piena di persone (i fluidi).

• La Griglia (Il Pavimento): La stanza è divisa in zone irregolari (come un mosaico fatto di sassi di forme diverse).
• Il Filo Elastico (La Fronte): C'è un gruppo di ballerini che tengono per mano un lungo nastro elastico colorato che separa due gruppi di persone.
• Il Nuovo Software: È come se avessimo insegnato a questi ballerini come muoversi perfettamente su quel pavimento irregolare, senza mai perdere il contatto con il nastro, anche se il pavimento cambia forma o si rompe.

🛠️ Cosa fa esattamente questo software?

Il paper descrive come questo software gestisce il "filo elastico" in 6 passi magici ad ogni istante di tempo:

1. Correzione del Volume (Il "Riempimento"): A volte, muovendosi, il filo elastico perde un po' di "aria" o ne guadagna troppo. Il software controlla: "Ehi, la bolla è diventata più piccola? Rimettila a posto!" e spinge delicatamente i punti del filo per recuperare il volume originale.
2. Rimodellamento (Il "Ritaglio"): Se il filo si allunga troppo o si accartoccia, il software lo "taglia e ricuce". Aggiunge nuovi nodi dove serve (per avere più dettagli) o toglie quelli dove non servono (per risparmiare energia), mantenendo sempre la forma liscia.
3. Calcolo della Tensione Superficiale (La "Pelle"): L'acqua ha una "pelle" che cerca di contrarsi. Il software calcola esattamente quanto questa pelle tira su ogni singolo triangolo del filo elastico.
4. Comunicazione (Il "Messaggero"): Il filo elastico (che si muove liberamente) deve dire al pavimento (la griglia del computer) quanto tira. Usa un metodo intelligente (chiamato RKPM) per distribuire queste forze in modo che il pavimento le senta tutte correttamente, anche se i sassi del mosaico sono di forme diverse.
5. Costruzione dell'Indicatore (La "Mappa"): Il software disegna una mappa digitale che dice al computer: "Qui c'è acqua, lì c'è olio". Lo fa risolvendo un'equazione matematica per assicurarsi che il passaggio sia netto e preciso.
6. Avanzamento (Il "Passo"): Infine, sposta tutto il filo elastico in base a come sta fluendo l'acqua in quel momento.

🏆 Perché è importante?

Il paper mostra che questo nuovo software è un campionato mondiale di precisione:

• Non si rompe: Quando una goccia viene stirata fino a diventare sottile come un filo, i metodi vecchi (sfocati) spesso la fanno "rompere" o sparire. Questo metodo la mantiene intatta.
• Funziona ovunque: Può gestire stanze con pavimenti irregolari (griglie non strutturate), cosa che i vecchi metodi "filo elastico" non sapevano fare.
• È Open Source: È gratuito e disponibile per tutti, così altri scienziati possono usarlo per studiare cose come le bolle di sapone, la combustione nei motori, o il flusso del sangue nelle vene.

In sintesi

cfdmfFTFoam è come aver dato ai ricercatori un pennarello indelebile e intelligente che può disegnare confini perfetti su qualsiasi superficie, anche quella più irregolare, permettendo di vedere i fluidi muoversi con una precisione che prima era impossibile da ottenere su computer potenti. È un passo avanti enorme per chi studia come si comportano le gocce, le bolle e i fluidi nel mondo reale.

Sommario tecnico

Titolo: cfdmfFTFoam: Un solver di tracciamento dell'interfaccia (Front-Tracking) per flussi multifase su griglie non strutturate generali in OpenFOAM

1. Il Problema

La previsione affidabile dei flussi multifase complessi (come la deformazione di bolle, collisioni di gocce, spruzzi e ebollizione) richiede metodi di simulazione a risoluzione completa (FRS - Fully Resolved Simulation). Esistono due categorie principali: il tracciamento dell'interfaccia (Lagrangiano) e la cattura dell'interfaccia (Euleriano, come VOF o Level-Set).
Sebbene i metodi di cattura (es. VOF in OpenFOAM) siano ampiamente disponibili e supportati da software open-source, i metodi di tracciamento dell'interfaccia (FTM - Front-Tracking Method) puri offrono vantaggi significativi in termini di accuratezza nella calcolo della curvatura, della normale e nella gestione della fisica di coalescenza e rottura senza errori di dissipazione numerica.
Tuttavia, esiste una carenza critica di software open-source per FTM puro che funzioni su griglie Euleriane non strutturate generali. Le implementazioni esistenti (come PARIS, FronTier++, TrioCFD) sono limitate a griglie cartesiane strutturate. L'unico solver recente basato su OpenFOAM (lentFoam) è un approccio ibrido (Level-Set/FTM) che non rappresenta l'interfaccia con una mesh Lagrangiana pura. Questo lavoro mira a colmare tale lacuna.

2. Metodologia

L'autore ha sviluppato cfdmfFTFoam, un pacchetto solver integrato in OpenFOAM v9, che implementa un metodo FTM puro su griglie non strutturate. La metodologia si basa sull'integrazione del codice Fortran Ftc3D (sviluppato da Gretar Tryggvason e Sadegh Dabiri) all'interno dell'architettura C++ di OpenFOAM.

I componenti chiave dell'algoritmo includono:

• Dual Mesh: Risoluzione delle equazioni di Navier-Stokes su una griglia Euleriana volumetrica (non strutturata) e tracciamento dell'interfaccia tramite una mesh superficiale Lagrangiana triangolare mobile.
• Comunicazione Front-Field: Implementazione del Metodo delle Particelle a Nucleo Riproduttivo (RKPM) per trasferire le forze (tensione superficiale) dalla mesh Lagrangiana alla griglia Euleriana e per costruire la funzione indicatrice, garantendo la conservazione del momento zero e primo.
• Algoritmi Sottostanti:
  • Correzione del Volume (VC2): Un algoritmo ad hoc per correggere la posizione dei vertici della frontiera e preservare il volume della goccia/bolla, compensando la non conservazione intrinseca del FTM.
  • Remeshing: Algoritmi di raffinamento, accoppiamento (coarsening) e smoothing (rimozione delle ondulazioni) per mantenere la qualità della mesh triangolare (es. TSUR3D, VCS).
  • Calcolo della Tensione Superficiale: Metodo diretto basato sugli elementi (direct element-based).
  • Avvezione: Utilizzo delle librerie Lagrangiane di OpenFOAM per il tracciamento dei vertici della frontiera.

3. Contributi Chiave

• Primo Solver FTM Puro Open-Source su Griglie Non Strutturate: cfdmfFTFoam è il primo strumento open-access che permette simulazioni FTM pure su griglie generali (es. poliedriche, tetraedriche, prismi), superando il limite delle griglie cartesiane strutturate.
• Integrazione Completa in OpenFOAM: Il codice è stato riscritto in C++ e integrato nativamente nell'ecosistema OpenFOAM, supportando il calcolo parallelo (decomposizione del dominio) e le funzionalità di mesh dinamica.
• Libreria "frontTracking": È stata sviluppata una libreria modulare che incapsula la logica del tracciamento, permettendo l'estensibilità e l'uso di diversi modelli per la correzione del volume, lo smoothing e la comunicazione.
• Validazione Estesa: Il solver è stato testato contro benchmark standard, dimostrando la superiorità rispetto ai metodi VOF tradizionali in termini di accuratezza geometrica e riduzione delle correnti spurie (parasitic currents).

4. Risultati

Il documento presenta una serie di benchmark per validare il solver:

• Deformazione 3D (Benchmark di Leveque): Il solver mostra un errore di avvezione dell'interfaccia dell'1%, contro il 41% e il 36% ottenuti rispettivamente dai solver VOF (MULES e isoAdvector) di OpenFOAM. Il metodo FTM evita la rottura non fisica dell'interfaccia durante grandi stiramenti.
• Goccia Stagnante: La differenza di pressione calcolata ha un errore dello 0.046% rispetto alla soluzione analitica, con correnti spurie molto basse ($C_{asc} \approx 2.1 \times 10^{-3}$), superando di gran lunga le prestazioni del solver VOF standard.
• Goccia in Flusso di Poiseuille: Il solver dimostra accuratezza su diverse topologie di griglia (esagonali, prismi esagonali, poliedriche, tetraedriche) e mostra una buona scalabilità parallela, con efficienza vicina a quella del solver standard interFoam.
• Bolla di Taylor e Bolla Libera: Il solver riproduce accuratamente la forma terminale e la velocità di risalita di bolle in regime di Stokes e ad alto numero di Reynolds. È stato notato che la scelta del raggio del nucleo nel metodo RKPM ($\alpha$) è critica: valori più piccoli migliorano l'accuratezza ma richiedono regolarizzazione di Tikhonov per evitare singolarità matriciali.

5. Significato e Impatto

Il lavoro di Amani rappresenta un passo fondamentale per la comunità CFD open-source.

• Accessibilità: Rende il potente metodo FTM, precedentemente confinato a codici proprietari o griglie semplici, accessibile per ricerche su geometrie complesse e flussi multifase avanzati.
• Flessibilità: La capacità di utilizzare griglie non strutturate permette di simulare problemi con geometrie complesse (es. scambiatori di calore, reattori chimici) che non possono essere modellati efficientemente con griglie cartesiane.
• Base per Futuri Sviluppi: Il codice fornisce una base solida per l'implementazione futura di dinamiche di linea di contatto, cambi di topologia (rottura/coalescenza esplicita) e trasporto di calore/massa interfacciale.
• Riproducibilità: Essendo open-source (licenza GPLv3) e ospitato su GitHub e CPC Library, favorisce la trasparenza scientifica e la collaborazione nel campo della simulazione multifase.

In sintesi, cfdmfFTFoam è uno strumento maturo e validato che estende le capacità di OpenFOAM verso simulazioni multifase ad alta fedeltà, combinando l'accuratezza geometrica del Front-Tracking con la flessibilità delle griglie non strutturate.