Exact expressions for the unresolved stress in a finite-volume based large-eddy simulation

Resumen Técnico: Expresiones Exactas para el Estrés No Resuelto en una Simulación de Remolinos de Gran Escala Basada en Volúmenes Finitos

Planteamiento del Problema
La simulación de remolinos de gran escala (LES) tiene como objetivo modelar flujos turbulentos resolviendo las escalas grandes mientras se modelan los efectos de las escalas sub-rejilla no resueltas. Tradicionalmente, el problema de clausura en LES se centra en el error de conmutación derivado de la no conmutación de los términos de advección no lineales y el filtrado espacial, lo que resulta en el tensor de esfuerzo residual (RST) clásico, $\tau_{ij} = \overline{u_i u_j} - \bar{u}_i \bar{u}_j$. Sin embargo, en simulaciones prácticas, la LES se resuelve utilizando esquemas de discretización numérica, comúnmente el Método de Volúmenes Finitos (FVM).

Un problema crítico es que el ancho del filtro $\Delta$ en LES suele ser del mismo orden de magnitud que el espaciamiento de la malla $h$ (es decir, $\Delta \approx h$). En tales regímenes, los errores de discretización —que surgen de los flujos numéricos y los operadores de divergencia discretos— son del mismo orden que los esfuerzos sub-filtro físicos. Los enfoques clásicos suelen definir el modelo de clausura basándose en ecuaciones continuas y luego discretizarlas, o tratan los errores de discretización como problemas de convergencia separados. Esta separación conduce a un desajuste entre los datos objetivo (el residuo real) y el modelo, particularmente cuando la propia discretización introduce efectos no locales o no simétricos significativos que los modelos de clausura estándar, simétricos y locales (como el de Smagorinsky), no pueden capturar.

Metodología
Los autores proponen un marco matemático unificado denominado LES-FVM, que integra el proceso de filtrado de LES directamente con la discretización FVM. La metodología central consiste en:

1. Formalismo Unificado: En lugar de aplicar el filtrado LES a las ecuaciones continuas y luego discretizarlas, los autores aplican el filtro FVM (un promedio de volumen sobre los volúmenes de control) a las ecuaciones filtradas por LES. Esto crea un sistema "doblemente filtrado" donde el ancho del filtro $\Delta$ y el espaciamiento de la malla $h$ se tratan como parámetros independientes.
2. Derivación de Residuos Exactos: Utilizando una propiedad de "intercambio de filtro discreto" (donde el operador de diferencia finita es equivalente a una derivada exacta filtrada), los autores derivan expresiones exactas para el flujo residual en las ecuaciones discretizadas.
3. Descomposición de Residuos: El flujo residual total $\tau_{\Delta,h}$ se descompone en tres componentes distintos:
   • Término Clásico: El error de conmutación estándar asociado con el doble filtro ($\Delta$ y $h$).
   • Término de Flujo: Errores que surgen de la diferencia entre el flujo exacto y el flujo numérico ($r$ frente a $r_h$).
   • Término de Divergencia: Errores que surgen de la diferencia entre la divergencia exacta y la divergencia discreta.
4. Extensión al Flujo Incompresible 3D: El marco se extiende a las ecuaciones de Navier-Stokes incompresibles en 3D. Esto requiere el manejo del término de presión mediante operadores de proyección para asegurar que el campo de velocidad permanezca libre de divergencia en el entorno discreto. Los autores derivan un RST proyectado que da cuenta de la interacción entre la proyección de presión y la discretización.

Contribuciones Clave

• RST Informado por la Discretización Exacta: El artículo deriva una expresión exacta para el tensor de esfuerzo no resuelto en una LES basada en FVM. A diferencia del RST clásico, esta nueva expresión incluye explícitamente contribuciones de los flujos numéricos y los operadores de divergencia discretos.
• No Simetría y No Localidad: Para flujos incompresibles, los autores demuestran que el RST propuesto de LES-FVM es inherentemente no simétrico y no local (debido a la proyección de presión). Esto contrasta fuertemente con los RST de LES clásicos, que típicamente se asumen simétricos y locales.
• Error A-Posteriori Cero: En experimentos "asistidos por DNS" (donde el residuo exacto se calcula a partir de datos de Simulación Numérica Directa e se inyecta en la ecuación LES), la expresión de RST propuesta produce un error a-posteriori de cero. En contraste, las expresiones de RST clásicas resultan en errores que se acumulan con el tiempo.
• Mejora del Modelado de Clausura: Cuando se utiliza el nuevo RST como datos objetivo para ajustar el coeficiente de un modelo de Smagorinsky, el coeficiente resultante es mayor que cuando se ajusta al RST clásico. Esto conduce a un rendimiento mejorado en términos de espectros de energía y error relativo, particularmente para relaciones de filtro-malla pequeñas ($\Delta/h$).

Resultos
Los autores validan su teoría utilizando dos casos de prueba: la ecuación de Burgers viscosa en 1D y la turbulencia isotrópica decreciente en 3D.

• Ecuación de Burgers 1D:

  • Análisis A-Priori: El estudio muestra que para relaciones de filtro comunes ($\Delta \le 4h$), los errores de discretización (términos de flujo y divergencia) constituyen una parte significativa (hasta el 70-80%) del residuo total. Ignorar estos términos conduce a una representación deficiente del residuo real.
  • Análisis A-Posteriori: El uso del residuo informado por la discretización exacta como término de clausura resulta en una reconstrucción "perfecta" de la solución de la DNS filtrada.
  • Modelo de Smagorinsky: El ajuste del coeficiente de Smagorinsky al nuevo RST produce un coeficiente más alto en comparación con el ajuste al RST clásico. Las simulaciones que utilizan este coeficiente informado muestran errores relativos menores y espectros de energía que coinciden mejor con la DNS filtrada, particularmente en relaciones $\Delta/h$ bajas.

• Ecuaciones de Navier-Stokes 3D:

  • La derivación confirma que el RST en 3D es no simétrico y no local.
  • Los experimentos asistidos por DNS en turbulencia decreciente 3D confirman que el marco propuesto captura correctamente el esfuerzo residual, validando la extensión de la teoría 1D al flujo incompresible completo.

Significación y Reivindicaciones
El artículo sostiene que la separación tradicional entre "modelar la física" (clausura LES) y "modelar la numérica" (discretización) es defectuosa cuando el ancho del filtro y el espaciamiento de la malla son comparables. Al unificar estos procesos en el marco LES-FVM, los autores proporcionan una definición matemáticamente rigurosa del esfuerzo no resuelto que es consistente con el esquema numérico.

La principal significación radica en la demostración de que los errores inducidos por la discretización no son meramente ruido numérico, sino que son partes integrales del esfuerzo residual que deben contabilizarse en el modelo de clausura. Los autores argumentan que utilizar el RST clásico como objetivo para modelos de clausura ajustados o basados en datos introduce una inconsistencia fundamental, lo que conduce a coeficientes de modelo subóptimos y errores acumulados. Su RST "informado por la discretización" proporciona los datos objetivo correctos, permitiendo el desarrollo de modelos de clausura más precisos y estables, particularmente para LES implícita o escenarios con mallas gruesas donde $\Delta \approx h$.

Los autores mantienen la modestia respecto a la aplicación inmediata de estos hallazgos en flujos industriales complejos, señalando que, si bien el nuevo RST mejora el desempeño del modelo de Smagorinsky simple, el pleno potencial de este marco se realizará probablemente con modelos de clausura más expresivos, no locales y no simétricos (por ejemplo, enfoques basados en aprendizaje automático) que puedan explotar plenamente la estructura del residuo derivado.