Ensemble algorithm for parametrized flow problems with energy stable open boundary conditions

Resumen Técnico: Algoritmo de Conjunto para Problemas de Flujo Parametrizados con Condiciones de Contorno Abiertas de Estabilidad Energética

Planteamiento del Problema
El artículo aborda el desafío computacional de resolver conjuntos (ensembles) de ecuaciones de Navier-Stokes incomprimibles sujetas a condiciones iniciales, fuerzas corporales y coeficientes de viscosidad ($\nu_j$) variables. En aplicaciones como el análisis de sensibilidad, la modelación de orden reducido y el pronóstico de conjuntos, el número de miembros del conjunto ($J$) puede ser grande. Los métodos semi-implícitos estándar requieren ensamblar y resolver un sistema lineal distinto para cada miembro $j$ en cada paso de tiempo debido a la dependencia de la matriz del sistema de $\nu_j$ y los términos no lineales. Esto resulta en costos computacionales prohibitivos. Aunque los enfoques totalmente explícitos permiten un único ensamblaje de sistema, imponen condiciones de paso de tiempo restrictivas, particularmente para flujos con números de Reynolds elevados. Los métodos de conjunto previos para viscosidad constante o variaciones específicas de viscosidad a menudo requerían supuestos restrictivos sobre los coeficientes de viscosidad (por ejemplo, diferencias relativas acotadas) para garantizar la estabilidad.

Metodología
Los autores proponen un nuevo método de cálculo de conjuntos que desacopla la matriz del sistema lineal del índice del miembro del conjunto $j$, permitiendo la solución de un único sistema lineal con múltiples lados derechos en cada paso de tiempo. El núcleo de la metodología consiste en:

1. Estrategia de División: El término convectivo no lineal se divide en un componente de la media del conjunto y un componente fluctuante. El componente medio se trata de forma semi-implícita, mientras que el componente fluctuante se trata de forma explícita.
2. Tratamiento de la Viscosidad: A diferencia de trabajos anteriores (por ejemplo, [8]) que requerían límites específicos en la variación de los coeficientes de viscosidad para asegurar la estabilidad, este artículo introduce un tratamiento diferente del término difusivo. Al definir una viscosidad máxima $\nu_\infty = \max_j \nu_j$ y dividir el término de difusión, los autores derivan un esquema donde la matriz del sistema lineal es independiente de $j$ sin imponer restricciones sobre la magnitud de las variaciones de la viscosidad.
3. Condiciones de Contorno Abiertas (OBC): El método se extiende a dominios con contornos abiertos ($\Gamma_N$). El contorno se particiona en regiones de flujo de salida (outflow) y de reflujo (backflow). Los autores utilizan una condición de contorno de estabilidad energética que involucra una escala de longitud característica $L$ y una función de tipo Heaviside para manejar la naturaleza convectiva de la salida.
4. Discretización: El artículo presenta esquemas de paso de tiempo de primer orden (Algoritmos 1–3) y de segundo orden (Algoritmos 4–6). La discretización espacial utiliza espacios de elementos finitos conformes (específicamente pares Taylor-Hood $P_2/P_1$) con formas trilineales sesgo-simetrizadas para asegurar la estabilidad energética.
5. Análisis de Estabilidad: Pruebas teóricas establecen la estabilidad energética para los esquemas de primer orden bajo condiciones CFL específicas que involucran fluctuaciones de velocidad. Para los esquemas de segundo orden, aunque no se proporciona una prueba de estabilidad rigurosa sin restricciones de viscosidad (debido a la complejidad del término $\nabla \partial_{tt} u$), los experimentos numéricos sugieren estabilidad bajo las mismas condiciones que los esquemas de primer orden.

Contribuciones Clave

• Estabilidad de Viscosidad Sin Restricciones: La principal contribución teórica es un esquema de conjunto para coeficientes de viscosidad variables que garantiza la estabilidad sin la restriciva suposición $|\nu_j - \nu|/\nu \leq \sqrt{\mu}$ encontrada en la literatura previa.
• Extensión a Contornos Abiertos: Los autores extienden los métodos de conjunto a problemas con contornos abiertos, derivando formulaciones débiles y probando la estabilidad energética para esquemas de primer orden bajo restricciones específicas de paso de tiempo.
• Eficiencia Algorítmica: Los algoritmos propuestos resultan en una única matriz de coeficientes por paso de tiempo, la cual es resoluble eficientemente utilizando solucionadores de múltiples lados derechos, reduciendo significamente el costo computacional en comparación con la resolución de $J$ sistemas independientes.

Resultos
Se realizaron experimentos numéricos utilizando el paquete FreeFem++:

1. Estudio de Convergencia: Una prueba de convergencia en un dominio con soluciones exactas confirmó las tasas de convergencia de segundo orden esperadas tanto para la velocidad como para la presión en los esquemas de segundo orden propuestos. La precisión fue comparable a la de las simulaciones secuenciales independientes.
2. Flujo alrededor de un Cilindro: Los algoritmos fueron probados en un problema de referencia que involucra el flujo alrededor de un cilindro con tres valores de viscosidad diferentes. Los resultados (arrastre, sustentación y caída de presión) coincidieron con los valores de referencia de la literatura utilizando condiciones de Dirichlet y de salida de "no hacer nada" (do-nothing). La condición de contorno abierta permitió que los vórtices salieran del dominio de forma natural, evitando los problemas de reflexión artificial asociados con los perfiles de salida parabólicos prescritos.
3. Flujo en Canal con Contracción: Se simuló un escenario de flujo complejo con dos salidas y viscosidades variables. Los resultados del conjunto coincidieron cualitativamente con las corridas de Simulación Numérica Directa (DNS) independientes, a pesar de que las simulaciones de conjunto se realizaron en una malla más gruesa estabilizada por un filtro no lineal adaptativo.

Significancia y Reivindicaciones
El artículo afirma ofrecer un marco robusto y eficiente para simulaciones de flujo parametrizadas. Al eliminar las suposiciones restrictivas sobre las variaciones de viscosidad requeridas por los métodos de conjunto anteriores, el algoritmo propuesto amplía la aplicabilidad de las simulaciones de conjuntos a problemas con parámetros físicos altamente variables. La extensión a condiciones de contorno abiertas con estabilidad energética demostrable es un paso significativo hacia simulaciones realistas de flujos no acotados. Los autores concluyen que, si bien los esquemas de conjunto funcionan bien cualitativa y cuantitativamente frente a las simulaciones independientes, se requiere más investigación para cuantificar plenamente las ventajas computacionales, particularmente respecto a las comparaciones de tiempo de CPU en mallas refinadas adaptativamente para flujos con números de Reynolds altos.