A family of second order, linear, unconditionally stable methods for the Cahn-Hilliard-Navier-Stokes equations

Resumen Técnico: Una familia de métodos lineales de segundo orden y estabilidad incondicional para las ecuaciones de Cahn-Hilliard-Navier-Stokes

Planteamiento del Problema
La simulación numérica precisa de la dinámica de fluidos multifásicos complejos, específicamente aquellos modelados por las ecuaciones de Cahn-Hilliard-Navier-Stokes (CHNS), presenta desafíos significativos debido a la rigidez matemática del operador de Cahn-Hilliard de cuarto orden y su acoplamiento no lineal con las ecuaciones de Navier-Stokes incompresibles. Los esquemas estándar de paso de tiempo explícito enfrentan restricciones severas de paso de tiempo (por ejemplo, $\Delta t \sim O(\Delta x^4)$) para mantener la estabilidad. Si bien los métodos totalmente implícitos o de división convexa (convex-splitting) ofrecen estabilidad de energía incondicional, a menudo requieren resolver sistemas no lineales acoplados y de gran escala en cada paso de tiempo, lo que aumenta el costo computacional. Existe la necesidad de algoritmos que combinen precisión temporal de segundo orden, estabilidad incondicional y solvabilidad lineal sin sacrificar las propiedades de conservación física del sistema.

Metodología
Los autores proponen una familia de métodos IMEX (implícito-explícito) lineales de segundo orden y estabilidad incondicional para flujos de dos fases con densidad igualada. El esquema semidiscreto integra tres componentes principales:

1. Extrapolación de Términos No Lineales: Los términos convectivos y de acoplamiento no lineales se tratan explícitamente utilizando cantidades extrapoladas (denotadas como $H_{n+\theta}$), evitando la necesidad de iteraciones no lineales.
2. Formulación de Variables Auxiliares: El término de energía libre no lineal se gestiona mediante una variable auxiliar $q$ (derivada del potencial de doble pozo $F(\phi)$), reformulando el problema para mantener la linealidad.
3. Regularización de Curvatura Temporal: Se introduce un mecanismo de estabilización controlado por un parámetro $\epsilon \geq 0$. Esto involucra un operador de interpolación específico $J^\epsilon_{n+\theta}$ que incorpora la curvatura temporal discreta de la solución, un concepto adaptado de trabajos previos sobre regularización de Navier-Stokes.

El algoritmo resultante (Algoritmo 2.1) requiere resolver únicamente sistemas lineales con coeficientes constantes en cada paso de tiempo. El esquema está parametrizado por $\theta \in (1/2, 1]$ y $\epsilon \geq 0$.

Contribuciones Clave y Resultados Teóricos

• Estabilidad Incondicional: El artículo establece una estimación de energía discreta rigurosa que demuestra que el esquema propuesto es incondicionalmente estable a largo plazo para $\theta \in (1/2, 1]$ y $\epsilon \geq 0$. La prueba utiliza matrices definidas positivas simétricas para definir normas discretas y demuestra que el esquema satisface una ley de disipación de energía discreta sin restricciones de paso de tiempo.
• Solvabilidad Lineal: A diferencia de los enfoques totalmente implícitos o de división convexa que pueden requerir resolvedores no lineales, este método produce sistemas lineales en cada paso de tiempo, reduciendo significativamente la complejidad computacional.
• Precisión de Segundo Orden: El análisis teórico y los experimentos numéricos confirman que el método mantiene una precisión temporal de segundo orden.

Resultados Numéricos
Los autores validan el método a través de una serie de computaciones de referencia:

• Análisis de Convergencia: Utilizando el Método de Soluciones Manufacturadas (MMS), el esquema demuestra una convergencia temporal aproximadamente de segundo orden para las variables de velocidad, presión y campo de fase a través de diferentes configuraciones de elementos finitos ($P2-P1-P2-P2$ y $P3-P2-P3-P3$).
• Robustez y Regularización: En el problema de la cavidad impulsada por tapa (lid-driven cavity) de dos fases, el estudio examina la interacción entre $\theta$ y $\epsilon$. Se observa que para valores de $\theta$ cercanos al límite de Crank-Nicolson (por ejemplo, $\theta=0.51$) con $\epsilon=0$, el esquema puede volverse inestable debido a oscilaciones temporales débilmente amortiguadas. La introducción de una regularización de curvatura positiva ($\epsilon > 0$) amortigua eficazmente estas oscilaciones, restaurando la acotación y la robustez incluso para pasos de tiempo mayores.
• Pruebas de Referencia Físicas:
  • Descomposición Espinodal: El método captura con precisión la disipación de energía y la conservación de masa en simulaciones de largo tiempo ($T=100$), mostrando una separación de fases y un crecimiento de grano (coarsening) efectivos.
  • Relajación de Gotas: El algoritmo resuelve con éxito el movimiento de la interfaz impulsado por la tensión superficial, transicionando de una gota cuadrada a una circular de equilibrio mientras preserva la masa.
  • Cavidad Impulsada por Tapa y la Inestabilidad de Rayleigh-Taylor: El esquema maneja la fuerte cizalladura, la deformación interfacial y las inestabilidades impulsadas por la densidad (incluyendo el roll-up de Kelvin-Helmholtz en regímenes de baja viscosidad) sin oscilaciones espurias visibles.

Significancia
El artículo afirma que la familia de métodos propuesta ofrece una alternativa práctica y eficiente para las discretizaciones de sistemas CHNS. Al desacoplar las no linealidades mediante extrapolación y variables auxiliares, y al estabilizar la discretización temporal mediante la regularización de la curvatura, el método logra un equilibrio entre precisión de segundo orden, estabilidad incondicional y costo computacional lineal. Los autores enfatizan que el esquema preserva las estructuras físicas principales del modelo, incluyendo la disipación de energía y la conservación de la masa, lo que lo hace adecuado para simular flujos interfaciales complejos en diversos regímenes.