A theoretical one-dimensional model for variable-density Rayleigh-Taylor turbulence

Este artículo revisa y amplía el modelo teórico de 1965 de Belen'kii y Fradkin para la turbulencia de Rayleigh-Taylor de densidad variable, demostrando que la ecuación de similitud completa captura con precisión las características del flujo no Boussinesq y que una solución simplificada corregida por masa aproxima eficazmente la compleja dinámica de mezcla regida por la competencia entre la difusión y la conservación de la masa.

Lo esencial

Imagina dos fluidos, uno pesado (como la miel) y otro ligero (como el aire), situados uno encima del otro. La gravedad desea que el pesado se hunda y el ligero ascienda, pero están atrapados en una batalla desordenada y agitada en la interfaz. Esto es la inestabilidad de Rayleigh-Taylor. A medida que se mezclan, forman una "sopa" turbulenta donde picos pesados se zambullen hacia abajo y burbujas ligeras flotan hacia arriba.

Durante décadas, los científicos han intentado predecir a qué velocidad crece esta capa de mezcla. La mayoría de las teorías modernas asumen que los fluidos tienen una densidad "casi" igual, utilizando una regla empírica simple. Sin embargo, este artículo revisa una teoría olvidada de hace 60 años, de 1965, por Belen'kii y Fradkin, que ofrece una forma diferente y más precisa de observar este caos, especialmente cuando la diferencia de densidad es enorme.

Aquí está el desglose de lo que hace el artículo, utilizando analogías simples:

1. La Receta Olvidada

Los autores encontraron una antigua "receta" (un modelo matemático) para cómo se mezclan estos fluidos. La receta original fue escrita en ruso, era un poco confusa de leer y tenía algunos errores tipográficos.

• Lo que hicieron: Limpian la receta, la tradujeron y la reescribieron utilizando un lenguaje moderno y claro.
• La Idea Central: En lugar de pensar en la mezcla como una explosión compleja en 3D, la trataron como un problema de difusión unidimensional. Imagina la capa de mezcla no como una tormenta caótica, sino como una única mancha que se extiende sobre un papel. Modelaron esta extensión de la "mancha" utilizando un concepto llamado difusividad turbulenta (qué tan rápido se extiende el caos).

2. La Regla del "Logaritmo" vs. la Regla "Lineal"

El gran descubrimiento en este artículo es sobre cómo crece la capa de mezcla con el tiempo.

• La Vieja Visión: La mayoría de los científicos pensaban que la tasa de crecimiento dependía de un número lineal llamado número de Atwood (que mide la diferencia entre el fluido pesado y el ligero). Si la diferencia se duplica, la velocidad de mezcla se duplica.
• La Nueva (Antigua) Visión: El modelo de 1965 sugiere que el crecimiento depende del logaritmo natural de la relación de densidades ($\ln R$).
  • La Analogía: Piensa en el número de Atwood como una línea recta en una gráfica. El logaritmo es como una curva que se aplana. El artículo argumenta que cuando la diferencia de densidad es enorme (como comparar plomo con aire), la mezcla no acelera linealmente; su tasa de crecimiento se ralentiza, siguiendo esta curva logarítmica. Esto coincide mejor con las simulaciones por computadora recientes que la vieja regla lineal.

3. La Asimetría "Pesada" y "Ligera"

Cuando los fluidos pesados y ligeros se mezclan, no se comportan de la misma manera.

• La Observación: El fluido pesado forma "picos" que se zambullen hacia abajo rápidamente, mientras que el fluido ligero forma "burbujas" que ascienden más lentamente.
• La Perspectiva del Artículo: El antiguo modelo de 1965 predice naturalmente esta asimetría sin necesidad de ajustes adicionales. Muestra que a medida que crece la diferencia de densidad, los "picos" se vuelven mucho más largos que las "burbujas".
• El Desplazamiento de Velocidad: El artículo también muestra que la velocidad de la mezcla se desplaza hacia el lado del fluido ligero.
  • La Analogía: Imagina un juego de tira y afloja donde un equipo es mucho más pesado. La cuerda no se mueve simplemente hacia el medio; todo el centro de la acción se desplaza hacia el equipo más ligero. El modelo captura perfectamente este "desplazamiento".

4. El Truco de la "Corrección de Masa"

El modelo original de 1965 tenía una versión simplificada que era fácil de resolver pero tenía un defecto: violaba la ley de conservación de la masa.

• El Problema: Si solo usas las matemáticas simples, es como un globo que gana o pierde aire mágicamente mientras se expande. La cantidad total de "cosa" (masa) no suma correctamente.
• La Solución: Los autores se dieron cuenta de que si simplemente desplazas todo el perfil de mezcla ligeramente hacia el lado del fluido ligero, las matemáticas de repente funcionan perfectamente.
  • La Analogía: Imagina una colina de arena perfectamente simétrica (el modelo simplificado). Se ve bien, pero si pesas la arena, le falta un poco. Si deslizas toda la colina unos centímetros hacia la izquierda, el peso se equilibra y de repente se ve exactamente como los datos desordenados del mundo real.
  • Este "desplazamiento" explica por qué los picos crecen más rápido que las burbujas: la difusión del "logaritmo de la densidad" es simétrica, pero la necesidad de conservar la masa obliga a toda la estructura a inclinarse hacia el lado ligero.

5. La Conclusión

El artículo concluye que este modelo simple y unidimensional de 1965 es en realidad una "mina de oro".

• Captura todos los comportamientos extraños y complejos de la mezcla de alta densidad (asimetría, desplazamiento de velocidades, crecimiento logarítmico) que los científicos modernos solo han confirmado recientemente con supercomputadoras.
• Sugiere que la física de esta turbulencia está gobernada por una competencia entre la difusión (extenderse) y la conservación de la masa (mantener la misma cantidad total de fluido).

En resumen: Los autores desenterraron una teoría antigua y polvorienta, la limpiaron y demostraron que explica las observaciones modernas de la mezcla de fluidos mejor que muchas teorías actuales. Probaron que un simple "desplazamiento" en las matemáticas corrige los errores del modelo antiguo y describe perfectamente por qué los fluidos pesados se zambullen más rápido de lo que los fluidos ligeros ascienden cuando son muy diferentes en densidad.

Resumen técnico

A continuación se presenta un resumen técnico detallado del artículo "Un modelo teórico unidimensional para la turbulencia de Rayleigh–Taylor de densidad variable" de Chian Yeh Goh y Guillaume Blanquart.

1. Planteamiento del Problema

La inestabilidad de Rayleigh–Taylor (RT) ocurre cuando un fluido pesado es acelerado hacia un fluido más ligero, lo que conduce a una mezcla turbulenta. Si bien está bien establecido que el crecimiento tardío de la altura de la capa de mezcla ($h$) escala como $h \sim g t^2$, la dependencia de la relación de densidades ($R = \rho_H / \rho_L$) sigue siendo objeto de debate en flujos de densidad variable (no Boussinesq).

• Consenso Actual: La mayoría de los estudios modernos interpretan los resultados utilizando la aproximación de Boussinesq, donde la tasa de crecimiento escala linealmente con el número de Atwood ($A = (R-1)/(R+1)$), es decir, $h \approx \alpha A g t^2$.
• La Discrepancia: Los datos experimentales y numéricos muestran que el crecimiento de las burbujas es relativamente insensible a $A$, mientras que el crecimiento de las espigas aumenta significativamente con $A$, lo que genera asimetrías en los perfiles de velocidad y densidad.
• El Vacío: Un trabajo seminal de 1965 de Belen'kii y Fradkin propuso un modelo de difusividad turbulenta que arrojaba una escala de $h \propto (\ln R) g t^2$. Sin embargo, este trabajo fue publicado en ruso, contenía errores tipográficos y fue en gran medida pasado por alto. Además, el análisis original se basaba en una Ecuación Diferencial Ordinaria (ODE) simplificada sin explorar completamente las propiedades de la ecuación de similitud completa ni compararla con datos modernos de alta fidelidad.

2. Metodología

Los autores revisan, aclaran y extienden el modelo de Belen'kii y Fradkin (1965) utilizando un marco riguroso unidimensional (1D).

Marco Matemático:

• Ecuaciones Gobernantes: Los autores comienzan con las ecuaciones de Navier-Stokes promediadas por Reynolds (RANS) para flujo de densidad variable. Derivan una ecuación de densidad media análoga a un problema de difusión 1D:
  $$\frac{\partial \bar{\rho}}{\partial t} = \frac{\partial}{\partial y} \left( D_t \frac{\partial \bar{\rho}}{\partial y} \right)$$
• Modelo de Difusividad Turbulenta: Derivan una forma específica para la difusividad turbulenta ($D_t$) basada en la teoría de la longitud de mezcla de Prandtl y argumentos de escalado de energía. El modelo postula:
  $$D_t^* = K \left( \frac{1}{\bar{\rho}} \frac{\partial \bar{\rho}}{\partial y} \right)^{1/2} g^{1/2} h_t^2$$
  donde $h_t$ es la escala de altura de la turbulencia.
• Solución Auto-similar: Introduciendo variables de similitud, la ecuación en derivadas parciales se reduce a una ODE no lineal.
  • ODE Completa: Una ODE no lineal de tercer orden (reducida a una ODE de primer orden para una variable transformada $x$) que captura la física completa.
  • ODE Simplificada: Una aproximación analítica derivada al despreciar un término de orden superior ($x^3$), válida para relaciones de densidad pequeñas ($R \lesssim 4$).

Estrategia de Análisis:

• Solución Numérica: La ODE completa se resuelve numéricamente para un amplio rango de números de Atwood ($A \in [0, 0.8]$).
• Solución Analítica: La ODE simplificada se resuelve analíticamente para recuperar la escala original de $\ln R$.
• Corrección de Masa: Los autores identifican que la ODE simplificada viola la conservación global de masa debido a las suposiciones de simetría. Proponen una "corrección de masa" (desplazamiento espacial de la solución) para alinear el modelo simplificado con la realidad física.
• Validación: Los perfiles teóricos (densidad, difusividad, tasas de crecimiento) se comparan con datos de Simulación Numérica Directa (DNS) y resultados experimentales de diversas fuentes bibliográficas.

3. Contribuciones Clave

1. Resurrección y Aclaración: El artículo proporciona una derivación moderna y limpia del modelo de 1965 de Belen'kii y Fradkin, corrigiendo la notación y aclarando suposiciones que anteriormente estaban oscurecidas.
2. Análisis de la ODE Completa: A diferencia del trabajo original, este artículo resuelve la ecuación de similitud completa (no solo la aproximación). Demuestra que el modelo completo captura inherentemente características complejas no Boussinesq sin ajustes ad hoc.
3. Mecanismo de Asimetría: Los autores proporcionan una explicación mecánica de la asimetría observada entre espigas y burbujas. Muestran que la variable natural para la mezcla turbulenta en flujos de densidad variable es $\ln \bar{\rho}$. La difusión de $\ln \bar{\rho}$ es simétrica, pero el mapeo exponencial de vuelta a la densidad ($\bar{\rho}$) combinado con la restricción de conservación de masa fuerza un desplazamiento espacial hacia el lado del fluido ligero.
4. Modelo Simplificado Corregido por Masa: Demuestran que la solución analítica simple (difusividad parabólica) puede reproducir con precisión la solución completa si se aplica una corrección global de masa (desplazamiento espacial). Esto ofrece una aproximación de forma cerrada y computacionalmente económica para el modelado de turbulencia.
5. Validación de la Ley de Escalamiento: El estudio valida la escala de $\ln R$ para la altura de la capa de mezcla, mostrando que actúa como una estimación "mediana" que se ajusta mejor a la dispersión de los datos experimentales y numéricos existentes que la escala lineal de $A$ en relaciones de densidad altas.

4. Resultados Clave

• Escalamiento de la Altura de la Capa de Mezcla: La altura total de la capa de mezcla escala como $h \propto (\ln R) g t^2$. Aunque se producen desviaciones en $A$ extremas, esta escala permanece consistente con la dispersión en los datos DNS y experimentales disponibles hasta $A \approx 0.84$.
• Asimetría Espiga vs. Burbuja:
  • El modelo predice que las alturas de las espigas ($h_s$) crecen más rápido que las alturas de las burbujas ($h_b$) a medida que aumenta $A$.
  • La relación $h_s/h_b$ aumenta con $A$, lo cual es consistente con las observaciones empíricas.
  • El parámetro de crecimiento $\alpha_s$ aumenta con $A$, mientras que $\alpha_b$ permanece relativamente constante (o se desvía menos), lo que explica por qué la aproximación de Boussinesq falla para las espigas en $A$ altos.
• Desplazamientos de Perfiles:
  • Velocidad/Difusividad: Los perfiles de difusividad (y velocidad) turbulenta se desplazan sistemáticamente hacia el lado del fluido ligero a medida que aumenta $A$.
  • Densidad: Los perfiles de fracción molar colapsan bien cuando se normalizan, pero los límites físicos (puntas de espiga/burbuja) se vuelven asimétricos ($\lambda_s > \lambda_b$).
• Comparación con Datos: Los perfiles teóricos para la difusividad normalizada y la fracción molar muestran un acuerdo razonable con los datos DNS (por ejemplo, de Goh et al., Livescu et al.), particularmente en el núcleo de la capa de mezcla. Las discrepancias se confinan principalmente a los bordes donde la difusión molecular (despreciada en el modelo) se vuelve relevante.

5. Significado

• Perspectiva Teórica: El artículo resuelve una ambigüedad de larga data en la turbulencia RT al mostrar que las dinámicas competitivas entre la difusión de $\ln \bar{\rho}$ (que es simétrica) y la conservación de masa (que requiere un desplazamiento) son los impulsores fundamentales de la asimetría no Boussinesq.
• Utilidad Práctica: La solución simplificada "corregida por masa" proporciona un conjunto de expresiones de forma cerrada para los perfiles de densidad y difusividad. Esto es altamente valioso para el desarrollo de modelos de turbulencia promediados por Reynolds (RANS), donde la eficiencia computacional es crítica.
• Reevaluación del Escalamiento: El trabajo desafía la dominancia del escalamiento lineal del número de Atwood ($A$) en favor del escalamiento logarítmico de la relación de densidades ($\ln R$) para flujos RT de densidad variable, sugiriendo que este último es una base teórica más robusta para aplicaciones de alta relación de densidad (por ejemplo, fusión por confinamiento inercial, astrofísica).

En resumen, Goh y Blanquart demuestran que un marco teórico 1D simple, propuesto originalmente en 1965, contiene la física esencial para explicar las observaciones modernas de la turbulencia RT de densidad variable, siempre que el modelo se resuelva completamente o se corrija para la conservación de masa.