Time integration as filtering: a space-time discretization-aware LES formulation

Resumen Técnico: La integración temporal como filtrado en la LES consciente de la discretización espacio-temporal

Planteamiento del problema
En la LES implícita (ILES), el filtro espacial suele identificarse con la propia malla computacional. Las formulaciones tradicionales de LES "conscientes de la discretización" abordan el problema del cierre de la turbulencia derivando un esfuerzo residual exacto basado en la propiedad de "intercambio de filtros" (filter-swap): se reconoce que la derivada central finita gruesa es una derivada espacial exacta filtrada mediante un filtro de tipo top-hat. Sin embargo, las formulaciones existentes tratan la integración temporal como un paso de post-procesamiento separado. En consecuencia, el error de la discretización temporal queda fuera de la garantía de esfuerzo residual exacto. A medida que los pasos de tiempo aumentan (números CFL más altos), los cierres estándar que solo consideran el espacio se degradan porque no logran contabilizar el error específico introducido por el integrador temporal, particularmente cuando se utilizan esquemas explícitos como Euler hacia adelante (Forward Euler).

Metodología
El artículo extiende el concepto de intercambio de filtros espacial a lo espacio-temporal al observar que la diferencia temporal de Euler hacia adelante es, en sí misma, una derivada temporal exacta filtrada (unilateral y de tipo top-hat).

1. Intercambio de filtros espacio-temporal: Al aplicar un filtro espacio-temporal $(\cdot)_{h,\tau} = (\cdot)_\tau \circ (\cdot)_h$ a la ley de conservación escalar, los autores derivan una ley discreta exacta donde tanto las derivadas espaciales como las temporales son reemplazadas por sus contrapartes filtradas.
2. Descomposición de flujo: El flujo espacio-temporal exacto se descompone en el flujo numérico más un flujo residual, $r_{h,\tau}$. Este residual se divide en cuatro términos de error distintos:
   • $r^{LES}_h$: El conmutador LES clásico (filtrado de malla frente a no linealidad).
   • $r^{num}_h$: El error de flujo numérico.
   • $r^{div}_h$: Una corrección para la divergencia discreta.
   • $r^{time}_\tau$: Un nuevo término que representa el error de cuadratura de flujo temporal (la diferencia entre el flujo promediado en el paso de tiempo y el flujo instantáneo).
3. Modelado: Los autores proponen un modelo de viscosidad de remolino consciente de la discretización espacio-temporal. Expanden el término de error de cuadratura temporal, $r^{time}_\tau$, demostrando que se comporta como un flujo difusivo proporcional a $\tau (f'(u))^2 \partial_x u$. Esto se identifica como una difusión de tipo Lax–Wendroff, la cual es no negativa por construcción y escala linealmente con el paso de tiempo $\tau$. El cierre final aumenta el modelo de Smagorinsky estándar con este término temporal.

Resultados clave
El marco de trabajo fue probado en la turbulencia de Burgers viscosa utilizando un esquema de volúmenes finitos de Euler hacia adelante con números CFL variables (0.06 a 0.89).

• Análisis a priori: En la descomposición del flujo residual, el término temporal ($r^{time}_\tau$) crece significamente con el número CFL, convirtiéndose en el término de error individual más grande en valores de CFL prácticos (por ejemplo, 0.59 del residual en un CFL de 0.89 para flujos centrales). Por el contrario, los términos espaciales mantienen proporciones fijas independientemente del paso de tiempo.
• Análisis a posteriori:
  • Flujo central: Los cierres estándar (sin modelo, Smagorinsky clásico y Smagorinsky de solo espacio) se degradan rápidamente a medida que el número CFL aumenta, con errores que suben de 0.1 a ~0.85. El cierre consciente del espacio-tiempo propuesto mantiene un error relativo bajo y constante (0.09) en todo el rango de números CFL.
  • Flujo hacia arriba (Upwind): El flujo numérico de tipo upwind es inherentemente sobre-disipativo. En este caso, el error temporal (que no es disipativo) cancela parcialmente la sobre-disipación del flujo upwind. Aunque el modelo espacio-temporal sigue siendo preciso, es ligeramente sobre-disipativo en pasos grandes porque ajusta los efectos de cancelación de forma independiente en lugar de explotar la cancelación fortuita.

Significado y afirmaciones
El artículo sostiene que reconocer la diferencia temporal de Euler hacia adelante como un filtro top-hat permite una descomposición de flujo discreta exacta en el espacio-tiempo. La contribución principal es la identificación del error de cuadratura de flujo temporal como un término de residual distinto y dominante en pasos de tiempo gruesos.

• Los autores demuestran que un cierre aumentado con el término de orden principal de este error (difusión Lax–Wendroff) permite que las simulaciones de Euler hacia adelante mantengan la precisión incluso en pasos de tiempo grandes donde los cierres de solo espacio fallan.
• El trabajo establece que el integrador temporal contribuye exactamente con un nuevo término al residual (el error de cuadratura) y ningún término de conmutación temporal separado.
• El artículo señala modestamente que gran parte del error de cuadratura temporal es no disipativo y, por lo tanto, invisible para los modelos estándar de viscosidad de remolino, sugiriendo que los cierres aprendidos por aprendizaje automático (machine-learned) entrenados sobre el objetivo exacto o sobre integradores implícitos/multietapa son extensiones naturales a futuro.