A hybrid Volume of Fluid Phase-Field method for Direct Numerical Simulations of soluble surfactant-laden interfacial flows

Questo lavoro presenta un metodo ibrido Volume-of-Fluid Phase-Field con raffinamento adattivo della griglia per simulazioni numeriche dirette di flussi contenenti tensioattivi solubili, che cattura accuratamente l'accoppiamento tra trasporto nel bulk e all'interfaccia per dimostrare come gli stress di Marangoni alterino significativamente la dinamica di risalita delle bolle in geometrie tridimensionali.

L'essenziale

Il Quadro Generale: Un Nuovo Modo per Simulare l'Acqua "Saponosa"

Immagina di osservare una bolla che sale attraverso un bicchiere d'acqua. Se l'acqua è pura, la bolla sale dritta verso l'alto, oscillando un po' ma muovendosi velocemente. Ma se aggiungi sapone (tensioattivi) all'acqua, la bolla si comporta diversamente. Potrebbe rallentare, oscillare di più o persino cambiare percorso.

Questo accade perché le molecole di sapone amano attaccarsi alla superficie della bolla. Mentre la bolla si muove, queste molecole vengono spinte in giro, creando una "tensione" disomogenea sulla pelle della bolla. Questa tensione disomogenea agisce come mani invisibili che spingono la bolla in direzioni diverse, modificando il modo in cui si muove.

Il Problema:
Simulare questo su un computer è incredibilmente difficile. È come cercare di filmare una bolla di sapone con una telecamera che ha due impostazioni in conflitto:

1. La Telecamera Nitida: Deve vedere il bordo della bolla come una linea sottile come un rasoio (per calcolare pressione e forma).
2. La Telecamera Sgranata: Deve vedere le molecole di sapone che si distribuiscono in modo uniforme su quel bordo (per calcolare come si muove il sapone).

La maggior parte dei metodi informatici ti costringe a scegliere un'impostazione della telecamera, rendendo la simulazione o fisicamente inaccurata o computazionalmente impossibile da eseguire per forme 3D complesse.

La Soluzione:
Gli autori di questo documento hanno costruito un metodo ibrido. Pensalo come una simulazione a "schermo diviso" che utilizza il meglio di entrambi i mondi simultaneamente:

• Il Bordo Nitido (Volume-of-Fluid): Usano un metodo che mantiene il bordo della bolla nitido e conserva perfettamente la quantità di liquido (come un contorno ad alta definizione).
• Il Sapone Liscio (Phase-Field): Usano un secondo strato "sfocato" che funge da autostrada liscia su cui le molecole di sapone possono viaggiare. Questo permette al sapone di muoversi naturalmente tra l'acqua e la superficie della bolla senza rimanere intrappolato o perdersi.

Come Funziona: L'Analogia del "Controllore del Traffico"

Per far funzionare questo, gli autori hanno creato un sistema di traffico digitale per le molecole di sapone:

1. L'Autostrada (L'Interfaccia): La superficie della bolla è un'autostrada affollata. Le molecole di sapone sono le auto.
2. Le Rampa di Immissione e Uscita (Adsorbimento/Desorbimento):
   • Adsorbimento: Le molecole di sapone dall'acqua (la massa) vogliono saltare sull'autostrada (la superficie della bolla).
   • Desorbimento: Le molecole di sapone si stancano e saltano giù dall'autostrada tornando nell'acqua.
   • Il nuovo metodo calcola esattamente quante auto salgono o scendono in ogni singolo istante, assicurandosi che nessuna auto scompaia o appaia dal nulla.
3. Il Ingorgo (Stress di Marangoni): Quando troppe auto di sapone si accumulano in un punto sulla bolla, quel punto diventa "appiccicoso" (alta tensione). La pelle della bolla cerca di allontanarsi da quel punto appiccicoso, creando una forza che rallenta la bolla o la fa oscillare. La simulazione cattura perfettamente questa lotta di trazione.

Cosa Hanno Testato (I "Test di Guida")

Prima di lasciare che la loro nuova auto guidasse sull'autostrada, l'hanno fatta passare attraverso una scuola guida con tre test specifici:

1. Il Test di Allungamento (Sfera in Espansione): Hanno gonfiato una bolla ricoperta di sapone. Hanno verificato se il sapone si distribuiva uniformemente mentre la bolla diventava più grande. La simulazione corrispondeva perfettamente alla matematica.
2. Il Test di Rotazione (Bolla Rotante): Hanno fatto ruotare una bolla con il sapone sopra. Hanno verificato se il sapone si muoveva correttamente lungo il cerchio senza fuoriuscire. Ancora una volta, la simulazione era perfetta.
3. Il Test di Scambio (Parete Piana): Hanno osservato il sapone muoversi dall'acqua a una parete piana e viceversa. Hanno testato tre scenari:
   • Solo salire: Il sapone si attacca? Sì.
   • Solo scendere: Il sapone se ne va? Sì.
   • Entrambi: Trova un equilibrio? Sì.

L'Evento Principale: La Bolla che Sale

Infine, hanno lasciato che il loro nuovo metodo simulasse una bolla che sale in un serbatoio 3D di acqua.

• La Bolla Pulita: È salita relativamente velocemente e dritta.
• La Bolla di Sapone "Insolubile": Il sapone era incollato alla superficie e non poteva staccarsi. Ha creato un forte "ingorgo" nella parte posteriore della bolla, rallentandola significativamente.
• La Bolla di Sapone "Solubile" (Quella Reale): È qui che il nuovo metodo brilla. Il sapone poteva saltare su e giù dalla bolla.
  • Se il sapone saltava giù facilmente (alto "desorbimento"), la bolla si comportava quasi come una bolla pulita.
  • Se il sapone saltava su facilmente (alto "adsorbimento"), la bolla si comportava come la versione con il sapone bloccato.
  • Nel mezzo, la bolla ha mostrato una danza complessa: ha rallentato, cambiato percorso e lasciato una "scia" di sapone nell'acqua dietro di sé mentre saliva.

Perché Questo È Importante (Secondo il Documento)

Gli autori affermano che questo metodo è robusto, scalabile e accurato.

• Robusto: Non si blocca quando la bolla assume forme strane o si spezza.
• Scalabile: Può essere eseguito su supercomputer per gestire simulazioni 3D enormi e complesse in modo efficiente.
• Accurato: Predice correttamente la velocità con cui le bolle salgono e come oscillano, corrispondendo alla fisica del mondo reale.

In breve: Hanno costruito un nuovo motore digitale che può finalmente simulare come si comportano le bolle di sapone nello spazio 3D, gestendo la danza complessa tra la forma della bolla e le molecole di sapone che si muovono su e giù dalla sua pelle, tutto senza perdere accuratezza o bloccare il computer.

Sommario tecnico

Sintesi Tecnica: Un Metodo Ibrido Volume-of-Fluid e Phase-Field per Flussi Interfacciali con Surfattanti Solubili

Enunciazione del Problema
La simulazione numerica accurata di flussi interfacciali che coinvolgono surfattanti solubili presenta una sfida formidabile a causa della natura multiscale e multiphysics del problema. A differenza dei surfattanti insolubili, le specie solubili sono governate da un sistema accoppiato che coinvolge convezione-diffusione nel bulk, trasporto interfacciale e processi di scambio cinetico (adsorbimento/desorbimento). Queste interazioni generano campi di tensione superficiale che variano spazialmente e temporalmente, portando a stress di Marangoni che alterano significativamente la topologia del flusso, la stabilità e il moto. I metodi numerici esistenti presentano limitazioni distinte:

• Metodi a Integrale/Elemento di Bordo: Limitati ai regimi di flusso di Stokes e con difficoltà nel gestire cambiamenti di topologia (rottura/coalescenza).
• Front-Tracking: Alta accuratezza ma implementazione complessa in 3D, che richiede un trattamento topologico esplicito e si scontra con problemi di bilanciamento del carico nel calcolo parallelo.
• Metodi Level-Set: Gestiscono naturalmente i cambiamenti di topologia ma spesso soffrono di problemi di conservazione della massa che richiedono riinizializzazione e correzioni post-hoc.
• Metodi Phase-Field: Gestiscono naturalmente i cambiamenti di topologia ma spesso si basano su ipotesi di equilibrio termodinamico che possono essere violate da forte convezione, portando a non conservazione della massa.
• Metodi Volume-of-Fluid (VoF): Eccellenti per la conservazione della massa e interfacce nette, ma tradizionalmente faticano a definire operatori differenziali superficiali necessari per il trasporto di surfattanti su un'interfaccia a spessore zero.

Esiste un'esigenza specifica di un metodo basato su VoF geometrico che supporti il Raffinamento Adattivo della Mesh (AMR), gestisca surfattanti solubili con scambio dinamico bulk-interfaccia e operi all'interno di un framework open-source scalabile.

Metodologia
Gli autori presentano un framework numerico ibrido implementato nel codice open-source Basilisk, che combina un metodo geometrico Volume-of-Fluid (VoF) con un approccio Phase-Field diffuso.

1. Rappresentazione dell'Interfaccia:

   • VoF ($c$): Un campo scalare che rappresenta la frazione volumetrica di un fluido, utilizzato per il tracciamento conservativo dell'interfaccia e la soluzione della quantità di moto. Garantisce la conservazione esatta della massa e supporta la ricostruzione geometrica.
   • Phase-Field ($\phi$): Una variabile ausiliaria che evolve tramite una formulazione Accurate Conservative Diffuse-Interface (ACDI) basata su Allen-Cahn. Funge da vettore liscio per il trasporto del surfattante, fornendo vettori normali regolarizzati, curvatura e gradienti superficiali senza la necessità di calcoli espliciti di derivate superficiali su un'interfaccia netta.
   • Accoppiamento: Il Phase-Field viene periodicamente riinizializzato dal campo VoF utilizzando una funzione a distanza firmata per mantenere la coerenza.

2. Equazioni Governative:

   • Moto del Fluido: Equazioni di Navier-Stokes per fluidi incomprimibili con densità e viscosità variabili spazialmente. Le forze di tensione superficiale sono incorporate come forze di volume: una forza capillare isotropa ($\sigma \kappa \mathbf{n} \delta_s$) e una forza tangenziale di Marangoni ($\nabla_s \sigma \delta_m$).
   • Trasporto del Surfattante: Vengono risolte due campi scalari: concentrazione nel bulk ($F$) e concentrazione interfacciale ($f$, una rappresentazione volumetrica di $\Gamma$).
     • L'equazione di trasporto interfacciale è regolarizzata per evitare derivate superficiali, convertendo la concentrazione superficiale $\Gamma$ in un campo volumetrico $f = \Gamma \delta_\phi$ localizzato all'interno dell'interfaccia diffusa.
     • Cinetica: Adsorbimento e desorbimento sono modellati tramite termini sorgente/pozzo regolarizzati ($J$) localizzati all'interfaccia, che collegano le concentrazioni nel bulk e all'interfaccia. Il modello supporta equazioni di stato generali (ad es. isoterma di Langmuir) per correlare la tensione superficiale alla concentrazione.

3. Discretizzazione Numerica:

   • Spaziale: Discretizzazione a volumi finiti su griglie quadtree/octree dinamicamente raffinate (AMR). La risoluzione è concentrata vicino all'interfaccia per risolvere i gradienti ripidi nella concentrazione e nel flusso.
   • Temporale: Viene utilizzata la suddivisione degli operatori. L'avvezione è trattata esplicitamente (schema geometrico non diviso Bell–Collela–Glaz), mentre diffusione, anti-diffusione e termini sorgente sono trattati implicitamente utilizzando un solver Poisson multigrid modificato.
   • Stabilità: Viene imposta una condizione Courant–Friedrichs–Lewy (CFL) basata su una velocità efficace che include sia contributi convettivi che anti-diffusivi per garantire la stabilità in presenza di forti gradienti interfacciali.

Contributi Chiave

• Framework Ibrido: Sviluppo di un metodo che mantiene i vantaggi di conservazione della massa e interfaccia netta del VoF, utilizzando una variabile Phase-Field esclusivamente come vettore liscio per il trasporto del surfattante, evitando le complessità geometriche delle derivate superficiali su interfacce lineari a tratti.
• Capacità per Surfattanti Solubili: Estensione del precedente lavoro VoF-Phase-Field (limitato a surfattanti insolubili) al caso completamente solubile, incorporando termini accoppiati di avvezione-diffusione nel bulk e scambio cinetico.
• Scalabilità e Adattività: Il framework è completamente adattivo (AMR) e supporta configurazioni 2D, assialsimmetriche e 3D con alta scalabilità parallela, consentendo la risoluzione di gradienti multiscale sia nei campi del bulk che in quelli interfacciali.
• Implementazione Open-Source: Il metodo è implementato in Basilisk, fornendo uno strumento riproducibile e accessibile alla comunità.

Risultati e Validazione
Il documento valida il metodo attraverso una serie di test di riferimento e casi applicativi:

1. Validazione senza Flusso:

   • Sfera in Espansione (3D): Verifica del recupero delle soluzioni analitiche per la diluizione di surfattanti insolubili durante l'espansione isotropa.
   • Interfaccia Rotante (2D): Conferma della convergenza del secondo ordine per l'avvezione-diffusione di surfattanti insolubili su un'interfaccia in movimento.
   • Adsorbimento/Desorbimento su Interfaccia Piana: Validazione del solver accoppiato bulk-interfaccia contro soluzioni analitiche per adsorbimento a flusso costante, desorbimento puro e cinetica completa di Langmuir. Il metodo cattura correttamente la transizione all'equilibrio e la formazione di strati di impoverimento/arricchimento nel bulk.

2. Bolla Ascensionale Assialsimmetrica:

   • Limite Insolubile: Dimostrazione della riduzione della velocità terminale di risalita dovuta agli stress di Marangoni, con risultati che convergono al limite di interfaccia pulita al diminuire della concentrazione di surfattante.
   • Dinamica Solubile: Variazione sistematica del numero di Biot (desorbimento) e del numero di Damköhler (adsorbimento).
     • Biot Alto (Desorbimento Rapido): Il sistema si avvicina al limite di interfaccia pulita poiché i surfattanti vengono rimossi rapidamente dall'interfaccia.
     • Damköhler Alto (Adsorbimento Rapido): Il sistema si avvicina al limite di superficie incomprimibile (insolubile) poiché l'interfaccia si satura rapidamente.
   • Struttura del Flusso: Le visualizzazioni hanno mostrato la formazione di "cappucci di stagnazione" di surfattante nella parte posteriore della bolla, che sopprimono la ricircolazione interna e modificano la struttura della scia.

3. Bolla Ascensionale Tridimensionale:

   • Le simulazioni di una bolla che risale in 3D con surfattanti solubili hanno catturato fenomeni complessi tra cui deformazione della forma, deviazioni della traiettoria laterale (moto elicoidale) e formazione di pennacchi di scia dovuti al desorbimento.
   • Il metodo ha recuperato con successo sia i limiti asintotici puliti che quelli insolubili in 3D, dimostrando robustezza nella gestione di pattern di flusso instazionari e geometrie deformabili.

Significato e Affermazioni
Il documento afferma di fornire un framework numerico robusto e versatile che colma il divario tra i modelli classici di surfattanti insolubili e formulazioni realistiche solubili. La sua importanza risiede in:

• Trattamento Unificato: Gestisce l'accoppiamento non lineare, le interfacce evolutive e le configurazioni completamente tridimensionali all'interno di un'unica piattaforma open-source.
• Fedeltà Fisica: Cattura l'interazione complessa tra idrodinamica, trasporto bulk/interfacciale e stress di Marangoni senza fare affidamento su modelli di chiusura della turbolenza o approssimazioni empiriche.
• Applicabilità Pratica: Il metodo è adatto per flussi a superficie libera che coinvolgono riconnessione dell'interfaccia (ad es. rottura di getti, onde che si infrangono) in presenza di surfattanti solubili rilevanti a livello industriale e biologico.

Gli autori notano modestamente che l'attuale framework è limitato a numeri di Péclet moderati a causa della diffusione numerica intrinseca nella discretizzazione. Identificano la risoluzione accurata dei gradienti interfacciali netti nel regime ad alto Péclet (dominato dalla convezione) come una sfida per il lavoro futuro, con l'obiettivo di sviluppare formulazioni che minimizzino la diffusione numerica.