Análisis Espectral de la Línea de Contacto en un Flujo Axisimétrico 3D Impulsado por Esfuerzos Superficiales de Marangoni

Abstract

Se investiga numéricamente la evolución de una gota cuando es depositada en el centro de un sustrato sólido, circular y plano cuya temperatura disminuye logarítmicamente con el radio. El gradiente térmico sobre el sustrato induce un gradiente de temperatura sobre la superficie del líquido, por lo que sobre esta se genera un esfuerzo de Marangoni. Dicho esfuerzo impulsa al líquido hacia la periferia del sustrato. Resultados experimentales publicados muestran que la línea de contacto avanza manteniendo inicialmente la simetría axial, pero a partir de cierto punto se desestabiliza y presenta ondulaciones que crecen con el tiempo. Con el objetivo de estudiar dicha inestabilidad, se obtienen soluciones numéricas tridimensionales de la ecuación que gobierna al fenómeno para diferentes condiciones iniciales, se detecta en cada caso la posición de la línea de contacto, y por medio de un análisis espectral se determina la evolución de cada modo normal azimutal, sus tasas de crecimiento y la posible interacción entre ellos. Para resolver numéricamente la ecuación se usó el método de elementos finitos y un esquema de diferenciación inversa con paso temporal adaptativo. Se utilizó un solver iterativo (Generalized Minimal Residuals - GMRES), preacondicionado por medio de un algoritmo multigrilla estándar MUMPS. Se utilizó un criterio de convergencia de 10−3 para el error relativo definido por una norma Euclidiana ponderada para dos iteraciones sucesivas. El dominio se discretizó con elementos triangulares no estructurados, con tamaños máximo y mínimo de 0,01 y 0,0008, respectivamente.