Droplet Spreading on Nanowall-Textured Hydrophobic Surfaces: A Nonlinear Reduced-Order Model

Resumen Técnico: Propagación de Gotas en Superficies Hidrofóbicas con Textura de Nanomuros

Planteamiento del Problema
El impacto de gotas en superficies no mojables implica una compleja conversión transitoria de la energía cinética traslacional en energía interfacial, movimiento interno y disipación. Si bien las leyes de escala clásicas (p. ej., $\beta_{max} \sim We^{1/4}$) resumen las tendencias de deformación máxima, no identifican de manera única cómo las texturas superficiales modifican los equilibrios de fuerza y energía subyacentes. Específicamente, las texturas superficiales introducen efectos contrapuestos: cambios en la geometría de mojado de equilibrio (estados Wenzel vs. Cassie–Baxter), barreras de energía que detienen la línea de contacto (anclaje o pinning) y la potencial penetración del líquido en la textura, lo que altera la disipación y la energía recuperable.

El entorno específico abordado es el impacto de gotas de agua sobre superficies de silicio lisas hidrofobizadas frente a superficies de silicio con textura de nanomuros, basadas en los datos experimentales de Yilbas et al. [11]. El desafío consiste en determinar el sistema radial de un solo grado de libertad (1-DOF) más simple capaz de reproducir la dependencia de la dirección y la condición de la propagación máxima sin recurrir a simulaciones multifásicas totalmente resueltas, que serían computacionalmente costosas (ya que requerirían resolver escalas de longitud que abarcan cuatro órdenes de magnitud).

Metodología
Los autores formulan un modelo de orden reducido no lineal de un solo grado de libertad (1-DOF) que actúa como un sustituto radial concentrado (lumped radial surrogate). El modelo no es una simulación de volumen de fluido o de campo de fase totalmente resuelta.

1. Dinámica Gobernante: El sistema rastrea una extensión radial efectiva $q(t)$, con un valor de referencia $q_0 = D_0/2$. La ecuación de movimiento es:
   $$m\ddot{q} = -k_{cap}(q - q_0)\left[1 + \kappa_s \left(\frac{q - q_0}{q_0}\right)^2\right] - c_{visc}\dot{q} - F_{pin}\tanh\left(\frac{\dot{q}}{v_s}\right)$$

   • Masa ($m$): Establecida igual a la masa total de la gota.
   • Rigidez Capilar ($k_{cap}$): Determinada por la escala de tiempo inercial-capilar, no mediante ajuste de parámetros.
   • Rigidez No Lineal: Un término cúbico ($\kappa_s$) proporciona un endurecimiento dependiente de la amplitud para evitar una complacencia poco realista ante grandes deformaciones.
   • Amortiguamiento Viscoso: Modelado con un coeficiente dependiente del número de Reynolds.
   • Resistencia Cinética: Un término regularizado se opone al movimiento, aproximándose a una magnitud $F_{pin}$ a altas velocidades.

2. Regla de Anclaje Híbrida: Dado que la fuerza cinética regularizada se anula a velocidad cero (lo que impide un anclaje estático), se impone un evento de movimiento-a-estancamiento explícito. La trayectoria termina cuando la velocidad cae por debajo de una banda de captura ($v_{stick}$) y la fuerza impulsora capilar es insuficiente para superar un umbral de capacidad estática ($\eta_s F_{pin}$).

3. Cierre de Anclaje: La magnitud de anclaje $F_{pin}$ se estima mediante un modelo de fuerza a escala de poro escalado por un recuento de encuentros de referencia fijo. Notablemente, bajo el supuesto de que el diámetro hidráulico es igual al espaciamiento de la textura ($D_H = 2a$), la magnitud de anclaje adimensional es independiente tanto del espaciamiento como del tamaño de la gota.

4. Calibración y Validación:

   • El coeficiente cúbico $\kappa_s$ se ajusta a un único punto de datos gráfico (superficie lisa, $V=20\,\mu$L, $h=20\,$mm).
   • La respuesta de la superficie texturizada se predice utilizando parámetros fijos, sin ajustes adicionales.
   • Análisis de Sensibilidad: Se realiza una sensibilidad determinista local sobre los parámetros de cierre ($\kappa_s, F_{pin}, \zeta_0, v_s$) y los ajustes numéricos para distinguir los controles físicos de los artefactos numéricos.

Contribuciones Clave

• Jerarquía de Modelos: El artículo establece una jerarquía de modelos (lineal sin anclaje, no lineal sin anclaje, no lineal con anclaje) para aislar los roles específicos de la rigidez no lineal y la resistencia de anclaje.
• Dinámica Híbrida Explícita: Introduce una regla matemática de movimiento-a-estancamiento para manejar el anclaje estático, reconociendo que la regularización suave por sí sola no puede mantener un estado anclado.
• Identificabilidad y Transparencia: El trabajo clasifica explícitamente los parámetros como valores de la literatura, entradas reconstruidas, coeficientes ajustados o cierres asumidos. Demuestra que el acuerdo en una condición no puede validar todos los cierres simultáneamente.
• Definición del Alcance: El modelo está rigurosamente limitado a la predicción del factor de propagación radial máxima ($\beta_{max}$). Excluye explícitamente las predicciones de rebote vertical, tiempos de contacto o transiciones de estado de mojado.

Resultados

• Caso Base ($V=20\,\mu$L, $h=20\,$mm):
  • El modelo lineal sin anclaje sobrepredice la propagación ($\beta_{max} = 3.955$).
  • El modelo no lineal sin anclaje (calibrado) coincide con el objetivo de la superficie lisa ($\beta_{max} = 2.253$).
  • El modelo no lineal texturizado predice una propagación reducida ($\beta_{max} = 1.810$), consistente con la tendencia cualitativa de reducción de propagación en superficies de nanomuros observada por Yilbas et al.
• Transferencia de Condición: Con los parámetros fijos, el modelo predice que la propagación máxima aumenta monótonamente con el número de Weber ($We$) reconstruido para ambas superficies. La superficie texturizada produce consistentemente un $\beta_{max}$ menor que la superficie lisa en los rangos de volumen y altura de liberación investigados.
• Sensibilidad:
  • La rigidez no lineal ($\kappa_s$) y la magnitud de anclaje ($F_{pin}$) se identifican como los controles principales de la propagación máxima.
  • La magnitud de anclaje es el control dominante específico de la superficie texturizada.
  • Los parámetros numéricos (tolerancia del solver, localización de eventos, velocidad de suavizado cinético) tienen efectos insignificantes en la propagación máxima predicha dentro de la precisión reportada.
  • Tiempo de Finalización: El tiempo en el que la trayectoria termina es altamente sensible al factor de capacidad estática ($\eta_s$) y a la magnitud de anclaje cerca de los umbrales de conmutación, lo que indica que el tiempo de finalización es un diagnóstico de clasificación de estado más que una predicción física robusta.

Significado y Reivindicaciones
El artículo afirma modestamente proporcionar un sustituto radial reproducible y de bajo costo computacional que define qué se puede y qué no se puede inferir a partir de la información disponible de la propagación máxima.

• Reproduce con éxito la tendencia cualitativa de que las texturas de nanomuros reducen la propagación máxima en comparación con las superficies lisas.
• Clarifica que la reducción en la propagación es impulsada por los efectos combinados del trabajo de anclaje hacia el exterior y la regla híbrida de capacidad estática.
• Los autores enfatizan que el modelo no predice el rebote vertical, los coeficientes de restitución normal o los tiempos de contacto físico. El "proxy de recuperación radial cero" para el caso texturizado es una convención de clasificación para la no recuperación radial, no una predicción de la altura de rebote experimental cero.
• El trabajo sirve como una línea base explícita para futuros modelos de dos coordenadas o multifásicos resueltos, resaltando las limitaciones de los enfoques de una sola coordenada para capturar transiciones complejas de estado de mojado y mecanismos de recuperación de energía.