Sharp-interface Simulations of Energetic Multiphase Flows with Large Density and Viscosity Ratios

Este artículo propone un marco de flujos de Región de Donante de Momento Sincronizados (SynDRoM, por sus siglas en inglés) combinado con un limitador de viscosidad para mejorar la robustez numérica y la fidelidad física de las simulaciones de interfaz nítida para flujos multifásicos energéticos con grandes relaciones de densidad y viscosidad.

Lo esencial

La visión general: Simulando un mar tormentoso

Imagina que estás intentando crear una simulación por computadora de una violenta tormenta oceánica. Quieres ver cómo chocan las olas, cómo el aire es succionado hacia el agua y cómo se forman las burbujas. Esto es complicado porque el agua es pesada y densa, mientras que el aire es ligero y tenue. En términos de física, tienen una diferencia masiva de "densidad".

Cuando las computadoras intentan simular esto, a menudo fallan o producen resultados extraños e imposibles (como que el agua de repente se convierta en un fantasma o que el aire atraviese el agua como una bala). Este artículo presenta un nuevo conjunto de reglas (algoritmos) para hacer que estas simulaciones sean estables, precisas y físicamente realistas, incluso cuando las olas chocan violentamente.

El problema: El "Fantasma" y el "Choque"

Los autores explican que los métodos antiguos para simular estos flujos tienen dos fallas principales:

1. El problema del "Fantasma" (Penetración de velocidad):
   Imagina un camión pesado (agua) y una pluma (aire) moviéndose uno al lado del otro. En las simulaciones antiguas, el "viento" de la pluma a veces soplaba el camión hacia atrás, o el camión empujaba la pluma a través de su propio cuerpo. Esto se llama "penetración de velocidad". Esto crea formas falsas y no físicas en el agua, como un "cuerno de diablo" que sobresale de la ola.

2. El problema del "Choque" (Picos de cantidad de movimiento):
   Para solucionar el problema del fantasma, los científicos probaron un nuevo método llamado CMOM (Momento de Masa-Cantidad de Movimiento Consistente). Es como llevar un registro estricto de cuánto "impulso" (cantidad de movimiento) tiene cada gota de agua. Sin embargo, este método tiene un efecto secundario. Cuando un poco de agua pesada se mueve hacia una celda llena de aire, las matemáticas se confunden. Es como dividir un número enorme por un número diminuto, lo que resulta en un pico de velocidad masivo e imposible. Esto crea "bloques de velocidad": bolsas de aire falsas que se mueven a velocidades supersónicas que no deberían existir.

La solución: El método "SynDRoM"

Los autores proponen una nueva solución llamada SynDRoM (Región de Donante Sincronizada del flujo de Momento). Así es como funciona, usando una analogía:

La analogía: La cinta transportadora en movimiento
Imagina una cinta transportadora que lleva cajas.

• La forma antigua: Cuentas las cajas (masa) y el peso de las cajas (momento) por separado. Si una caja se mueve, podrías accidentalmente contar su peso en un lugar donde la caja aún no ha llegado. Esto causa el "choque" o el "pico" de velocidad.
• La forma SynDRoM: Este método actúa como un equipo sincronizado. Antes de mover el peso, observas exactamente qué parte de la cinta transportadora está trayendo ese peso.
  • Pregunta: "Si estoy moviendo este trozo específico de aire, ¿exactamente qué trozo de momento está unido a él?"
  • Asegura que el momento solo se mueva si la masa está realmente allí para transportarlo.
  • El resultado: No más picos de velocidad falsos. El aire permanece lento y el agua permanece pesada, tal como en la vida real. La simulación se mantiene fluida y no "explota".

El segundo problema: La viscosidad "resbaladiza"

El artículo también aborda un segundo problema: la Viscosidad (qué tan espeso o pegajoso es un fluido).

• El problema: El agua es pegajosa; el aire es resbaladizo. Cuando se mezclan en un límite agudo (como una ola rompiendo), la computadora intenta adivinar la "pegajosidad" en medio. Si adivina mal, las matemáticas se vuelcen inestables, como intentar equilibrar un lápiz sobre su punta.
• La solución: Los autores introducen un Limitador de Viscosidad.
  • La analogía: Imagina una señal de límite de velocidad. Incluso si las matemáticas intentan calcular una "pegajosidad" que haría que el fluido se moviera de forma imposiblemente rápida (inestable), el limitador dice: "No, no puedes ir más rápido que la velocidad del fluido más delgado aquí". Limita el cálculo para evitar que la simulación colapse, sin cambiar la física real del agua o el aire.

La prueba: ¿Funciona?

Los autores probaron sus nuevas reglas de tres maneras:

1. La rotura de la presa (Dam Break): Simularon una pared de agua colapsando.
   • Métodos antiguos: El agua se veía distorsionada con picos falsos.
   • SynDRoM: El agua chocó naturalmente y el aire no fue succionado hacia el agua de formas extrañas.
2. La inestabilidad de Kelvin-Helmholtz: Esto es cuando el viento sopla sobre el agua, creando olas que ruedan (como las nubes).
   • Resultado: La simulación mostró correctamente cómo las olas rodaban y crecían, sin que la computadora añadiera energía falsa o amortiguara las olas. Demostró que el método respeta las leyes de la física.
3. La ola que rompe: Simularon una ola masiva y diagonal rompiendo.
   • Resultado: La ola rompió, salpicó y creó espuma tal como un océano real. La energía total del sistema se mantuvo equilibrada (no desapareció mágicamente ni explotó). Incluso cuando añadieron "pegajosidad" (viscosidad), la simulación se mantuvo estable.

Conclusión

Este artículo presenta un nuevo "policía de tráfico" para las simulaciones por computadora de agua y aire.

• Evita que el aire atraviese el agua como un fantasma.
• Evita que el agua cree picos de velocidad imposibles.
• Mantiene los cálculos de "pegajosidad" sin que las matemáticas se rompan.

Al sincronizar exactamente qué se está moviendo con dónde se está moviendo, los autores han creado una herramienta de simulación mucho más robusta y confiable para estudiar eventos oceánicos violentos, como los que los ingenieros navales necesitan comprender para el diseño de barcos.

Resumen técnico

Resumen Técnico: Simulaciones de Interfaz Aguda de Flujos Multifásicos Energéticos con Grandes Relaciones de Densidad y Viscosidad

Planteamiento del Problema
La simulación de flujos multifásicos energéticos con relaciones de densidad extremas (por ejemplo, sistemas aire-agua con relaciones $\sim 10^3$ a $10^4$) sigue siendo un desafío significativo en la hidrodinámica naval. Las formulaciones tradicionales basadas en la velocidad suelen sufrir divergencia numérica debido al desalineamiento de los gradientes de presión y las discontinuidades de densidad, lo que provoca transferencia de momento espuria y explosiones de energía. Aunque los marcos de transporte de Momento-Masa Consistente (CMOM) mejoran la robustez al imponer la conservación del momento y la conservación de la energía semi-discreta, introducen nuevas dificultades cuando se combinan con métodos de interfaz aguda. Específicamente, el CMOM puede generar choques de momento severos, resultando en oscilaciones y sobreimpulsos de velocidad no físicos (violando la condición de Variación Total Limitada, o TVD) cuando las interfaces agudas cruzan los límites de la malla. Además, la estimación inadecuada de la viscosidad cerca de las interfaces, particularmente al dividir la viscosidad dinámica por la densidad en mallas escalonadas (staggered grids), puede introducir inestabilidades numéricas en pasos de tiempo finitos, socavando la robustez del solver.

Metodología
El artículo propone un conjunto de esquemas numéricos que preservan la física, diseñados para conciliar la fidelidad física con la robustez numérica para flujos con altas relaciones de densidad.

1. Flujo de Momento de Región Donante Sincronizado (SynDRoM):
   Para abordar las oscilaciones de velocidad inherentes al CMOM estándar, los autores introducen SynDRoM. Este método impone la monotonicidad del campo de velocidad transportado sincronizando el flujo de momento con el flujo de masa.

   • Mecanismo: A diferencia de los limitadores de pendiente existentes que se centran únicamente en la reconstrucción de la velocidad, SynDRoM reconstruye un par de distribución de densidad-velocidad dentro de una celda para satisfacer las restricciones de promedio de celda y los límites TVD. Posteriormente, evoluciona el momento utilizando un concepto de "región donante" determinado estrictamente por el flujo de masa. La velocidad advectada se calcula como un promedio ponderado sobre esta región donante específica, asegurando que la velocidad actualizada permanezca acotada.
   • Implementación: El método opera dentro de un marco de división direccional consistente con la advección de Volumen de Fracción (cVOF) conservativa. Utiliza un esquema de integración temporal de segundo orden de Runge-Kutta (RK2) ponderado por densidad para manejar la evolución no lineal de la velocidad respecto a la densidad y el momento.

2. Limitador de Viscosidad Cinemática Acotado:
   Para abordar las inestabilidades derivadas de los saltos de viscosidad en las interfaces, el artículo introduce un limitador de viscosidad.

   • Mecanismo: Los autores identifican que la interpolación estándar (aritmética o armónica) de la viscosidad dinámica en los vértices de la malla escalonada, seguida de la división por la densidad, puede producir una viscosidad cinemática efectiva no acotada, violando las restricciones de CFL viscoso. El limitador propuesto restringe la viscosidad dinámica efectiva utilizada en los términos de diferencia cruzada mediante la densidad mínima de las celdas de momento circundantes y un límite superior basado en la ocupación de fase para la viscosidad cinemática. Esto asegura que la viscosidad cinemática efectiva permanezca dentro de un rango que preserva la estabilidad.

Contribuciones Clave

• Algoritmo SynDRoM: Un método generalizado para CMOM que elimina las oscilaciones espurias de la velocidad sin sacrificar la conservación del momento o la nitidez de la interfaz; tiene en cuenta explícitamente la distribución del momento y la correspondiente región donante basada en el flujo de masa, evitando la necesidad de una reconstrucción completa de la interfaz en mallas escalonadas.
• Estabilización de la Viscosidad: Un limitador de viscosidad simple pero efectivo que evita la divergencia numérica en flujos con altas relaciones de viscosidad al acotar la viscosidad cinemática efectiva, preservando así la estabilidad sin comprometer la fidelidad física.
• Consistencia de División Direccional: Una formulación que asegura que las actualizaciones de momento sean consistentes con la advección de masa de división direccional, evitando flujos ficticios que podrían violar la acotación de la velocidad.

Resultados
Los métodos propuestos fueron validados mediante una serie de casos de prueba:

• Rotura de Presa 2D (2D Dam Break): La formulación SynDRoM resolvió con éxito la deformación de la interfaz y eliminó las estructuras de velocidad espurias (como las distorsiones de "cuernos de diablo" vistas en formulaciones de velocidad y los "bloques de velocidad" en el CMOM original) manteniendo la acotación de la energía.
• Advección de Gotícula 1D: En una prueba cuasi-unidimensional con una alta relación de densidad ($10^{-3}$), SynDRoM mantuvo la monotonicidad de la velocidad y las propiedades TVD, mientras que el CMOM original y las variantes de promedio de Favre produjeron sobreimpulsos no físicos.
• Inestabilidad de Kelvin-Helmholtz: El esquema capturó correctamente el crecimiento lineal de las perturbaciones y su subsiguiente evolución no lineal (emparejamiento de vórtices, enrollamiento) sin amortiguación artificial ni explosión de energía. El momento horizontal se conservó según la tolerancia del solver, y las regiones de alta velocidad se confinaron correctamente a la fase de fluido ligero.
• Flujo de Poiseuille Estratificado: El limitador de viscosidad previno con éxito la divergencia numérica observada en el promedio aritmético estándar y corrigió los gradientes de velocidad excesivos vistos en el promedio armónico, produciendo soluciones estables y físicamente precisas.
• Rompedor de Ola Estacionaria 3D: Las simulaciones de rompedores de olas estacionarias diagonales demostraron la capacidad del solver para manejar la extrema no linealidad, la formación de chorros y la entalpía de aire. El método mantuvo la energía mecánica acotada incluso en ausencia de viscosidad explícita, donde la disipación de energía ocurrió naturalmente a través del mecanismo iLES (Simulación de Grandes Edificios Implícita) a medida que el flujo se volvía irregular. Los estudios de convergencia de malla confirmaron tendencias consistentes en la evolución de la energía y las tasas de disipación.

Significancia
El artículo sostiene que las modificaciones propuestas a los esquemas de advección y viscosidad mejoran significamente la fidelidad y robustez de las simulaciones de flujos multifásicos, particularmente para flujos de superficie libre energéticos encontrados en la hidrodinámica naval. Al asegurar que los métodos CMOM satisfagan la condición TVD cerca de las celdas de interfaz vaciándose y al acotar la viscosidad cinemática, el enfoque elimina artefactos numéricos que podrían desencadenar eventos de ruptura falsos. Los autores enfatizan que los algoritmos SynDRoM y el limitador de viscosidad son ortogonales al método específico de captura de interfaz utilizado (demostrado aquí con VOF geométrico), lo que hace que el enfoque sea ampliamente aplicable. El trabajo proporciona una vía para simular flujos multifásicos complejos con altas relaciones de densidad con consistencia física, conservación del momento y acotación empírica de la energía, abordando desafíos de larga data en el campo.