A generalized, quasi-universal subgrid model for macroscopic moving contact line flows

Stel je een druppel water voor die op een plat oppervlak landt. Terwijl deze uitspreidt, beweegt de rand waar de vloeistof, de lucht en het vaste oppervlak samenkomen niet simpelweg zoals een autoband over het asfalt glijdt. In plaats daarvan gedraagt deze grens, bekend als een contactlijn, zich op een manier die de standaardregels van de vloeistofmechanica tart die worden gebruikt om de meeste stromende vloeistoffen te beschrijven. Decennialang hebben wetenschappers geprobeerd deze beweging op computers te simuleren, omdat de fysica aan de uiterste rand van de druppel zich afspeelt op schalen die zo minuscuul zijn — tot op de grootte van individuele moleculen — dat het proberen te berekenen van elk detail meer rekenkracht zou vereisen dan er bestaat. Om dit te omzeilen, hebben onderzoekers gebruikgemaakt van kortere wegen, of "subgrid-modellen", om te gokken wat er op die microscopische rand gebeurt. Echter, deze kortere wegen hebben een groot gebrek: de resultaten die ze produceren veranderen vaak afhankelijk van hoe fijn de computer de ruimte in een rooster verdeelt, wat de simulaties onbetrouwbaar maakt voor voorspellingen in de echte wereld.

Een team van onderzoekers aan het Indian Institute of Science Education and Research in Bhopal heeft een nieuwe manier ontwikkeld om dit probleem aan te pakken. Ze hebben een model gecreëerd dat computers in staat stelt om bewegende contactlijnen nauwkeurig te simuleren zonder de onmogelijke-om-te-meten details van de microscopische wereld te hoeven kennen. Hun aanpak overbrugt de kloof tussen de kleine, onzichtbare fysica aan de rand van een druppel en het grotere, zichtbare gedrag van de vloeistof. Door dit te doen, hebben ze simulaties geproduceerd die consistent en nauwkeurig blijven, ongeacht de grootte van het computerraster, een cruciale stap naar het voorspellen van hoe vloeistoffen zich gedragen in alles van inkjetprinters tot olieherwinning, zonder dat het model voor elke nieuwe situatie aangepast hoeft te worden.

De kern van het probleem ligt in een conflict tussen hoe vloeistoffen bewegen en hoe ze vaste oppervlakken raken. In de standaard vloeistofmechanica stelt een regel genaamd de "no-slip condition" dat een vloeistof die een vast oppervlak raakt, volledig moet stoppen met bewegen bij die grens. Hoक dit werkt goed voor de meeste stromingen, het breekt bij de contactlijn van een bewegende druppel, wat een wiskundige singulariteit creëert waarbij de krachten oneindig worden. Om dit op te lossen, hebben wetenschappers historisch gezien aangenomen dat de vloeistof een klein beetje glijdt aan het oppervlak, maar de afstand van dit glijden is zo klein dat het onmogelijk direct te meten is. Gevolgelijk moesten eerdere computermodellen deze microscopische waarden raden of vertrouwen op passende parameters om overeen te komen met experimentele gegevens. Dit betekende dat als een onderzoeker een ander type vloeistof of een andere snelheid wilde simuleren, hij vaak opnieuw moest beginnen en het model opnieuw moest kalibreren, wat het vermogen van het model om nieuwe scenario's te voorspellen ondermijnde.

De onderzoekers in deze studie realiseerden zich dat de inconsistentie in eerdere simulaties voortkwam uit een misverstand over hoe de hoek van de druppel verandert naarmate men zich van de contactlijn verwijdert. Ze maakten gebruik van een algemeen aanvaarde theorie uit de jaren 1980 die beschrijft hoe het oppervlak van de vloeistof vervormt door de competitie tussen oppervlaktespanning en de interne wrijving van de vloeistof. Deze theorie laat zien dat de hoek die je ziet afhangt van hoe ver je van de rand verwijderd bent. De onderzoekers combineerden dit theoretische begrip met bestaande empirische formules — regels afgeleid van het observeren van echte druppels — om een nieuw, verenigd model te creëren. In plaats van de computer te dwingen de microscopische glijlengte of de exacte hoek aan de wand te raden, berekent hun model de juiste hoek te gebruiken op basis van de grootte van de cellen van het computerraster.

Om hun idee te testen, voerde het team een reeks computersimulaties uit van een druppel die zich verspreidt op een platte plaat. Ze vergeleken hun nieuwe methode met oudere technieken met verschillende roostergroottes, variërend van grof tot zeer fijn. Bij de oudere methoden veranderde de roostergrootte de resultaten aanzienlijk, waardoor de gesimuleerde druppel met verschillende snelheden uitspreidde of andere vormen aannam. Met het nieuwe model bleven de resultaten echter vrijwel identiek over alle roostergroottes heen. De gesimuleerde druppel gedroeg zich consistent, of de computer nu een grof of een fijn raster gebruikte. Dit toonde aan dat het model erin geslaagd was om de kunstmatige afhankelijkheid van het raster, een fenomeen dat bekend staat als roosterafhankelijkheid, te elimineren.

De studie onderzocht ook hoe dit nieuwe model presteert onder verschillende omstandigheden, zoals verschillende stroomsnelheden en verschillende soorten vloeistoffen. Ze testten scenario's waarin de vloeistof zich verspreidde (voortschrijdend) en waarin deze zich terugtrok (terugwijkend). In elk geval produceerde het nieuwe model stabiele, roosteronafhankelijke resultaten. Een belangrijke bevinding was dat het model geen "passende parameters" vereiste die afkomstig waren van specifieke experimenten om te werken. Zodra de onderzoekers de schaal hadden vastgesteld waarop de oorspronkelijke empirische formules waren afgeleid, kon het model het gedrag van de stroming voorspellen zonder verdere aanpassingen. Dit is een aanzienlijke afwijking van eerdere methoden, die vaak een unieke kalibratie vereisten voor elke nieuwe set omstandigheden.

De onderzoekers onderzochten ook waarom sommige oudere simulaties in het verleden goed leken te werken, zelfs zonder deze nieuwe correctie. Ze ontdekten dat in bepaalde situaties, met name bij gas-vloeistofsystemen, de fout die door de roostergrootte werd geïntroduceerd van nature klein genoeg was om niet opgemerkt te worden. Bovendien, als het computerraster toevallig overeenkwam met de specifieke schaal waarop de oorspronkelijke experimentele gegevens waren verzameld, zouden de resultaten toevallig overeenkomen met de werkelijkheid. Dit was echter een kwestie van geluk eerder dan een robuuste oplossing. Het nieuwe model verwijdert dit element van toeval en biedt een betrouwbaar kader dat werkt over een breed scala aan fysieke omstandigheden, inclusief verschillende dichtheden en viscositeiten.

Een van de meest veelbelovende aspecten van dit werk is de potentie voor voorspellende kracht. Omdat het model niet vertrouwt op het kennen van de microscopische glijlengte of de exacte contacthoek aan de wand — grootheden die berucht moeilijk te meten zijn — kan het worden toegepast op nieuwe problemen waar deze waarden onbekend zijn. De onderzoekers toonden aan dat door simpelweg de schaal te kennen waarop de oorspronkelijke empirische observaties zijn gedaan, het model universeel kan worden toegepast. Dit betekent dat ingenieurs en wetenschappers deze simulaties kunnen gebruiken om betere systemen te ontwerpen voor het coaten van oppervlakken, het beheren van vloeistofstroming in poreuze gesteenten of het begrijpen van biologische processen in de longen, met de zekerheid dat de resultaten geen artefacten zijn van het computerraster.

De studie concludeert dat hoewel de eerste benadering van hun model goed werkt voor de meeste hoeken, meer complexe scenario's met zeer grote contacthoeken een iets gedetailleerdere wiskundige correctie kunnen vereisen. Echter, de fundamentele aanpak van het koppelen van de microscopische theorie aan macroscopische observaties is succesvol gebleken. Door de contactlijn niet te behandelen als een vast punt met mysterieuze eigenschappen, maar als een dynamisch kenmerk dat reageert op de schaal van observatie, hebben de onderzoekers een instrument geleverd dat zowel nauwkeurig als breed toepasbaar is. Dit werk vertegenwoordigt een verschuiving van louter het aanpassen van gegevens naar het begrijpen van de onderliggende fysica, wat een pad biedt naar simulaties die werkelijk kunnen voorspellen hoe vloeistoffen zich in de echte wereld zullen gedragen.