A Hybrid Finite-Volume Reconstruction Framework for Efficient High-Order Shock-Capturing on Unstructured Meshes

Resumen Técnico: Un Marco de Reconstrucción de Volumen Finito Híbrido para la Captura de Choques de Alto Orden Eficiente en Mallas No Estructuradas

Planteamiento del Problema
La Dinámica de Fluidos Computacional (CFD) para flujos compresibles enfrenta un compromiso persistente entre la precisión y el costo computacional, particularmente en mallas no estructuradas requeridas para geometrías industriales complejas. Los métodos de bajo orden son robustos y no oscilatorios cerca de los choques, pero sufren de una disipación numérica excesiva, fallando al resolver estructuras de flujo de escala fina en regiones suaves. Por el contrario, los esquemas de alto orden no oscilatorios, como los esquemas de WENO (Weighted Essentially Non-Oscillatory) y sus variantes, ofrecen una precisión superior en regiones suaves pero incurren en una carga computacional significativa y pueden tener dificultades con la robustez en flujos dominados por choques fuertes. Los enfoques híbridos existentes suelen depender de la detección a posteriori (por ejemplo, el paradigma MOOD), que puede ser computacionalmente costoso debido a la recomputación iterativa de celdas "problemáticas", o de detectores a priori que son demasiado sensibles o demasiado restrictivos, lo que conduce a una disipación o oscilaciones innecesarias.

Metodología
Los autores proponen un marco de reconstrucción híbrido de orden arbitrario y multidimensional implementado dentro del solver de volumen finito de código abierto UCNS3D. La innovación central es una estrategia de detección a priori de un solo paso que clasifica el campo de flujo local en tres regiones distintas, asignando un método de reconstrucción óptimo a cada una:

1. Regiones Suaves: Se emplea una reconstrucción lineal de alto orden. Esto evita el costo computacional de los esquemas no lineales donde no ofrecen beneficio de precisión y previene los problemas de degeneración de orden que a veces se ven en los esquemas WENO tradicionales en extremos suaves.
2. Regiones Débilmente No Suaves: Se aplica el esquema Compact Weighted Essentially Non-Oscillatory (CWENOZ). Este maneja regiones con gradientes fuertes o turbulencia en desarrollo donde se necesita precisión de alto orden pero las discontinuidades aún no se han formado completamente.
3. Regiones Discontinuas: Se activa un esquema MUSCL de segundo orden (con limitadores de pendiente) cerca de choques y discontinuidades reales para garantizar la estabilidad y el comportamiento no oscilatorio.

Componentes Clave:

• Detector de Admisibilidad Numérica (NAD) Rediseñado: A diferencia del enfoque iterativo a posteriori del paradigma MOOD, este NAD utiliza un Principio de Máximo Discreto (DMP) relajado para clasificar las celdas en un solo paso. Introduce parámetros auxiliares ($\delta_m$ y $\delta_w$) para definir márgenes para las regiones MUSCL y CWENOZ, respectivamente. Esto permite una sintonización flexible para equilibrar la robustez y la resolución.
• Criterio de Desactivación: Para evitar la activación falsa en regiones de campo lejano suaves, se utiliza un criterio de desactivación similar al limitador de Venkatakrishnan, asegurando que los esquemas lineales de alto orden permanezcan activos donde el flujo es genuinamente suave.
• Reconstrucción: El marco utiliza un método de reconstrucción de mínimos cuadrados con un algoritmo de Estructura de Estarcido Compacto (SBC), que soporta órdenes polinómicos arbitrarios ($P3, P5, P7$, etc.) en mallas de elementos mixtos (triángulos, cuadriláteros, tetraedros, hexaedros).
• Integración de Flujo y Tiempo: El esquema emplea resolvedores de Riemann HLLC o HLL para flujos invariantes, mínimos cuadrados restringidos para gradientes viscosos, y un integrador de tiempo de Runge-Kutta de Estabilidad Fuerte Preservada (SSP) de 4.º orden.

Contribuciones Clave

• Estrategia Híbrida: La integración de esquemas lineales, CWENOZ y MUSCL en un marco unificado que minimiza el uso de reconstrucciones no lineales costosas sin sacrificar la estabilidad.
• Eficiencia: Al restringir las reconstrucciones costosas de CWENOZ y MUSCL solo a las regiones necesarias, el método reduce significativamente la carga computacional en comparación con los esquemas no lineales uniformes de alto orden.
• Flexibilidad: Los parámetros del NAD permiten que el marco opere como un híbrido de tres esquemas, un híbrido de dos esquemas (por ejemplo, Lineal-CWENOZ o MUSCL-CWENOZ), o incluso un enfoque de un solo esquema, adaptándose a la física del flujo y a los requisitos de resolución específicos.
• Implementación: Implementación completa en el solver de código abierto UCNS3D, con soporte para reconstrucciones de orden arbitrario en mallas no estructuradas.

Resultados
El marco fue validado mediante una suite exhaustiva de pruebas de referencia 2D y 3D:

• Evolución de Vórtice 2D: En flujo suave, el método híbrido produjo resultados idénticos al esquema lineal de alto orden, confirmando que no activa erróneamente los limitadores no lineales. Logró la misma precisión que CWENOZ pero con un tiempo de CPU significativamente menor (reduciendo el costo casi a la mitad en comparación con CWENOZ).
• Problema de Shu-Osher: El método capturó exitosamente las interacciones de choque-onda de entropía. El estudio destacó que la elección del limitador de pendiente (Barth & Jespersen vs. MOGE) impacta significativamente el rendimiento del detector, siendo el limitador de Barth & Jespersen el que proporcionó el mejor equilibrio entre robustez y precisión para este montaje híbrido específico.
• Reflexión de Mach Doble y Escalón de Cara Frontal: El método resolvió interacciones de choque complejas y capas de cizalladura. Aunque las simulaciones de mayor orden (P7) mostraron un poco más de disipación que el CWENOZ puro al resolver pequeños vórtices, esto se atribuyó a la activación del esquema MUSCL y podría mitigarse relajando los criterios híbridos (aumentando la región CWENOZ).
• Problema de Riemann 2D: El método híbrido demostró una excelente robustez en diversas configuraciones. La configuración por defecto mejoró la eficiencia computacional en un 41.67% sobre CWENOZ manteniendo capacidades de captura de choque comparables.
• Vórtice Taylor-Green Viscoso Supersónico (3D): En una simulación de flujo turbulento ($M=1.25, Re=1600$), el esquema híbrido de 6.º orden resolvió un rango más amplio de números de onda y coincidió mejor con las soluciones de referencia (TENO y DG) que el esquema CWENOZ6 uniforme. Crucialmente, logró una aceleración de 2.5× en comparación con el esquema CWENOZ6 uniforme en el mismo hardware.
• Flujo Mach 3 sobre un Cilindro: El método capturó el choque de proa y las inestabilidades de la estela de manera robusta, con el esquema lineal de alto orden dominando el dominio y los esquemas no lineales activándose solo cerca de los choques. Esto resultó en una reducción de tres veces el costo computacional por iteración en comparación con una discretización uniforme de CWENO6.

Significado y Reivindicaciones
El artículo afirma que este marco híbrido acerca la precisión de alto orden a las simulaciones de CFD prácticas de escala industrial. Al combinar el costo computacional reducido de los esquemas lineales en regiones suaves con la robustez de los esquemas no lineales cerca de los choques, el método ofrece una alternativa fiable a los esquemas no lineales uniformes de alto orden. Los autores enfatizan que el marco mantiene los órdenes de precisión diseñados en regiones suaves mientras garantiza soluciones no oscilatorias cerca de las discontinuidades. La capacidad de lograr aceleraciones significativas (hasta 2.5× en turbulencia 3D) sin comprometer la estabilidad o la precisión sugiere que este enfoque es un camino viable para simulaciones de alta fidelidad de flujos compresibles complejos en mallas no estructuradas.