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

Technische Samenvatting: Deformeerbare lichamen in 3D viskeuze stroming via de vorticiteit-stroomvector-formulering

Probleemstelling
Het simuleren van de interactie tussen driedimensionale (3D) incompressibele viskeuze stromingen en deformeerbare of elastische obstakels (zoals vesikels, druppels en cellen) vormt aanzienlijke uitdagingen. Bestaande methoden, waaronder de Lattice Boltzmann-methode (LBM), grensvolg-technieken en deeltjesgebaseerde benaderingen zoals Smoothed Dissipative Particle Dynamics (SDPD), kampen vaak met een hoge computationele complexiteit, ingewikkelde numerieke implementaties of moeilijkheden bij het hanteren van complexe interface-mechanica. Hoewel 2D-simulaties met vorticiteit-stroomfunctie-formuleringen gekoppeld aan faseveldmodellen robuust zijn gebleken, ontbrak er een vergelijkbaar, lichtgewicht en klassiek gefundeerd kader voor 3D deformeerbare interfaces die worden beheerst door Canham–Helfrich (membraanbuiging) of Cahn–Hilliard (oppervlaktespanning) energieën. Specifiek is er een tekort in de literatuur met betrekking tot vorticiteitsgebaseerde formuleringen voor 3D-interfaces die worden beheerst door deze energieën.

Methodologie
De auteurs introduceren een nieuwe numerieke formulering gebaseerd op de vorticiteit-stroomvector-benadering, specifiek afgestemd op stromingen met een laag Reynoldsgetal. De methode generaliseert de 2D stroomfunctie naar een 3D vectorpotentiaal ($\vec{\psi}$), waarbij het snelheidveld wordt gedefinieerd als $\vec{v} = \nabla \times \vec{\psi}$.

De kern van de methodologie bestaat uit het oplossen van een gekoppeld systeem van vergelijkingen voor drie primaire observabelen:

1. Interface Chemisch Potentieel ($\mu$): Afgeleid als de functionele afgeleide van een gekozen vrije energie ($F$) met betrekking tot de faseveld-ordeparameter ($\phi$). De auteurs demonstreren het kader met behulp van zowel de Canham–Helfrich-energie (voor membraanbuiging) als de Cahn–Hilliard-energie (voor oppervlaktespanning).
2. Fluïdum Vorticiteit ($\vec{\omega}$): Gedefinieerd als $\vec{\omega} = \nabla \times \vec{v}$.
3. Stroomvector ($\vec{\psi}$): De vectorpotentiaal die voldoet aan de gauge-conditie $\nabla \cdot \vec{\psi} = 0$.

De beheersende vergelijkingen reduceren tot een reeks gekoppelde Poisson-problemen:

• De evolutie van het faseveld $\phi$ wordt beheerst door een advectie-diffusievergelijking die de chemische potentiaal en het snelheidveld omvat.
• De vorticiteitsvergelijking wordt afgeleid van de curl van de Stokes-vergelijking (of Navier-Stokes voor inertiële uitbreidingen), wat resulteert in een Poisson-vergelijking waarbij de bronterm afhankelijk is van het kruisproduct van de fasegradiënt en de gradiënt van de chemische potentiaal: $\nabla^2 \vec{\omega} = \frac{1}{\eta} (\nabla \phi) \times (\nabla \mu)$.
• De stroomvector wordt verkregen door $\nabla^2 \vec{\psi} = -\vec{\omega}$ op te lossen.

Numerieke Implementatie
Het systeem wordt opgelost met behulp van een eindverschillenmethode op een vast, niet-geconformeerd kubisch rooster. De interface wordt behandeld als een diffusiegebied via de faseveld-benadering, waardoor de noodzaak voor complexe roosterdeformatie of immersed boundary-technieken wordt geëlimineerd. De tijdsontwikkeling maakt gebruik van een forward Euler-integratie voor het faseveld en een iteratieve solver voor de Poisson-vergelijkingen, gebruikmakend van een gestaffeld (sequentieel gesegregeerd) schema. De code is geïmplementeerd in Python en is beschikbaar als een open-source repository.

Belangrijkste Bijdragen

• Vereenvoudiging van 3D-solvers: De formulering transformeert de fluïdumsolver in een reeks gekoppelde Poisson-problemen, waardoor directe berekening van druk wordt vermeden en strikte incompressibiliteit door constructie wordt gegarandeerd.
• Verenigd Kader voor Interface-energieën: De methode integreert naadloos verschillende vrije-energierelaties (Canham–Helfrich voor buiging en Cahn–Hilliard voor oppervlaktespanning) binnen dezelfde numerieke structuur, wat eenvoudige aanpassing van membraanmechanica mogelijk maakt.
• Diagnostische Observabelen: De auteurs stellen nieuwe observabelen voor op basis van de gekwadrateerde fluctuaties van de vorticiteits- en stroomvectorvelden ten opzichte van de ongeperturbreerde stroming. Deze metrieken blijken zeer gevoelig voor vormdeformaties en morfologische transities.

Resultaten
De methode werd gevalideerd door simulaties van enkelvoudige vesikels en druppels in Newtoniaanse Poiseuille- en Couette-stromingen:

• Poiseuille-stroming: De simulaties herstelden succesvol canonieke axiale vormen voor rode bloedcellen (RBC's), inclus[ief] "parachute" en "slipper" morfologieën, consistent met de literatuur. De diagnostische observabelen volgden de relaxatiedynamica effectief, waarbij duidelijke pieken werden waargenomen die overeenkomen met vormveranderingen (bijv. het verdwijnen van een deuk/dimple) en een plateau bereikten zodra het object een stationaire toestand bereikte.
• Couette-stroming: Het kader legde laterale migratie vast, een niet-inertieel liftfenomeen waarbij vesikels loodrecht op de stromingsrichting migreren. De simulaties toonden aan dat de evenwichtspositie en migratiedynamiek systematisch afhangen van het viscositeitscontrast ($\zeta$) tussen de interne en externe fluïda. Hogere viscositeitscontrasten leidden tot meer uitgesproken drift richting de bewegende wand.
• Stressanalyse: De methode berekende succesvol de in-plane afschuifspanningverdelingen, waarbij werd aangetoond dat hogere afschuifsnelheden leiden tot grotere stressheterogeniteit en morfologische instabiliteiten kunnen veroorzaken.

Betekenis en Claims
Het artikel beweert dat deze formulering een lichtgewicht, flexibel en computationeel efficiënt alternatief biedt voor bestaande mesoscopische of hoog-complexe solvers voor 3D low-Reynolds hydrodynamica. Door te vertrouren op klassieke vloeistofmechanica principes in plaats van kinetische theorie (zoals bij LBM), vermijdt de methode de noodzaak voor consistentiebewijzen met de continuümmechanica.

De auteurs benadrukken dat hoewel het huidige werk zich richt op single-body dynamica in Newtoniaanse fluïda, het kader ontworpen is voor uitbreidbaarheid. Het kan van nature lichaamskrachten, inertie-effecten, visco-elastische media en multi-body interacties accommoderen. Het vermogen om 3D in-plane afschuifdynamiek te vangen — een kenmerk dat vaak verloren gaat in 2D-benaderingen — positioneert deze methode als een waardevol instrument voor het bestuderen van cellulaire biomechanica en zachte materie fysica in beperkte stromingen. Het open-source karakter van de code is bedoeld om benchmarking en verdere ontwikkeling in dit opkomende gebied te stimuleren.