A 3D sharp and conservative VOF method for modeling the contact line dynamics with hysteresis on complex boundaries

Deze paper introduceert een scherp en behoudend 3D VOF-numerieke methode binnen een ingebedde grensvlakframework die, door een aangepast advectie- en reconstructieschema voor gemengde cellen en een innovatieve contacthoek-impositietechniek met hysteresis, contactlijndynamica op complexe oppervlakken nauwkeurig simuleert zonder CFL-tijdstapbeperkingen.

De kern

Stel je voor dat je een druppel water ziet die over een oppervlak rolt. Soms rolt hij soepel over een gladde tafel, maar wat gebeurt er als die tafel vol zit met oneffenheden, gaatjes, of als hij schuin staat? En wat als die druppel vastplakt aan de rand van een gaatje?

Dit is precies het probleem dat de auteurs van dit wetenschappelijke artikel proberen op te lossen. Ze hebben een nieuwe, slimme computermethode ontwikkeld om te simuleren hoe vloeistoffen (zoals water of inkt) bewegen over complexe, onregelmatige oppervlakken.

Hier is een uitleg in gewone taal, met een paar creatieve vergelijkingen:

1. Het Probleem: De "Vastzittende" Druppel

In de echte wereld zijn oppervlakken zelden perfect vlak. Denk aan een regenjas met plooien, een insectenpootje met haren, of een printer met een ingewikkeld patroon. Als je een computermodel maakt om te zien hoe een druppel zich gedraagt op zo'n oppervlak, krijg je vaak twee grote problemen:

• De "Gaten" in de rekensom: Computers werken vaak met een rooster van kleine kubusjes (net als een 3D-puzzel). Als een oppervlak schuin door een kubusje loopt, ontstaat er een klein, onregelmatig stukje vloeistof. De computer raakt hierdoor vaak de vloeistof kwijt of telt er teveel bij, alsof je water uit een lekke emmer haalt.
• De "Kleine Steentjes" in de weg: Soms zijn die onregelmatige stukjes zo klein dat de computer extreem langzaam moet rekenen om ze nauwkeurig te volgen. Het is alsof je een auto moet besturen, maar elke keer als je over een klein steentje rijdt, moet je volledig stoppen om te kijken of je er niet overheen rijdt.

2. De Oplossing: Een Slimme "Puzzel" Methode

De auteurs hebben een nieuwe manier bedacht om deze problemen op te lossen, bestaande uit twee hoofdstukken:

Deel A: De Onuitputtelijke Emmer (Behoud van Volume)

Stel je voor dat je vloeistof verplaatst door een doolhof van muren. In de oude methoden lekte er soms water uit de emmer als je langs een schuine muur liep.

• De nieuwe truc: Ze hebben een systeem bedacht dat elke druppel water die de emmer verlaat, precies meet. Als er door een klein gaatje in de muur (een "cut cell") iets te veel of te weinig water lijkt te verdwijnen, wordt dit direct gecompenseerd.
• De "Verdelings-strategie": Als er een heel klein stukje water vastzit in een hoekje dat de computer bijna niet kan zien, gooien ze dat niet weg. In plaats daarvan verdelen ze dat beetje water netjes over de buren. Zo blijft de totale hoeveelheid water exact hetzelfde, en hoeft de computer niet te wachten op die kleine steentjes. De auto rijdt nu soepel door, zelfs over de oneffenheden.

Deel B: De Slimme Hoekmeter (Contacthoek)

Wanneer een druppel een oppervlak raakt, vormt hij een hoek. Op een vlakke tafel is dit makkelijk te meten. Maar op een bolle of gekartelde oppervlakte? Dat is lastig.

• De oude manier: De computer probeerde de hoek te raden door een rechte lijn te trekken. Dit werkte goed op vlakke dingen, maar op ronde of geknikte oppervlakken gaf dit een scheef beeld, alsof je een bolle spiegel gebruikt om je gezicht te meten.
• De nieuwe manier: Ze gebruiken een "parabool" (een boogvormige lijn) in plaats van een rechte lijn. Ze kijken niet alleen naar het punt waar de druppel raakt, maar ook naar de vorm van de druppel eromheen. Ze passen de vorm van de druppel aan alsof ze een stuk klei vormen die perfect tegen de onregelmatige muur aan past. Hierdoor weet de computer precies hoe de druppel moet liggen, zelfs als het oppervlak eruitziet als een berg met pieken en dalen.

3. Het "Kleef-effect" (Hysteresis)

Soms plakt een druppel even vast voordat hij weer gaat bewegen (net als een druppel regen op een raam die pas wegrolt als het raam schuin genoeg staat). De nieuwe methode kan dit "plakken" en "loslaten" ook simuleren. Het houdt rekening met het feit dat de druppel soms even vastzit aan de oneffenheden voordat hij verder rolt.

Waarom is dit belangrijk?

Deze methode is als een superkracht voor ingenieurs en wetenschappers:

• Inktjetprinters: Om druppels inkt perfect op papier te plaatsen, zelfs op ruw papier.
• Medische technologie: Om te begrijpen hoe bloed of medicijnen door kleine, kromme buisjes stromen.
• Waterzuivering: Om te zien hoe water door complexe filters stroomt.

Kortom: De auteurs hebben een computerprogramma gebouwd dat vloeistoffen op complexe oppervlakken kan simuleren zonder dat er water "verdwijnt" in de rekensom en zonder dat de computer vastloopt op kleine details. Het is alsof ze een perfecte, onuitputtelijke emmer hebben bedacht die over elke denkbare berg en dal kan rollen, terwijl de vorm van de druppel altijd perfect wordt bewaard.

Technische samenvatting

Hier is een gedetailleerde technische samenvatting van het artikel "A 3D sharp and conservative VOF method for modeling the contact line dynamics with hysteresis on complex boundaries", gepubliceerd in het Journal of Computational Physics (2026).

1. Probleemstelling

De dynamiek van het contactlijn (de lijn waar vloeistof, gas en vast materiaal samenkomen) speelt een cruciale rol in talloze engineering- en wetenschappelijke processen, zoals inkjet-printen, microfluidica en water-elektrolyse. Bestaande numerieke methoden voor het simuleren van deze fenomenen stuiten op twee fundamentele uitdagingen wanneer ze worden toegepast op complexe, niet-uitgelijnde geometrieën (zoals gebogen of onregelmatige oppervlakken):

1. Massa- en volumebehoud in gemengde cellen: In Cartesische roosters snijden complexe vaste lichamen de roostercellen, waardoor "gemengde cellen" ontstaan die zowel vloeistof, gas als vast materiaal bevatten. Bestaande VOF (Volume-of-Fluid) advectieschema's leiden vaak tot lokale massa-ongenauwkeurigheden en niet-fysische accumulatie van vloeistof in vaste gebieden, vooral wanneer de vloeistofvolume-fractie zeer klein is (kleine "cut cells"). Dit leidt ook tot extreme tijdstapbeperkingen (CFL-constraints).
2. Accurate oplegging van het contacthoek-voorwaarde: Het nauwkeurig bepalen van de kromming en het opleggen van de contacthoek op onregelmatige oppervlakken is moeilijk. Bestaande methoden (zoals lineaire fitting of ghost-cell extrapolatie) falen vaak bij extreme contacthoeken of gekromde oppervlakken, wat leidt tot asymmetrische druppelverspreiding en onnauwkeurige oppervlaktespanningskrachten.

2. Methodologie

De auteurs hebben een nieuw 3D-numeriek kader ontwikkeld binnen de open-source solver Basilisk, dat een gekoppelde geometrische VOF-methode en een Embedded Boundary Method (EBM) combineert. De kern van de methode ligt in het scherp oplossen van zowel het vloeistof/gas-interface ($\Gamma_l$) als het vaste grensvlak ($\Gamma_s$) binnen een Cartesisch rooster.

Belangrijkste methodologische componenten:

• Geometrische VOF-advectie in gemengde cellen:

  • Een aangepast flood-fill algoritme wordt gebruikt om het vloeistof/gas-interface te reconstrueren binnen de onregelmatige polyhedrale sub-volumes van gemengde cellen.
  • Om lokale massa-ongenauwkeurigheden te corrigeren, wordt een gecorrigeerde advectiewijdte ($u_n^*$) geïntroduceerd. Deze compenseert voor het volume dat wordt geblokkeerd door het vaste materiaal, waardoor de flux over celvlakken exact behouden blijft.
  • Om de strenge tijdstapbeperkingen veroorzaakt door kleine cut cells te omzeilen, wordt een flux-herverdelingsstrategie (redistribution advection) ontwikkeld. Deze strategie verplaatst overtollig of tekortkomend volume naar aangrenzende cellen. Hierdoor wordt de CFL-beperking volledig verwijderd, terwijl zowel lokaal als globaal volumebehoud wordt gegarandeerd.

• 3D Contacthoek-voorwaarde met Paraboloïde Fitting:

  • In plaats van lineaire fitting, wordt een hoogtefunctie (Height Function - HF) methode gebruikt die is uitgebreid met een paraboloïde fitting strategie.
  • Pre-fitting: Eerst wordt een voorlopige paraboloïde gefit op basis van HF-punten in het vloeistofgebied om de lokale interface-geometrie te schatten.
  • Gecombineerde fitting: Vervolgens wordt een finale paraboloïde gefit die de contacthoek-voorwaarde expliciet integreert op het snijpunt van het contactlijn en het vaste oppervlak. Dit gebeurt door een niet-lineair systeem op te lossen dat de HF-data en de geometrische contacthoek-voorwaarde ($\cos \theta = \mathbf{n}_l \cdot \mathbf{n}_s$) combineert.
  • Verticale HF: Voor extreme contacthoeken of vlakken die bijna evenwijdig lopen aan de roosterassen, wordt een "verticale" hoogtefunctie-methode geïntroduceerd om de robuustheid te waarborgen.
  • Hysterese: Een iteratief bisection-algoritme wordt gebruikt om contacthoek-hysterese (verschil tussen vooruitgaande en terugtrekkende hoek) te modelleren, waarbij de contactlijn vastgepind blijft totdat de lokale hoek de drempelwaarden bereikt.

3. Belangrijkste Bijdragen

1. Eerste volledig geometrische en conservatieve 3D VOF-scheme: Voor het eerst wordt een VOF-methode gepresenteerd die strikt lokaal en globaal massabehoud garandeert op willekeurig complexe 3D-oppervlakken, zonder de tijdstapbeperkingen van traditionele cut-cell methoden.
2. Robuuste 3D Contacthoek-implementatie: De ontwikkeling van een nieuwe paraboloïde-fitting methode die contacthoeken nauwkeurig oplegt op onregelmatige, gebogen oppervlakken, zelfs bij extreme hoeken. Dit lost het probleem op van asymmetrie en onnauwkeurigheid bij bestaande lineaire methoden.
3. Integratie van Hysterese: Een numerieke strategie voor het modelleren van contactlijn-pinning en hysterese in 3D, gebaseerd op volumebehoudsprincipes.

4. Resultaten en Validatie

De methode is gevalideerd via een reeks uitdagende benchmarktests:

• Volumebehoud: In een test met een roterende vloeistofschaal rond een vaste bol (met complexe cut cells) toonde de gecorrigeerde methode lokale fouten aan op het niveau van machineprecisie ($L_\infty \approx 10^{-15}$), terwijl traditionele methoden grote lokale fouten ($L_\infty \approx 10^{-2}$) vertoonden.
• Druppelverspreiding op vlakke en bolle oppervlakken: De methode produceerde evenwichtscontactlijnen die perfect overeenkwamen met theoretische voorspellingen, ongeacht de oriëntatie van het oppervlak of de diepte van de embedding. In tegenstelling tot lineaire fitting, behield de paraboloïde-methode de cirkelvormige symmetrie van de contactlijn.
• Schuifgedreven beweging met hysterese: Simulaties van druppels die glijden onder schuifstroom toonden correcte pinning/depining dynamiek en eindvormen die overeenkwamen met bestaande literatuur.
• Complexe 3D Geometrieën: De methode slaagde erin druppelimpact op platen met gaten (ronde en scherpe randen), spontane klimbeweging op kegels en verspreiding op roosterachtige structuren nauwkeurig te simuleren. De methode bleef stabiel bij scherpe geometrische overgangen en hoge krommingen.

5. Betekenis

Deze studie vertegenwoordigt een significante doorbraak in de numerieke simulatie van multiphase flows. De voorgestelde methode overwint de beperkingen van bestaande scherpe-interface methoden door:

• Efficiëntie: Het elimineren van de tijdstapbeperkingen die vaak 3D-simulaties met complexe geometrieën onpraktisch maken.
• Nauwkeurigheid: Het bieden van een geometrisch consistente manier om contacthoeken op willekeurige 3D-oppervlakken op te leggen, wat essentieel is voor het modelleren van natting en afstoting.
• Toepasbaarheid: Het openen van de deur voor het simuleren van realistische, complexe wettingsproblemen in engineering (zoals coatingprocessen, microfluidica en materiaalwetenschap) met een hoge mate van betrouwbaarheid en conservatisme.

De methode is geïmplementeerd in de open-source solver Basilisk en biedt een krachtig hulpmiddel voor fundamenteel onderzoek en praktische toepassingen op het gebied van vloeistof-dynamica.