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

Dit artikel introduceert een bevochtigingsrandvoorwaardelijke conditie voor gebogen vaste oppervlakken in de kleurgradiënt lattice Boltzmann-methode door ordeparameters op ghost-nodes bij te werken, een schema dat op GPU-hardware is gevalideerd om effectief grote dichtheids- en viscositeitscontrasten te verwerken terwijl het spurious currents minimaliseert en zowel statische als dynamische contactlijngedragingen nauwkeurig reproduceert.

De kern

Stel je voor dat je probeert te simuleren hoe een druppel water zich gedraagt wanneer deze een gebogen oppervlak raakt, zoals een regendruppel die op een blad landt of een bubbel die langs een gebogen glas glijdt. Om dit op een computer te doen, gebruiken wetenschappers een methode genaamd de Lattice Boltzmann Method. Denk aan deze methode als een gigantisch, onzichtbaar rooster van kleine tegels die het computerscherm bedekken. Elke tegel bevat een beetje "vloeistofinformatie", en de computer werkt deze tegels stap voor stap bij om te zien hoe de vloeistof beweegt.

Het lastige deel is de randvoorwaarde—specifiek hoe de vloeistof zich gedraagt wanneer deze een vaste wand raakt. In de echte wereld stopt water niet abrupt bij een wand; het vormt een specifieke hoek (een contacthoek) afhankelijk van of het oppervlak nat is (zoals schoon glas) of droog (zoals een gewaxte auto).

Het Probleel: De "Geest" in de Machine

In de computersimulatie is de vaste wand geen gladde lijn; het is grillig omdat het gemaakt is van vierkante roostertegels. Om de wiskunde te laten werken, moet de computer weten wat de vloeistof doet binnen de vaste wand, ook al is daar geen vloeistof aanwezig. Deze denkbeeldige plekken binnen de wand worden "ghost nodes" (geestknopen) genoemd.

Eerdere methoden om deze ghost nodes aan te sturen, hadden bepaalde gebreken:

• Ze creëerden soms valse "geeststromingen" (spurious velocities), waarbij de vloeistof leek te bewegen zonder dat er een kracht aan te pas kwam.
• Ze hadden moeite met gebogen oppervlakken, omdat ze vaak leken te zijn ontworpen voor alleen platte wanden.
• Ze vereisten soms speciale, ingewikkelde wiskunde om een neutrale hoek te verwerken (waar het water noch uitspreidt, noch een druppel vormt).

De Oplossing: Een Nieuwe Regel voor de Geesten

De auteurs van dit artikel introduceerden een nieuwe, simpelere regel voor deze ghost nodes.

De Analogie: Stel je voor dat de vloeistof een "gemoedstoestand" heeft (vertegenwoordigd door een kleur, van 0 voor gas tot 1 voor vloeistof). In de echte wereld verandert deze gemoedstoestand geleidelijk van gas naar vloeistof terwijl je het oppervlak kruist.

• Oude Methode: Het was alsover proberen de gemoedstoestand van een persoon die achter een muur staat te raden door een willekeurige gok te schreeuwen.
• Nieuwe Methode: De auteurs realiseerden zich dat als je de "gemoedstoestand" van de persoon die net buiten de muur staat (in de vloeistof) kent, je die vloeiende curve van de gemoedstoestand wiskundig door de muur heen naar de ghost node kunt doortrekken. Ze vragen simpelweg: "Als de vloeistof hier een hoek van 45 graden wil vormen, wat moet de gemoedstoestand van de ghost node dan zijn om dat te laten gebeuren?"

Deze nieuwe regel is als een naadloze brug. Het breidt de natuurlijke vorm van de vloeistofdruppel uit tot vlak aan, en zelfs een klein stukje in, de vaste wand, waardoor de hoek die de druppel met de wand maakt precies is wat de wetenschapper heeft gevraagd.

Wat Ze Testten

Om te bewijzen dat hun nieuwe regel werkt, draaiden ze verschillende simulaties op een zeer krachtige computerchip (een NVIDIA A100 GPU):

1. De Statische Druppel: Ze plaatsten een waterdruppel op een platte plaat en een gebogen cilinder. Ze controleerden of de druppel op exact de hoek terechtkwam die ze hadden gevraagd.
   • Resultaat: Hun nieuwe regel was nauwkeuriger dan de vorige beste methode, vooral wanneer de hoek zeer scherp was (zoals een beading druppel) of zeer plat (zoals een uitspreidende druppel).
2. Het Zwevende Deeltje: Ze simuleerden een cilinder die drijft op de grens tussen olie en water.
   • Resultaat: Hun methode berekende de positie van de waterlijn nauwkeuriger dan voorheen.
3. De Vallende Druppel: Ze simuleerden een vallende druppel die een cilinder raakt, waarbij ze keken hoe de druppel spatte en zich verspreidde.
   • Resultaat: De druppel gedroeg zich realistisch, en de nieuwe regel veroorzaakte geen vreemde, valse bewegingen in de vloeistof.

Belangrijkste Conclusies

• Nauwkeurigheid: De nieuwe methode gaat veel beter om met gebogen oppervlakken, waardoor de hoek van de vloeistof correct blijft, of de wand nu plat of rond is.
• Stabiliteit: Het creëert zeer weinig "valse ruis" (spurious currents) in de simulatie, waardoor de vloeistof er natuurlijker uitziet.
• Eenvoud: Het vermijdt de noodzaak voor speciale, ingewikkelde wiskunde wanneer de contacthoek precies 90 graden (neutraal) is, wat een hoofdpijndossier was voor eerdere methoden.
• Snelheid: Door gebruik te maken van moderne computerchips (GPU's) en een specifieke programmeerstijl, maakten ze de simulaties zeer snel. Ze ontdekten dat het gebruik van een iets minder precieze getalindeling (single precision) de computer twee keer zo snel liet draaien zonder de resultaten voor de meeste tests te verslechteren.

Kortom, de auteurs hebben een betere "spelregel" gebouwd voor hoe computersimulaties de rand waar vloeistof een vaste wand raakt afhandelen, waardoor digitale druppels meer en meer lijken op echte druppels, zelfs op gebogen oppervlakken.

Technische samenvatting

Technische Samenvatting: Kleurgradiënt Lattice Boltzmann-modellering van Bevochtingsrandvoorwaarden op Gebogen Vaste Grensvlakken

Probleemstelling
Nauwkeurige modellering van bevochtingsverschijnselen is cruciaal voor het simuleren van technologisch relevante processen zoals stroming in poreuze media, microfluïdica en het drogen van colloïdale films. In Lattice Boltzmann-methoden (LBM) die gebruikmaken van diffuse interfaces, vormt de implementatie van bevochtingsrandvoorwaarden op gebogen vaste grensvlakken een aanzienlijke uitdaging. Bestaande benaderingen kampen vaak met beperkingen, waaronder de introductie van kunstmatige massa-bronnen of -putten, het genereren van spurieuze (parasitaire) snelheden nabij interfaces, beperkingen tot planaire geometrieën, en inconsistenties bij het gelijktijdig afdwingen van no-slip en bevochtigingscondities. Bov�ür maken veel methoden speciale behandeling vereist voor neutrale contacthoeken ($\theta = 90^\circ$) of vertrouwen ze op lineaire extrapolaties die fouten nabij vaste grenzen introduceren.

Methodologie
De auteurs introduceren een bevochtingsrandvoorwaarde die specifiek is ontworpen voor gebogen vaste grensvlakken binnen een tweedimensionaal kleurgradiënt LBM-raamwerk. De kern van de methodologie betreft het updaten van de ordeparameter (kleurveld, $\phi$) op "ghost nodes" die zich binnen de vaste fase bevinden.

• Ghost Node Update Regel: In plaats van een constante waarde voor te schrijven of lineaire extrapolatie te gebruiken, breidt dit schema het evenwichtige kleurprofiel uit in de vaste fase. De waarde van het kleurveld op een ghost node ($x_G$) wordt berekend op basis van het kleurveld op een overeenkomstig punt in het vloeistofdomein ($x_F$), gelegen op een afstand $\Delta x$ langs de vaste normaal.
• Wiskundige Formulering: Uitgaande van de aanname dat de kleurverdeling nabij evenwicht is, wordt de afgeleide van de ordeparameter in de richting van de vaste normaal afgeleid uit het evenwichtsprofiel en de voorgeschreven contacthoek $\theta$. Integratie van deze relatie levert een niet-lineaire update regel op:
  $$\phi_G = \frac{\phi_F}{\phi_F + (1 - \phi_F) \exp(\varepsilon)}$$
  waarbij $\varepsilon = -4\Delta x \cos \theta / W$, en $W$ de interfacebreedte is.
• Belangrijkste Kenmerken: Deze formulering zorgt ervoor dat het kleurveld van de ghost node begrensd blijft tussen 0 en 1, wat kunstmatige massaoverdracht voorkomt. Cruciaal is dat de regel geen speciale behandeling vereist wanneer $\theta = 90^\circ$ (waarbij $\varepsilon = 0$), in tegenstelling tot eerdere schema's. De methode is geïmplementeerd met een JAX-gebaseerd raamwerk dat draait op NVIDIA A100 GPU's, waarbij gebruik wordt gemaakt van kernel fusion en automatische differentiatie voor hoge prestaties.

Belangrijkste Bijdragen en Resultaten
Het voorgestelde schema werd gevalideerd tegen analytische oplossingen en het benchmark-schema van Fakhari et al. [20] over verschillende testgevallen:

1. Statische druppel op een vlakke plaat: Het schema produceerde consistent lagere fouten in gemeten contacthoeken vergeleken met Fakhari et al. over een bereik van $30^\circ$ tot $150^\circ$. Het genereerde ook kleinere spurieuze stromingen ($5,74 \times 10^{-7}$ tot $7 \times 10^{-7}$ rooster-eenheden) vergeleken met het benchmark-schema.
2. Druppel op een vaste cilinder: Voor een druppel die op een gebogen cilinder rust, toonde de voorgestelde methode een verbeterde nauwkeurigheid in het voorspellen van de apexhoogte van de druppel ($h_{max}$), met name bij extreme contacthoeken ($30^\circ$ en $150^\circ$), waarbij de fouten werden verminderd van $\sim 0,76\%$ naar $\sim 0,35\%$ vergeleken met het benchmark-schema.
3. Vaste deeltje op het vloeistof-gas interface: In simulaties van een periodieke reeks cilinders die gevangen is in een interface, verminderde het voorgestelde schema de fouten in verticale interfaceverplaatsing tot $\sim 9\%$ bij extreme hoeken, vergeleken met $\sim 14\%$ voor het benchmark-schema.
4. Dynamische druppelimpact: Het schema simuleerde succesvol de impact van een druppel op een cilinder bij hoge Reynolds- en Bond-getallen, waarbij complexe dynamische gedragingen werden vastgelegd, zoals het spreiden van druppels op hydrofiele oppervlakken en het splitsen op hydrofobe oppervlakken. De resultaten kwamen goed overeen met analytische schalingwetten voor filmdikte.
5. Computationele Prestaties: De GPU-implementatie toonde aan dat single-precision (FP32) berekeningen minstens twee keer zo snel waren als double-precision (FP64) berekeningen, terwijl de nauwkeurigheid behouden bleef voor statische en matige dynamische gevallen. Echter, FP32 werd instabiel voor de hoog-Reynolds/Weber-getal druppelimpacttest, waardoor FP64 noodzakelijk werd voor deze scenario's met hoge energie.

Betekenis en Claims
Het artikel claimt dat de voorgestelde bevochtigingsrandvoorwaarde een robuust en nauwkeurig alternatief biedt voor het simuleren van onmengbare twee-fasenstromingen op gebogen vaste grensvlakken. Door het evenwichtsprofiel via een specifieke ghost node update regel in de vaste fase uit te breiden, kan de methode:

• De capaciteit van het model om grote dichtheids- en viscositeitscontrasten aan te kunnen, behouden.
• Relatief kleine spurieuze stromingen produceren (orders van grootte lager dan populaire pseudopotentiaal-modellen).
• De numerieke singulariteiten vermijden die geassocieerd worden met neutrale contacthoeken in andere methoden.
• Een superieure nauwkeurigheid demonstreren ten opzichte van bestaande schema's (specifiek Fakhari et al.) voor zowel statische als dynamische contactlijnen op gebogen geometrieën.

De auteurs merken op dat hoewel de huidige implementatie beperkt is tot single-GPU uitvoering en twee dimensies, de onderliggende fysica direct toepasbaar is op phase-field en free-energy modellen en uitgebreid kan worden naar drie dimensies. Het werk vormt een fundament voor toekomstige simulaties waarbij vaste deeltjes op een vloeistof-gas interface betrokken zijn, met name wanneer dit gekoppeld wordt aan vooruitgang in berekeningen van capillaire krachthysteresis-schema's.