Droplets sitting on thin elastic sheets: A study with the boundary element method

Technische Samenvatting: Druppels op dunne elastische vellen: Een studie met de Boundary Element Method

Probleemstelling
Het artikel onderzoekt de elasto-capillaire interactie tussen een vloeistofdruppel en een dun elastisch vel. In tegenstelling tot rigide substraten vervormen flexibele vaste stoffen onder invloed van capillaire krachten, wat leidt tot complexe evenwichtsconfiguraties waarbij de druppel in het vel zinkt en een lensachtige structuur vormt. Hoewel eerdere studies rimpelpatronen, durotaxis (beweging naar zachtere regio's) en anisotrope vormen onder spanning hebben verkend, was er behoefte aan een verenigd numeriek kader om de evenwichtsconfiguraties van druppels op elastische vellen onder verschillende randvoorwaarden en spanningsprotocollen te simuleren. Specifiek beogen de auteurs te kwantificeren hoe de dikte van het vel, isotrope spanning en uniaxiale rek de morfologie, stressverdeling en optische eigenschappen van de druppel beïnvloeden.

Methodologie
De auteurs breiden hun eerder ontwikkelde Boundary Element Method (BEM) uit om de gekoppelde dynamica van een druppel en een elastisch vel te simuleren. De kern van de methodologie omvat:

1. Dynamische Vergelijkingen: Het systeem wordt gemodelleerd met overgedempte dynamica waarbij de snelheidvelden van het druppeloppervlak ($\dot{R}_0$) en het elastische vel ($\dot{R}_{sh}$) gekoppeld zijn. De leidende vergelijking is een lineair systeem $G\dot{R} = -\nabla_R F$, waarbij $G$ de wrijvingsmatrix is en $F$ de totale vrije energie.
2. Vrije Energie Formulering:
   • Druppel ($F_0$): Omvat vloeistof-gas oppervlaktespanning en het werk verricht tegen het substraat-vloeistof interface, gestuurd door de Young-contacthoek. Een lijnspanningsterm is opgenomen maar gesteld op nul omdat deze verwaarloosbaar is voor de overwogen elastische vellen.
   • Vel ($F_{sh}$): Gemodelleerd met behulp van het Skalak-model voor in-plane afschuiving en oppervlakte-dilatatie (hyperelasticiteit) en het Helfrich-model voor out-of-plane buiging. De elastische moduli (afschuiving $\kappa_S$, oppervlakte $\kappa_A$ en buiging $\kappa_b$) zijn afgeleid van bulk materiaaleigenschappen (Youngs modulus $E$ en Poisson-ratio $\nu$) en de dikte van het vel $d$.
   • Restricties ($F_{cons}$): Afgedwongen via Lagrange-multiplicatoren om een constant druppelvolume te handhaven en de geometrische connectiviteit tussen de druppel- en de velmeshes bij het interface te waarborgen.
3. Wrijvingsmatrix: De matrix $G$ incorporeert viskeuze afschuifspanningen (afgeleid van de Oseen-tensor), contactlijnwrijving (Moffatt's resultaat) en vloeistof-substraatwrijving gestuurd door de Navier slip-conditie.
4. Numerieke Implementatie: Het systeem wordt gediscretiseerd met behulp van driehoekige meshes gegenereerd door gmsh. De auteurs gebruiken een gewogen interpolatieschema (met behulp van Surrogates.jl) om krachten en snelheden tussen de verschillende meshes van de druppel en het vel in kaart te brengen. Simulaties worden gedimensieloos gemaakt met behulp van de contactlijnstraal $R_0$ en karakteristieke tijdschalen afgeleid van de Cox-Voinov wet.

Belangrijkste Resultaten
De studie past de uitgebreide BEM toe op drie specifieke protocollen:

1. Vastgeklemd Elastisch Vel:

   • De evenwichtsconfiguratie van de druppel is een lens die steeds asymmetrischer wordt ten opzichte van het contactlijnvlak naarmate de dikte van het vel toeneemt.
   • Dikteafhankelijkheid: Naarmate de dikte ($d$) toeneemt, neemt het onderste deel van de druppel ($h_l$) af terwijl het bovenste deel ($h_u$) toeneemt, hoewel de totale hoogte relatief constant blijft.
   • Contacthoeken: De microscopische contacthoek bij de contactlijn komt overeen met de Young-hoek, wat het model valideert. Echter, de macroscopische contacthoek (bepaald door de algehele lensvorm) verschilt aanzienlijk van de Young-hoek en is afhankelijk van de dikte van het vel.
   • Ontkoppeling: De vorm van de onderste helft van de druppel wordt primair bepaald door de elasticiteit van het vel en is bijna onafhankelijk van de evenwichtscontacthoek ($\theta_{eq}$), terwijl $\theta_{eq}$ primair de bovenste helft beïnvloedt.

2. Isotroop Gestrekt Vel:

   • Het toepassen van isotrope spanning tilt de druppel op, waardoor de vervorming van het vel afneemt totdat het weer vlak wordt bij hoge rek ($\epsilon \approx 0.1$).
   • Stressprofielen: De auteurs kwantificeren radiale profielen van azimutale ($T_{\phi\phi}$) en radiale ($T_{rr}$) elastische spanningen. Spanningen nemen toe met de rek, en de contactlijn bevindt zich nabij een buigpunt in de stresscurves.
   • Optische Afstemming: De lensvormige druppel fungeert als een vloeistofoptische lens. De brandpuntsafstand, berekend via de lensmakerformule, kan reversibel met ongeveer 15% worden afgestemd door de toegepaste spanning te variëren. Dit bereik is vergelijkbaar met het afstemmingsbereik van lenzen in moderne smartphones.

3. Uniaxiaal Gestrekt Vel:

   • Het strekken van het vel langs één richting verlengt de druppel. De mate van verlenging hangt af van de randvoorwaarden van de laterale zijden (vastgeklemd, gecontroleerde inwaartse beweging, of vrij).
   • Morfologie: Vrije laterale zijden resulteren in de meest uitgesproken verlenging en significante verticale verplaatsingen, wat leidt tot opwaartse plooien en inkepingen in het vel. De contactlijn wordt niet-planair, verhoogd in de laterale richting en verlaagd in de strekrichting.
   • Anisotropie: De aspectratio van de contactlijn-voetafdruk en de anisotropie van de macroscopische contacthoek worden gekwantificeerd. De verlenging wordt toegeschreven aan het feit dat de druppel verschillende effectieve stijfheden van het vel ervaart langs de strekrichting en de laterale richting.

Betekenis en Claims
Het artikel beweert dat de uitgebreide BEM een robuust numeriek instrument biedt voor het simuleren van elasto-capillaire systemen, waarbij het de interactie tussen vloeistofdynamica en elastische vervormingen succesvol vastlegt. De primaire bijdragen zijn:

• Methodologisch: Een geformaliseerde BEM die de oppervlaktedynamica van een druppel koppelt aan de elastische vervormingen van een vel, inclusclusief in-plane en buigenergieën.
• Fysisch Inzicht: De demonstratie dat de dikte van het vel een cruciale parameter is die de asymmetrie van de druppel controleert, en dat de microscopische contacthoek consistent blijft met de Wet van Young ondanks macroscopische vervormingen.
• Toepassing: De studie toont expliciet aan dat het aanbrengen van spanning op een elastisch vel de brandpuntsafstand van vloeistoflenzen dynamisch en reversibel kan afstemmen. Dit biedt een potentieel mechanisme voor het controleren van optische eigenschappen in micro-optofluïdische apparaten.
• Toekomstige Richtingen: De auteurs suggereren dat de methode deuren opent voor het bestuderen van rimpelpatronen, durotaxis in spanning gradiënten en elasto-capillaire interacties tussen meerdere druppels, waarbij zij opmerken dat uniaxiale spanning kan leiden tot anisotrope interactiepotentialen en complexe dynamica.