Eddy thermal diffusivity model and mean temperature profiles in turbulent vertical convection

Este artículo propone un modelo de difusividad térmica turbulenta dependiente del espacio que permite derivar perfiles analíticos de temperatura media en convección natural vertical, los cuales muestran dos funciones de escalado universales en las regiones interna y externa que concuerdan bien con datos de simulación numérica directa.

Lo esencial

¡Hola! Imagina que tienes dos paredes verticales infinitas, como los muros de un edificio muy alto. Una pared está muy caliente (como un radiador encendido) y la otra muy fría (como un bloque de hielo). Entre ellas hay un fluido, como el aire o el agua.

El problema que estudia este artículo es: ¿Cómo se mueve el calor y el fluido en ese espacio?

Cuando el aire caliente sube y el frío baja, se crea un movimiento caótico y turbulento, como una olla de agua hirviendo pero vertical. Los científicos quieren predecir exactamente cómo cambia la temperatura desde la pared caliente hasta la fría.

Aquí te explico lo que hicieron estos investigadores (Ng y Ching) usando analogías sencillas:

1. El problema de la "caja negra"

Antes de este estudio, los científicos tenían dos formas de ver este problema:

• La teoría vieja: Decía que la temperatura cambiaba de una manera muy específica (como una curva suave) cerca de las paredes y otra manera en el centro. Pero, al compararlo con simulaciones por computadora muy precisas (que son como "laboratorios virtuales"), se veía que la teoría vieja fallaba, especialmente cuando el fluido era más viscoso (como la miel) o menos (como el agua).
• El misterio: Nadie tenía una fórmula matemática simple que funcionara para todos los tipos de fluidos y temperaturas.

2. La nueva idea: Un "mapa de tráfico" para el calor

Los autores proponen una nueva forma de ver el problema. Imagina que el calor no viaja solo, sino que es transportado por remolinos de fluido (como pequeños remolinos de viento). A esto lo llaman "difusividad turbulenta".

Su gran invento es un modelo de tres capas para describir cómo viaja el calor, como si fuera una carretera con tres zonas distintas:

• Zona 1: La "Zona de Cuidado" (Junto a la pared caliente/fría).
  Aquí, el fluido está pegado a la pared y no se mueve mucho. Es como caminar por un pasillo muy estrecho y lleno de gente. El calor se mueve lento y depende mucho de las propiedades del fluido (como su "viscosidad" o resistencia). Los autores dicen que aquí el calor se comporta de una manera matemática muy específica (una curva cúbica).
• Zona 2: La "Zona de Transición" (El medio).
  Es como una autopista de entrada donde el tráfico empieza a acelerar. Aquí conectan la zona lenta con la zona rápida.
• Zona 3: La "Zona de Alta Velocidad" (El centro).
  En el medio, entre las dos paredes, el fluido se mueve muy rápido y de forma muy caótica. Aquí, el calor viaja a toda velocidad, como un coche de Fórmula 1 en una recta. En esta zona, las propiedades moleculares del fluido importan menos; lo que cuenta es la turbulencia.

3. La "Receta Universal"

Lo genial de su trabajo es que, al usar este modelo de tres capas, descubrieron que todos los datos (desde el aire hasta aceites muy espesos) encajan en dos fórmulas universales:

1. Una fórmula para la pared (Zona 1).
2. Otra fórmula para el centro (Zona 3).

Es como si, sin importar si estás conduciendo un camión o una moto, en la ciudad (cerca de la pared) sigues las mismas reglas de tráfico lentas, y en la autopista (el centro) sigues las reglas de alta velocidad. Sus fórmulas capturan perfectamente ambas situaciones.

4. ¿Por qué es importante?

Antes, si querías diseñar un sistema de ventilación para un rascacielos o entender cómo el hielo del Ártico se derrite en el océano, tenías que adivinar o usar modelos que fallaban.

• La analogía del "GPS": Los modelos anteriores eran como un GPS antiguo que te decía "gira a la derecha" cuando en realidad había un muro. El nuevo modelo de Ng y Ching es un GPS moderno y preciso que te dice exactamente cómo va la temperatura en cada metro del camino.
• Validación: Probaron su teoría contra millones de datos de simulaciones por computadora y funcionó perfecto. Además, demostraron que los modelos anteriores (que sugerían curvas logarítmicas o raíces cúbicas simples) no eran correctos para todos los casos.

En resumen

Estos científicos crearon un mapa matemático mejor para entender cómo se mueve el calor en fluidos turbulentos entre dos paredes.

• Dividieron el problema en una zona lenta (cerca de la pared) y una zona rápida (en el centro).
• Encontraron dos reglas universales que explican la temperatura en todo el sistema, sin importar qué tan caliente o viscoso sea el fluido.
• Su modelo es más preciso que los anteriores, lo que ayudará a ingenieros a diseñar edificios más eficientes y a climatólogos a entender mejor el cambio climático en los polos.

Es un trabajo elegante que toma un fenómeno caótico y turbulento y le pone orden con una fórmula matemática que funciona para casi todo.

Resumen técnico

A continuación presento un resumen técnico detallado del artículo "Eddy thermal diffusivity model and mean temperature profiles in turbulent vertical convection" (Modelo de difusividad térmica turbulenta y perfiles de temperatura media en convección vertical turbulenta), publicado en el Journal of Fluid Mechanics.

1. Planteamiento del Problema

El estudio se centra en la convección natural turbulenta de un fluido confinado entre dos paredes verticales infinitas mantenidas a diferentes temperaturas ($T_h$ y $T_c$). Este flujo es fundamental en aplicaciones de ingeniería (ventilación de edificios) y fenómenos naturales (interacción hielo-océano).

• Parámetros de control: El estado del flujo está determinado por el número de Rayleigh ($Ra$) y el número de Prandtl ($Pr$).
• Desafío actual: Aunque la convección laminar está bien entendida mediante ecuaciones de capa límite, la convección turbulenta carece de una comprensión completa.
• Limitaciones de modelos previos:
  • Modelos anteriores (George & Capp, 1979; Hölling & Herwig, 2005) propusieron perfiles de temperatura con dependencia cúbica inversa o logarítmica en una capa de superposición, pero estos resultados no se han validado para un rango general de $Pr$ y a menudo violan las condiciones de contorno en las paredes.
  • Estudios recientes (Li et al., 2023) utilizaron modelos de flujo de calor turbulento que no cumplen con las condiciones de anulación de derivadas en la pared ($x=0$), lo que los hace físicamente incorrectos cerca de los límites.

2. Metodología

Los autores proponen un modelo de difusividad térmica turbulenta (o eddy diffusivity) dependiente del espacio, denotado como $K(x)$, para cerrar las ecuaciones de transporte de momento y energía promediadas en el tiempo.

• Ecuaciones gobernantes: Se utilizan las ecuaciones de Navier-Stokes bajo la aproximación de Oberbeck-Boussinesq. Tras promediar en el tiempo, se obtienen ecuaciones para la velocidad media vertical $w(x)$ y la temperatura media $T(x)$.
• Modelo de Difusividad ($K(x)$):
  • Se propone un modelo de tres capas para $K(x)$ basado en el comportamiento físico esperado:
    1. Región interior (cerca de la pared): Se modela con una función cúbica ($Ay^3$) para satisfacer las condiciones de contorno donde el flujo de calor turbulento y sus derivadas deben anularse en la pared.
    2. Región intermedia: Se utiliza una función lineal para conectar suavemente las regiones interior y exterior.
    3. Región exterior (cerca de la línea central): Se modela con una función cuadrática simétrica que alcanza un máximo en el centro ($x=H/2$), aprovechando la simetría del problema.
  • El modelo depende de dos parámetros independientes: $A$ (relacionado con la región interior) y $C_m$ (relacionado con la región exterior).
• Derivación Analítica:
  • Integrando la ecuación de energía con el modelo de $K(x)$, los autores derivan expresiones analíticas para el perfil de temperatura media $T(x)$ y el número de Nusselt ($Nu$).
  • Se identifican escalas de temperatura y longitud características para las regiones interior ($T_i, l_i$) y exterior ($T_o, H$).

3. Contribuciones Clave

1. Modelo Cerrado y Universal: Se presenta un modelo analítico cerrado que describe los perfiles de temperatura media mediante dos funciones de escala universales ($F_i$ para la región interior y $F_o$ para la exterior), válidas para diferentes números de Rayleigh y Prandtl.
2. Resolución de Inconsistencias: A diferencia de modelos previos, este enfoque respeta estrictamente las condiciones de contorno en las paredes (anulación de $u'T'$ y sus derivadas), evitando errores físicos cerca de $x=0$.
3. Validación Exhaustiva: Los resultados se validan utilizando datos de Simulación Numérica Directa (DNS) de alta fidelidad (Howland et al., 2022) cubriendo un amplio rango de parámetros: $1 \le Pr \le 100$ y $10^6 \le Ra \le 10^9$.

4. Resultados Principales

• Perfiles de Temperatura:
  • Región Interior: El perfil de temperatura se describe mediante la función $F_i(z)$, que difiere de una función lineal, indicando que el flujo de calor turbulento es significativo incluso cerca de la pared. La escala de temperatura $T_i$ depende de $Pr$y$Ra$ de una manera que incluye la viscosidad cinemática, a diferencia de modelos anteriores que ignoraban este efecto.
  • Región Exterior: El perfil se describe mediante una función logarítmica $F_o(z)$, válida en la zona cercana a la línea central.
• Escalas de Velocidad Características: Los parámetros del modelo ($A$ y $C_m$) determinan las velocidades características de transferencia de calor en las regiones interior ($V_i$) y exterior ($V_o$). Se demuestra que $V_i$ está relacionada con la velocidad de fricción $u_\tau$ mediante una función dependiente de $Pr$.
• Comportamiento del Número de Nusselt: El modelo predice correctamente que $Nu \sim Ra^{1/3}$ en el límite de alto $Ra$, consistente con teorías previas y datos DNS.
• Comparación con la Literatura:
  • El modelo supera a las propuestas de George & Capp (1979) y Hölling & Herwig (2005), que fallan al colapsar los datos para valores altos de $Pr$.
  • Se demuestra que el modelo de Li et al. (2023), que predice una dependencia de raíz cúbica inversa en la región exterior, no coincide con los datos DNS, mientras que la función logarítmica propuesta en este trabajo sí lo hace.

5. Significado e Impacto

Este trabajo proporciona una comprensión teórica más robusta de la convección vertical turbulenta. Al ofrecer un modelo analítico que:

1. Es consistente con las condiciones físicas en las fronteras.
2. Captura la dependencia correcta con el número de Prandtl.
3. Se ajusta excelentemente a datos DNS de alta resolución.

El estudio establece un nuevo estándar para la modelización de perfiles de temperatura en convección natural, con implicaciones directas para la mejora de modelos de ventilación, diseño térmico de edificios y la comprensión de procesos geofísicos en interfaces verticales. La identificación de dos funciones de escala universales simplifica la predicción de perfiles térmicos sin necesidad de simulaciones costosas para cada configuración de parámetros.