Evolution of passive scalar mixing layers in stratified and unstratified homogeneous turbulence

Las simulaciones de grandes remolinos de alta resolución revelan que, si bien la mezcla de escalares pasivos en turbulencia homogénea estratificada se comporta de manera similar a los casos no estratificados en la dirección transversal, la estratificación inhibe severamente la mezcla vertical al suprimir la agitación a gran escala, lo que conduce a límites de crecimiento e intensidades de fluctuación distintos que informan enfoques de modelado específicos para el flujo escalar.

Lo esencial

Imagina que estás observando cómo una gota de tinta se dispersa a través de un vaso de agua. Si remueves el agua suavemente, la tinta se distribuye uniformemente. Así es como los científicos suelen pensar en la mezcla en fluidos: la turbulencia actúa como una cuchara gigante, removiendo todo hasta que se vuelve uniforme.

Pero, ¿qué sucede si ese agua está estratificada? Imagina que el agua en el fondo es pesada y salada, mientras que el agua en la parte superior es ligera y dulce. Esto se llama estratificación. En el mundo real, esto ocurre en el océano (donde el agua profunda es más densa) y en la atmósfera (donde el aire se vuelve más delgado y ligero a medida que subes).

Este artículo es un experimento informático de alta tecnología que se pregunta: ¿Cómo cambia esta estratificación la forma en que una "mancha" (como la contaminación o el humo) se dispersa a través de un fluido turbulento?

Aquí está la historia de lo que descubrieron, desglosada en conceptos simples.

El Escenario: Dos Tipos de "Manchas"

Los investigadores crearon un fluido virtual que inicialmente giraba de manera caótica (turbulento). Luego, introdujeron dos "manchas" diferentes (escalares pasivos) para ver cómo se comportaban:

1. La Mancha Horizontal: Una capa de tinta extendida de lado a lado, como una hoja de papel plana flotando en el agua.
2. La Mancha Vertical: Una capa de tinta extendida de arriba a abajo, como un muro vertical de color.

Realizaron dos simulaciones: una en agua normal (no estratificada) y otra en agua "estratificada" (estratificada).

El Gran Descubrimiento: La Diferencia "Arriba-Abajo" vs. "De Lado a Lado"

1. La Mancha de Lado a Lado (Capa Transversal)

Qué sucedió: Cuando la mancha se extendía de lado a lado, la estratificación no impidió su dispersión. De hecho, se dispersó ligeramente más rápido en el agua estratificada que en el agua normal.
La Analogía: Imagina una multitud de personas corriendo en un pasillo. Si el suelo es perfectamente plano (no estratificado), corren en todas direcciones. Si el suelo tiene una ligera, invisible pendiente (estratificada), aún corren de lado perfectamente bien, quizás incluso con un poco más de energía. La "tinta" se extiende ampliamente en ambos casos.
El Truco: Aunque la dispersión general fue similar, la "tinta" en el agua estratificada fue más "picuda". En lugar de un gradiente suave, tenía bordes más afilados y dentados. Fue más "intermitente", lo que significa que había bolsillos de tinta pura y bolsillos de agua pura con menos un terreno medio suave.

2. La Mancha Arriba-Abajo (Capa Vertical)

Qué sucedió: Aquí es donde ocurrió la magia (y la restricción). En el agua normal, la mancha vertical se extendía arriba y abajo fácilmente, al igual que la de lado a lado. Pero en el agua estratificada, la dispersión casi se detuvo por completo.
La Analogía: Imagina intentar remover un batido espeso con una cuchara. Si intentas mover la cuchara arriba y abajo, las capas del batido te resisten. Las cosas pesadas quieren quedarse en el fondo, y las cosas ligeras quieren quedarse en la parte superior. El movimiento de "remover" se aplasta.
El Resultado: La mancha vertical creció un poco al principio, pero luego chocó contra un "techo". No pudo hacerse más ancha porque las capas estables del fluido actuaron como una tapa, impidiendo que la turbulencia mezclara las cosas verticalmente. El fluido aún podía girar de lado a lado, pero no podía mezclarse arriba y abajo.

¿Por qué importa esto? (El "Por qué" detrás del "Qué")

Los investigadores descubrieron que en la dirección vertical, el fluido se comporta como un resorte. Una vez que la turbulencia intenta empujar una capa pesada hacia arriba o una capa ligera hacia abajo, la gravedad la devuelve. Esto detiene el movimiento de "remover".

Sin embargo, el fluido aún puede girar de lado a lado. Por lo tanto, la "longitud vertical" de la turbulencia se bloquea en un tamaño específico (determinado por la fuerza de la gravedad y las capas del fluido), y la mancha no puede crecer más allá de ese tamaño.

La "Receta" para la Predicción

El artículo también intentó crear una "receta" matemática simple para predecir cómo se dispersarían estas manchas sin necesidad de una supercomputadora.

• Si conoces la forma de la mancha: Puedes usar una fórmula simple de un solo número para predecir qué tan rápido se dispersa de lado a lado. Funciona muy bien.
• Si no conoces la forma: Tienes que adivinar la forma (asumiendo que parece una curva suave). Si haces esto, necesitas una fórmula de dos números. Esto funciona muy bien después de que la mancha ha tenido tiempo de asentarse en un ritmo con el fluido que gira.

La Conclusión

• Mezcla de lado a lado: Las capas estables (como en el océano profundo o la atmósfera superior) no detienen la mezcla de lado a lado; incluso podrían hacerla un poco más intensa y dentada.
• Mezcla vertical: Las capas estables actúan como un freno. Detienen que el fluido se mezcle arriba y abajo casi por completo.
• La distinción entre "Remover" y "Mezclar": El fluido aún puede "remover" (moverse) de lado a lado, pero no puede "mezclar" (combinar) verticalmente porque las capas resisten ser intercambiadas.

Los autores señalan que su experimento utilizó un tipo específico de propiedad del fluido (número de Prandtl de 0.7). Advierten que si el fluido fuera "más espeso" o tuviera propiedades diferentes (número de Prandtl > 1), los resultados podrían cambiar debido a un efecto "inverso" donde la mezcla crea su propia flotabilidad. Pero para las condiciones que probaron, la regla de "de lado a lado es libre, arriba-abajo está bloqueado" se mantiene verdadera.

Resumen técnico

Resumen Técnico: Evolución de Capas de Mezcla de Escalares Pasivos en Turbulencia Homogénea Estratificada y No Estratificada

Planteamiento del Problema
La dispersión de partículas y escalares en ambientes estratificados estables, como la troposfera superior, la estratosfera y las capas mixtas oceánicas, representa un desafío crítico en las ciencias atmosféricas y oceánicas. Si bien los efectos de la estratificación estable sobre el campo de densidad y la decadencia de la turbulencia están bien documentados, la dispersión de escalares pasivos (por ejemplo, aerosoles, contaminantes, microplásticos) sigue siendo menos comprendida. Específicamente, no está claro cómo evoluciona una capa de mezcla de escalar pasivo (PSML) cuando se introduce en turbulencia homogénea isotrópica (HIT) que posteriormente decae bajo la influencia de un gradiente de densidad estabilizante. La pregunta central es cómo la introducción de escalas de flotabilidad y anisotropía altera la evolución auto-similar de las estadísticas del escalar observada en HIT no estratificada.

Metodología
Los autores emplearon simulaciones numéricas de grandes remolinos (LES) de alta resolución para investigar la decadencia de la turbulencia homogénea tanto no estratificada como estratificada estable. Las simulaciones utilizaron las ecuaciones de Navier-Stokes no divergentes con la aproximación de Boussinesq no hidrostática.

• Configuración del Flujo: Se realizaron dos simulaciones. La primera fue HIT no estratificada. La segunda fue idéntica en las condiciones iniciales, pero incluía un gradiente de densidad estabilizante uniforme ($g \neq 0$), lo que provocó que el flujo transitara hacia una turbulencia fuertemente estratificada después de aproximadamente un periodo de flotabilidad.
• Escalares Pasivos: Se inicializaron dos capas de mezcla de escalares pasivos: una con un gradiente medio en la dirección vertical ($\xi_z$) y otra en la dirección transversal ($\xi_y$). Estas se modelaron como fracciones de mezcla ($0 \le \xi \le 1$).
• Detalles Numéricos: Para minimizar los efectos difusivos que causan desviaciones de la auto-similitud, las simulaciones utilizaron hiperviscosidad e hiperdifusividad con coeficientes iguales. El número de Prandtl molecular se estableció en $Pr = 0.7$. Los números de Reynolds fueron altos ($Re_h > 15,000$ y $Re_\lambda > 450$ para el caso no estratificado), asegurando la existencia de rangos inercial-convectivos e inercial-difusivos. La resolución fue suficiente para resolver la escala de longitud de flotabilidad ($L_b/\Delta \gg 1$).
• Análisis: El estudio rastreó la evolución de las anchuras de los perfiles del escalar, la varianza, momentos de orden superior (asimetría, curtosis), flujos y eficiencia de mezcla. También evaluó modelos de difusividad turbulenta para el flujo de escalar transversal bajo dos condiciones: cuando el perfil medio del escalar es conocido y cuando debe asumirse.

Resultados Clave

1. Mezcla Transversal ($\xi_y$): En la dirección perpendicular a la gravedad, la evolución del escalar pasivo en el caso estratificado es ampliamente similar al caso no estratificado, aunque se dispersa ligeramente más rápido.

   • Las fluctuaciones del escalar son más intensas en el caso estratificado, con una varianza pico aproximadamente un 15% mayor.
   • La interfaz turbulenta/no turbulenta es más intermitente en el caso estratificado.
   • Una vez que el campo escalar alcanza un cuasi-equilibrio con el campo de velocidad, los perfiles de flujo normalizados para ambos casos, estratificado y no estratificado, colapsan, indicando mecanismos de transporte subyacentes similares a pesar de la estratificación.

2. Mezcla Vertical ($\xi_z$): La estratificación inhibe drásticamente la mezcla vertical.

   • La anchura de la capa vertical ($\delta_z$) crece inicialmente, pero deja de expandirse rápidamente después de aproximadamente un periodo de flotabilidad ($t^* \approx 3$).
   • Esta cesación del crecimiento se atribuye a la supresión de la agitación vertical a gran escala por fuerzas de flotabilidad.
   • La anchura final de la capa es proporcional a la escala de longitud integral vertical de la velocidad horizontal ($L_v$), que está restringida para mantener un número de Froude vertical del orden de uno ($Fr_v \sim O(1)$).
   • Mientras que el perfil medio muestra un crecimiento monótono debido a la difusión molecular, el análisis de flujos revela un fenómeno de "flujo inverso" donde la agitación barre ocasionalmente los filamentos del escalar de vuelta hacia su lado original, consistente con un comportamiento similar a ondas internas.

3. Eficiencia de Mezcla: La mezcla acumulada (disipación de la varianza del escalar) en el caso estratificado eventualmente supera a la del caso no estratificado en tiempos tardíos. Esto se atribuye al hecho de que la mezcla sigue a la agitación, y la tasa de agitación horizontal en el flujo estratificado es proporcional a la escala de longitud horizontal, la cual crece más que en el caso no estratificado.

4. Modelado del Flujo:

   • Perfil Conocido: Si el perfil medio del escalar es conocido, un modelo de difusividad turbulenta de una sola constante ($F_y = -c_1 \delta_y v' \partial \xi / \partial y$) predice efectivamente el flujo para ambos casos, estratificado y no estratificado, con $c_1 \approx 0.4$.
   • Perfil Asumido: Si la forma del perfil (asumida como una función de error) y la escala de longitud del escalar deben estimarse a partir de estadísticas de velocidad, se requiere un modelo de dos constantes. Este modelo funciona bien en tiempos tardíos cuando el escalar está en cuasi-equilibrio con el campo de velocidad, permitiendo que la escala de longitud del escalar se escale a partir de la tasa de disipación de energía cinética.

Significado y Afirmaciones
El artículo afirma proporcionar una caracterización detallada de la mezcla de escalares pasivos en turbulencia estratificada en decadencia, cerrando la brecha entre las predicciones teóricas de difusión anisotrópica y las observaciones numéricas.

• Anisotropía: El estudio confirma que la estratificación crea una dicotomía marcada en la mezcla: la mezcla horizontal permanece activa e incluso se ve potenciada, mientras que la mezcla vertical se suprime efectivamente una vez que el flujo se vuelve fuertemente estratificado.
• Utilidad de Modelado: Los autores demuestran que los enfoques estándar de difusividad turbulenta siguen siendo viables para flujos estratificados, siempre que el escalar esté en cuasi-equilibrio con el campo de velocidad. Destacan que asumir un perfil conocido simplifica significativamente los requisitos de modelado en comparación con asumir una forma de perfil.
• Limitaciones y Trabajo Futuro: Los autores notan modestamente que sus resultados son específicos para un número de Prandtl de 0.7. Declaran explícitamente que la literatura reciente sugiere que números de Prandtl más altos ($Pr > 1$) inducen flujos de flotabilidad inversos a pequeñas escalas, lo que podría alterar las tasas de decadencia del flujo y, en consecuencia, las tasas de mezcla de escalares pasivos. Por lo tanto, los hallazgos actuales no se aplican directamente a escenarios de alto número de Prandtl sin una investigación adicional.

El estudio concluye que, si bien la estratificación altera fundamentalmente la dinámica vertical de una capa de mezcla, la mezcla transversal conserva muchas características de la turbulencia no estratificada, aunque con mayor intermitencia e intensidad de fluctuación.