A hybrid Volume of Fluid Phase-Field method for Direct Numerical Simulations of soluble surfactant-laden interfacial flows

Este trabajo presenta un método híbrido de Volumen de Fluido-Campo de Fase con refinamiento adaptativo de malla para simulaciones numéricas directas de flujos con surfactantes solubles, que captura con precisión el acoplamiento entre el transporte en el volumen y en la interfaz para demostrar cómo las tensiones de Marangoni alteran significativamente la dinámica de ascenso de burbujas en geometrías tridimensionales.

Lo esencial

El Panorama General: Una Nueva Forma de Simular Agua "Jabonosa"

Imagina que estás observando una burbuja subir a través de un vaso de agua. Si el agua es pura, la burbuja dispara hacia arriba en línea recta, tambaleándose un poco pero moviéndose rápido. Pero si añades jabón (tensioactivos) al agua, la burbuja se comporta de manera diferente. Podría ralentizarse, tambalearse más, o incluso cambiar su trayectoria.

Esto sucede porque las moléculas de jabón aman pegarse a la superficie de la burbuja. A medida que la burbuja se mueve, estas moléculas son empujadas, creando una "tensión" desigual en la piel de la burbuja. Esta tensión desigual actúa como manos invisibles que empujan la burbuja en diferentes direcciones, cambiando cómo se mueve.

El Problema:
Simular esto en una computadora es increíblemente difícil. Es como intentar filmar una burbuja de jabón con una cámara que tiene dos configuraciones conflictivas:

1. La Cámara Nítida: Necesita ver el borde de la burbuja como una línea tan fina como una navaja (para calcular la presión y la forma).
2. La Cámara Borrosa: Necesita ver las moléculas de jabón extendiéndose suavemente sobre ese borde (para calcular cómo se mueve el jabón).

La mayoría de los métodos informáticos te obligan a elegir una configuración de cámara, haciendo que la simulación sea físicamente inexacta o computacionalmente imposible de ejecutar para formas 3D complejas.

La Solución:
Los autores de este artículo construyeron un método híbrido. Piensa en ello como una simulación de "pantalla dividida" que utiliza lo mejor de ambos mundos simultáneamente:

• El Borde Nítido (Volumen de Fluido): Utilizan un método que mantiene el borde de la burbuja nítido y conserva perfectamente la cantidad de líquido (como un contorno de alta definición).
• El Jabón Suave (Campo de Fase): Utilizan una segunda capa "difusa" que actúa como una autopista suave para que las moléculas de jabón viajen. Esto permite que el jabón se mueva naturalmente entre el agua y la superficie de la burbuja sin quedarse atascado ni perderse.

Cómo Funciona: La Analogía del "Controlador de Tráfico"

Para hacer que esto funcione, los autores crearon un sistema de tráfico digital para las moléculas de jabón:

1. La Autopista (La Interfaz): La superficie de la burbuja es una autopista concurrida. Las moléculas de jabón son los coches.
2. Las Ramps de Entrada y Salida (Adsorción/Desorción):
   • Adsorción: Las moléculas de jabón del agua (el volumen) quieren saltar a la autopista (la superficie de la burbuja).
   • Desorción: Las moléculas de jabón se cansan y saltan de la autopista de vuelta al agua.
   • El nuevo método calcula exactamente cuántos coches saltan dentro o fuera en cada instante, asegurando que ningún coche desaparezca o aparezca de la nada.
3. El Embotellamiento (Efecto Marangoni): Cuando demasiados coches de jabón se amontonan en un solo punto de la burbuja, ese punto se vuelve "pegajoso" (alta tensión). La piel de la burbuja intenta alejarse de ese punto pegajoso, creando una fuerza que ralentiza la burbuja o la hace tambalearse. La simulación captura esta lucha de tirones perfectamente.

Lo Que Probaron (Los "Exámenes de Conducción")

Antes de dejar que su nuevo coche condujera por la autopista, lo hicieron pasar por una escuela de conducción con tres pruebas específicas:

1. La Prueba de Estiramiento (Esfera en Expansión): Inflaron una burbuja cubierta de jabón. Verificaron si el jabón se extendía uniformemente a medida que la burbuja se hacía más grande. La simulación coincidió perfectamente con las matemáticas.
2. La Prueba de Giro (Burbuja Giratoria): Hicieron girar una burbuja con jabón sobre ella. Verificaron si el jabón se movía correctamente alrededor del círculo sin filtrarse. Nuevamente, la simulación fue exacta.
3. La Prueba de Intercambio (Pared Plana): Observaron cómo el jabón se movía desde el agua hacia una pared plana y viceversa. Probaron tres escenarios:
   • Solo saltando hacia arriba: ¿Se adhiere el jabón? Sí.
   • Solo saltando hacia abajo: ¿Se va el jabón? Sí.
   • Ambos: ¿Encuentra un equilibrio? Sí.

El Evento Principal: La Burbuja que Sube

Finalmente, dejaron que su nuevo método simulara una burbuja subiendo en un tanque tridimensional de agua.

• La Burbuja Limpia: Subió relativamente rápido y en línea recta.
• La Burbuja de Jabón "Insoluble": El jabón estaba pegado a la superficie y no podía salir. Creó un fuerte "embotellamiento" en la parte trasera de la burbuja, ralentizándola significativamente.
• La Burbuja de Jabón "Soluble" (La Realidad): Aquí es donde brilla el nuevo método. El jabón podía saltar dentro y fuera de la burbuja.
  • Si el jabón saltaba fuera fácilmente (alta "desorción"), la burbuja se comportaba casi como una limpia.
  • Si el jabón saltaba dentro fácilmente (alta "adsorción"), la burbuja se comportaba como la versión de jabón atascado.
  • En el medio, la burbuja mostró un baile complejo: se ralentizó, cambió su trayectoria y dejó un "rastro" de jabón en el agua detrás de ella mientras subía.

Por Qué Esto Importa (Según el Artículo)

Los autores afirman que este método es robusto, escalable y preciso.

• Robusto: No se bloquea cuando la burbuja adquiere formas extrañas o se rompe.
• Escalable: Puede ejecutarse en superordenadores para manejar simulaciones 3D enormes y complejas de manera eficiente.
• Preciso: Predice correctamente a qué velocidad suben las burbujas y cómo se retuercen, coincidiendo con la física del mundo real.

En resumen: Construyeron un nuevo motor digital que finalmente puede simular cómo se comportan las burbujas de jabón en el espacio 3D, manejando la danza complicada entre la forma de la burbuja y las moléculas de jabón moviéndose dentro y fuera de su piel, todo sin perder precisión ni bloquear la computadora.

Resumen técnico

Resumen Técnico: Un Método Híbrido de Volumen de Fluido y Campo de Fase para Flujos Interfaciales con Surfactantes Solubles

Planteamiento del Problema
La simulación numérica precisa de flujos interfaciales que involucran surfactantes solubles presenta un desafío formidable debido a la naturaleza multiescala y multipísica del problema. A diferencia de los surfactantes insolubles, las especies solubles están gobernadas por un sistema acoplado que involucra convección-difusión en el volumen, transporte interfacial y procesos de intercambio cinético (adsorción/desorción). Estas interacciones generan campos de tensión superficial que varían espacial y temporalmente, dando lugar a tensiones de Marangoni que alteran significativamente la topología, la estabilidad y el movimiento del flujo. Los métodos numéricos existentes enfrentan limitaciones distintas:

• Métodos de Integral/Elemento de Frontera: Restringidos a regímenes de flujo de Stokes y tienen dificultades con cambios de topología (ruptura/coalescencia).
• Seguimiento de Frente (Front-Tracking): Alta precisión pero implementación compleja en 3D, requiriendo tratamiento topológico explícito y enfrentando problemas de equilibrio de carga en computación paralela.
• Métodos de Nivel de Conjunto (Level-Set): Manejan naturalmente los cambios de topología pero a menudo sufren problemas de conservación de masa que requieren reinitialización y correcciones posteriores.
• Métodos de Campo de Fase: Manejan naturalmente los cambios de topología pero a menudo dependen de suposiciones de equilibrio termodinámico que pueden ser violadas por convección fuerte, llevando a la no conservación de masa.
• Métodos de Volumen de Fluido (VoF): Excelentes para la conservación de masa e interfaces nítidas, pero tradicionalmente tienen dificultades para definir operadores diferenciales superficiales requeridos para el transporte de surfactantes en una interfaz de espesor cero.

Existe una necesidad específica de un método basado en VoF geométrico que soporte Refinamiento de Malla Adaptativo (AMR), maneje surfactantes solubles con intercambio dinámico volumen-interfaz y opere dentro de un marco escalable y de código abierto.

Metodología
Los autores presentan un marco numérico híbrido implementado en el código de código abierto Basilisk, que combina un método geométrico de Volumen de Fluido (VoF) con un enfoque difuso de Campo de Fase.

1. Representación de la Interfaz:

   • VoF ($c$): Un campo escalar que representa la fracción volumétrica de un fluido, utilizado para el seguimiento conservativo de la interfaz y la solución de momento. Garantiza la conservación exacta de masa y soporta la reconstrucción geométrica.
   • Campo de Fase ($\phi$): Una variable auxiliar que evoluciona mediante una formulación de Interfaz Difusa Conservativa Precisa basada en Allen-Cahn (ACDI). Sirve como un portador suave para el transporte de surfactantes, proporcionando vectores normales regularizados, curvatura y gradientes superficiales sin la necesidad de cálculos explícitos de derivadas superficiales en una interfaz nítida.
   • Acoplamiento: El Campo de Fase se reinitializa periódicamente desde el campo VoF utilizando una función de distancia con signo para mantener la consistencia.

2. Ecuaciones Gobernantes:

   • Movimiento del Fluido: Ecuaciones de Navier-Stokes incompresibles con densidad y viscosidad que varían espacialmente. Las fuerzas de tensión superficial se incorporan como fuerzas corporales volumétricas: una fuerza capilar isotrópica ($\sigma \kappa \mathbf{n} \delta_s$) y una fuerza tangencial de Marangoni ($\nabla_s \sigma \delta_m$).
   • Transporte de Surfactantes: Se resuelven dos campos escalares: concentración en el volumen ($F$) y concentración interfacial ($f$, una representación volumétrica de $\Gamma$).
     • La ecuación de transporte interfacial se regulariza para evitar derivadas superficiales, convirtiendo la concentración superficial $\Gamma$ en un campo volumétrico $f = \Gamma \delta_\phi$ localizado dentro de la interfaz difusa.
     • Cinética: La adsorción y desorción se modelan mediante términos fuente/sumidero regularizados ($J$) localizados en la interfaz, vinculando las concentraciones en el volumen y en la interfaz. El modelo soporta ecuaciones de estado generales (por ejemplo, isoterma de Langmuir) para relacionar la tensión superficial con la concentración.

3. Discretización Numérica:

   • Espacial: Discretización de volumen finito en mallas de cuadrícula/octree dinámicamente refinadas (AMR). La resolución se concentra cerca de la interfaz para resolver gradientes pronunciados en concentración y flujo.
   • Temporal: Se utiliza descomposición de operadores. La advección se trata explícitamente (esquema geométrico no dividido Bell–Collela–Glaz), mientras que la difusión, anti-difusión y términos fuente se tratan implícitamente utilizando un solver Poisson multigrilla modificado.
   • Estabilidad: Se impone una condición de Courant–Friedrichs–Lewy (CFL) basada en una velocidad efectiva que incluye tanto contribuciones convectivas como anti-difusivas para garantizar la estabilidad en presencia de gradientes interfaciales fuertes.

Contribuciones Clave

• Marco Híbrido: Desarrollo de un método que retiene las ventajas de conservación de masa e interfaz nítida del VoF mientras utiliza una variable de Campo de Fase únicamente como un portador suave para el transporte de surfactantes, evitando las complejidades geométricas de las derivadas superficiales en interfaces lineales por partes.
• Capacidad para Surfactantes Solubles: Extensión del trabajo previo VoF-Campo de Fase (limitado a surfactantes insolubles) al caso completamente soluble, incorporando términos acoplados de advección-difusión en el volumen e intercambio cinético.
• Escalabilidad y Adaptabilidad: El marco es totalmente adaptativo (AMR) y soporta configuraciones 2D, axisimétricas y 3D con alta escalabilidad paralela, permitiendo la resolución de gradientes multiescala tanto en campos volumétricos como interfaciales.
• Implementación de Código Abierto: El método está implementado en Basilisk, proporcionando una herramienta reproducible y accesible para la comunidad.

Resultados y Validación
El artículo valida el método mediante una serie de pruebas de referencia y casos de aplicación:

1. Validación sin Flujo:

   • Esfera en Expansión (3D): Verificó la recuperación de soluciones analíticas para la dilución de surfactantes insolubles durante la expansión isotrópica.
   • Interfaz Giratoria (2D): Confirmó convergencia de segundo orden para la advección-difusión de surfactantes insolubles en una interfaz en movimiento.
   • Adsorción/Desorción en Interfaz Plana: Validó el solver acoplado volumen-interfaz contra soluciones analíticas para adsorción de flujo constante, desorción pura y cinética completa de Langmuir. El método captura correctamente la transición al equilibrio y la formación de capas de agotamiento/richez en el volumen.

2. Burbuja Ascendente Axisimétrica:

   • Límite Insoluble: Demostró la reducción de la velocidad terminal de ascenso debido a las tensiones de Marangoni, con resultados convergiendo al límite de interfaz limpia a medida que disminuye la concentración de surfactante.
   • Dinámica Soluble: Varió sistemáticamente el número de Biot (desorción) y el número de Damköhler (adsorción).
     • Alto Biot (Desorción Rápida): El sistema se acerca al límite de interfaz limpia a medida que los surfactantes son removidos rápidamente de la interfaz.
     • Alto Damköhler (Adsorción Rápida): El sistema se acerca al límite de superficie incompresible (insoluble) a medida que la interfaz se satura rápidamente.
   • Estructura del Flujo: Las visualizaciones mostraron la formación de "tapas de estancamiento" de surfactante en la parte trasera de la burbuja, las cuales suprimen la recirculación interna y modifican la estructura de la estela.

3. Burbuja Ascendente Tridimensional:

   • Las simulaciones de una burbuja ascendiendo en 3D con surfactantes solubles capturaron fenómenos complejos que incluyen deformación de forma, desviaciones de trayectoria lateral (movimiento helicoidal) y la formación de estelas traseras debido a la desorción.
   • El método recuperó con éxito tanto los límites asintóticos limpios como insolubles en 3D, demostrando robustez al manejar patrones de flujo inestables y geometrías deformables.

Significado y Afirmaciones
El artículo afirma proporcionar un marco numérico robusto y versátil que cierra la brecha entre los modelos clásicos de surfactantes insolubles y formulaciones realistas solubles. Su significado radica en:

• Tratamiento Unificado: Maneja el acoplamiento no lineal, interfaces en evolución y configuraciones completamente tridimensionales dentro de una sola plataforma de código abierto.
• Fidelidad Física: Captura la compleja interacción entre hidrodinámica, transporte volumen/interfaz y tensiones de Marangoni sin depender de modelos de cierre de turbulencia o aproximaciones empíricas.
• Aplicabilidad Práctica: El método es adecuado para flujos de superficie libre que involucran reconexión de interfaz (por ejemplo, ruptura de chorros, olas rompientes) en presencia de surfactantes solubles relevantes industrial y biológicamente.

Los autores notan modestamente que el marco actual está limitado a números de Péclet moderados debido a la difusión numérica inherente en la discretización. Identifican la resolución precisa de gradientes interfaciales agudos en el régimen de alto Péclet (dominado por convección) como un desafío para el trabajo futuro, con el objetivo de desarrollar formulaciones que minimicen la difusión numérica.