Turbulence Modelling of Mixing Layers under Anisotropic Strain

Este estudio investiga el impacto de las tasas de deformación anisotrópicas en las capas de mezcla turbulentas utilizando el modelo K-L RANS, demostrando que un cierre de deformación transversal mejora la precisión predictiva frente al enfoque isotrópico predeterminado y sugiriendo modificaciones correspondientes para los modelos K-$\epsilon$ y K-$\omega$.

Lo esencial

Imagina que tienes un tazón con dos líquidos de diferentes colores, como aceite y agua, sobre una mesa. Si sacudes la mesa, el límite entre ellos se vuelve desordenado y comienza a mezclarse. Esto es similar a lo que ocurre en una capa de mezcla turbulenta en física: dos fluidos de diferentes densidades son empujados juntos, creando un caos revuelto y giratorio.

Este artículo trata sobre comprender qué sucede cuando no solo sacudes la mesa, sino que también estiras o comprimes toda la habitación donde ocurre la mezcla.

Aquí tienes un desglose de la historia del artículo, utilizando analogías sencillas:

1. El Escenario: Estirar la Habitación

En muchos escenarios del mundo real, como una estrella explotando (supernova) o una bomba de fusión nuclear siendo comprimida, el espacio donde los fluidos se mezclan no está simplemente quieto. El espacio mismo se está expandiendo o contrayendo.

• La Analogía: Imagina que la capa de mezcla es un trozo de masa siendo amasada. Por lo general, los científicos estudian cómo se mezcla la masa cuando solo la empujas alrededor. Pero en este artículo, los autores preguntan: "¿Qué sucede si, mientras amasas, alguien también tira de la mesa donde está la masa, estirándola longitudinalmente o comprimiéndola transversalmente?"
• El Problema: El "estiramiento" (deformación) no es el mismo en todas las direcciones. Si estiras una banda elástica, se vuelve más larga en una dirección pero más delgada en las otras. Esto se llama deformación anisotrópica. La mayoría de los modelos informáticos utilizados para predecir estas mezclas asumen que el estiramiento es el mismo en todas las direcciones (como inflar un globo perfecto), lo cual no coincide con la realidad.

2. La Herramienta: El Modelo "K-L"

Para predecir cómo se mezclan los fluidos, los autores utilizan un programa informático llamado el modelo de turbulencia K-L.

• La Analogía: Piensa en este modelo como un libro de recetas para predecir el caos. Tiene dos ingredientes principales que rastrea:
  1. Cuánta energía hay en los remolinos (Energía Cinética Turbulenta).
  2. Qué tan grandes son los remolinos (Escala de Longitud Turbulenta).
• El modelo intenta adivinar qué tan grandes se volverán los remolinos a medida que los fluidos se mezclan. La parte complicada es una regla en la receta llamada el término de "compresión volumétrica". Esta regla le dice al modelo cómo cambia el tamaño de los remolinos cuando toda la habitación está siendo comprimida o estirada.

3. El Experimento: Probar Tres Reglas Diferentes

Los autores ejecutaron simulaciones informáticas para ver qué "regla" para la compresión volumétrica funcionaba mejor cuando la habitación se estiraba en direcciones específicas. Probaron tres versiones de la receta:

1. La Regla "Promedio": Asume que el estiramiento es el mismo en todas las direcciones (la configuración predeterminada).
2. La Regla "Longitudinal": Asume que el tamaño de los remolinos cambia basándose únicamente en cuánto se estira la habitación a lo largo de la dirección de mezcla.
3. La Regla "Transversal": Asume que el tamaño de los remolinos cambia basándose en cuánto se estira la habitación a través de la dirección de mezcla (perpendicular al flujo).

4. Los Resultados: Gana la Regla "Transversal"

Los autores compararon sus predicciones informáticas con simulaciones de alta resolución y muy detalladas (que actúan como una referencia "perfecta").

• El Hallazgo: La regla predeterminada "Promedio" fue aceptable, pero no excelente. La regla "Longitudinal" en realidad hizo que las predicciones fueran peores.
• El Ganador: La Regla "Transversal" (usando la deformación transversal) fue la más precisa.
• ¿Por qué? Los autores explican que cuando estiras una capa de mezcla, los grandes "torbellinos" (remolinos) se comportan de manera diferente dependiendo de la dirección. Resulta que el tamaño de estos remolinos es más sensible a cómo cambia el espacio transversalmente (de lado) que a cómo cambia longitudinalmente. Al usar el estiramiento lateral para ajustar el tamaño de los remolinos en la receta, el modelo predijo el ancho de la mezcla y la energía con mucha mayor precisión.

5. El Panorama General: Una Nueva Receta de "Tres Partes"

El artículo también examinó cómo simplificar estas ecuaciones complejas en un modelo de "Flotabilidad-Resistencia" (una forma más sencilla de pensar en la mezcla).

• Se dieron cuenta de que el "ancho de la mezcla" y el "tamaño de los remolinos" en realidad reaccionan a fuerzas diferentes. El ancho se estira con la tracción longitudinal, pero el tamaño del remolino reacciona a la compresión transversal.
• La Conclusión: Para obtener la mejor predicción, necesitas un modelo que trate estas dos cosas por separado. En lugar de una regla para todo, necesitas un modelo de tres partes que evolucione el ancho y el tamaño del remolino de forma independiente.

Resumen

En resumen, este artículo trata sobre corregir un modelo informático utilizado para predecir cómo se mezclan los fluidos cuando el espacio que los rodea se está distorsionando. Los autores descubrieron que la forma estándar de calcular cómo los "remolinos" se encogen o crecen era incorrecta para estas condiciones específicas. Al cambiar la regla para observar cómo se estira el espacio transversalmente en lugar de simplemente promediarlo, hicieron que el modelo fuera mucho más preciso. Esto ayuda a los científicos a comprender mejor eventos complejos como explosiones estelares o experimentos de energía de fusión, donde los fluidos se comprimen y estiran constantemente de manera desigual.

Resumen técnico

A continuación se presenta un resumen técnico detallado del artículo "Modelado de turbulencia de capas de mezcla bajo deformación anisotrópica" de Pascoe, Groom y Thornber.

1. Planteamiento del Problema

Las capas de mezcla turbulentas, como las generadas por las inestabilidades de Rayleigh-Taylor (RTI) y Richtmyer-Meshkov (RMI), son fundamentales para la física de alta densidad de energía, incluyendo la fusión por confinamiento inercial (ICF) y fenómenos astrofísicos como las supernovas. En estos escenarios, el flujo a menudo experimenta tasas de deformación anisotrópicas debido al movimiento radial en geometrías convergentes (implosiones) o divergentes (explosiones).

Los modelos actuales de turbulencia de Navier-Stokes promediados de Reynolds (RANS), particularmente el ampliamente utilizado modelo K-L (un modelo de dos ecuaciones que resuelve la energía cinética turbulenta $K$ y la escala de longitud turbulenta $L$), suelen depender de supuestos de cierre derivados para deformación isotrópica o flujos incompresibles. El cierre estándar para el término de compresión volumétrica en la ecuación de transporte de la escala de longitud asume que la escala de longitud evoluciona basándose en la tasa de deformación isotrópica media (divergencia de la velocidad dividida por 3). Sin embargo, en escenarios anisotrópicos realistas (por ejemplo, implosiones esféricas), las tasas de deformación en las direcciones radial (axial) y circunferencial (transversal) difieren significativamente. El artículo aborda la falta de definición teórica y modelado preciso de cómo la compresión volumétrica afecta la escala de longitud turbulenta bajo estas condiciones anisotrópicas.

2. Metodología

Los autores realizaron una investigación sistemática utilizando simulaciones RANS unidimensionales para aislar los efectos de la dirección de la deformación en la evolución de la capa de mezcla.

• Marco Numérico: Las simulaciones utilizaron el modelo de turbulencia K-L implementado en el código FLAMENCO. El modelo resuelve ecuaciones de transporte para la energía cinética turbulenta promediada de Favre ($\tilde{K}$) y la escala de longitud turbulenta ($L$).
• Casos de Prueba: El estudio se basó en el caso de inestabilidad de Richtmyer-Meshkov del "grupo $\theta$ a escala de cuarto" (configuración de pesado a ligero, número de Atwood $A_t=0.5$, número de Mach del choque $M_a \approx 1.84$).
• Aplicación de Deformación: Para aislar los efectos de la anisotropía, se aplicaron tasas de deformación normal uniformes al dominio de simulación en dos direcciones distintas:
  • Deformación Axial ($S_A$): Deformación aplicada en la dirección de la mezcla (normal a la interfaz).
  • Deformación Transversal ($S_T$): Deformación aplicada en las direcciones paralelas a la interfaz (plano homogéneo).
  • Se probaron perfiles de velocidad constante y de tasa de deformación constante.
• Variaciones de Cierre: Se compararon tres cierres diferentes para el término de compresión volumétrica ($\bar{\rho} L S_\phi$) en la ecuación de la escala de longitud:
  1. Cierre Isotrópico (Predeterminado): Escala $L$ con la tasa de deformación media ($S_I = \nabla \cdot \mathbf{u} / 3$).
  2. Cierre Axial: Escala $L$ con la tasa de deformación axial ($S_A$).
  3. Cierre Transversal: Escala $L$ con la tasa de deformación transversal ($S_T$).
• Validación: Los resultados RANS se compararon con datos de simulación de grandes remolinos implícita de alta fidelidad (ILES) de estudios anteriores (Pascoe et al.), los cuales sirvieron como la verdad fundamental.

3. Contribuciones Clave

• Identificación de Deficiencias en el Cierre: El artículo demuestra que el supuesto estándar de cierre isotrópico falla al predecir con precisión la evolución de las capas de mezcla turbulentas bajo deformación anisotrópica, particularmente en lo que respecta a la energía cinética turbulenta (TKE) y el ancho integral.
• Propuesta de Cierre Transversal: Los autores proponen y validan una nueva estrategia de cierre donde la evolución de la escala de longitud turbulenta bajo compresión volumétrica es impulsada por la tasa de deformación transversal ($S_T$) en lugar de la deformación isotrópica media.
• Derivación de un Modelo de Flotación-Resistencia de Tres Ecuaciones: Mediante un análisis de similitud, los autores derivaron un modelo simplificado de flotación-resistencia. Descubrieron que una sola escala de longitud no puede capturar simultáneamente la escala del ancho integral (que sigue la deformación axial) y la escala de longitud turbulenta (que sigue la deformación transversal). En consecuencia, desarrollaron un modelo de tres ecuaciones que evoluciona el ancho integral y la escala de longitud turbulenta por separado.
• Generalización a Otros Modelos: El estudio extiende los hallazgos a otros modelos populares de dos ecuaciones (K-$\epsilon$ y K-$\omega$), derivando modificaciones equivalentes para sus términos de compresión volumétrica basadas en el cierre de deformación transversal.

4. Resultados

• Ancho Integral ($W$):
  • Bajo deformación axial, el ancho de la capa de mezcla se estira o comprime directamente por el gradiente de velocidad del flujo medio. Todos los cierres capturaron esta tendencia, pero el cierre transversal proporcionó la predicción más precisa de la tasa de crecimiento turbulento, minimizando el Error Porcentual Absoluto Medio (MAPE) en un factor de cuatro en comparación con el cierre isotrópico.
  • Bajo deformación transversal, el cierre isotrópico mostró casi ninguna variación en el ancho (cancelando incorrectamente los efectos de producción de cizalladura y compresión), mientras que el cierre transversal predijo correctamente el ligero aumento en la tasa de crecimiento observado en los datos ILES.
• Energía Cinética Turbulenta (TKE):
  • El cierre transversal superó significativamente a los otros métodos en la predicción de la TKE.
  • El cierre axial tuvo el peor desempeño, a menudo prediciendo el signo incorrecto del cambio para la TKE bajo expansión.
  • El cierre isotrópico sobreestimó la influencia de la deformación en la TKE.
• Análisis de Similitud:
  • El análisis reveló que bajo deformación anisotrópica, la relación entre el pico de la escala de longitud turbulenta y el ancho integral ($\beta_W$) no es constante, violando el supuesto estándar de similitud utilizado en muchas calibraciones K-L.
  • Esto hizo necesario el desarrollo del modelo de flotación-resistencia de tres ecuaciones, donde el ancho integral escala con la deformación axial ($S_A$) y la escala de longitud turbulenta escala con la deformación transversal ($S_T$).
• Equivalencia de Modelos: Las modificaciones derivadas para el modelo K-L se mapearon exitosamente a los modelos K-$\epsilon$ y K-$\omega$, sugiriendo que el término de compresión volumétrica en estos modelos también debe ajustarse para tener en cuenta la deformación transversal en lugar de solo la divergencia media.

5. Significado

Este trabajo proporciona un avance crítico en el modelado de la mezcla turbulenta compresible en aplicaciones de alta densidad de energía.

• Capacidad Predictiva Mejorada: Al corregir el cierre de compresión volumétrica para tener en cuenta la anisotropía, el modelo K-L (y sus equivalentes) puede ahora predecir con mayor precisión el crecimiento de la capa de mezcla y la TKE en escenarios de implosión/explosión, que son centrales para la Fusión por Confinamiento Inercial (ICF) y las simulaciones astrofísicas (por ejemplo, supernovas).
• Perspectiva Física: El estudio aclara los roles físicos distintos de las deformaciones axial y transversal: la deformación axial impulsa principalmente el estiramiento geométrico de la capa de mezcla, mientras que la deformación transversal gobierna la evolución de los remolinos turbulentos (escala de longitud) y la disipación.
• Aplicación Práctica: Las modificaciones propuestas son relativamente simples de implementar en códigos RANS existentes, pero ofrecen mejoras sustanciales en la precisión sin el costo computacional de LES o DNS. El modelo de flotación-resistencia de tres ecuaciones derivado ofrece una herramienta computacionalmente eficiente para aplicaciones de ingeniería donde las simulaciones RANS completas son demasiado costosas.

En conclusión, el artículo argumenta que para entornos de deformación anisotrópica, la escala de longitud turbulenta debe modelarse evolucionando con la tasa de deformación transversal, un hallazgo que desafía la suposición de larga data de escalado isotrópico en los cierres de turbulencia estándar.