Consistent closure modeling in large eddy simulations by direct approximation of the filtered advection term

Este artículo propone un modelo de cierre consistente para la simulación de grandes remolinos (LES) mediante la aproximación directa del término de advección filtrada, lo que elimina inconsistencias conceptuales y mejora la precisión de los espectros de energía cinética y las correlaciones de velocidad en comparación con los métodos LES clásicos.

Lo esencial

Imagina que estás intentando predecir cómo se moverá el agua en un río muy turbulento. El agua tiene remolinos gigantes, remolinos medianos y remolinos diminutos que giran a toda velocidad.

Para los científicos, simular esto en una computadora es un desafío enorme. Si intentas calcular cada gota y cada micro-remolino, la computadora se desbordaría. Así que, en lugar de eso, usan una técnica llamada Simulación de Grandes Remolinos (LES).

El Problema: La "Lupa" imperfecta

Imagina que tienes una lupa (un filtro) para mirar el río.

• Lo que quieres ver: Los grandes remolinos (los que puedes ver a simple vista).
• Lo que ignoras: Los micro-remolinos (los que son demasiado pequeños para tu lupa).

En la física tradicional (la que se ha usado durante décadas), hay un error conceptual en cómo se usa esta lupa. Los científicos calculaban el movimiento de los grandes remolinos multiplicando sus velocidades entre sí. Pero, matemáticamente, cuando multiplicas dos cosas grandes, a veces creas "fantasmas" de cosas muy pequeñas (ondas de alta frecuencia) que tu lupa no debería ver.

La analogía de la foto borrosa:
Es como si intentaras tomar una foto de un paisaje borroso (los grandes remolinos), pero al procesar la foto, el software de edición añadiera accidentalmente ruido de estática y píxeles de colores extraños (los micro-remolinos fantasma).

Para arreglar esto, los científicos actuales tienen que usar "parches" o "filtros de emergencia" (llamados limitadores o estabilizadores) para borrar manualmente esos píxeles extraños. El problema es que estos parches hacen que el resultado dependa de qué tan buena sea la cámara (la malla de la computadora) y no de la realidad física. Si cambias la cámara, cambian los resultados. ¡No es justo!

La Solución: La "Receta Exacta"

Este artículo propone una forma nueva y más limpia de hacer las cosas. En lugar de multiplicar las velocidades y luego intentar arreglar los errores, los autores (Max Hausmann y Berend van Wachem) crearon una receta matemática exacta para calcular cómo se mueven los grandes remolinos, sin generar esos "fantasmas" de micro-remolinos.

La analogía de la escalera:
Imagina que quieres describir el movimiento del agua con una precisión perfecta.

1. El método antiguo: Era como intentar subir una escalera saltando de dos en dos escalones. A veces te caías (errores) y tenías que usar una red de seguridad (los parches) para no caer al vacío.
2. El nuevo método: Es como construir una escalera infinita donde cada peldaño es un poco más fino. Los autores dicen: "No necesitamos la escalera infinita completa. Solo necesitamos los primeros 3 o 4 peldaños para estar casi perfectamente seguros".

Esta "receta" (una serie matemática) se asegura de que solo veamos lo que nuestra lupa permite ver. No crea micro-remolinos falsos.

¿Qué descubrieron?

Los autores probaron su nueva receta en dos escenarios:

1. Turbulencia que se desvanece: Como el agua que deja de moverse después de una tormenta.
2. Flujo de corte: Como el viento que sopla fuerte en una dirección y lento en otra.

Los resultados fueron sorprendentes:

• Sin parches: Su método funcionaba perfectamente sin necesidad de esos "parches" de emergencia. La computadora no se confundía.
• Más realista: La forma en que se movía el agua en la simulación se parecía mucho más a la realidad (como si vieras una foto HD en lugar de una pixelada).
• Convergencia: Si hacías la simulación con una computadora más potente (más píxeles), el resultado no cambiaba drásticamente; se volvía simplemente más preciso. Esto significa que el método es sólido y confiable.

En resumen

Este artículo es como decir: "Dejemos de intentar arreglar los errores de nuestra lupa con cinta adhesiva. En su lugar, vamos a usar una lupa que no hace errores desde el principio."

Han encontrado una manera de simular el caos de la turbulencia que es matemáticamente consistente, más precisa y que no depende de trucos numéricos para funcionar. Es un paso gigante hacia predicciones del clima, diseño de aviones y motores que sean mucho más fiables y seguros.

Resumen técnico

Aquí presento un resumen técnico detallado del artículo "Consistent closure modeling in large eddy simulations by direct approximation of the filtered advection term" (Modelado consistente de cierre en simulaciones de grandes remolinos mediante aproximación directa del término de advección filtrada), escrito por Max Hausmann y Berend van Wachem.

1. El Problema: Inconsistencia Conceptual en la LES Clásica

El artículo identifica una inconsistencia conceptual fundamental y a menudo ignorada en el marco de las Simulaciones de Grandes Remolinos (LES, por sus siglas en inglés).

• La Inconsistencia: En la LES clásica, el término de advección filtrada ($\overline{u_i u_j}$) se aproxima mediante el producto de las velocidades filtradas ($\overline{u}_i \overline{u}_j$) más un tensor de esfuerzos de subescala modelado ($\tau_{ij}$).
• El Dilema de los Números de Onda: Aunque la velocidad filtrada $\overline{u}_i$ está limitada en banda (contiene solo números de onda hasta un corte $k_c$ definido por el filtro), el producto de dos funciones filtradas ($\overline{u}_i \overline{u}_j$) genera matemáticamente componentes de frecuencia más altos (hasta $2k_c$).
• Consecuencias:
  • El término de advección clásico introduce números de onda que exceden la definición de una ecuación filtrada, violando la premisa de que la solución debe contener solo escalas resolubles.
  • Esto obliga a los simuladores a utilizar limitadores de flujo, términos de estabilización o desaliasing para amortiguar estos números de onda no resueltos.
  • La solución se vuelve dependiente de la malla y del esquema numérico, violando el principio fundamental de que la solución debe converger al refinar la malla.
  • En la LES implícita, se confía exclusivamente en el error de discretización, lo que impide una validación rigurosa con DNS filtrada explícitamente (FDNS).

2. Metodología: Aproximación Directa del Término de Advección Filtrada

Los autores proponen una alternativa consistente: aproximar directamente el término de advección filtrada ($\overline{u_i u_j}$) en lugar de descomponerlo en un producto de velocidades filtradas más un modelo de estrés.

• Derivación Matemática:

  • Se asume un filtro Gaussiano con desviación estándar $\sigma$.
  • Utilizando una expansión en serie de Taylor del término de advección en función del ancho del filtro al cuadrado ($\sigma^2$), derivan una expresión exacta infinita.
  • La expansión se basa en la relación diferencial del filtro Gaussiano: $\frac{\partial \Phi}{\partial (\sigma^2)} = \frac{1}{2} \nabla^2 \Phi$.
  • Al multiplicar las series de Taylor de dos cantidades de flujo genéricas, obtienen una expresión para $\overline{u_i u_j}$ que es una suma de términos que involucran derivadas espaciales de las velocidades filtradas.

• Formulación de los Términos:

  • Orden 0: $\overline{u}_i \overline{u}_j$ (equivalente a la LES clásica).
  • Orden 2: Incluye términos de segundo orden en $\sigma$ con derivadas espaciales de segundo orden.
  • Orden 4: Incluye términos de cuarto orden en $\sigma$.
  • Propiedad Clave: Cada término en esta expansión es explícitamente filtrado. Esto garantiza que el término de advección aproximado esté limitado en banda y decaiga tan rápido como el núcleo del filtro en el espacio espectral, evitando la introducción de números de onda inconsistentes.

• Validación A Priori:

  • Se utilizó un campo de velocidad de Turbulencia Isotrópica Homogénea (HIT) obtenido de DNS.
  • Se compararon las aproximaciones de orden 0, 2 y 4 con el término de advección filtrado explícitamente.
  • Se analizó la correlación y la transferencia de energía (espectro de energía cinética).

3. Contribuciones Clave

1. Expresión Exacta en Serie Infinita: Derivación de una expresión matemática exacta para el término de advección filtrada basada en una expansión en potencias del ancho del filtro, que contiene únicamente cantidades filtradas.
2. Consistencia Espectral: Demostración de que truncar la serie después de unos pocos términos (orden 2 o 4) produce un modelo que respeta estrictamente el límite de ancho de banda de la definición de la LES, eliminando la necesidad de limitadores de flujo o estabilización artificial.
3. Mejora en la Transferencia de Energía: El modelo propuesto captura con mayor precisión la transferencia bidireccional de energía (incluyendo el backscatter) entre las escalas filtradas y las subescala, superando a los modelos de viscosidad turbulenta clásicos (como Smagorinsky) que a menudo fallan en la distribución espectral.
4. Independencia de la Malla: El método converge bajo refinamiento de malla sin depender de errores de discretización para disipar energía, a diferencia de la LES clásica o la LES implícita.

4. Resultados

Los autores realizaron estudios a posteriori en dos casos de prueba:

A. Turbulencia en Decaimiento (Decaying Turbulence)

• Evolución Temporal: La LES con la aproximación de orden 2 y 4 predice la evolución de la energía cinética y la disipación viscosa con una precisión muy superior a la LES clásica (que decae demasiado lento).
• Espectro de Energía: La LES clásica con viscosidad turbulenta predice demasiada energía en los números de onda altos (insuficiente disipación). La aproximación de orden 4, especialmente cuando se combina con un término de viscosidad de Smagorinsky, reproduce el espectro de energía de la FDNS con gran precisión en todo el rango de números de onda.
• Estructuras del Flujo: La LES clásica produce estructuras de flujo elongadas (típicas de modelos de viscosidad turbulenta). La aproximación propuesta genera estructuras más isotrópicas y realistas, similares a las de la FDNS.
• Convergencia: Al refinar la malla (de $64^3$ a $128^3$), el espectro de energía de la LES propuesta permanece esencialmente invariante, confirmando la independencia de la malla. La LES clásica muestra cambios significativos en el espectro con el refinamiento.
• Invarianza Galileana: Aunque la aproximación no es exactamente invariante Galileana (solo lo es hasta el orden de la aproximación), los estudios muestran que los errores introducidos al cambiar el marco de referencia son mínimos y de la magnitud del error de aproximación mismo, sin relevancia práctica significativa.

B. Flujo de Corte Turbulento (Turbulent Shear Flow)

• En un flujo de corte inhomogéneo y anisotrópico, la LES clásica sobreestimó significativamente las correlaciones de velocidad ($\langle u'u' \rangle$) y falló en mantener la turbulencia en mallas muy gruesas.
• La aproximación de orden 4 predijo con precisión tanto el perfil de velocidad media como las correlaciones de las fluctuaciones, incluso en mallas más gruesas donde la LES clásica se relaminarizaba.

5. Significado e Impacto

Este trabajo representa un avance conceptual importante en la modelado de turbulencia:

• Resolución de una Paradoja: Resuelve la paradoja de que la LES clásica requiere intervenciones numéricas (limitadores) para funcionar, lo que la hace dependiente de la malla y menos interpretable físicamente.
• Modelado Físico vs. Numérico: Cambia el enfoque de "corregir" errores numéricos a "modelar" correctamente la física de la advección filtrada.
• Robustez: El método es numéricamente robusto, no requiere esquemas de advección especiales (como limitadores de flujo) y permite el uso de esquemas centrales puros.
• Futuro: Proporciona una base para desarrollar LES más confiables, interpretables y precisas, especialmente en configuraciones donde la resolución es limitada o donde la consistencia física es crítica.

En conclusión, los autores demuestran que aproximar directamente el término de advección filtrada mediante una expansión en serie de potencias del ancho del filtro ofrece una alternativa superior a la descomposición clásica, logrando una consistencia espectral, una mejor predicción de la transferencia de energía y una convergencia robusta bajo refinamiento de malla.