Deformable bodies in a 3-dimensional viscous flow: Vorticity-Stream vector formulation

Sintesi Tecnica: Corpi Deformabili in Flusso Viscoso 3D tramite la Formulazione Vorticità-Vettore Stream

Definizione del Problema
Simulare l'interazione tra flussi viscosi incompressibili tridimensionali (3D) e ostacoli deformabili o elastici (come vescicole, goccioline e cellule) presenta sfide significative. I metodi esistenti, inclusi i Metodi Lattice Boltzmann (LBM), le tecniche di integrale di bordo e gli approcci basati su particelle come la Smoothed Dissipative Particle Dynamics (SDPD), soffrono spesso di un'elevata complessità computazionale, implementazioni numeriche intricate o difficoltà nel gestire meccaniche interfacciali complesse. Sebbene le simulazioni 2D che utilizzano formulazioni vorticità-funzione di stream accoppiate con modelli di campo di fase si siano dimostrate robuste, è mancato un quadro comparabile, leggero e classicamente fondato per interfacce deformabili 3D a bassi numeri di Reynolds. Nello specifico, esiste una lacuna nella letteratura riguardante le formulazioni basate sulla vorticità per interfacce 3D governate da energie di tipo Canham–Helfrich (flessione di membrana) o Cahn–Hildland (tensione superficiale).

Metodologia
Gli autori introducono una nuova formulazione numerica basata sull'approccio vorticità-vettore stream, specificamente adattata per flussi a basso numero di Reynolds. Il metodo generalizza la funzione di stream 2D in un potenziale vettoriale 3D ($\vec{\psi}$), dove il campo di velocità è definito come $\vec{v} = \nabla \times \vec{\psi}$.

Il nucleo della metodologia consiste nel risolvere un sistema accoppiato di equazioni per tre osservabili primarie:

1. Potenziale Chimico dell'Interfaccia ($\mu$): Derivato dalla derivata funzionale di un'energia libera scelta ($F$) rispetto al parametro d'ordine del campo di fase ($\phi$). Gli autori dimostrano il framework utilizzando sia l'energia di Canham–Helfrich (per la flessione della membrana) che l'energia di Cahn–Hilliard (per la tensione superficiale).
2. Vorticità del Fluido ($\vec{\omega}$): Definita come $\vec{\omega} = \nabla \times \vec{v}$.
3. Vettore Stream ($\vec{\psi}$): Il potenziale vettoriale che soddisfa la condizione di gauge $\nabla \cdot \vec{\psi} = 0$.

Le equazioni governanti si riducono a un insieme di problemi di Poisson accoppiati:

• L'evoluzione del campo di fase $\phi$ è governata da un'equazione di advezione-diffusione che coinvolge il potenziale chimico e il campo di velocità.
• L'equazione della vorticità è derivata dal rotore dell'equazione di Stokes (o di Navier-Stokes per estensioni inerziali), risultando in un'equazione di Poisson dove il termine sorgente dipende dal prodotto vettoriale del gradiente di fase e del gradiente del potenziale chimico: $\nabla^2 \vec{\omega} = \frac{1}{\eta} (\nabla \phi) \times (\nabla \mu)$.
• Il vettore stream viene ottenuto risolvendo $\nabla^2 \vec{\psi} = -\vec{\omega}$.

Implementazione Numerica
Il sistema viene risolto utilizzando un metodo a differenze finite su una griglia cubica fissa e non adattata (unfitted). L'interfaccia è trattata come una regione diffusiva tramite l'approccio field-phase, eliminando la necessità di complessa deformazione della mesh o tecniche di metodo del bordo immerso (immersed boundary). L'evoluzione temporale impiega un'integrazione Euler forward per il campo di fase e un risolutore iterativo per le equazioni di Poisson, utilizzando uno schema scalato (sequenziale segregato). Il codice è implementato in Python ed è disponibile come repository open-source.

Contributi Chiave

• Semplificazione dei Solutori 3D: La formulazione trasforma il solutore del fluido in un insieme di problemi di Poisson accoppiati, evitando il calcolo diretto della pressione e garantendo per costruzione una rigorosa incompressibilità.
• Framework Unificato per Energie Interfacciali: Il metodo integra perfettamente diversi modelli di energia libera (Canham–Helfrich per la flessione e Cahn–Hilliard per la tensione superficiale) all'interno della stessa struttura numerica, permettendo una facile personalizzazione della meccanica delle membrane.
• Osservabili Diagnostiche: Gli autori propongono nuove osservabili basate sulle fluttuazioni al quadrato dei campi di vorticità e del vettore stream rispetto al flusso non perturbato. Queste metriche si sono dimostrate altamente sensibili alle deformazioni di forma e alle transizioni morfologiche.

Risultati
Il metodo è stato validato attraverso simulazioni di singole vescicole e goccioline in flussi di Poiseuille e Couette newtoniani:

• Flusso di Poiseuille: Le simulazioni hanno recuperato con successo le forme assialsimmetriche canoniche per i globuli rossi (RBC), inclusi i morfotipi "paracadute" e "slipper", coerentemente con la letteratura. Le osservabili diagnostiche hanno tracciato efficacemente la dinamica di rilassamento, mostrando picchi distinti corrispondenti ai cambiamenti di forma (ad esempio, la scomparsa di un incavo) e stabilizzandosi una volta che l'oggetto ha raggiunto uno stato stazionario.
• Flusso di Couette: Il framework ha catturato la migrazione laterale, un fenomeno di lift non inerziale in cui le vescicole migrano perpendicolarmente alla direzione del flusso. Le simulazioni hanno dimostrato che la posizione di equilibrio e la dinamica di migrazione dipendono sistematicamente dal contrasto di viscosità ($\zeta$) tra i fluidi interno ed esterno. Contrasti di viscosità più elevati hanno portato a una deriva più pronunciata verso la parete mobile.
• Analisi dello Stress: Il metodo ha calcolato con successo le distribuzioni di sforzo di taglio nel piano, rivelando che tassi di taglio più elevati inducono una maggiore eterogeneità dello stress e possono portare a instabilità morfologiche.

Significatività e Rivendicazioni
L'articolo sostiene che questa formulazione offre un'alternativa leggera, flessibile ed efficiente dal punto di vista computazionale ai solutori mesoscopici o ad alta complessità per l'idrodinamica 3D a basso numero di Reynolds. Basandosi sui principi classici della meccanica dei fluidi piuttosto che sulla teoria cinetica (come nell'LBM), il metodo evita la necessità di prove di coerenza con la meccanica del continuo.

Gli autori sottolineano che, sebbene l'attuale lavoro si concentri sulla dinamica di corpo singolo in fluidi newtoniani, il framework è progettato per essere estensibile. Può naturalmente accogliere forze di corpo, effetti inerziali, mezzi viscoelastici e interazioni multi-corpo. La capacità di catturare la dinamica di taglio nel piano 3D — una caratteristica spesso persa nelle approssimazioni 2D — posiziona questo metodo come uno strumento prezioso per lo studio della biomeccanica cellulare e della fisica della materia soffice. La natura open-source del codice è intesa per favorire il benchmarking e lo sviluppo ulteriore in questo ambito emergente.