cfdmfFTFoam: A front-tracking solver for multiphase flows on general unstructured grids in OpenFOAM

Este artículo presenta cfdmfFTFoam, un nuevo solver de código abierto integrado en OpenFOAM que implementa el Método de Rastreo de Frentes (FTM) en mallas no estructuradas generales, superando las limitaciones de software existente y facilitando futuras investigaciones en flujos multifásicos.

Lo esencial

¡Claro que sí! Imagina que quieres estudiar cómo se comportan dos líquidos que no se mezclan, como el aceite y el agua, o cómo una burbuja de aire sube por un vaso de agua. Para hacer esto en una computadora, los científicos necesitan "dibujar" la frontera exacta donde se tocan estos dos mundos.

Este paper presenta una nueva herramienta llamada cfdmfFTFoam. Aquí te explico de qué se trata usando analogías sencillas:

1. El Problema: Dibujar la frontera en un laberinto

Imagina que tienes un líquido (como el agua) y quieres simular cómo se mueve una gota de aceite dentro de él.

• La vieja forma (VOF): Es como intentar dibujar la forma de la gota usando solo los bloques de un rompecabezas cuadrado. Si la gota es redonda, los bloques cuadrados no encajan bien y la gota parece pixelada o "escalonada". Además, con el tiempo, la gota puede perder su forma o romperse de forma falsa porque el dibujo es impreciso.
• La nueva forma (FTM - Front Tracking): En lugar de usar bloques, imagina que la gota es una burbuja de jabón real hecha de una malla de triángulos flexibles que flota dentro del agua. Esta burbuja se mueve, se estira y se deforma, pero siempre mantiene su forma perfecta y nítida.

2. La Innovación: Un puente entre dos mundos

El problema es que la burbuja (la malla triangular) se mueve libremente, pero el agua (el fluido) se calcula en una red fija de celdas (como una rejilla de pesca).

• El desafío: ¿Cómo le dices a la red de pesca qué está pasando en la burbuja de jabón que está flotando encima?
• La solución del paper: Los autores crearon un "traductor" o un puente. Han tomado un código antiguo y muy bueno (llamado Ftc3D) y lo han integrado en un programa muy popular de ingeniería llamado OpenFOAM.
• La analogía: Piensa en OpenFOAM como un GPS muy avanzado que puede manejar cualquier tipo de terreno (caminos rectos, curvas, montañas, o lo que sea, es decir, "mallas no estructuradas"). Antes, el código de la burbuja de jabón solo funcionaba en terrenos cuadrados y simples. Ahora, gracias a este nuevo programa, la burbuja puede viajar por cualquier tipo de terreno complejo sin perderse.

3. ¿Qué hace exactamente este nuevo programa?

El programa hace cuatro cosas principales para mantener la simulación realista:

1. El "Gimnasio" de la burbuja (Corrección de volumen): A veces, la burbuja de jabón digital se encoge o se expande un poco por errores de cálculo. El programa tiene un algoritmo que actúa como un entrenador personal: si la burbuja pierde volumen, le da un pequeño "empujón" a sus esquinas para que vuelva a tener el tamaño correcto. ¡Nunca deja que la burbuja desaparezca!
2. El "Peluquero" (Remallado): Cuando la burbuja se estira mucho, sus triángulos se vuelven muy largos y feos (como un pelo enredado). El programa corta esos triángulos largos y crea otros más pequeños y ordenados, o une los muy pequeños. Así, la malla siempre se ve bien y no se rompe.
3. El "Abogado" (Tensión superficial): Calcula con mucha precisión la fuerza que intenta encoger la burbuja (como la piel de la burbuja de jabón). Lo hace de una manera directa y precisa, evitando errores que otros programas cometen.
4. El "Mensajero" (Comunicación): Usa un método matemático sofisticado (RKPM) para pasar la información de la burbuja a la red de agua y viceversa, asegurando que las fuerzas se transmitan suavemente sin causar saltos bruscos en la simulación.

4. ¿Por qué es importante?

Antes, si querías simular algo complejo (como una burbuja subiendo por un tubo curvo o gotas chocando en un terreno irregular), tenías que usar programas muy limitados o métodos menos precisos.

• La ventaja: Este nuevo programa permite simular estos fenómenos con mucha más precisión y en cualquier tipo de geometría (tubos curvos, máquinas complejas, etc.).
• El resultado: Es como pasar de ver un video de baja resolución (pixelado) a ver una película en 4K. Las gotas no se rompen falsamente, las burbujas mantienen su forma y los científicos pueden confiar más en los resultados.

En resumen

cfdmfFTFoam es un nuevo "superpoder" para los ingenieros y científicos. Les permite simular cómo interactúan líquidos y gases con una precisión increíble, incluso en escenarios muy complicados, manteniendo la "piel" de las gotas y burbujas siempre nítida y real, sin importar por dónde viajen. Es una herramienta de código abierto, lo que significa que cualquiera puede usarla y mejorarla para descubrir cosas nuevas en el mundo de los fluidos.

Resumen técnico

Resumen Técnico: cfdmfFTFoam

1. Planteamiento del Problema
La simulación de flujos multifásicos con fenómenos complejos (como deformación de burbujas, colisión de gotas, ruptura de sprays y ebullición) requiere métodos de simulación totalmente resueltos (FRS). Aunque existen métodos de captura de interfaz (como VOF o Level-Set) ampliamente disponibles en software de código abierto como OpenFOAM, los métodos de seguimiento de interfaz (Front-Tracking Method - FTM) puros son escasos en el ámbito de código abierto.

• Limitaciones actuales: La mayoría de los códigos FTM existentes (como PARIS, FronTier++, TrioCFD) están limitados a grids cartesianos estructurados.
• El vacío: No existía una herramienta de código abierto que implementara un FTM puro en grids no estructurados generales, lo cual es crucial para geometrías complejas y simulaciones paralelas eficientes. Además, las soluciones híbridas existentes (como lentFoam) no son FTM puros, ya que carecen de una representación lagrangiana de la superficie y calculan fuerzas basándose en campos eulerianos.

2. Metodología
El trabajo presenta cfdmfFTFoam, un nuevo solver integrado en el framework de OpenFOAM (versión 9) que combina el código FTM clásico Ftc3D (desarrollado por Tryggvason y Dabiri) con la infraestructura de grids no estructurados de OpenFOAM.

• Enfoque Híbrido Lagrangiano-Euleriano:
  • Euleriano: Se resuelven las ecuaciones de Navier-Stokes en un grid volumétrico no estructurado (hexaédrico, tetraédrico, polihedral, etc.).
  • Lagrangiano: La interfaz entre fases se representa mediante una malla superficial triangular (el "frente") que se mueve independientemente del grid subyacente.
• Algoritmos Clave Implementados:
  1. Comunicación Front-Campo (RKPM): Para transferir datos (como la tensión superficial) desde la malla lagrangiana a las celdas eulerianas, se implementó el Método de Partículas con Núcleo Reproductor (RKPM). Esto permite una distribución precisa de fuerzas en grids no estructurados, conservando momentos de orden cero y uno.
  2. Corrección de Volumen (VC2): Dado que el FTM no conserva masa/volumen por defecto, se implementó un algoritmo que corrige la posición de los vértices del frente resolviendo una ecuación cúbica para minimizar el error de volumen.
  3. Remallado (Remeshing): Se incluyen algoritmos de refinamiento, coarsening (agrupamiento) y suavizado (TSUR3D, VCS) para mantener la calidad de la malla triangular y controlar la relación de aspecto de los elementos.
  4. Cálculo de Tensión Superficial: Se utiliza un método directo basado en elementos para calcular las fuerzas interfaciales con alta precisión.
  5. Construcción de Función Indicadora: Se ofrecen dos métodos: resolución de la Ecuación de Poisson (PE) y Transformación del Punto Más Cercano (CPT), siendo PE preferido por su suavidad en grids no estructurados.
  6. Advección del Frente: Se aprovechan las herramientas de seguimiento de partículas de OpenFOAM para mover los vértices del frente según la velocidad del fluido interpolada.

3. Contribuciones Clave

• Primera implementación FTM pura en grids no estructurados: Llena un vacío significativo en el software CFD de código abierto, permitiendo simulaciones en geometrías complejas.
• Integración robusta: Logra la comunicación entre la malla lagrangiana y el grid euleriano no estructurado en arquitecturas paralelas (descomposición de dominio), un desafío técnico mayor debido a la topología variable.
• Biblioteca modular: El paquete incluye un solver principal (cfdmfFTFoam), una biblioteca (frontTracking) y una serie de tutoriales, facilitando su adopción por la comunidad.
• Validación exhaustiva: Se demuestra la capacidad del solver para manejar grandes deformaciones, altos números de Reynolds y altas relaciones de densidad/viscosidad.

4. Resultados y Validación
El solver fue validado mediante varios benchmarks estándar:

• Deformación 3D (Leveque): El solver demostró una precisión superior a los métodos VOF (MULES e isoAdvector) en la preservación de la interfaz, evitando rupturas no físicas y logrando un error de advección del 1% (vs. 36-41% en VOF). También mostró una convergencia de segundo orden y una excelente conservación de volumen con el algoritmo VC2.
• Gota Estancada: Se evaluó la precisión de la tensión superficial. El solver logró un error de presión de 0.046% y corrientes parasitas muy bajas ($2.1 \times 10^{-3}$), superando significativamente al solver estándar interFoam (error de 12.6%).
• Gota en Flujo de Poiseuille: Se validó en múltiples topologías de celdas (hexaédricas, prismas, poliedros, tetraedros), mostrando concordancia con soluciones analíticas. La escalabilidad paralela fue comparable a la del solver estándar interFoam.
• Ascenso de Gota Taylor: Se reprodujeron resultados experimentales de una burbuja Taylor en un tubo, demostrando la capacidad de manejar la formación de cavidades y la dinámica compleja de la estela.
• Burbuja Libre Ascendente (Alta Inercia): En un caso con altos números de Eötvös y Morton, el solver predijo correctamente la velocidad terminal y la forma general, aunque subestimó ligeramente la deformación axial debido a inconsistencias en la discretización de momento en OpenFOAM para relaciones de densidad extremas ($10^4$).

5. Significado e Impacto
El desarrollo de cfdmfFTFoam es un avance significativo para la investigación en dinámica de fluidos computacional (CFD) multifásica:

• Accesibilidad: Pone a disposición de la comunidad científica un solver FTM de alta fidelidad, gratuito y de código abierto, anteriormente limitado a grupos de investigación privados o códigos comerciales.
• Flexibilidad Geométrica: Permite aplicar la alta precisión del FTM (sin disipación de interfaz, cálculo exacto de curvatura) a problemas industriales con geometrías complejas que requieren grids no estructurados.
• Base para Futuras Mejoras: El código sirve como plataforma para desarrollar futuras estrategias híbridas, dinámicas de líneas de contacto, transporte de calor/masa interfacial y algoritmos de cambio topológico.

En conclusión, el artículo presenta una solución robusta y validada que supera las limitaciones de los métodos de captura de interfaz tradicionales en términos de precisión geométrica, al mismo tiempo que supera las limitaciones de los métodos FTM existentes en cuanto a la flexibilidad de la malla.