A HHO formulation for variable density incompressible flows where the density is purely advected

Este artículo presenta una formulación de métodos híbridos de alto orden (HHO) para flujos incompresibles de densidad variable que garantiza la conservación exacta del volumen y la puridad de la advección de la densidad, ofreciendo ventajas como robustez frente a la presión, conservación de límites de densidad y eficiencia computacional, las cuales se validan mediante pruebas de convergencia y simulaciones de la inestabilidad de Rayleigh-Taylor.

Lo esencial

Imagina que estás en una cocina, preparando una ensalada de dos ingredientes que no se mezclan: aceite y vinagre. Si agitas la botella, ves cómo el aceite (más ligero) sube y el vinagre (más pesado) baja, creando remolinos y formas extrañas. Este es un ejemplo de lo que los físicos llaman "flujo de fluidos de densidad variable".

Ahora, imagina que quieres predecir exactamente cómo se moverá esa ensalada en una computadora. El problema es que las computadoras no "ven" el fluido como un todo suave; lo dividen en millones de trocitos pequeños (como un mosaico) y tratan de adivinar qué hace cada trocito. Si los cálculos no son perfectos, la computadora puede cometer errores: puede hacer que el aceite desaparezca, que el vinagre se convierta en agua, o que la presión se vuelva loca, arruinando toda la simulación.

¿Qué propone este artículo?

Los autores, Lorenzo Botti y Francesco Carlo Massa, han creado una nueva "receta" matemática (un algoritmo) llamada HHO (Hybrid High-Order) combinada con un reloj muy preciso llamado ESDIRK. Su objetivo es simular estos fluidos mezclados (como el aceite y el vinagre, o incluso magma y roca) de una manera increíblemente precisa y eficiente.

Aquí tienes los puntos clave explicados con analogías sencillas:

1. El Mosaico Inteligente (Discretización Híbrida)

La mayoría de los métodos antiguos miran solo el centro de cada trocito del mosaico. Este nuevo método es como un detective que mira tanto el centro de la habitación como las paredes.

• La analogía: Imagina que quieres saber cómo se mueve el agua en una piscina. Los métodos viejos solo miran el agua en el medio. Este nuevo método también vigila las orillas. Al hacerlo, puede reconstruir el movimiento del agua de forma mucho más realista, asegurándose de que el agua no "se filtre" mágicamente a través de las paredes (conservación exacta del volumen).

2. El Guardiana de la Presión (Robustez a la presión)

En la física de fluidos, la presión actúa como un árbitro que mantiene el equilibrio. A veces, el árbitro se equivoca y el juego se arruina.

• La analogía: En este nuevo método, el "árbitro" (la presión) hace su trabajo, pero si comete un pequeño error, no afecta a los jugadores (la velocidad y la densidad). Es como si los jugadores fueran tan buenos que, aunque el árbitro silbe mal, ellos siguen jugando perfectamente. Esto es lo que llaman "robustez a la presión".

3. El Vinagre que no se desvanece (Advección pura)

Cuando mezclas dos fluidos, la densidad (qué tan pesado es cada trozo) debe viajar con el fluido, como un pasajero en un autobús.

• La analogía: En simulaciones antiguas, a veces el "pasajero" (la densidad) se caía del autobús o se desvanecía en el aire. Este nuevo método asegura que la densidad sea transportada exactamente donde va el fluido, sin perder ni una gota. Es como si el autobús tuviera un sistema de seguridad que impide que nadie se caiga.

4. El Reloj de Alta Precisión (ESDIRK)

Para simular el movimiento en el tiempo, necesitan un reloj.

• La analogía: Imagina que tomas fotos de una carrera de coches. Si tomas una foto cada segundo, te pierdes los detalles. Si tomas una foto cada milisegundo, es mucho mejor, pero cuesta mucho trabajo. Este método usa un "reloj" (ESDIRK) que es capaz de tomar fotos muy precisas en intervalos de tiempo complejos, permitiendo ver detalles rápidos sin que la computadora se vuelva loca.

5. El Truco del Espacio (Condensación Estática)

Hacer estos cálculos requiere mucha memoria, como intentar guardar una biblioteca entera en un teléfono móvil.

• La analogía: Los autores usan un truco de magia llamado "condensación estática". Imagina que tienes que resolver un rompecabezas gigante. En lugar de guardar todas las piezas sueltas, agrupan las piezas del centro de cada cuadro y solo guardan las piezas de los bordes. Esto reduce el tamaño del rompecabezas en la memoria de la computadora, haciéndolo mucho más rápido y ligero, sin perder precisión.

6. La Prueba de Fuego: La Inestabilidad de Rayleigh-Taylor

Para probar su receta, simularon un fenómeno famoso: cuando un fluido pesado cae sobre uno ligero (como aceite sobre agua, pero al revés), creando remolinos complejos.

• El resultado:
  • Con baja diferencia de densidad (Aceite y agua), su método funcionó perfecto, incluso usando polinomios de alto grado (como usar lentes de aumento muy potentes) en mallas gruesas, ahorrando tiempo de cálculo.
  • Con alta diferencia de densidad (Aceite y mercurio, por ejemplo), el método se mantuvo estable y no permitió que la densidad se volviera negativa (algo imposible en la realidad), lo cual es crucial para evitar que la simulación explote.

En resumen

Este artículo presenta una nueva herramienta matemática para simular fluidos mezclados que es:

1. Precisa: No pierde ni una gota de fluido ni de densidad.
2. Robusta: No se rompe aunque la presión sea difícil de calcular.
3. Eficiente: Usa menos memoria y tiempo de computadora gracias a trucos matemáticos inteligentes.

Es como pasar de usar un mapa de papel arrugado para navegar un océano, a usar un GPS de última generación que no solo te dice dónde estás, sino que predice las olas y te asegura que nunca te perderás, incluso en la tormenta más fuerte.

Resumen técnico

Aquí presento un resumen técnico detallado del artículo en español, estructurado según los puntos solicitados:

Resumen Técnico: Formulación HHO para Flujos Incompresibles de Densidad Variable

Título: Una formulación HHO para flujos incompresibles de densidad variable donde la densidad es puramente advectada.
Autores: Lorenzo Botti y Francesco Carlo Massa (Universidad de Bérgamo, Italia).

1. Planteamiento del Problema

El trabajo aborda la simulación numérica de mezclas de fluidos inmiscibles e incompresibles con densidad variable. El desafío principal en este tipo de problemas es garantizar la conservación exacta del volumen y asegurar que la densidad se transporte únicamente por advección (sin difusión numérica artificial ni fuentes espurias), manteniendo la estabilidad en regímenes dominados por la convección y con altos números de Reynolds.

El problema se modela mediante las ecuaciones de Navier-Stokes incompresibles con densidad variable, donde la ecuación de conservación de masa se formula como la derivada temporal material de la densidad ($\partial_t \rho + \nabla \cdot (\rho \mathbf{u}) = 0$), aprovechando la condición de incompresibilidad ($\nabla \cdot \mathbf{u} = 0$).

2. Metodología

Los autores proponen una extensión de los métodos Hybrid High-Order (HHO) para manejar densidad variable, combinada con esquemas de integración temporal de alto orden.

• Discretización Espacial (HHO):

  • Se utilizan espacios híbridos para velocidad, densidad y presión, definidos por polinomios en las celdas y en las caras de la malla.
  • Conformidad H(div): La elección de un espacio de presión híbrido y la reconstrucción adecuada permiten obtener una aproximación de velocidad que es $H(\text{div})$-conforme y puntualmente divergente libre. Esto es crucial para garantizar la conservación de volumen celda por celda.
  • Tratamiento de la Advección: Se emplean formulaciones antisimétricas (skew-symmetric) para los términos convectivos de momento y masa, complementadas con estabilizaciones de tipo upwind para garantizar robustez en regímenes de alta convección.
  • Dos variantes: Se presentan dos esquemas, HHO-$\pi^k$ y HHO-$\pi^{k+1}$, que difieren en el grado de proyección $L^2$ aplicado a la velocidad advectada en los términos convectivos. La variante $\pi^k$ es demostrablemente estable, mientras que $\pi^{k+1}$ ofrece mayor precisión en regímenes dominados por la convección.
  • Robustez de Presión: El esquema es "pressure-robust", lo que significa que los errores de velocidad y densidad son insensibles a la precisión de la aproximación de la presión; las fuerzas de cuerpo irrotacionales afectan solo al campo de presión.

• Discretización Temporal (ESDIRK):

  • Se utilizan esquemas ESDIRK (Explicit Singly Diagonal Implicit Runge-Kutta) de alto orden, diseñados específicamente para Ecuaciones Diferenciales Algebraicas (DAEs) de índice 2.
  • Esto es necesario porque el sistema discretizado espacialmente implica variables primitivas ($\mathbf{u}, \rho$) en lugar de conservativas, creando una estructura DAE.
  • Se incluyen estabilizaciones temporales para la densidad para garantizar la robustez cuando la velocidad normal en una cara se anula.

• Resolución Algebraica:

  • Se emplea la condensación estática para eliminar las incógnitas internas de las celdas, reduciendo el sistema global a solo las incógnitas de las caras. Esto disminuye drásticamente la huella de memoria y el costo computacional, especialmente a altos órdenes polinómicos ($k$).
  • Se utiliza un precondicionador p-multigrid para acelerar la convergencia de los solucionadores lineales iterativos (FGMRES) dentro del método de Newton.

3. Contribuciones Clave

• Conservación Exacta: Garantiza la conservación de volumen celda por celda hasta la precisión de la máquina, lo que asegura que la densidad sea una cantidad puramente advectada.
• Robustez y Estabilidad: Demuestra estabilidad discreta mediante análisis de energía, incluso en regímenes dominados por la convección y con condiciones de contorno complejas.
• Eficiencia Computacional: La combinación de condensación estática y precondicionamiento p-multigrid permite simulaciones de alto orden con un costo computacional y de memoria significativamente reducidos en comparación con métodos tradicionales de alto orden.
• Preservación de Límites: En órdenes bajos ($k=0$), el esquema preserva los límites de la densidad, evitando densidades negativas en flujos con altas relaciones de densidad (Atwood alto).

4. Resultados Numéricos

Los autores validan el método mediante tres casos de prueba:

1. Caso Kovasznay (Flujo Estacionario):

   • Se verifica la tasa de convergencia espacial. Se observan tasas óptimas de convergencia ($k+1$ para presión, $k+2$ para velocidad en norma $L^2$).
   • La densidad se mantiene constante (igual a 1) con errores de orden de la precisión de la máquina, confirmando la advección pura.
   • La variante HHO-$\pi^{k+1}$ muestra una precisión superior a $k=0$ en comparación con HHO-$\pi^k$.

2. Solución Manufacturada de Guermond y Quartapelle (Flujo Transitorio):

   • Se validan las tasas de convergencia temporal y espacial.
   • Se confirman las tasas de convergencia temporales esperadas para esquemas ESDIRK de 3 y 6 etapas (2º y 4º orden, respectivamente).
   • Se demuestra nuevamente la robustez de la presión: el error de velocidad no se degrada debido a errores en la presión.

3. Inestabilidad de Rayleigh-Taylor (RTI):

   • Bajo número de Atwood (At=0.5): Se simula a números de Reynolds 1000 y 5000. Se compara una malla gruesa con alto orden ($k=6$) contra una malla fina con bajo orden ($k=1$). Los resultados muestran un excelente acuerdo, demostrando que el alto orden permite capturar la física con menos grados de libertad, reduciendo el tiempo de cálculo.
   • Alto número de Atwood (At=0.75): Se simula con $k=0$ en una malla fina. El esquema preserva los límites de densidad (evitando densidades negativas) y captura la evolución del penacho principal, aunque con menor resolución de las características más finas de la interfaz debido al bajo orden polinómico.

5. Significado e Impacto

Este trabajo representa un avance significativo en la simulación de flujos multifásicos incompresibles. La formulación propuesta combina la exactitud matemática (conservación estricta, conformidad $H(\text{div})$) con la eficiencia computacional (condensación estática, p-multigrid).

La capacidad de simular mezclas de fluidos inmiscibles con alta precisión en espacio y tiempo, manteniendo la estabilidad en regímenes de alta convección y altas relaciones de densidad, abre nuevas posibilidades para el estudio de fenómenos complejos como la inestabilidad de Rayleigh-Taylor y, potencialmente, flujos multifásicos más complejos (como se menciona en la conclusión, la extensión futura al sistema Cahn-Hilliard-Navier-Stokes). La metodología demuestra que es posible lograr soluciones de alto orden con un costo computacional manejable, superando las limitaciones de memoria y estabilidad de métodos anteriores.