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

Dit artikel presenteert een hybride Volume-of-Fluid-Phase-Field-methode met adaptief meshverfijning voor directe numerieke simulaties van stromingen met oplosbare surfactanten, die de koppeling tussen bulk- en interfaciaal transport nauwkeurig vastlegt om aan te tonen hoe Marangoni-spanningen de stijgdynamica van bellen in driedimensionale geometrieën aanzienlijk veranderen.

De kern

Het Grote Plaatje: Een Nieuwe Manier om "Zeepachtig" Water te Simuleren

Stel je voor dat je een bel ziet opstijgen door een glas water. Als het water zuiver is, schiet de bel recht omhoog, met een beetje wiebelen maar snel bewegend. Maar als je zeep (oppervlakte-actieve stoffen) aan het water toevoegt, gedraagt de bel zich anders. Hij kan vertragen, meer wiebelen, of zelfs zijn pad veranderen.

Dit gebeurt omdat zeepmoleculen graag aan het oppervlak van de bel plakken. Terwijl de bel beweegt, worden deze moleculen rondgeduwd, waardoor er ongelijke "spanning" ontstaat op de huid van de bel. Deze ongelijke spanning werkt als onzichtbare handen die de bel in verschillende richtingen duwen, waardoor de manier waarop hij beweegt verandert.

Het Probleem:
Dit op een computer simuleren is ongelooflijk moeilijk. Het is alsof je probeert een zeepbel te filmen met een camera die twee tegenstrijdige instellingen heeft:

1. De Scherpe Camera: Moet de rand van de bel zien als een mesdunne lijn (om druk en vorm te berekenen).
2. De Vage Camera: Moet zien hoe de zeepmoleculen zich glad over die rand verspreiden (om te berekenen hoe de zeep beweegt).

De meeste computermethoden dwingen je om één camera-instelling te kiezen, waardoor de simulatie ofwel fysiek onnauwkeurig is, ofwel computertechnisch onmogelijk om uit te voeren voor complexe 3D-vormen.

De Oplossing:
De auteurs van dit artikel hebben een hybride methode ontwikkeld. Denk hierbij aan een "split-screen" simulatie die tegelijkertijd het beste van beide werelden gebruikt:

• De Scherpe Rand (Volume-of-Fluid): Ze gebruiken een methode die de rand van de bel scherp houdt en de hoeveelheid vloeistof perfect behoudt (zoals een hoog-resolutie omtrek).
• De Gladde Zeep (Phase-Field): Ze gebruiken een tweede, "wazige" laag die fungeert als een gladde snelweg voor de zeepmoleculen om op te reizen. Hierdoor kan de zeep natuurlijk bewegen tussen het water en het oppervlak van de bel zonder vast te komen zitten of verloren te gaan.

Hoe Het Werkt: De "Verkeersleider" Analogie

Om dit werkend te krijgen, hebben de auteurs een digitaal verkeerssysteem voor de zeepmoleculen gecreëerd:

1. De Snelweg (Het Interface): Het oppervlak van de bel is een drukke snelweg. De zeepmoleculen zijn auto's.
2. De Op- en Afritten (Adsorptie/Desorptie):
   • Adsorptie: Zeepmoleculen uit het water (de bulk) willen de snelweg op (het oppervlak van de bel) springen.
   • Desorptie: Zeepmoleculen worden moe en springen van de snelweg af terug het water in.
   • De nieuwe methode berekent precies hoeveel auto's op elk moment op- of afstappen, zodat er geen auto's verdwijnen of uit het niets verschijnen.
3. De File (Marangoni-spanning): Wanneer te veel zeepauto's zich op één plek op de bel ophopen, wordt die plek "plakkerig" (hoge spanning). De huid van de bel probeert zich van die plakkerige plek af te trekken, waardoor een kracht ontstaat die de bel vertraagt of laat wiebelen. De simulatie vangt dit trek-krachtspel perfect op.

Wat Ze Testten (De "Rijproeven")

Voordat ze hun nieuwe auto op de snelweg lieten rijden, hebben ze hem door een rijschool gestuurd met drie specifieke tests:

1. De Rektest (Uitdijende Bol): Ze bliezen een met zeep bedekte bel op. Ze controleerden of de zeep zich gelijkmatig verspreidde naarmate de bel groter werd. De simulatie kwam perfect overeen met de wiskunde.
2. De Draaitest (Roterende Bel): Ze draaiden een bel met zeep erop. Ze controleerden of de zeep correct rond de cirkel bewoog zonder te lekken. Ook hier was de simulatie spot-on.
3. De Uitwisseltest (Vlakke Muur): Ze keken hoe zeep zich verplaatste van het water naar een vlakke muur en terug. Ze testten drie scenario's:
   • Alleen opstappen: Plakt de zeep? Ja.
   • Alleen afstappen: Verlaat de zeep? Ja.
   • Beide: Vindt het een balans? Ja.

Het Hoofdevenement: De Opstijgende Bel

Tot slot lieten ze hun nieuwe methode een bel simuleren die opstijgt in een 3D-tank met water.

• De Schone Bel: Hij steeg relatief snel en recht omhoog.
• De "Onoplosbare" Zeepbel: De zeep zat vast aan het oppervlak en kon niet weg. Dit creëerde een sterke "file" aan de achterkant van de bel, waardoor deze aanzienlijk vertraagde.
• De "Oplosbare" Zeepbel (De Echte Kwestie): Hier blinkt de nieuwe methode uit. De zeep kon op en van de bel springen.
  • Als de zeep makkelijk afstapte (hoge "desorptie"), gedroeg de bel zich bijna als een schone bel.
  • Als de zeep makkelijk opstapte (hoge "adsorptie"), gedroeg de bel zich als de versie met vastzittende zeep.
  • In het midden toonde de bel een complexe dans: hij vertraagde, veranderde zijn pad en liet een "spoor" van zeep in het water achter zich naarmate hij opstijg.

Waarom Dit Belangrijk Is (Volgens Het Artikel)

De auteurs beweren dat deze methode robuust, schaalbaar en nauwkeurig is.

• Robuust: Het crasht niet wanneer de bel vreemde vormen krijgt of uit elkaar valt.
• Schaalbaar: Het kan draaien op supercomputers om enorme, complexe 3D-simulaties efficiënt te verwerken.
• Nauwkeurig: Het voorspelt correct hoe snel bellen opstijgen en hoe ze wiebelen, en komt overeen met de fysica uit de echte wereld.

Kortom: Ze hebben een nieuwe digitale motor gebouwd die eindelijk kan simuleren hoe zeepbellen zich in de 3D-ruimte gedragen, waarbij ze de lastige dans tussen de vorm van de bel en de zeepmoleculen die op en af hun huid bewegen, afhandelen, allemaal zonder nauwkeurigheid te verliezen of de computer te laten crashen.

Technische samenvatting

Technische Samenvatting: Een Hybride Volume-of-Fluid en Faseveld-methode voor Oplosbare Surfactant-bevattende Interfaciale Stromingen

Probleemstelling
De nauwkeurige numerieke simulatie van interfaciale stromingen waarbij oplosbare surfactanten betrokken zijn, vormt een formidable uitdaging vanwege de multischaal- en multiphysica-aard van het probleem. In tegenstelling tot onoplosbare surfactanten worden oplosbare soorten beheerst door een gekoppeld systeem dat bulkconvectie-diffusie, interfaciaal transport en kinetische uitwisselingsprocessen (adsorptie/desorptie) omvat. Deze interacties genereren ruimtelijk en tijdelijk variërende oppervlaktespanningsvelden, wat leidt tot Marangoni-spanningen die de stromingstopologie, stabiliteit en beweging aanzienlijk veranderen. Bestaande numerieke methoden staan voor duidelijke beperkingen:

• Grensintegraal/-elementmethoden: Beperkt tot Stokes-stromingsregimes en hebben moeite met topologieveranderingen (breuk/samenvoeging).
• Front-tracking: Hoge nauwkeurigheid maar complexe implementatie in 3D, vereist expliciete topologische behandeling en ondervindt lastbalanceringsproblemen in parallelle computing.
• Level-setmethoden: Handhaven topologieveranderingen op natuurlijke wijze maar lijden vaak aan massa- behoudsproblemen die herinitialisatie en correcties achteraf vereisen.
• Faseveldmethoden: Handhaven topologieveranderingen op natuurlijke wijze maar vertrouwen vaak op thermodynamische evenwichtsveronderstellingen die kunnen worden geschonden door sterke convectie, wat leidt tot niet-behoud van massa.
• Volume-of-Fluid (VoF)-methoden: Uitstekend voor massabehoud en scherpe interfaces maar hebben traditioneel moeite met het definiëren van oppervlakedifferentiaaloperatoren die vereist zijn voor surfactanttransport op een interface met nul dikte.

Er is een specifieke behoefte aan een geometrische VoF-gebaseerde methode die Adaptief Netwerkverfijning (AMR) ondersteunt, oplosbare surfactanten met dynamische bulk-interface-uitwisseling verwerkt, en opereert binnen een schaalbaar, open-source kader.

Methodologie
De auteurs presenteren een hybride numeriek kader geïmplementeerd in de open-source code Basilisk, dat een geometrische Volume-of-Fluid (VoF)-methode combineert met een diffuse Faseveld-benadering.

1. Interface Representatie:

   • VoF ($c$): Een scalair veld dat het volumedeel van één vloeistof vertegenwoordigt, gebruikt voor conservatief interface-tracking en momentumoplossing. Het zorgt voor exact massabehoud en ondersteunt geometrische reconstructie.
   • Faseveld ($\phi$): Een hulpvariabele die evolueert via een op Allen-Cahn gebaseerde Accurate Conservatieve Diffuse-Interface (ACDI)-formulering. Het dient als een gladde drager voor surfactanttransport, waarbij geregulariseerde normaalvectoren, kromming en oppervlaktegradiënten worden geleverd zonder de noodzaak van expliciete oppervlakteafgeleidenberekeningen op een scherpe interface.
   • Koppeling: Het Faseveld wordt periodiek herinitialiseerd vanuit het VoF-veld met behulp van een functie met tekenafstand om consistentie te handhaven.

2. Beheersende Vergelijkingen:

   • Vloeistofbeweging: Incompressibele Navier-Stokes-vergelijkingen met ruimtelijk variërende dichtheid en viscositeit. Oppervlaktespanningskrachten worden opgenomen als volumetrische lichaamskrachten: een isotrope capillaire kracht ($\sigma \kappa \mathbf{n} \delta_s$) en een tangentiële Marangoni-kracht ($\nabla_s \sigma \delta_m$).
   • Surfactanttransport: Twee scalair velden worden opgelost: bulkconcentratie ($F$) en interfaciale concentratie ($f$, een volumetrische representatie van $\Gamma$).
     • De interfaciale transportvergelijking wordt geregulariseerd om oppervlakteafgeleiden te vermijden, waarbij de oppervlakteconcentratie $\Gamma$ wordt omgezet in een volumetrisch veld $f = \Gamma \delta_\phi$ dat gelokaliseerd is binnen de diffuse interface.
     • Kinetiek: Adsorptie en desorptie worden gemodelleerd via geregulariseerde bron/zink-termen ($J$) die gelokaliseerd zijn aan de interface, waardoor bulk- en interfaciale concentraties met elkaar verbonden zijn. Het model ondersteunt algemene toestandsvergelijkingen (bijv. Langmuir-isotherm) om oppervlaktespanning te relateren aan concentratie.

3. Numerieke Discretisatie:

   • Ruimtelijk: Finite-volume-discretisatie op dynamisch verfijnde quadtree/octree-roosters (AMR). De resolutie is geconcentreerd nabij de interface om steile gradiënten in concentratie en stroming op te lossen.
   • Tijdelijk: Operator-splitting wordt gebruikt. Advektie wordt expliciet behandeld (geometrisch ongesplitst Bell–Collela–Glaz-schema), terwijl diffusie, anti-diffusie en bron-termen impliciet worden behandeld met een gewijzigde multigrid Poisson-oplosser.
   • Stabiliteit: Een Courant–Friedrichs–Lewy (CFL)-voorwaarde wordt afgedwongen op basis van een effectieve snelheid die zowel convectieve als anti-diffusieve bijdragen omvat om stabiliteit te garanderen in aanwezigheid van sterke interface-gradiënten.

Belangrijkste Bijdragen

• Hybride Kader: Ontwikkeling van een methode die de massabehoudende en scherpe-interface-voordelen van VoF behoudt, terwijl een Faseveld-variabele uitsluitend wordt gebruikt als een gladde drager voor surfactanttransport, waardoor de geometrische complexiteiten van oppervlakteafgeleiden op stuksgewijs lineaire interfaces worden vermeden.
• Vermogen voor Oplosbare Surfactanten: Uitbreiding van eerdere VoF-Faseveld-werk (beperkt tot onoplosbare surfactanten) naar het volledig oplosbare geval, met inbegrip van gekoppelde bulk-advektie-diffusie en kinetische uitwisselingstermen.
• Schaalbaarheid en Adaptiviteit: Het kader is volledig adaptief (AMR) en ondersteunt 2D, axissymmetrische en 3D-configuraties met hoge parallelle schaalbaarheid, waardoor de oplossing van multischaal-gradiënten in zowel bulk- als interfaciale velden mogelijk wordt.
• Open-source Implementatie: De methode is geïmplementeerd in Basilisk, en biedt een reproduceerbaar en toegankelijk instrument voor de gemeenschap.

Resultaten en Validatie
Het artikel valideert de methode door middel van een reeks benchmarktesten en toepassingsgevallen:

1. Validatie zonder Stroming:

   • Uitdijende Bol (3D): Geverifieerd het herstel van analytische oplossingen voor verdunning van onoplosbare surfactanten tijdens isotrope uitdijing.
   • Roterende Interface (2D): Bevestigde tweede-orde convergentie voor advektie-diffusie van onoplosbare surfactanten op een bewegende interface.
   • Adsorptie/Desorptie op Vlakke Interface: De gekoppelde bulk-interface-oplosser is gevalideerd tegen analytische oplossingen voor constante flux-adsorptie, pure desorptie en volledige Langmuir-kinetiek. De methode vangt correct de overgang naar evenwicht en de vorming van uitputtings/verrijkinglagen in de bulk.

2. Axissymmetrisch Opstijgend Bellen:

   • Onoplosbare Limiet: Gemonstreerd de reductie van de terminale stijgsnelheid door Marangoni-spanningen, waarbij de resultaten convergeren naar de schone-interface-limiet naarmate de surfactantconcentratie afneemt.
   • Oplosbare Dynamica: Systematisch gevarieerd met het Biot-getal (desorptie) en het Damköhler-getal (adsorptie).
     • Hoog Biot (Snelle Desorptie): Het systeem nadert de schone-interface-limiet terwijl surfactanten snel van de interface worden verwijderd.
     • Hoog Damköhler (Snelle Adsorptie): Het systeem nadert de oppervlakte-oncompressibele (onoplosbare) limiet terwijl de interface snel verzadigt.
   • Stromingsstructuur: Visualisaties toonden de vorming van "stagnatiekappen" van surfactant aan de achterzijde van de bel, die interne recirculatie onderdrukken en de wake-structuur modificeren.

3. Driedimensionaal Opstijgend Bellen:

   • Simulaties van een bel die opstijgt in 3D met oplosbare surfactanten vingen complexe fenomenen op, waaronder vormdeformatie, zijwaartse trajectafwijkingen (helicoïdale beweging) en de vorming van achterwaartse pluimen door desorptie.
   • De methode slaagde erin zowel schone als onoplosbare asymptotische limieten in 3D te herstellen, wat robuustheid demonstreert bij het hanteren van onstabiele stromingspatronen en vervormende geometrieën.

Betekenis en Claims
Het artikel claimt een robust en veelzijdig numeriek kader te bieden dat de kloof overbrugt tussen klassieke modellen voor onoplosbare surfactanten en realistische oplosbare formuleringen. De betekenis hiervan ligt in:

• Unificatie: Het behandelt niet-lineaire koppeling, evoluerende interfaces en volledig driedimensionale configuraties binnen één open-source platform.
• Fysische Fideliteit: Het vangt de complexe wisselwerking tussen hydrodynamica, bulk/interfaciaal transport en Marangoni-spanningen zonder te vertrouwen op turbulentie-sluitingsmodellen of empirische benaderingen.
• Praktische Toepasbaarheid: De methode is geschikt voor vrije-vloeistofstromen waarbij interface-herconnectie optreedt (bijv. straalbreuk, brekende golven) in aanwezigheid van industriële en biologisch relevante oplosbare surfactanten.

De auteurs merken bescheiden op dat het huidige kader beperkt is tot matige Péclet-getallen vanwege inherente numerieke diffusie in de discretisatie. Zij identificeren de nauwkeurige oplossing van scherpe interfaciale gradiënten in het hoge Péclet-regime (convectie-gedomineerd) als een uitdaging voor toekomstig werk, met als doel formuleringen te ontwikkelen die numerieke diffusie minimaliseren.