A reduced model for droplet dynamics with interfacial viscosity

Este artículo extiende el modelo de Maffettone-Minale para incorporar las viscosidades de cizalladura e dilatación interfaciales para describir la deformación de gotas en flujo de cizalladura y valida la precisión del modelo extendido resultante frente a simulaciones numéricas completamente resueltas a través de diversos números de Capilaridad y de Boussinesq.

Lo esencial

Imagina una diminuta gota de aceite flotando en una corriente de agua. Si el agua empieza a fluir más rápido, la gota se aplasta y se estira, pasando de ser una esfera perfecta a una forma ovalada. Los científicos han intentado predecir durante mucho tiempo con exactitud cuánto se estira esta gota usando las matemáticas.

Durante décadas, utilizaron una "receta" famosa (llamada modelo de Maffettone–Minale) que funciona bien para gotas limpias. Pero en el mundo real, las gotas suelen tener una "piel" hecha de jabón, proteínas u otras moléculas. Esta piel no es solo un límite; tiene su propio grosor y viscosidad, conocida como viscosidad interfacial. Piensa en esto como si la gota llevara puesto un suéter pegajoso y elástico.

Este artículo presenta una nueva receta mejorada (el modelo EMM o Maffettone–Minale Extendido) que tiene en cuenta este suéter pegajoso. Así es como los autores lo desglosaron:

1. Los dos tipos de "pegajosidad"

Los autores se dieron cuenta de que la piel de la gota resiste el movimiento de dos maneras diferentes, y necesitaban medir ambas:

• Viscosidad de cizalladura (El efecto "banda elástica"): Imagina que intentas deslizar tu mano sobre la superficie de la gota. Si la piel es "viscosa por cizalladura", resiste ese movimiento de deslizamiento, como si arrastraras tu mano a través de la miel.
• Viscosidad dilatacional (El efecto "respiración"): Imagina que la gota intenta expandir o encoger su área superficial (como un globo inflándose). Si la piel es "viscosa por dilatación", resiste ese estiramiento o encogimiento, como una tela tensa y rígida que no quiere expandirse.

El artículo utiliza números especiales (llamados números de Boussinesq) para medir qué tan fuertes son estas dos resistencias en comparación con el grosor de la gota.

2. La nueva receta (El modelo EMM)

Los autores tomaron la vieja y simple receta matemática y añadieron nuevos ingredientes para manejar estos dos tipos de pegajosidad.

• El objetivo: Querían saber: ¿Hasta qué punto podemos estirar esta nueva receta antes de que deje de funcionar?
• El método: No se limitaron a adivinar. Construyeron una simulación computacional superdetallada (como una película de alta definición de la gota) que resolvía cada pequeña regla de la física desde cero. Esto sirvió como la "verdad".
• La prueba: Ejecutaron la nueva receta EMM junto con la simulación superdetallada. Compararon los resultados para ver si la receta simple coincidía con la película compleja.

3. Lo que encontraron

Los resultados fueron sorprendentes y específicos:

• Cuando el "suéter" es uniforme: Si la piel de la gota resiste el deslizamiento y el estiramiento por igual (un suéter equilibrado), la nueva receta funciona increíblemente bien, incluso cuando la gota se está estirando bastante. Predice con exactitud qué tan rápido se estira la gota y cuánto tiempo tarda en establecer su forma final.
• Cuando el "suéter" está desequilibrado: Si la piel es muy buena resistiendo el deslizamiento pero mala resistiendo el estiramiento (o viceversa), la receta simple empieza a volverse un poco imprecisa. Sigue funcionando para flujos suaves, pero si el flujo se vuelve demasiado fuerte, la receta pierde precisión.
• El efecto de "ralentización": El hallazgo más interesante fue sobre el tiempo. Cuando la gota tiene ambos tipos de pegajosidad al mismo tiempo, tarda mucho más en cambiar de forma. Es como si la gota se quedara "atascada" en su propia piel. Los autores descubrieron que su nueva receta captura perfectamente este efecto de "cámara lenta".
• El punto de ruptura: Si la gota tiene casi ninguna resistencia al deslizamiento (pero una alta resistencia al estiramiento), se estira tanto que eventualmente se rompe. La nueva receta predice correctamente que esto sucede antes en estas condiciones específicas.

4. La conclusión fundamental

Los autores crearon con éxito una herramienta simple, rápida y fiable para predecir cómo se comportan las gotas con "pieles pegajosas" en líquidos en movimiento.

• Por qué es importante: Ahorra a los científicos el tener que ejecutar simulaciones computacionales masivas y lentas para cada uno de los problemas.
• El matiz: La herramienta es muy precisa para estiramientos pequeños o medianos, especialmente cuando la resistencia de la piel está equilibrada. Si el flujo es extremadamente violento o la resistencia de la piel está muy desequilibrada, la herramienta empieza a perder su precisión y se necesitan las simulaciones computacionales de gran potencia.

En resumen, actualizaron la "calculadora de gotas" para manejar pieles pegajosas, demostrando que funciona de maravilla para la mayoría de los escenarios cotidianos, mientras marcan claramente los límites donde podría necesitar un poco de ayuda.

Resumen técnico

Resumen Técnico: Un Modelo Reducido para la Dinámica de Gotas con Viscosidad Interfacial

Planteamiento del Problema
Comprender la deformación de las gotas en flujos hidrodinámicos es un tema central en la reología. Si bien el modelo clásico de Maffettone–Minale (MM) proporciona una descripción fenomenológica exitosa para gotas limpias (aquellas sin viscosidad interfacial) en flujo de cizalla, no da cuenta de las complejidades introducidas por la viscosidad interfacial. En muchos escenarios prácticos, como cuando las interfaces están recubiertas con surfactantes, polímeros o proteínas, las viscosidades de cizalla interfacial ($\mu_s$) y de dilatación ($\mu_d$) se vuelven significativas. Estas viscosidades introducen una fricción adicional que altera el equilibrio de tensiones entre los fluidos voluminosos, influyendo en la deformación, la relajación y los umbrales de ruptura. El desafío radica en desarrollar un modelo de orden reducido que incorpore estos efectos interfaciales manteniendo la eficiencia computacional y la tractabilidad analítica del marco MM original.

Metodología
Los autores proponen un modelo de Maffettone–Minale Extendido (EMM), el cual generaliza la ecuación MM original para la evolución de un tensor morfológico $S$ (que describe la forma de la gota como un elipsoide).

• Extensión Teórica: El modelo EMM incorpora las viscosidades interfaciales a través de los números de Boussinesq, $Bqs = \mu_s/\mu R$ y $Bqd = \mu_d/\mu R$. Basándose en cálculos perturbativos para interfaces viscosas (específicamente la extensión de Narsimhan al trabajo de Rallison), los autores derivan expresiones generalizadas para los coeficientes del modelo $f_1$ y $f_2$. Estos coeficientes, que anteriormente eran funciones únicamente de la relación de viscosidad $\lambda$ para gotas limpias, ahora dependen de $\lambda$, $Bqs$y$Bqd$.
• Validación Numérica: Para cuantificar la precisión y el rango de validez del modelo EMM, los autores realizan una comparación sistemática contra Simulaciones de Resolución Completa (FRS). Las FRS utilizan un método de lattice Boltzmann con frontera inmersa para resolver la dinámica acoplada de los fluidos interno y externo, así como la interfaz deformable, sin asumir deformaciones pequeñas.
• Análisis de Error: El estudio evalúa el error relativo medio $\langle \epsilon \rangle_T$ entre la deformación dependiente del tiempo $D(t)$ predicha por el modelo EMM y los resultados de las FRS a través de un amplio espacio de parámetros definido por el número capilar ($Ca$) y los dos números de Boussinesq.

Contribuciones Clave

1. Derivación de Coeficientes Generalizados: El artículo proporciona expresiones analíticas explícitas para los coeficientes $f_1^{(EMM)}$ y $f_2^{(EMM)}$ del EMM que dan cuenta tanto de las viscosidades de cizalla como de las de dilatación interfacial. Estas expresiones se reducen a los coeficientes originales de MM cuando las viscosidades interfaciales son cero.
2. Caracterización de la Dinámica Transitoria: El estudio revela que la viscosidad interfacial afecta el tiempo de deformación transitoria de manera no lineal. Mientras que una gota con viscosidad predominantemente de cizalla o predominantemente de dilatación se comporta de manera similar a una gota limpia respecto al tiempo de relajación, una gota con valores comparables de ambas viscosidades experimenta una ralentización dramática en la dinámica de deformación.
3. Mapeo del Dominio de Validez: A través de extensas pruebas numéricas, los autores delinean las regiones específicas en el espacio de parámetros $(Ca, Bqs, Bqd)$ donde el modelo EMM permanece cuantitativamente preciso.

Resultos

• Deformación en Estado Estacionario: El modelo EMM predice con precisión la deformación en estado estacionario $D_s$ para un amplio rango de parámetros. El modelo captura la tendencia contraintuitiva donde el aumento de la viscosidad de dilatación (con viscosidad de cizalla cero) aumenta ligeramente la deformación, mientras que el aumento de la viscosidad de cizalla (o tener viscosidades de cizalla y de dilatación iguales) reduce la deformación.
• Comportamiento Transitorio: El modelo reproduce con éxito la transición de la relajación monotónica a las respuestas oscilatorias observadas en las FRS, particularmente cuando $Bqs = Bqd$. El tiempo de deformación característico aumenta bruscamente solo cuando ambas viscosidades interfaciales son apreciables.
• Límites de Precisión:
  • Para viscosidades interfaciales pequeñas ($Bqs, Bqd \ll 1$) y números capilares pequeños ($Ca \le 0.2$), el modelo EMM reproduce los datos de las FRS con errores típicamente inferiores a unos pocos puntos porcentuales.
  • La precisión del modelo se degrada más rápidamente en configuraciones de viscosidad fuertemente asimétricas (donde un número de Boussinesq domina) en comparación con el caso equiviscoso ($Bqs = Bqd$).
  • A medida que $Ca$ aumenta (por ejemplo, $Ca = 0.5$), los acoplamientos no lineales en la tensión interfacial, que el límite de pequeña deformación del modelo EMM aproxima, conducen a mayores discrepancias. Con un $Ca$ alto y un $Bqs$ bajo, el modelo se acerca al umbral de ruptura, donde las FRS indican inestabilidad numérica debido a grandes deformaciones.

Significancia
El artículo establece al modelo EMM como una herramienta robusta y computacionalmente económica para predecir la deformación de gotas en presencia de viscosidad interfacial. Demuestra que el modelo mantiene una conexión clara con las teorías clásicas de pequeña deformación al tiempo que ofrece un control independiente sobre los efectos de cizalla y de dilatación. El trabajo proporciona un mapa cuantitativo de la validez del modelo, destacando que, si bien es altamente efectivo para deformaciones moderadas y relaciones de viscosidad específicas, las simulaciones de resolución completa siguen siendo indispensables para deformaciones grandes o condiciones interfaciales fuertemente asimétricas. Los autores sugieren que este marco podría servir como base para futuras extensiones a configuraciones de flujo más complejas e interfaces cargadas con surfactantes, aunque se señala que tales generalizaciones son direcciones futuras en lugar de logros actuales.