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

Sintesi Tecnica: Analisi della Stabilità Lineare di un Film Liquido Verticale su un Substrato in Movimento

Enunciato del Problema
Questo studio investiga la stabilità lineare di un film liquido verticale formato su un substrato che si muove contro la gravità, una configurazione centrale nel processo industriale di dip-coating. In questo processo, un film liquido viene estratto da un bagno e la sua stabilità determina l'uniformità del rivestimento finale. Sebbene la letteratura precedente abbia affrontato la stabilità tramite espansioni asintotiche di onde lunghe (valide solo per numeri di Reynolds bassi e lunghezze d'onda elevate) o modelli di strato limite integrale, è mancata un'analisi completa del problema degli autovalori di Orr-Sommerfeld (OS) per flussi di film controllati e non controllati attraverso un ampio intervallo di parametri. Nello specifico, il documento affronta i meccanismi di instabilità e la soglia tra instabilità assoluta e convettiva per liquidi che spaziano in un ampio intervallo di numeri di Kapitza ($Ka$), da fluidi a basso $Ka$ come l'olio di mais ($Ka \approx 4$) a metalli liquidi ad alto $Ka$ come lo zinco ($Ka \approx 11525$).

Metodologia
Gli autori formulano il problema utilizzando le equazioni di Navier-Stokes 2D incomprimibili, adimensionalizzate utilizzando una velocità di riferimento ($U_p$) e una scala di lunghezza basata sul bilancio viscosità-gravità. Lo stato base è definito come un film piatto con una relazione non monotona tra lo spessore del film e la portata, producendo rami di film sottili ($\hat{h} < 1$) e spessi ($\hat{h} > 1$).

Per analizzare la stabilità, gli autori:

1. Linearizzano le equazioni governanti attorno allo stato base, derivando il problema degli autovalori di Orr-Sommerfeld.
2. Risolvono il problema degli autovalori numericamente utilizzando il metodo spettrale Chebyshev-Tau con $N=20$ e $N=80$ polinomi per garantire l'accuratezza e filtrare gli autovalori spuri.
3. Calcolano le curve neutre e i tassi di crescita ($\omega_i$) attraverso gli spazi del numero d'onda ($k$) e del numero di Reynolds ($Re$).
4. Investigano i meccanismi di instabilità attraverso espansioni asintotiche di onde lunghe, analisi della distribuzione della vorticità e un bilancio dettagliato dell'energia cinetica delle perturbazioni (tenendo conto di energia cinetica, energia potenziale della tensione superficiale, lavoro di taglio, sforzo di Reynolds e dissipazione viscosa).
5. Determinano la soglia Assoluta/Convettiva (AC) utilizzando il metodo di Briggs. Ciò ha comportato la mappatura dello spazio del numero d'onda complesso nello spazio del numero di frequenza complesso per identificare punti di sella che soddisfano il criterio di causalità (collisione), distinguendo così i flussi instabili convettivi (disturbi trasportati via) dai flussi instabili assoluti (disturbi che crescono localmente).

Risultati Chiave

• Stabilità e Curve Neutre: Lo studio rivela che i liquidi con basso $Ka$ (ad es. olio di mais) possiedono un'area di stabilità più piccola rispetto ai liquidi ad alto $Ka$ (ad es. zinco). Per i film spessi ($\hat{h} > 1$), le curve neutre per diversi fluidi convergono, suggerendo un meccanismo di stabilizzazione dominato dal bilancio tra forze viscose e gravitazionali dove gli effetti della tensione superficiale diminuiscono.
• Meccanismi di Instabilità:
  • Vorticità e Sfasamento: Per i film sottili, l'instabilità è guidata dall'ampiezza della vorticità con uno sfasamento minimo. Per i film spessi, l'ampiezza della vorticità diminuisce mentre lo sfasamento aumenta.
  • Ruolo della Tensione Superficiale: La tensione superficiale presenta un ruolo duplice, sia stabilizzante che destabilizzante, particolarmente per alti $Ka$. Per le onde lunghe, essa curva le linee di vorticità vicino al substrato, riducendo il flusso sotto le creste (stabilizzante). Per le onde corte, favorisce la produzione di vorticità all'interfaccia e crea regioni di vorticità intensa vicino al substrato (destabilizzante).
  • Bilancio Energetico: L'analisi mostra che la tensione superficiale contribuisce sia alla produzione che alla dissipazione dell'energia della perturbazione a seconda del $Ka$. Per l'alto $Ka$ (zinco), la tensione superficiale induce una significativa produzione di vorticità sia alla superficie libera che alla parete, portando a intensi sforzi di Reynolds ed effetti dissipativi. Per il basso $Ka$ (olio di mais), gli effetti di taglio dominano i termini viscosi.
• Soglia Assoluta/Convettiva:
  • Esiste una finestra di instabilità assoluta nello spazio dei parametri $\hat{h}-Re$, delimitata da soglie inferiori e superiori.
  • Soluzione di Derjaguin ($\hat{h}=1$): Il documento trova che per $Ka < 17$, la soluzione di Derjaguin è sempre instabile convettiva indipendentemente dal numero di Reynolds.
  • Soluzione Landau-Levich-Derjaguin (LLD): La soluzione LLD ($\hat{h} \approx 0.945 Re^{1/9} Ka^{-1/6}$) si rivela essere sempre instabile convettiva per zinco, acqua e acqua-glicerina. Per l'olio di mais, essa interseca la finestra di instabilità assoluta solo ad alti numeri di Reynolds ($100 < Re < 500$), un regime considerato irrealistico per le tipiche condizioni sperimentali a causa della lunghezza del dominio richiesta.
  • Scaling: Le branche inferiore e superiore della finestra di instabilità assoluta possono essere riscalate per far convergere i dati per diversi fluidi ($Ka > 195$). Viene proposta una specifica relazione di legge di potenza ($\check{h}/\check{h}_{min} = (Re/Re_{min})^{4/9}$) per approssimare la soglia AC inferiore.
  • Biforcazione: Per l'acqua-glicerina, è stata identificata una soglia di biforcazione intorno a $Re \approx 40$, dove appare una seconda branca di soluzione, indipendente da $Ka$ e governata dal bilancio gravità-inerzia.

Significato e Rivendicazioni
Gli autori affermano che questo lavoro fornisce la prima analisi approfondita del problema degli autovalori di Orr-Sommerfeld per i flussi di dip-coating attraverso un ampio intervallo di numeri $Ka$, colmando una lacuna lasciata dai precedenti modelli asintotici e integrali. Lo studio chiarisce i meccanismi fisici che guidano l'instabilità, specificamente il complesso e duplice ruolo della tensione superficiale nella generazione di vorticità e nel bilancio energetico.

Fondamentalmente, il documento stabilisce l'esistenza e la posizione della finestra di instabilità assoluta nello spazio dei parametri. Ciò è significativo per la progettazione delle leggi di controllo nei processi industriali di rivestimento, poiché la natura dell'instabilità (assoluta vs convettiva) determina se le azioni di controllo debbano influenzare l'intero film o solo regioni specifiche. I risultati suggeriscono che per i fluidi a basso $Ka$, la soluzione standard di Derjaguin è intrinsecamente stabile contro l'instabilità assoluta, mentre per i fluidi ad alto $Ka$, il panorama della stabilità è più complesso. I risultati sono destinati a servire come benchmark per la validazione di modelli a ordine ridotto e per la comprensione delle proprietà fondamentali di stabilità dei film liquidi in movimento.