Linear stability analysis of a vertical liquid film over a moving substrate

Technische Samenvatting: Lineaire Stabiliteitsanalyse van een Verticale Vloeistofilm over een Bewegend Substraat

Probleemstelling
Deze studie onderzoekt de lineaire stabiliteit van een verticale vloeistofilm die gevormd wordt op een substraat dat tegen de zwaartekracht in beweegt, een configuratie die centraal staat in het industriële dip-coating proces. In dit proces wordt een vloeistoffilm uit een bad getrokken, en de stabiliteit ervan bepaalt de uniformiteit van de uiteindelijke coating. Hoewel eerdere literatuur de stabiliteit heeft behandeld via asymptotische langgolven-expansies (geldig voor alleen lage Reynoldsgetallen en lange golflengten) of integrale grenslaagmodellen, ontbrak een uitgebreide analyse van het volledige Orr-Sommerfeld (OS) eigenwaardeprobleem voor gecontroleerde en ongecontroleerde filmstromingen over een breed scala aan parameters. Specifiek adresseert dit artikel de instabiliteitsmechanismen en de drempel tussen absolute en convectieve instabiliteit voor vloeistoffen die een breed bereik aan Kapitza-getallen ($Ka$) beslaan, van laag-$Ka$ vloeistoffen zoals maïsolie ($Ka \approx 4$) tot hoog-$Ka$ vloeibare metalen zoals zink ($Ka \approx 11525$).

Methodologie
De auteurs formuleren het probleem met behulp van de incompressibele 2D Navier-Stokes vergelijkingen, gedimensioneerd met behulp van een referentiesnelheid ($U_p$) en een viscositeits-zwaartekrachtbalans lengteschaal. De basistoestand wordt gedefinieerd als een vlakke film met een niet-monotone relatie tussen filmdikte en debiet, wat resulteert in dunne ($\hat{h} < 1$) en dikke ($\hat{h} > 1$) filmtakken.

Om de stabiliteit te analyseren, hebben de auteurs:

1. De leidende vergelijkingen gelineariseerd rond de basistoestand, waarbij het Orr-Sommerfeld eigenwaardeprobleem is afgeleid.
2. Het eigenwaardeprobleem numeriek opgelost met behruik van de Chebyshev-Tau spectrale methode met $N=20$ en $N=80$ polynomen om nauwkeurigheid te garanderen en spookeigenwaarden te filteren.
3. Neutrale curven berekend en groeisnelheden ($\omega_i$) over de golfgetal- ($k$) en Reynoldsgetal- ($Re$) ruimtes.
4. Instabiliteitsmechanismen onderzocht via langgolven asymptotische expansies, vorticiteitsverdeling-analyse en een gedetailleerde kinetische energiebalans van de perturbaties (rekening houdend met kinetische energie, oppervlaktespanning potentiele energie, afschuivingsarbeid, Reynolds-spanning en viskeuze dissipatie).
5. De Absolute/Convectieve (AC) drempel bepaald met behulp van Briggs' methode. Dit omvatte het mappen van de complexe golfgetalruimte naar de complexe frequentieruimte om zadelpunten te identificeren die voldoen aan het causaliteitscriterium (botsingscriterium), waardoor onderscheid kon worden gemaakt tussen convectief instabiele stromingen (waarbij verstoringen worden weggevoerd) en absoluut instabiele stromingen (waarbij verstoringen lokaal groeien).

Belangrijkste Resultaten

• Stabiliteit en Neutrale Curven: De studie laat zien dat vloeistoffen met een lage $Ka$ (bijv. maïsolie) een kleiner stabiliteitsgebied bezitten vergeleken met hoog-$Ka$ vloeistoffen (bijv. zink). Voor dikke films ($\hat{h} > 1$) convergeren de neutrale curven voor verschillende vloeistoffen, wat suggereert dat er een stabiliserend mechanisme aanwezig is dat gedomineerd wordt door de balans tussen viskeuze en gravitationele krachten, waarbij de effecten van oppervlaktespanning afnemen.
• Instabiliteitsmechanismen:
  • Vorticiteit en Faseverschuiving: Voor dunne films wordt instabiliteit gedreven door de amplitude van de vorticiteit met minimale faseverschuiving. Voor dikke films neemt de vorticiteitsamplitude af terwijl de faseverschuiving toeneemt.
  • Rol van Oppervlaktespanning: Oppervlaktespanning vertoont een dubbelvoudige stabiliserende en destabiliserende rol, met name voor hoge $Ka$. Voor lange golven buigt oppervlaktespanning de vorticiteitslijnen nabij het substraat af, wat de stroming onder de kammen vermindert (stabiliserend). Voor korte golven bevordert het de productie van vorticiteit bij het grensvlak en creëert het intense vorticiteitsregio's nabij het substraat (destabiliserend).
  • Energiebalans: De analyse toont aan dat oppervlaktespanning bijdraagt aan zowel de productie als de dissipatie van perturbatie-energie, afhankelijk van de $Ka$. Voor hoge $Ka$ (zink) induceert oppervlaktespanning significante vorticiteitsproductie bij zowel het vrije oppervlak als de wand, wat leidt tot intense Reynolds-spanningen en dissipatieve effecten. Voor lage $Ka$ (maïsolie) domineren afschuifeffecten de viskeuze termen.
• Absolute/Convectieve Drempel:
  • Er bestaat een venster van absolute instabiliteit in de $\hat{h}-Re$ parameterrumte, begrensd door lagere en hogere drempels.
  • Derjaguin's Oplossing ($\hat{h}=1$): Het artikel stelt vast dat voor $Ka < 17$ de Derjaguin-oplossing altijd convectief instabiel is, ongeacht het Reynoldsgetal.
  • Landau-Levich-Derjaguin (LLD) Oplossing: De LLD-oplossing ($\hat{h} \approx 0.945 Re^{1/9} Ka^{-1/6}$) wordt gevonden als altijd convectief instabiel voor zink, water en water-glycerine. Voor maïsolie snijdt deze de absolute instabiliteitsvenster alleen bij hoge Reynoldsgetallen ($100 < Re < 500$), een regime dat als onrealistisch wordt beschouwd voor typische experimentele condities vanwege de vereiste domeinlengte.
  • Schaling: De lagere en hogere takken van de absolute instabiliteitsvenster kunnen worden geschaald om data voor verschillende vloeistoffen te laten samenvallen ($Ka > 195$). Een specifieke machtswet-relatie ($\check{h}/\check{h}_{min} = (Re/Re_{min})^{4/9}$) wordt voorgesteld om de lagere AC-drempel te benaderen.
  • Bifurcatie: Voor water-glycerine is een bifurcatiepunt geïdentificeerd rond $Re \approx 40$, waar een tweede oplossingstak verschijnt, onafhankelijk van $Ka$ en beheerst door de balans tussen zwaartekracht en traagheid.

Significantie en Claims
De auteurs claimen dat dit werk de eerste grondige analyse biedt van het volledige Orr-Sommerfeld eigenwaardeprobleem voor dip-coating stromingen over een breed bereik van $Ka$-getallen, waarmee een gat wordt opgevuld dat is achtergelaten door eerdere asymptotische en integrale modellen. De studie verheldert de fysieke mechanismen die instabiliteit drijven, specif gezien de complexe, tweeledige rol van oppervlaktespanning in de generatie van vorticiteit en de energiebalans.

Cruciaal is dat het papier de existentie en locatie van de absolute instabiliteitsvenster in de parameterrumte vaststelt. Dit is significant voor het ontwerp van controlewetten in industriële coatingprocessen, aangezien de aard van de instabiliteit (absoluut versus convectief) bepaalt of controleacties de gehele film of slechts specifieke regio's moeten beïnvloeden. De bevindingen suggereren dat voor laag-$Ka$ vloeistoffen de standaard Derjaguin-oplossing inherent stabiel is tegen absolute instabiliteit, terwijl voor hoog-$Ka$ vloeistoffen het stabiliteitslandschap complexer is. De resultaten zijn bedoeld om te dienen als benchmark voor het valideren van reduced-order modellen en het begrijpen van de fundamentele stabiliteitseigenschappen van bewegende vloeistoffilms.