Droplets sitting on thin elastic sheets: A study with the boundary element method

Resumen Técnico: Gotas sobre Láminas Elásticas Delgadas: Un Estudio con el Método de Elementos de Contorno

Planteamiento del Problema
El artículo investiga la interacción elasto-capilar entre una gota de líquido y una lámina elástica delgada. A diferencia de los sustratos rígidos, los sólidos flexibles se deforman bajo la influencia de las fuerzas capilares, lo que da lugar a formas de equilibrio complejas donde la gota se hunde en la lámina, formando una estructura similar a una lente. Si bien estudios previos han explorado patrones de arrugas, la durotaxis (movimiento hacia regiones más blandas) y formas anisotrópicas bajo tensión, se necesitaba un marco numérico unificado para simular las formas de equilibrio de las gotas sobre láminas elásticas bajo diversas condiciones de contorno y protocolos de tensión. Específicamente, los autores pretenden cuantificar cómo el espesor de la lámina, la tensión isotrópica y el estiramiento uniaxial influyen en la morfología de la gota, la distribución de tensiones y las propiedades ópticas.

Metodología
Los autores extienden su Método de Elementos de Contorno (BEM, por sus siglas en inglés) desarrollado previamente para simular la dinámica acoplada de una gota y una lámina elástica. El núcleo de la metodología consiste en:

1. Ecuaciones Dinámicas: El sistema se modela mediante una dinámica sobreamortiguada donde los campos de velocidad de la superficie de la gota ($\dot{R}_0$) y de la lámina elástica ($\dot{R}_{sh}$) están acoplados. La ecuación gobernante es un sistema lineal $G\dot{R} = -\nabla_R F$, donde $G$ es la matriz de fricción y $F$ es la energía libre total.
2. Formulación de la Energía Libre:
   • Gota ($F_0$): Incluye la tensión superficial líquido-gas y el trabajo realizado contra la interfaz sustrato-líquido, regido por el ángulo de contacto de Young. Se incluye un término de tensión de línea, pero se establece en cero al ser despreciable para las láminas elásticas consideradas.
   • Lámina ($F_{sh}$): Modelada utilizando el modelo de Skalak para el cizallamiento en el plano y la dilatación de área (hiperelasticidad) y el modelo de Helfrich para la flexión fuera del plano. Los módulos elásticos (cizalladura $\kappa_S$, área $\kappa_A$ y flexión $\kappa_b$) se derivan de las propiedades del material en volumen (módulo de Young $E$ y razón de Poisson $\nu$) y el espesor de la lámina $d$.
   • Restricciones ($F_{cons}$): Aplicadas mediante multiplicadores de Lagrange para mantener un volumen de gota constante y asegurar la conectividad geométrica entre las mallas de la gota y la lámina en la interfaz.
3. Matriz de Fricción: La matriz $G$ incorpora tensiones de cizalladura viscosa (derivadas del tensor de Oseen), fricción de la línea de contacto (resultado de Moffatt) y fricción líquido-sustrato gobernada por la condición de deslizamiento de Navier.
4. Implementación Numérica: El sistema se discretiza utilizando mallas triangulares generadas con gmresh. Los autores emplean un esquema de interpolación ponderada (usando Surrogates.jl) para mapear fuerzas y velocidades entre las distintas mallas de la gota y la lámina. Las simulaciones se adimensionalizan utilizando el radio de la línea de contacto $R_0$ y escalas de tiempo características derivadas de la ley de Cox-Voinov.

Resultos Clave
El estudio aplica el BEM extendido a tres protocolos específicos:

1. Lámina Elástica Sujeta (Clamped):

   • La forma de equilibrio de la gota es una lente que se vuelve cada vez más asimétrica con respecto al plano de la línea de contacto a medida que aumenta el espesor de la lámina.
   • Dependencia del Espesor: A medida que el espesor ($d$) aumenta, la parte inferior de la gota ($h_l$) disminuye mientras que la parte superior ($h_u$) aumenta, aunque la altura total permanece relativamente constante.
   • Ángulos de Contacto: El ángulo de contacto microscópico en la línea de contacto concuerda con el ángulo de Young, validando el modelo. Sin embargo, el ángulo de contacto macroscópico (determinado por la forma general de la lente) difiere significativamente del ángulo de Young y depende del espesor de la lámina.
   • Desacoplamiento: La forma de la mitad inferior de la gota está determinada principalmente por la elasticidad de la lámina y es casi independiente del ángulo de contacto de equilibrio ($\theta_{eq}$), mientras que $\theta_{eq}$ influye principalmente en la mitad superior.

2. Lámina Estirada Isotrópicamente:

   • La aplicación de tensión isotrópica eleva la gota, reduciendo la deformación de la lámina hasta que esta vuelve a ser plana en extensiones altas ($\epsilon \approx 0.1$).
   • Perfiles de Tensión: Los autores cuantifican perfiles radiales de tensiones elásticas azimutales ($T_{\phi\phi}$) y radiales ($T_{rr}$). Las tensiones aumentan con la extensión, y la línea de contacto se localiza cerca de un punto de inflexión en las curvas de tensión.
   • Sintonización Óptica: La gota con forma de lente actúa como una lente óptica líquida. La distancia focal, calculada mediante la ecuación del fabricante de lentes, puede sintonizarse reversiblemente en aproximadamente un 15% variando la tensión aplicada. Este rango es comparable al rango de sintonización de las lentes en los teléfonos inteligentes modernos.

3. Lámina Estirada Uniaxialmente:

   • El estiramiento de la lámina a lo largo de una dirección elonga la gota. El grado de elongación depende de las condiciones de contorno de los bordes laterales (sujetos, movimiento controlado hacia adentro o libres).
   • Morfología: Los bordes laterales libres resultan en la elongación más pronunciada y desplazamientos verticales significativos, creando pliegues ascendentes y hoyuelos en la lámina. La línea de contacto se vuelve no planar, elevada en la dirección lateral y deprimida en la dirección de estiramiento.
   • Anisotropía: Se cuantifica la relación de aspecto de la huella de la línea de contacto y la anisotropía del ángulo de contacto macroscópico. La elongación se atribuye a que la gota experimenta diferentes rigideces efectivas de la lámina a lo largo de las direcciones de estiramiento y lateral.

Significación y Reivindicaciones
El artículo afirma que el BEM extendido proporciona una herramienta numérica robusta para simular sistemas elasto-capilares, capturando con éxito la interacción entre la dinámica de fluidos y las deformaciones elásticas. Las principales contribuciones son:

• Metodológicas: Un BEM formalizado que acopla la dinámica de la superficie de la gota con las deformaciones de la lámina elástica, incluyendo energías de plano y de flexión.
• Perspectiva Física: Demostración de que el espesor de la lámina es un parámetro crítico que controla la asimetría de la gota y que el ángulo de contacto microscópico permanece consistente con la ley de Young a pesar de las deformaciones macroscópicas.
• Aplicación: El estudio demuestra explícitamente que la aplicación de tensión a una lámina elástica permite la sintonización dinámica y reversible de la distancia focal de las lentes líquidas. Esto ofrece un mecanismo potencial para controlar las propiedades ópticas en dispositivos micro-optofluídicos.
• Direcciones Futuras: Los autores sugieren que el método abre vías para estudiar patrones de arrugas, durotaxis en gradientes de tensión e interacciones elastocapilares entre múltiples gotas, señalando que la tensión uniaxial podría conducir a potenciales de interacción anisotrópicos y dinámicas complejas.