Color-gradient lattice Boltzmann modeling of wetting boundary condition on curved solid boundaries

Este artículo introduce una condición de contorno de mojabilidad para superficies sólidas curvas en el método de lattice Boltzmann de gradiente de color mediante la actualización de los parámetros de orden en nodos fantasma, un esquema validado en hardware GPU para manejar eficazmente grandes contrastes de densidad y viscosidad mientras minimiza las corrientes espurias y reproduce con precisión tanto los comportamientos de la línea de contacto estática como los dinámicos.

Lo esencial

Imagina que estás intentando simular cómo se comporta una gota de agua cuando golpea una superficie curva, como una gota de lluvia aterrizando en una hoja o una burbuja deslizándose por un vidrio curvo. Para hacer esto en una computadora, los científicos utilizan un método llamado Método de Lattice Boltzmann. Piensa en este método como una cuadrícula gigante e invisible de pequeñas teselas que cubren la pantalla de la computadora. Cada tesela contiene un poco de información de "fluido", y la computadora actualiza estas teselas paso a paso para ver cómo se mueve el fluido.

La parte difícil es la condición de contorno —específicamente, cómo se comporta el fluido cuando toca una pared sólida—. En el mundo real, el agua no se detiene en seco ante una pared; forma un ángulo específico (llamado ángulo de contacto) dependiendo de si la superficie es húmeda (como un vidrio limpio) o seca (como un coche encerado).

El Problema: El "Fantasma" en la Máquina

En la simulación por computadora, la pared sólida no es una línea suave; es dentada porque está hecha de teselas de una cuadrícula cuadrada. Para que las matemáticas funcionen, la computadora necesita saber qué está haciendo el fluido dentro de la pared sólida, a pesar de que no hay fluido allí. Estos puntos imaginarios dentro de la pared se llaman "nodos fantasma" (ghost nodes).

Los métodos anteriores para decirle a estos nodos fantasma qué hacer tenían algunos defectos:

• A veces creaban "corrientes fantasma" (velocidades espurias) donde el fluido parecía moverse por sí solo sin ninguna fuerza.
• Tenían dificultades con las superficies curvas, actuando a menudo como si solo estuvieran diseñados para paredes planas.
• A veces requerían matemáticas especiales y complicadas solo para manejar un ángulo neutro (donde el agua ni se extiende ni se encoge).

La Solución: Una Nueva Regla para los Fantasmas

Los autores de este artículo introdujeron una regla nueva y más simple para estos nodos fantasma.

La Analogía: Imagina que el fluido tiene un "estado de ánimo" (representado por un color, de 0 para gas a 1 para líquido). En el mundo real, este estado de ánimo cambia suavemente de gas a líquido a medida que cruzas la superficie.

• Método Antiguo: Era como intentar adivinar el estado de ánimo de una persona que está detrás de una pared gritando una suposición al azar.
• Nuevo Método: Los autores se dieron cuenta de que, si conoces el "estado de ánimo" de la persona que está justo fuera de la pared (en el fluido), puedes extender matemáticamente esa curva suave del estado de ánimo a través de la pared hacia el nodo fantasma. Simplemente preguntan: "Si el fluido quiere formar un ángulo de 45 grados aquí, ¿cuál debe ser el estado de ánimo del nodo fantasma para que eso suceda?".

Esta nueva regla es como un puente sin costuras. Extiende la forma natural de la gota de fluido justo hasta y ligeramente dentro de la pared sólida, asegurando que el ángulo que la gota hace con la pared sea exactamente el que el científico solicitó.

Lo que Probaron

Para demostrar que su nueva regla funciona, realizaron varias simulaciones en un chip de computadora muy potente (un NVIDIA A100 GPU):

1. La Gota Estática: Colocaron una gota de agua sobre una placa plana y un cilindro curvo. Verificaron si la gota se asentaba en el ángulo exacto que pidieron.
   • Resultado: Su nueva regla fue más precisa que el mejor método anterior, especialmente cuando el ángulo era muy agudo (como una gota que se encoge) o muy plano (como una gota que se extiende).
2. La Partícula Flotante: Simularon un cilindro flotando en el límite entre aceite y agua.
   • Resultado: Su método calculó la posición de la línea de agua con mayor precisión que antes.
3. La Gota en Caída: Simularon una gota cayendo y golpeando un cilindro, observando cómo salpicaba y se extendía.
   • Resultado: La gota se comportó de manera realista, y la nueva regla no causó movimientos extraños o falsos en el fluido.

Conclusiones Clave

• Precisión: El nuevo método maneja mucho mejor las superficies curvas, manteniendo el ángulo del fluido correcto ya sea que la pared sea plana o redonda.
• Estabilidad: Crea muy poco "ruido falso" (corrientes espurias) en la simulación, lo que significa que el fluido se ve más natural.
• Simplicidad: Evita la necesidad de matemáticas especiales y complicadas cuando el ángulo de contacto es exactamente de 90 grados (neutro), lo cual era un dolor de cabeza para los métodos anteriores.
• Velocidad: Al utilizar chips de computadora modernos (GPUs) y un estilo de programación específico, hicieron que las simulaciones corrieran muy rápido. Descubrieron que usar un formato numérico ligeramente menos preciso (precisión simple) hacía que la computadora funcionara el doble de rápido sin arruinar los resultados para la mayoría de las pruebas.

En resumen, los autores construyeron un mejor "libro de reglas" para cómo las simulaciones por computadora manejan el borde donde el líquido se encuentra con una pared sólida, haciendo que las gotas digitales se vean y actúen más como las reales, incluso en superficies curvas.

Resumen técnico

Resumen Técnico: Modelado de Lattice Boltzmann de Gradiente de Color para Condiciones de Contorno de Humectación en Fronteras Sólidas Curvas

Planteamiento del Problema
El modelado preciso de los fenómenos de humectación es crítico para simular procesos tecnológicamente relevantes como el flujo en medios porosos, la microfluídica y el secado de películas coloidales. En los marcos de la metodología de Lattice Boltzmann (LBM) que utilizan interfaces difusas, la implementación de condiciones de contorno de humectación en fronteras sólidas curvas presenta desafíos significativos. Los enfoques existentes suelen sufrir limitaciones que incluyen la introducción de fuentes o sumideros artificiales de masa, la generación de velocidades espurias (parásitas) cerca de las interfaces, restricciones a geometrías planas e inconsistencias al imponer simultáneamente las condiciones de no deslizamiento (no-slip) y de humectación. Además, muchos métodos requieren un tratamiento especial para ángulos de contacto neutros ($\theta = 90^\circ$) o dependen de extrapolaciones lineales que introducen errores cerca de las fronteras sólidas.

Metodología
Los autores presentan una condición de contorno de humectación diseñada específicamente para fronteras sólidas curvas dentro de un marco de LBM de gradiente de color en dos dimensiones. El núcleo de la metodología consiste en actualizar el parámetro de orden (campo de color, $\phi$) en los "nodos fantasma" (ghost nodes) ubicados dentro de la fase sólida.

• Regla de Actualización de Nodos Fantasma: En lugar de prescribir un valor constante o utilizar la extrapolación lineal, el esquema extiende el perfil de equilibrio del color hacia la fase sólida. El valor del campo de color en un nodo fantasma ($x_G$) se calcula basándose en el campo de color en un punto correspondiente en el dominio de fluido ($x_F$) ubicado a una distancia $\Delta x$ a lo largo de la normal sólida.
• Formulación Matemática: Asumiendo que la distribución de color está cerca del equilibrio, la derivada del parámetro de orden en la dirección de la normal sólida se deriva del perfil de equilibrio y del ángulo de contacto $\theta$ prescrito. La integración de esta relación produce una regla de actualización no lineal:
  $$\phi_G = \frac{\phi_F}{\phi_F + (1 - \phi_F) \exp(\varepsilon)}$$
  donde $\varepsilon = -4\Delta x \cos \theta / W$, y $W$ es el ancho de la interfaz.
• Características Clave: Esta formulación asegura que el campo de color del nodo fantasma permanezca acotado entre 0 y 1, evitando la transferencia de masa artificial. Crucialmente, la regla no requiere un tratamiento especial cuando $\theta = 90^\circ$ (donde $\varepsilon = 0$), a diferencia de otros esquemas. El método se implementa utilizando un marco basado en JAX que se ejecuta en GPUs NVIDIA A100, utilizando fusión de kernels (kernel fusion) y diferenciación automática para un alto rendimiento.

Contribuciones Clave y Resultados
El esquema propuesto fue validado contra soluciones analíticas y el esquema de referencia de Fakhari et al. [20] a través de varios casos de prueba:

1. Gota Estática sobre una Placa Plana: El esquema produjo consistentemente errores menores en los ángulos de contacto medidos en comparación con Fakhari et al. en un rango de $30^\circ$ a $150^\circ$. También generó corrientes espurias más pequeñas ($5.74 \times 10^{-7}$ a $7 \times 10^{-7}$ unidades de red) en comparación con el esquema de referencia.
2. Gota sobre un Cilindro Sólido: Para una gota situada sobre un cilindro curvo, el método propuesto demostró una mejora en la precisión al predecir la altura del ápice de la gota ($h_{max}$), particularmente en ángulos de contacto extremos ($30^\circ$ y $150^\circ$), donde los errores se redujeron de $\sim 0.76\%$ a $\sim 0.35\%$ en comparación con el esquema de referencia.
3. Partícula Fija en la Interfaz Líquido-Gas: En simulaciones de un arreglo periódico de cilindros atrapados en una interfaz, el esquema propuesto redujo los errores en el desplazamiento vertical de la interfaz a $\sim 9\%$ en ángulos extremos, frente al $\sim 14\%$ del esquema de referencia.
4. Impacto de Gota Dinámica: El esquema simuló con éxito el impacto de una gota sobre un cilindro en números de Reynolds y Bond altos, capturando comportos complejos como el esparcimiento de la gota en superficies hidrofílicas y la división en superficies hidrofóbicas. Los resultados se alinearon bien con las leyes de escala analíticas para el espesor de la película.
5. Rendimiento Computacional: La implementación en GPU mostró que los cálculos de precisión simple (FP32) fueron al menos dos veces más rápidos que los de precisión doble (FP64), manteniendo la precisión para casos estáticos y dinámicos moderados. Sin embargo, el FP32 se volvió inestable para la prueba de impacto de gota con números de Reynolds/Weber elevados, necesitando FP64 para estos escenarios de alta energía.

Significancia y Reivindicaciones
El artículo sostiene que la condición de contorno de humectación propuesta ofrece una alternativa robusta y precisa para simular flujos bifásicos inmiscibles en fronteras sólidas curvas. Al extender el perfil de equilibrio hacia la fase sólida mediante una regla específica de actualización de nodos fantasma, el método:

• Mantiene la capacidad del modelo para manejar grandes contrastes de densidad y viscosidad.
• Produce corrientes espurias relativamente pequeñas (órdenes de magnitud menores que los modelos de pseudopotencial populares).
• Evita las singularidades numéricas asociadas con los ángulos de contacto neutros presentes en otros métodos.
• Demuestra una precisión superior sobre los esquemas existentes (específicamente Fakhari et al.) tanto para líneas de contacto estáticas como dinámicas en geometrías curvas.

Los autores señalan que, si bien la implementación actual está restringida a la ejecución en una sola GPU y en dos dimensiones, la física subyacente es directamente aplicable a los modelos de campo de fase (phase-field) y de energía libre, y puede extenderse a tres dimensiones. El trabajo proporciona una base para futuras simulaciones que involucren partículas sólidas en interfaces líquido-gas, particularmente cuando se acoplen con avances en los cálculos de la fuerza capilar de la línea de contacto y esquemas de histéresis.