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

Resumen Técnico: Análisis de Estabilidad Lineal de una Película Líquida Vertical sobre un Sustrato en Movimiento

Planteamiento del Problema
Este estudio investiga la estabilidad lineal de una película líquida vertical formada sobre un sustrato que se mueve contra la gravedad, una configuración central en el proceso industrial de recubrimiento por inmersión (dip-coating). En este proceso, una película líquida es extraída de un baño, y su estabilidad determina la uniformidad del recubrimiento final. Si bien la literatura previa ha abordado la estabilidad mediante expansiones asintóticas de onda larga (válidas solo para números de Reynolds pequeños y longitudes de onda largas) o modelos integrales de capa límite, ha faltado un análisis exhaustivo del problema de autovalores completo de Orr-Sommerfeld (OS) para flujos de película controlados y no controlados a través de un amplio rango de parámetros. Específicamente, el artículo aborda los mecanismos de inestabilidad y el umbral entre la inestabilidad absoluta y convectiva para líquidos que abarcan un amplio rango de números de Kapitza ($Ka$), desde fluidos de bajo $Ka$ como el aceite de maíz ($Ka \approx 4$) hasta metales líquidos de alto $Ka$ como el zinc ($Ka \approx 11525$).

Metodología
Los autores formulan el problema utilizando las ecuaciones de Navier-Stokes 2D incompresibles, adimensionalizadas utilizando una velocidad de referencia ($U_p$) y una escala de longitud de balance viscosidad-gravedad. El estado base se define como una película plana con una relación no monotónica entre el espesor de la película y el caudal, lo que genera ramas de película delgada ($\hat{h} < 1$) y gruesa ($\hat{h} > 1$).

Para analizar la estabilidad, los autores:

1. Linealizaron las ecuaciones gobernantes alrededor del estado base, derivando el problema de autovalores de Orr-Sommerfeld.
2. Resolvieron el problema de autovalores numéricamente utilizando el método espectral Chebyshev-Tau con $N=20$ y $N=80$ polinomios para asegurar la precisión y filtrar autovalores espurios.
3. Computaron curvas neutras y tasas de crecimiento ($\omega_i$) a través de los espacios de número de onda ($k$) y número de Reynolds ($Re$).
4. Investigaron los mecanismos de inestabilidad mediante expansiones asintóticas de onda larga, análisis de la distribución de vorticidad y un balance detallado de la energía cinética de las perturbaciones (considerando la energía cinética, la energía potencial de tensión superficial, el trabajo de cizalla, el esfuerzo de Reynolds y la disipación viscosa).
5. Determinaron el umbral Absoluto/Convectivo (AC) mediante el método de Briggs. Esto implicó el mapeo del espacio complejo del número de onda al espacio complejo de la frecuencia para identificar puntos de silla que satisfacen el criterio de causalidad (colisión), distinguiendo así entre flujos inestables convectivamente (perturbaciones arrastradas) y flujos absolutamente inestables (perturbaciones que crecen localmente).

Resultos Clave

• Estabilidad y Curvas Neutras: El estudio revela que los líquidos con bajo $Ka$ (por ejemplo, aceite de maíz) poseen un área de estabilidad más pequeña en comparación con los líquidos de alto $Ka$ (por ejemplo, zinc). Para películas gruesas ($\hat{h} > 1$), las curvas neutras para diferentes fluidos convergen, lo que sugiere un mecanismo estabilizador dominado por el balance de las fuerzas viscosas y gravitatorias donde los efectos de la tensión superficial disminuyen.
• Mecanismos de Inestabilidad:
  • Vorticidad y Desfase: Para películas delgadas, la inestabilidad es impulsada por la amplitud de la vorticidad con un desfase mínimo. Para películas gruesas, la amplitud de la vorticidad disminuye mientras que el desfase aumenta.
  • Papel de la Tensión Superficial: La tensión superficial exhibe un papel dual, tanto estabilizador como desestabilizador, particularmente para alto $Ka$. Para ondas largas, curva las líneas de vorticidad cerca del sustrato, reduciendo el flujo bajo las crestas (estabilizador). Para ondas cortas, fomenta la producción de vorticidad en la interfaz y crea regiones de vorticidad intensa cerca del sustrato (desestabilizador).
  • Balance de Energía: El análisis muestra que la tensión superficial contribuye tanto a la producción como a la disipación de la energía de la perturbación dependiendo del $Ka$. Para alto $Ka$ (zinc), la tensión superficial induce una producción significativa de vorticidad tanto en la superficie libre como en la pared, lo que conduce a intensos esfuerzos de Reynolds y efectos disipativos. Para bajo $Ka$ (aceite de maíz), los efectos de cizalla dominan los términos viscosos.
• Umbral Absoluto/Convectivo:
  • Existe una ventana de inestabilidad absoluta en el espacio de parámetros $\hat{h}-Re$, delimitada por umbrales inferiores y superiores.
  • Solución de Derjaguin ($\hat{h}=1$): El artículo encuentra que para $Ka < 17$, la solución de Derjaguin es siempre convectivamente inestable independientemente del número de Reynolds.
  • Solución de Landau-Levich-Derjaguin (LLD): Se encuentra que la solución LLD ($\hat{h} \approx 0.945 Re^{1/9} Ka^{-1/6}$) es siempre convectivamente inestable para el zinc, el agua y el agua-glicerol. Para el aceite de maíz, esta intersecta la ventana de inestabilidad absoluta solo en números de Reynolds altos ($100 < Re < 500$), un régimen considerado poco realista para condiciones experimentales típicas debido a la longitud de dominio requerida.
  • Escalamiento: Las ramas inferior y superior de la ventana de inestabilidad absoluta pueden reescalarse para colapsar los datos de diferentes fluidos ($Ka > 195$). Se propone una relación de ley de potencia específica ($\check{h}/\check{h}_{min} = (Re/Re_{min})^{4/9}$) para aproximar el umbral AC inferior.
  • Bifurcación: Para el agua-glicerol, se identifica un punto de bifurcación alrededor de $Re \approx 40$, donde aparece una segunda rama de solución, independiente de $Ka$ y gobernada por el balance gravedad-inercia.

Significancia y Reivindicaciones
Los autores afirman que este trabajo proporciona el primer análisis exhaustivo del problema de autovalores completo de Orr-Sommerfeld para flujos de recubrimiento por inmersión a través de un amplio rango de números $Ka$, abordando un vacío dejado por modelos asintóticos e integrales previos. El estudio clarifica los mecanismos físicos que impulsan la inestabilidad, específicamente el papel complejo y dual de la tensión superficial en la generación de vorticidad y el balance de energía.

Crucialmente, el artículo establece la existencia y ubicación de la ventana de inestabilidad absoluta en el espacio de parámetros. Esto es significativo para el diseño de leyes de control en procesos industriales de recubrimiento, ya que la naturaleza de la inestabilidad (absoluta vs. convectiva) dicta si las acciones de control deben afectar a toda la película o solo a regiones específicas. Los hallazgos sugieren que para fluidos de bajo $Ka$, la solución estándar de Derjaguin es inherentemente estable frente a la inestabilidad absoluta, mientras que para fluidos de alto $Ka$, el panorama de estabilidad es más complejo. Los resultados están destinados a servir como un referente para validar modelos de orden reducido y comprender las propiedades fundamentales de estabilidad de las películas líquidas en movimiento.