Internal flow and concentration in neighbouring evaporating binary droplets and rivulets

Este estudio numérico y teórico analiza cómo la proximidad entre gotas o riachuelos binarios adyacentes induce un efecto de apantallamiento que reduce la tasa de evaporación y rompe la simetría de los campos de flujo y concentración, revelando que la posición del punto de estancamiento interfacial depende del ángulo de contacto, la distancia entre ellos y, en el caso de las gotas, de los efectos de Marangoni debido al flujo azimutal adicional.

Lo esencial

¡Hola! Imagina que estás en una cocina y has dejado caer dos gotas de un líquido especial (una mezcla de agua y un alcohol) muy cerca una de la otra sobre una mesa. Con el tiempo, el agua se evapora, pero el alcohol no. Lo que hace este artículo es estudiar qué pasa "por dentro" de esas gotas mientras se secan, especialmente cuando están tan juntas que se "molestan" entre sí.

Aquí te explico los hallazgos principales usando analogías sencillas:

1. El problema de las "gotas vecinas"

Imagina que tienes dos personas gritando en una habitación vacía. Si están lejos, cada una se oye perfectamente. Pero si se ponen muy cerca, sus voces se mezclan y se "taponean" entre sí; el sonido de una llega más débil a la otra.

En el mundo de las gotas, pasa algo similar con el vapor. Cuando dos gotas están cerca, el vapor que sale de una se acumula en el espacio entre ellas. Esto crea una especie de "nube de humedad" local que hace que sea más difícil para la otra gota evaporarse. A esto los científicos lo llaman efecto de escudo.

2. El baile interior: El efecto "Marangoni"

Ahora, imagina que dentro de esas gotas hay un baile. Como el agua se evapora más rápido que el alcohol, la composición de la gota cambia. Donde hay más alcohol, la "tensión superficial" (la piel de la gota) es más fuerte; donde hay más agua, es más débil.

Esto crea una fuerza invisible que arrastra el líquido desde las zonas débiles hacia las fuertes. Es como si alguien tirara de una manta desde un lado: el líquido se mueve para compensar. A este movimiento lo llamamos flujo de Marangoni. En una sola gota, este baile suele ser simétrico (se mueve igual hacia todos lados).

3. El caos cuando hay dos gotas

Aquí es donde se pone interesante. Cuando las dos gotas están juntas:

• La evaporación no es igual en todos lados: la parte de la gota que mira a su vecina se evapora más lento (por el "escudo" de vapor).
• Esto rompe la simetría del baile. El líquido ya no gira uniformemente. Se crea un desequilibrio.

Los autores descubrieron que este desequilibrio se puede medir mirando un punto específico dentro de la gota llamado punto de estancamiento. Imagina un río que fluye hacia un lado y luego gira; el punto de estancamiento es el lugar exacto donde el agua se detiene antes de girar. En las gotas vecinas, este punto se mueve hacia el centro, pero no se queda quieto; depende de qué tan inclinada esté la gota y qué tan lejos esté su vecina.

4. La diferencia entre "Gotas" y "Rivulets" (Canales)

El estudio comparó dos formas:

• Gotas (Droplets): Son como pequeñas montañas redondas.
• Rivulets: Son como canales largos y estrechos (imagina una gota estirada como una laguna).

El hallazgo sorprendente:

• En los canales (rivulets), el punto de estancamiento se mueve dependiendo de la forma de la gota y la distancia al vecino, pero no le importa qué tan fuerte sea la fuerza del baile (el número de Marangoni). Es como si el canal fuera tan estrecho que el líquido solo tiene una dirección para ir: adelante o atrás.
• En las gotas redondas, el punto de estancamiento sí depende de la fuerza del baile. ¿Por qué? Porque las gotas redondas tienen una "tercera dimensión" de movimiento (pueden girar en círculo, como un remolino). Esta libertad extra hace que la mezcla sea más compleja y que el punto de estancamiento reaccione a la intensidad de la fuerza.

5. ¿Por qué importa esto?

Esto no es solo teoría. Piensa en:

• Impresoras de inyección de tinta: Si las gotas de tinta se secan de forma desigual, la imagen sale mal.
• Rociado de pesticidas: Si las gotas de veneno no se distribuyen bien, las plagas no mueren.
• Refrigeración por spray: Necesitas que el líquido se evapore uniformemente para enfriar bien.

En resumen

El artículo nos dice que cuando dos gotas de una mezcla se secan juntas, se crean corrientes internas desordenadas porque se "taponean" entre sí.

• Si son canales largos, el desorden depende de la distancia y la forma, pero es predecible.
• Si son gotas redondas, el desorden es más complicado porque pueden girar en todas direcciones, y la fuerza interna del líquido juega un papel crucial.

Los científicos han creado modelos matemáticos (como recetas simplificadas) para predecir exactamente dónde se moverá ese punto de estancamiento, lo cual ayuda a diseñar mejores tecnologías de impresión y rociado en el futuro. ¡Es como aprender a predecir el clima dentro de una gota de agua!

Resumen técnico

Aquí presento un resumen técnico detallado del artículo en español, estructurado según los puntos solicitados:

Título: Flujo interno y concentración en gotas y riachuelos binarios adyacentes en evaporación

1. Planteamiento del Problema

El estudio aborda la evaporación de gotas o riachuelos (rivulets) adyacentes de mezclas binarias, un fenómeno común en aplicaciones como la impresión por inyección de tinta, el enfriamiento por pulverización y la entrega de pesticidas.

• El desafío principal: Cuando múltiples cuerpos líquidos están en proximidad, se produce un efecto de blindaje (shielding effect). La concentración de vapor en la fase gaseosa aumenta cerca de los cuerpos líquidos, lo que reduce la tasa de evaporación individual y altera la dinámica del flujo interno.
• Complejidad física: A diferencia de las gotas puras (donde domina el efecto "mancha de café" o coffee-stain), las mezclas binarias (ej. agua y 1,2-hexanediol) experimentan evaporación selectiva. Esto genera gradientes de tensión superficial que inducen flujos de Marangoni. La interacción entre el blindaje de vapor, la asimetría inducida por la proximidad y los flujos de Marangoni rompe la simetría del campo de flujo y de concentración, haciendo que los modelos axisimétricos tradicionales sean insuficientes.
• Limitación computacional: Resolver el problema completo tridimensional (3D) para gotas adyacentes es computacionalmente prohibitivo debido a la necesidad de alta resolución para capturar gradientes de concentración y tasas de evaporación.

2. Metodología

Los autores emplean una combinación de simulaciones numéricas directas y modelos teóricos simplificados:

• Simulaciones Transientes (2D y 3D):

  • Se utiliza un método de elementos finitos Lagrangiano-Euleriano de interfaz nítida (sharp-interface) implementado en el paquete pyoomph.
  • Se resuelven las ecuaciones completas de Navier-Stokes, advección-difusión y difusión de vapor para riachuelos (2D) como punto de partida, y se extrapola a gotas (3D).
  • Se asume que solo un componente es volátil, la evaporación es lenta y los efectos térmicos son despreciables frente a los efectos composicionales.

• Modelo Cuasi-Estacionario:

  • Dado que la evaporación es lenta y el flujo de Marangoni induce una mezcla suficiente, se asume que el sistema es cuasi-estacionario en cada instante.
  • Se linealizan las propiedades alrededor de la composición media y se desprecia el flujo de "mancha de café" (asumiendo números de Marangoni altos).

• Modelo de Lubricación:

  • Para reducir la complejidad, se aplica la aproximación de lubricación (válido para ángulos de contacto bajos) a las ecuaciones cuasi-estacionarias.
  • Se incorpora la dispersión de Taylor-Aris para modelar la mezcla vertical de la composición.
  • Se utilizan expresiones analíticas para el flujo de evaporación local, considerando el efecto de blindaje entre los cuerpos adyacentes.

• Parámetros de Control:

  • Ángulo de contacto ($\theta$).
  • Distancia entre centros ($b$).
  • Número de Marangoni ($Ma$), que cuantifica la fuerza de los gradientes de tensión superficial.

3. Contribuciones Clave

1. Desarrollo de modelos simplificados: Se valida un modelo cuasi-estacionario y una versión de lubricación que capturan la física esencial de la evaporación asimétrica sin el costo computacional de las simulaciones transientes 3D completas.
2. Identificación de la asimetría: Se cuantifica la asimetría del flujo mediante la posición del punto de estancamiento ($x_0$) en la interfaz o cerca del sustrato, donde dos vórtices de flujo se encuentran.
3. Diferenciación entre Geometrías: Se establece una distinción fundamental en el comportamiento de los riachuelos (2D) frente a las gotas (3D) respecto a la dependencia del número de Marangoni.

4. Resultados Principales

• Validación del Modelo Cuasi-Estacionario:

  • Existe un acuerdo perfecto entre las simulaciones transientes completas y el modelo cuasi-estacionario para riachuelos, siempre que el número de Marangoni sea suficientemente alto para superar el transporte por el efecto "mancha de café".

• Comportamiento en Riachuelos (2D):

  • Independencia de $Ma$: El desplazamiento del punto de estancamiento ($x_0$) en riachuelos es independiente del número de Marangoni. La asimetría está controlada únicamente por el ángulo de contacto ($\theta$) y la distancia entre riachuelos ($b$).
  • Efecto del Ángulo y Distancia: A medida que aumenta el ángulo de contacto o disminuye la distancia, el punto de estancamiento se desplaza hacia el centro, pero la asimetría persiste incluso a grandes distancias.
  • Solución Analítica: Se derivó una expresión analítica para el desplazamiento $x_0$ en riachuelos que confirma su independencia de $Ma$ y su dependencia de la geometría y el flujo de evaporación.

• Comportamiento en Gotas (3D):

  • Dependencia de $Ma$: A diferencia de los riachuelos, en las gotas el punto de estancamiento sí depende del número de Marangoni.
  • Mecanismo: La presencia de un grado de libertad adicional en la dirección azimutal permite flujos de Marangoni más complejos. Esto genera una evolución no lineal de los gradientes de concentración, lo que hace que la simetría del flujo sea sensible a la intensidad de las fuerzas de Marangoni.
  • Gradientes Verticales: A altos valores de $Ma$, se observan grandes gradientes verticales en la velocidad, haciendo que la posición del punto de estancamiento varíe significativamente entre la interfaz y el sustrato.

• Limitaciones del Modelo de Lubricación:

  • El modelo de lubricación predice bien la tendencia general, pero subestima el desplazamiento del punto de estancamiento en ángulos de contacto altos. Esto se debe principalmente a que las expresiones analíticas actuales para el flujo de evaporación local asumen ángulos de contacto nulos, ignorando la dependencia real con $\theta$.

5. Significado e Impacto

• Avance Teórico: El trabajo proporciona una comprensión fundamental de cómo la proximidad altera la dinámica interna de fluidos multicomponente, desafiando la suposición de simetría axial en sistemas de gotas múltiples.
• Herramienta de Diseño: Los modelos simplificados (cuasi-estacionario y de lubricación) ofrecen herramientas computacionalmente eficientes para predecir el comportamiento de evaporación en aplicaciones industriales (impresión, recubrimientos) donde el control de la deposición y la mezcla es crítico.
• Insight Físico: La distinción crucial entre la independencia de $Ma$ en 2D (riachuelos) y su fuerte dependencia en 3D (gotas) revela la importancia de la dimensionalidad y la libertad azimutal en la dinámica de fluidos confinados.
• Futuro: El estudio señala la necesidad urgente de desarrollar expresiones analíticas para el flujo de evaporación local en gotas con ángulos de contacto no nulos para mejorar la precisión de los modelos de lubricación.

En resumen, el artículo demuestra que la asimetría en la evaporación de gotas binarias adyacentes es un fenómeno gobernado por una interacción compleja entre la geometría, la distancia y los flujos de Marangoni, con comportamientos cualitativamente diferentes dependiendo de si el sistema es bidimensional (riachuelos) o tridimensional (gotas).