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

Resumo Técnico: Análise de Estabilidade Linear de um Filme Líquido Vertical sobre um Substrato em Movimento

Enunciado do Problema
Este estudo investiga a estabilidade linear de um filme líquido vertical formado sobre um substrato que se move contra a gravidade, uma configuração central para o processo industrial de dip-coating (revestimento por imersão). Neste processo, um filme líquido é extraído de um banho, e sua estabilidade determina a uniformidade do revestimento final. Embora a literatura anterior tenha abordado a estabilidade por meio de expansões de ondas longas assintóticas (válidas apenas para números de Reynolds baixos e comprimentos de onda longos) ou modelos de camada limite integral, carecia-se de uma análise abrangente do problema de autovalor de Orr-Sommerfeld completo para fluxos de filmes controlados e não controlados através de uma ampla gama de parâmetros. Especificamente, o artigo aborda os mecanismos de instabilidade e o limiar entre a instabilidade absoluta e convectiva para líquidos que abrangem uma vasta gama de números de Kapitza ($Ka$), desde fluidos de baixo $Ka$, como o óleo de milho ($Ka \approx 4$), até metais líquidos de alto $Ka$, como o zinco ($Ka \approx 11525$).

Metodologia
Os autores formulam o problema utilizando as equações de Navier-Stokes 2D incompressíveis, adimensionalizadas usando uma velocidade de referência ($U_p$) e uma escala de comprimento de balanço viscoso-gravitacional. O estado base é definido como um filme plano com uma relação não monotônica entre a espessura do filme e a taxa de fluxo, resultando em ramos de filmes finos ($\hat{h} < 1$) e espessos ($\hat{h} > 1$).

Para analisar a estabilidade, os autores:

1. Linearizaram as equações governantes em torno do estado base, derivando o problema de autovalor de Orr-Sommerfeld.
2. Resolveram o problema de autovalor numericamente usando o método espectral Chebyshev-Tau com $N=20$ e $N=80$ polinômios para garantir a precisão e filtrar autovalores espúrios.
3. Computaram curvas neutras e taxas de crescimento ($\omega_i$) nos espaços de número de onda ($k$) e número de Reynolds ($Re$).
4. Investigaram os mecanismos de instabilidade por meio de expansões assintóticas de ondas longas, análise de distribuição de vorticidade e um balanço detalhado de energia cinética das perturbações (contabilizando energia cinética, energia potencial de tensão superficial, trabalho de cisalhamento, tensão de Reynolds e dissipação viscosa).
5. Determinaram o limiar Absoluto/Convectivo (AC) usando o método de Briggs. Isso envolveu o mapeamento do espaço de número de onda complexo para o espaço de frequência complexa para identificar pontos de sela que satisfaçam o critério de causalidade (colisão), distinguindo assim fluxos instáveis convectivamente (perturbações varridas para longe) de fluxos absolutamente instáveis (perturbações crescendo localmente).

Principais Resultados

• Estabilidade e Curvas Neutras: O estudo revela que líquidos com baixo $Ka$ (ex: óleo de milho) possuem uma área de estabilidade menor comparada a líquidos de alto $Ka$ (ex: zinco). Para filmes espessos ($\hat{h} > 1$), as curvas neutras para diferentes fluidos convergem, sugerindo um mecanismo estabilizador dominado pelo balanço entre forças viscosas e gravitacionais, onde os efeitos da tensão superficial diminuem.
• Mecanismos de Instabilidade:
  • Vorticidade e Deslocamento de Fase: Para filmes finos, a instabilidade é impulsionada pela amplitude da vorticidade com deslocamento de fase mínimo. Para filmes espessos, a amplitude da vorticidade diminui enquanto o deslocamento de fase aumenta.
  • Papel da Tensão Superficial: A tensão superficial exibe um papel duplo, tanto estabilizador quanto desestabilizador, particularmente para alto $Ka$. Para ondas longas, ela curva as linhas de vorticidade próximo ao substrato, reduzindo o fluxo sob as cristas (estabilizando). Para ondas curtas, ela promove a produção de vorticidade na interface e cria regiões de vorticidade intensa próximo ao substrato (desestabilizando).
  • Balanço de Energia: A análise mostra que a tensão superficial contribui tanto para a produção quanto para a dissipação da energia da perturbação dependendo do $Ka$. Para alto $Ka$ (zinco), a tensão superficial induz produção significativa de vorticidade tanto na superfície livre quanto na parede, levando a tensões de Reynolds intensas e efeitos dissipativos. Para baixo $Ka$ (óleo de milho), os efeitos de cisalhamento dominam os termos viscosos.
• Limiar Absoluto/Convectivo:
  • Uma janela de instabilidade absoluta existe no espaço de parâmetros $\hat{h}-Re$, delimitada por limiares inferiores e superiores.
  • Solução de Derjaguin ($\hat{h}=1$): O artigo encontra que, para $Ka < 17$, a solução de Derjaguin é sempre instável convectivamente, independentemente do número de Reynolds.
  • Solução de Landau-Levich-Derjaguin (LLD): A solução LLD ($\hat{h} \approx 0.945 Re^{1/9} Ka^{-1/6}$) é encontrada como sempre instável convectivamente para zinco, água e água-glicerina. Para o óleo de milho, ela intercepta a janela de instabilidade absoluta apenas em altos números de Reynolds ($100 < Re < 500$), um regime considerado irrealista para condições experimentais típicas devido ao comprimento de domínio exigido.
  • Escalonamento: Os ramos inferior e superior da janela de instabilidade absoluta podem ser reescalonados para colapsar dados de diferentes fluidos ($Ka > 195$). Uma relação específica de lei de potência ($\check{h}/\check{h}_{min} = (Re/Re_{min})^{4/9}$) é proposta para aproximar o limiar AC inferior.
  • Bifurcação: Para água-glicerina, um ponto de bifurcação é identificado em torno de $Re \approx 40$, onde um segundo ramo de solução aparece, independente de $Ka$ e governado pelo balanço gravidade-inércia.

Significância e Alegações
Os autores afirmam que este trabalho fornece a primeira análise completa do problema de autovalor de Orr-Sommerfeld para fluxos de dip-coating através de uma ampla gama de números de $Ka$, abordando uma lacuna deixada por modelos anteriores assintóticos e integrais. O estudo esclarece os mecanismos físicos que impulsionam a instabilidade, especificamente o papel complexo e duplo da tensão superficial na geração de vorticidade e no balanço de energia.

Crucialmente, o artigo estabelece a existência e a localização da janela de instabilidade absoluta no espaço de parâmetros. Isso é significativo para o design de leis de controle em processos industriais de revestimento, pois a natureza da instabilidade (absoluta vs. convectiva) dita se as ações de controle devem afetar todo o filme ou apenas regiões específicas. Os resultados sugerem que, para fluidos de baixo $Ka$, a solução padrão de Derjaguin é inerentemente estável contra a instabilidade absoluta, ao passo que para fluidos de alto $Ka$, o cenário de estabilidade é mais complexo. Os resultados pretendem servir como um referencial para validar modelos de ordem reduzida e compreender as propriedades fundamentais de estabilidade de filmes líquidos em movimento.