Color-gradient lattice Boltzmann modeling of wetting boundary condition on curved solid boundaries

Questo articolo introduce una condizione al contorno di bagnabilità per superfici solide curve nel metodo lattice Boltzmann a gradiente di colore aggiornando i parametri d'ordine sui nodi fantasma, uno schema validato su hardware GPU per gestire efficacemente grandi contrasti di densità e viscosità minimizzando le correnti spurie e riproducendo accuratamente sia i comportamenti della linea di contatto statica che di quella dinamica.

L'essenziale

Immagina di cercare di simulare il comportamento di una goccia d'acqua quando colpisce una superficie curva, come una goccia di pioggia che atterra su una foglia o una bolla che scivola lungo un vetro curvo. Per farlo su un computer, gli scienziati usano un metodo chiamato Metodo Lattice Boltzmann. Pensa a questo metodo come a una gigantesca griglia invisibile di minuscole piastrelle che copre lo schermo del computer. Ogni piastrella contiene un po' di informazione sul "fluido", e il computer aggiorna queste piastrelle passo dopo passo per vedere come si muove il fluido.

La parte complicata è la condizione al contorno — specificamente, come il fluido si comporta quando tocca una parete solida. Nel mondo reale, l'acqua non si ferma semplicemente davanti a una parete; forma un angolo specifico (chiamato angolo di contatto) a seconda che la superficie sia bagnata (come un vetro pulito) o asciutta (come un'auto cerata).

Il Problema: Il "Fantasma" nella Macchina

Nella simulazione al computer, la parete solida non è una linea liscia; è frastagliata perché è fatta di piastrelle di una griglia quadrata. Per far funzionare la matematica, il computer deve sapere cosa sta facendo il fluido dentro la parete solida, anche se lì non c'è fluido. Questi punti immaginari all'interno della parete sono chiamati "nodi fantasma" (ghost nodes).

I metodi precedenti per dire a questi nodi fantasma cosa fare presentavano alcuni difetti:

• A volte creavano "correnti fantasma" (velocità spurie) dove il fluido sembrava muoversi da solo senza alcuna forza.
• Facevano fatica con le superfici curve, agendo spesso come se fossero stati progettati solo per pareti piatte.
• A volte richiedevano una matematica speciale e complicata solo per gestire un angolo neutro (dove l'acqua non si diffonde né si raggruppa).

La Soluzione: Una Nuova Regola per i Fantasmi

Gli autori di questo articolo hanno introdotto una nuova regola più semplice per questi nodi fantasma.

L'Analogia: Immagina che il fluido abbia un "umore" (rappresentato da un colore, da 0 per il gas a 1 per il liquido). Nel mondo reale, questo umore cambia gradualmente da gas a liquido mentre attraversi la superficie.

• Vecchio Metodo: Era come cercare di indovinare l'umore di una persona che sta dietro un muro gridando un tentativo casuale.
• Nuovo Metodo: Gli autori hanno capito che se conosci l' "umore" della persona che si trova appena fuori dal muro (nel fluido), puoi matematicamente estendere quella curva dell'umore fluida attraverso il muro fino al nodo fantasma. Chiedono semplicemente: "Se il fluido vuole formare un angolo di 45 gradi qui, quale deve essere l'umore del nodo fantasma per far sì che ciò accada?"

Questa nuova regola è come un ponte senza soluzione di continuità. Estende la forma naturale della goccia di fluido fin quasi dentro la parete solida, assicurando che l'angolo che la goccia forma con la parete sia esattamente quello richiesto dallo scienziato.

Cosa hanno testato

Per dimostrare che la loro nuova regola funziona, hanno eseguito diverse simulazioni su un chip per computer molto potente (un NVIDIA A100 GPU):

1. La Goccia Statica: Hanno posto una goccia d'acqua su una piastra piatta e su un cilindro curvo. Hanno controllato se la goccia si assestava esattamente sull'angolo richiesto.
   • Risultato: La loro nuova regola era più accurata rispetto al precedente miglior metodo, specialmente quando l'angolo era molto acuto (come una goccia che si raggruppa) o molto piatto (come una goccia che si diffonde).
2. La Particella Galleggiante: Hanno simulato un cilindro che galleggia al confine tra olio e acqua.
   • Risultato: Il loro metodo calcolava la posizione della linea dell'acqua in modo più accurato rispetto al passato.
3. La Goccia in Caduta: Hanno simulato una goccia che cade e colpisce un cilindro, osservando come schizza e si diffonde.
   • Risultato: La goccia si è comportata in modo realistico e la nuova regola non ha causato strani movimenti falsi nel fluido.

Punti Chiave

• Accuratezza: Il nuovo metodo gestisce le superfici curve molto meglio dei vecchi metodi, mantenendo corretto l'angolo del fluido sia che la parete sia piatta sia che sia rotonda.
• Stabilità: Crea pochissimo "rumore falso" (correnti spurie), il che significa che il fluido appare più naturale.
• Semplicità: Evita la necessità di una matematica speciale e complicata quando l'angolo di contatto è esattamente di 90 gradi (neutro), il che era un problema per i metodi precedenti.
• Velocità: Utilizzando chip per computer moderni (GPU) e uno stile di programmazione specifico, hanno reso le simulazioni molto veloci. Hanno scoperto che l'uso di un formato numerico leggermente meno preciso (precisione singola) rendeva il computer due volte più veloce senza rovinare i risultati per la maggior parte dei test.

In breve, gli autori hanno costruito un migliore "libro di regole" per come le simulazioni al computer gestiscono il punto di incontro tra un liquido e una parete solida, rendendo le gocce digitali simili e comportarsi come quelle reali, anche su superfici curve.

Sommario tecnico

Sintesi Tecnica: Modellazione LBM a Gradiente di Colore per Condizioni al Contorno di Bagnabilità su Confini Solidi Curvi

Definizione del Problema
La modellazione accurata dei fenomeni di bagnabilità è fondamentale per simulare processi tecnologicamente rilevanti come il flusso in mezzi porosi, la microfluidica e l'essiccazione di film colloidali. Nei framework del metodo Lattice Boltzmann (LBM) che utilizzano interfacce diffuse, l'implementazione di condizioni al contorno di bagnabilità su confini solidi curvi presenta sfide significative. Gli approcci esistenti soffrono spesso di limitazioni che includono l'introduzione di sorgenti o pozzi di massa artificiali, la generazione di velocità spurie (parassite) vicino alle interfacce, restrizioni a geometrie planari e incongruenze nel soddisfare simultaneamente le condizioni di no-slip e di bagnabilità. Inoltre, molti metodi richiedono un trattamento speciale per angoli di contatto neutri ($\theta = 90^\circ$) o si affidano a estrapolazioni lineari che introducono errori vicino ai confini solidi.

Metodologia
Gli autori introducono una condizione al contorno di bagnabilità progettata specificamente per confini solidi curvi all'interno di un framework LBM a gradiente di colore bidimensionale. Il nucleo della metodologia consiste nell'aggiornare il parametro d'ordine (campo di colore, $\phi$) sui "nodi fantasma" (ghost nodes) situati all'interno della fase solida.

• Regola di Aggiornamento dei Nodi Fantasma: Invece di prescrivere un valore costante o utilizzare l'estrapolazione lineare, lo schema estende il profilo di equilibrio del colore nella fase solida. Il valore del campo di colore in un nodo fantasma ($x_G$) viene calcolato sulla base del campo di colore in un punto corrispondente nel dominio fluido ($x_F$) situato a una distanza $\Delta x$ lungo la normale al solido.
• Formulazione Matematica: Assumendo che la distribuzione del colore sia vicina all'equilibrio, la derivata del parametro d'ordine nella direzione della normale al solido viene derivata dal profilo di equilibrio e dall'angolo di contatto prescritto $\theta$. L'integrazione di questa relazione produce una regola di aggiornamento non lineare:
  $$\phi_G = \frac{\phi_F}{\phi_F + (1 - \phi_F) \exp(\varepsilon)}$$
  dove $\varepsilon = -4\Delta x \cos \theta / W$, e $W$ è la larghezza dell'interfaccia.
• Caratteristiche Chiave: Questa formulazione garantisce che il campo di colore del nodo fantasma rimanga compreso tra 0 e 1, prevenendo il trasferimento di massa artificiale. Fondamentalmente, la regola non richiede trattamenti speciali quando $\theta = 90^\circ$ (dove $\varepsilon = 0$), a differenza di altri schemi. Il metodo è implementato utilizzando un framework basato su JAX in esecuzione su GPU NVIDIA A100, sfruttando la fusione dei kernel (kernel fusion) e l'differenziazione automatica per alte prestazioni.

Contributi Chiavi e Risultati
Lo schema proposto è stato validato rispetto a soluzioni analitiche e allo schema di riferimento di Fakhari et al. [20] attraverso diversi casi di test:

1. Goccia Statica su una Piastra Piatta: Lo schema ha prodotto costantemente errori minori negli angoli di contatto misurati rispetto a Fakhari et al. in un intervallo da $30^\circ$ a $150^\circ$. Ha inoltre generato correnti spurie minori ($5.74 \times 10^{-7}$ a $7 \times 10^{-7}$ unità di reticolo) rispetto allo schema di riferimento.
2. Goccia su un Cilindro Solido: Per una goccia posta su un cilindro curvo, il metodo proposto ha dimostrato una migliore accuratezza nella previsione dell'altezza dell'apice della goccia ($h_{max}$), in particolare per angoli di contatto estremi ($30^\circ$ e $150^\circ$), dove gli errori sono stati ridotti da $\sim 0.76\%$ a $\sim 0.35\%$ rispetto allo schema di riferimento.
3. Particella Fissa all'Interfaccia Liquido-Gas: Nelle simulazioni di un array periodico di cilindri intrappolati a un'interfaccia, lo schema proposto ha ridotto gli errori nello spostamento verticale dell'interfaccia a circa il $9\%$ per angoli estremi, rispetto al $\sim 14\%$ dello schema di riferimento.
4. Impatto di Goccia Dinamica: Lo schema ha simulato con successo l'impatto di una goccia su un cilindro ad alti numeri di Reynolds e Bond, catturando comportamenti dinamici complessi come la diffusione della goccia su superfici idrofile e la scissione su superfici idrofobiche. I risultati sono in linea con le leggi di scala analitiche per lo spessore del film.
5. Prestazioni Computazionali: L'implementazione su GPU ha mostrato che i calcoli a singola precisione (FP32) erano almeno due volte più veloci di quelli a doppia precisione (FP64), mantenendo l'accuratezza per i casi statici e dinamici moderati. Tuttavia, l'FP32 è diventato instabile per il test di impatto della goccia ad alto numero di Reynolds/Weber, rendendo necessario l'uso di FP64 per tali scenari ad alta energia.

Significatività e Rivendicazioni
Il documento sostiene che il proposto schema di bagnabilità offre un'alternativa robusta e accurata per la simulazione di flussi bifase immiscibili su confini solidi curvi. Estendendo il profilo di equilibrio nella fase solida tramite una specifica regola di aggiornamento dei nodi fantasma, il metodo:

• Mantiene la capacità del modello di gestire grandi contrasti di densità e viscosità.
• Produce correnti spurie relativamente piccole (ordini di grandezza inferiori rispetto ai popolari modelli pseudopotenziali).
• Evita le singolarità numeriche associate agli angoli di contatto neutri presenti in altri metodi.
• Dimostra una precisione superiore rispetto agli schemi esistenti (specificamente Fakhari et al.) sia per linee di contatto statiche che dinamiche su geometrie curve.

Gli autori osservano che, sebbene l'attuale implementazione sia limitata all'esecuzione su singola GPU e in due dimensioni, la fisica sottostante è direttamente applicabile ai modelli di campo di fase (phase-field) e di energia libera e può essere estesa a tre dimensioni. Il lavoro fornisce una base per future simulazioni coinvolgenti particelle solide all'interfaccia liquido-gas, in particolare quando accoppiate con avanzamenti nei calcoli delle forze capillari della linea di contatto e degli schemi di isteresi.