Two pathways to diapycnal mixing in strongly stratified flows with no initial vertical shear

Este estudio combina la teoría lineal y las simulaciones numéricas directas para revelar que en flujos fuertemente estratificados sin cizalladura vertical inicial, las inestabilidades de cizalladura horizontal impulsan inevitablemente la mezcla diapicnal a través de dos vías distintas —ya sea mediante la emergencia directa de cizalladura vertical o a través de una evolución no lineal hacia vórtices columnares—, las cuales finalmente desencadenan inestabilidades de Kelvin-Helmholtz a pequeña escala, pero producen diferentes eficiencias de mezcla debido a su excitación de distintas escalas verticales.

Lo esencial

Imagina el océano o la atmósfera como un pastel gigante y estratificado. Las capas están hechas de fluidos con diferentes densidades (como diferentes sabores de pastel) y no les gusta mezclarse fácilmente. Normalmente, los científicos estudian qué sucede cuando empujas estas capas lateralmente y hacia arriba o hacia abajo al mismo tiempo (cizalladura vertical). Pero este artículo hace una pregunta diferente: ¿Qué sucede si solo empujas las capas lateralmente, sin un movimiento inicial de arriba hacia abajo, en un entorno fuertemente estratificado?

Los investigadores descubrieron que, incluso si comienzas con un flujo horizontal perfectamente plano, la naturaleza tiene dos "recetas" o vías distintas para crear caos (turbulencia) y mezclar las capas. La vía que la naturaleza elige depende enteramente de cómo "siembres" el experimento; esencialmente, qué pequeño empujoncito le des al fluido al principio.

Aquí está el desglose de las dos vías utilizando analogías simples:

La configuración: El río tranquilo

Imagina un río ancho y tranquilo fluyendo horizontalmente. El agua está estratificada en capas (como una pila de panqueques). Al principio, el flujo es suave y bidimensional (solo se mueve de izquierda a derecha, no de arriba hacia abajo).

Vía 1: El efecto de la "Habitación abarrotada" (La ruta directa)

Cómo comienza: Le das al río un pequeño empujón aleatorio en todas partes a la vez (como lanzar un puñado de confeti al aire).
Qué sucede:

1. La ondulación: Debido a las capas, el fluido no solo ondula de izquierda a derecha; inmediatamente empieza a ondular hacia arriba y hacia abajo en muchos tamaños diferentes al mismo tiempo. Piensa en esto como una multitud de personas en una habitación que intentan moverse todas a la vez, creando un forcejeo caótico y multidireccional.
2. La cizalladura: Estas ondulaciones crean fuertes corrientes verticales (cizalladura) muy rápidamente.
3. La ruptura: Estas corrientes verticales se vuelven tan fuertes que se rompen, creando pequeños vórtices violentos (como diminutos torbellinos). Esta es la inestabilidad de "Kelvin-Helmholtz", que se ve como las olas que se rompen cuando el viento sopla sobre el agua.
   El resultado: La mezcla ocurre de manera eficiente. Debido a que la energía se distribuye en muchas escalas de ondulación, la "fricción" (disipación viscosa) es menor, lo que hace que el proceso de mezcla sea relativamente eficiente.

Vía la 2: El efecto de la "Danza sincronizada" (La ruta indirecta)

Cómo comienza: Le das al río un empujón muy específico y organizado (como un director de orquesta agitando una batuta para que todos se muevan con un patrón específico).
Qué sucede:

1. El vórtice: En lugar de ondulaciones caóticas, el fluido se organiza en largas columnas verticales de agua giratoria (como gigantescos tornados parados en el río). Durante mucho tiempo, el flujo permanece perfectamente bidimensional, compuesto solo por estas grandes columnas giratorias.
2. El bamboleo: Eventualmente, estas gigantescas columnas se vuelven inestables. Comienzan a bambolearse de una manera muy específica y de alta frecuencia. Los investigadores llaman a esto una "inestabilidad hiperbólica". Imagina un trompo que comienza a bambolearse violentamente justo antes de caerse.
3. La ruptura: Este bamboleo violento crea capas verticales de cizalladura muy delgadas y agudas. Estas capas delgadas luego se rompen en pequeños vórtices violentos, al igual que en la Vía 1.
   El resultado: La mezcla ocurre, pero es menos eficiente. ¿Por qué? Porque esta vía crea capas extremadamente delgadas y agudas. Requiere mucha energía (fricción) crear y romper estas capas tan pequeñas y afiladas. Es como intentar cortar un bloque grueso de queso con un cuchillo sin filo (Vía 1) frente a una hoja de afeitar (Vía 2); la hoja de afeitar crea un corte mucho más nítido y con mayor gasto de energía.

La gran conclusión

El artículo demuestra que la cizalladura vertical (movimiento de arriba hacia abajo) no necesita estar presente desde el principio. Es un subproducto inevitable de la cizalladura horizontal en fluidos fuertemente estratificados, siempre que el fluido sea lo suficientemente espeso (alto número de Reynolds).

• Si comienzas con ruido aleatorio: Obtienes la Vía 1 (Mezcla directa y eficiente).
• Si comienzas con un patrón específico: Obtienes la Vía 2 (Mezcla indirecta y menos eficiente).

Los investigadores utilizaron potentes simulaciones por computadora para demostrar que estas dos vías son reales, distintas, y que la "receta" que eliges al principio determina cuánta energía se desperdicia como calor frente a cuánta se utiliza para realmente mezclar las capas.

En resumen: Incluso en un fluido estratificado perfectamente tranquilo, un empuje horizontal eventualmente creará caos vertical. Pero dependiendo de cómo comiences el empuje, ese caos se verá diferente y mezclará las capas con diferentes niveles de eficiencia.

Resumen técnico

Planteamiento del Problema
La estabilidad de flujos horizontales a gran escala y fuertemente estratificados ante perturbaciones tridimensionales (3D) es un mecanismo crítico para la transferencia descendente de energía y el transporte diapicnal en sistemas geofísicos y astrofísicos. Mientras que los flujos estratificados con cizalla vertical están bien estudiados, sus contrapartes con cizalla horizontal (sin cizalla vertical inicial) han recibido menos atención. Una cuestión central persiste: ¿cómo transitan los flujos de cizalla horizontal, inicialmente bidimensionales (2D) y fuertmente estratificados ($Fr \ll 1$), hacia la turbulencia y generan cizalla vertical? La literatura previa sugiere dos vías distintas: una que involucra el crecimiento directo de modos primarios 3D que conducen a inestabilidades de Kelvin-Helmholtz (KH), y otra que involucra la evolución no lineal hacia vórtices columnares seguidos de inestabilidades secundarias. Sin embargo, no ha quedado claro si estas vías son mutuamente excluyentes debido a diferencias en las condiciones iniciales (por ejemplo, tangente hiperbólica frente a perfiles sinusoidales) o en la formulación del problema (Problema de Valor Inicial frente a forzado por cuerpo), o si representan rutas alternativas disponibles dentro del mismo sistema físico. Además, el impacto de la vía seleccionada en la eficiencia de la mezcla diapicnal no ha sido cuantificado sistemáticamente.

Metodología
Los autores combinan la teoría de estabilidad lineal con Simulaciones Numéricas Directas (DNS) para investigar estas transiciones en un dominio triperiodico. El estudio se centra en un flujo de Kolmogorov plano, horizontal y verticalmente invariante ($\bar{u}(y) = \sin(y)\mathbf{e}_x$) en un fluido linealmente estratificado con frecuencia de flotabilidad constante $N$. Las simulaciones se realizan a números de Reynolds ($Re = 10,000$) y Péclet ($Pe = 10,000$) elevados, con números de Froude bajos ($Fr = 0.15$y$Fr = 0.1$), lo que resulta en números de Reynolds de flotabilidad ($Re_b = Fr^2 Re$) de 225 y 100, respectivamente.

Para aislar el efecto de las condiciones iniciales en la selección de la vía, se realizan dos conjuntos de DNS:

1. DNS1 (Vía P1): Inicializado con ruido blanco a escala de malla, permitiendo que todos los modos inestables de la inestabilidad de cizalla horizontal primaria crezcan.
2. DNS2 (Vía P2): Inicializado con ruido a escala de malla más una "semilla" específica para el modo propio 2D ($k_z=0$) de más rápido crecimiento de la inestabilidad primaria.

Adicionalmente, los autores realizan un análisis de estabilidad lineal de un modelo analítico "congelado" del flujo columnar que emerge en la DNS2. Este modelo aproxima el flujo 2D evolucionado no linealmente como un arreglo de vórtices de Taylor-Green superpuesto con flujos sinusoidales, permitiendo el cálculo de los modos propios y las tasas de crecimiento de las perturbaciones 3D secundarias.

Contribuciones Clave y Resultados

1. Confirmación de Dos Vías en un Mismo Configuración:
   El estudio demuestra que ambas vías hacia la turbulencia pueden ocurrir dentro de la misma configuración del modelo, distinguidas únicamente por el espectro de perturbación inicial:

• Vía P1 (Crecimiento 3D Directo): Cuando se inicializa con ruido blanco, los modos propios de $k_z \neq 0$ de la inestabilidad de cizalla horizontal primaria crecen directamente. Estos modos generan cizalla vertical en un amplio rango de escalas, la cual se vuelve posteriormente inestable a inestabilidades de KH de pequeña escala cuando el $Re_b$ es suficientemente grande. Esta vía se alinea con las observaciones en flujos de Kolmogorov forzados por cuerpo.
• Vía P2 (Evolución Columnar): Cuando se siembra el modo $k_z=0$, este domina el flujo, impulsando el fondo hacia un flujo columnar (vortical) 2D de larga duración y dependiente del tiempo. Este estado 2D es posteriormente inestable a modos 3D secundarios. Los autores identifican estos modos secundarios como inestabilidades hiperbólicas que surgen de los puntos de estancamiento hiperbólicos en el arreglo de vórtices. Estas inestabilidades poseen altos números de onda verticales ($k_z \sim Fr^{-1}$) e impulsan la emergencia de cizalla vertical, la cual eventualmente desencadena inestabilidades de KH terciarias.

2. Caracterización de la Inestabilidad Secundaria en la Vía P2:
   A través del análisis de estabilidad lineal del flujo columnar congelado, los autores muestran que la inestabilidad secundaria es una inestabilidad hiperbólica de alto $k_z$. La tasa de crecimiento de esta inestabilidad escala con la tasa de deformación local de los puntos hiperbólicos ($\Delta_h$) y es aproximadamente $\lambda \approx 0.46 \Delta_h$, siendo ampliamente independiente de la fuerza de la estratificación para amplitudes grandes. La inestabilidad se manifiesta como "ondulaciones" en el campo de velocidad vertical centradas alrededor de los puntos hiperbólicos, distintas de la inestabilidad de "zigzag" observada tradicionalmente en pares de vórtices aislados.

3. Impacto en la Eficiencia de Mezcla:
   El estudio cuantifica la mezcla diapicnal y la eficiencia de mezcla ($\eta$) para ambas vías.

• Disipación de Flotabilidad ($\chi$): La tasa máxima de disipación de la varianza de la flotabilidad es comparable en ambas vías, ya que está controlada por movimientos de pequeña escala ($k_z \gtrsim Fr^{-1}$) en ambos casos.
• Disipación Viscosa ($\epsilon$): La Vía P2 exhibe una disipación viscosa significativamente mayor que la Vía P1. Esto se debe a que la inestabilidad hiperbólica en P2 excita un espectro estrecho de alto $k_z$, lo que conduce a una disipación intensa en capas de cizalla vertical delgadas. En contraste, la P1 excita un espectro amplio de escalas verticales, distribuyendo la disipación de manera más extensa.
• Eficiencia de Mezcla: Consecuentemente, la eficiencia de mezcla máxima es significativamente mayor en la Vía P1 ($\eta_{max} \approx 0.25$) que en la Vía P2 ($\eta_{max} \approx 0.165$). La diferencia radica en que la P2 disipa más energía cinética en relación con la cantidad de flotabilidad mezclada.

Significancia y Reivindicaciones
El artículo sostiene que la emergencia de la cizalla vertical que impulsa las inestabilidades de KH de pequeña escala es un subproducto inevitable de las inestabilidades de cizalla horizontal en flujos fuertemente estratificados a un $Re_b$ suficientemente grande, siempre que el flujo posea un punto de inflexión. La selección entre las dos vías está determinada por las condiciones iniciales: específicamente, la relación entre la amplitud inicial del modo $k_z=0$ respecto a los modos $k_z \neq 0$.

Los autores sugieren que la prevalencia de la Vía P2 en la literatura previa (específicamente en estudios de capas de cizalla de tangente hiperbólica) puede ser un artefacto de las prácticas de inicialización donde a menudo se añade una semilla $k_z=0$ para acelerar los cálculos, suprimiendo así el crecimiento directo de los modos 3D (Vía P1). Por el contrario, la ausencia de la P1 en algunos contextos puede deberse a límites de estratificación más débiles donde el rango de modos inestables $k_z \neq 0$ es limitado.

En última instancia, este trabajo establece que las inestabilidades de cizalla horizontal son actores clave en la generación de turbulencia estratificada y la mezcla diapicnal, incluso en ausencia de cizalla vertical preexistente. Destaca que la ruta específica hacia la turbulencia (crecimiento 3D directo frente a la transición 2D-a-3D) altera fundamentalmente el espectro de disipación de energía y, por consiguiente, la eficiencia de la mezcla, un factor que debe considerarse al parametrizar la turbulencia en modelos climáticos y geofísicos. Los autores señalan que, aunque sus resultados son robustos para flujos idealizados, la interacción con la rotación y los campos de ondas ambientales en entornos geofísicos realistas sigue siendo una cuestión abierta.