Numerical Solution of Natural Convection on a Heated Vertical Plate at Low Diffusivity Ratios

Resumen Técnico: Solución Numérica de la Convección Natural en una Placa Vertical Calefactada con Bajas Relaciones de Difusividad

Planteamiento del Problema
Este estudio investiga el flujo de convección natural transitorio a lo largo de una placa vertical isotérmica bajo condiciones de bajas relaciones de difusividad, centrándose específicamente en fluidos donde la difusividad térmica excede significativamente la difusividad de cantidad de movimiento (números de Prandtl bajos, $Pr \ll 1$). Si bien la convección natural sobre placas verticales es un referente bien establecido, la evolución transitoria y las estructuras de la capa límite en regímenes de bajo $Pr$ —característicos de metales líquidos y sales fundidas— difieren sustancialmente de los fluidos de $Pr$ moderado o alto (como el aire o los aceites). El artículo aborda un vacío en la literatura respecto al comportamiento transitorio específico de las capas límite térmicas y los mecanismos de transferencia de calor en estas condiciones altamente difusivas, donde la difusión térmica domina el transporte de cantidad de movimiento impulsado por la flotabilidad.

Metodología
La investigación emplea un enfoque numérico basado en el Método de Diferencias Finitas (MDF) para resolver las ecuaciones adimensionales gobernantes derivadas de la aproximación de Boussinesq. El sistema de ecuaciones incluye la ecuación de continuidad, las ecuaciones de cantidad de movimiento y la ecuación de energía.

• Discretización: El dominio espacial utiliza una malla cartesiana uniforme estructurada. Los términos difusivos se aproximan mediante diferencias finitas centrales de segundo orden para preservar la precisión espacial, mientras que los términos convectivos no lineales se discretizan utilizando un esquema de upwind (corriente arriba) de primer orden para suprimir oscilaciones numéricas espurias en regiones con gradientes pronunciados.
• Integración Temporal: Se utiliza un esquema de Euler explícito para el avance temporal. Para garantizar la estabilidad numérica, el paso de tiempo ($\Delta\tau$) se adhiere estrictamente a la condición de Courant–Friedrichs–Lewy (CFL), la cual es particularmente restrictiva para fluidos de bajo $Pr$ debido a la difusividad térmica mejorada.
• Validación: La metodología numérica se valida mediante un estudio de independencia de malla (convergiendo en una malla de $100 \times 100$) y comparando los resultados con correlaciones clásicas establecidas y datos recientes de la literatura para la convección de bajo $Pr$. Las simulaciones cubren un rango de número de Prandtl de $0.01$a$0.71$, con números de Reynolds ($Re=10^2$) y Rayleigh ($Ra=10^5$) fijos.

Resultos Clave
Las simulaciones numéricas revelan comportamientos físicos distintos en regímenes de bajo $Pr$ en comparación con la convección convencional:

1. Estructura de la Capa Límite: A medida que el número de Prandtl disminuye, la capa límite térmica se vuelve significativamente más gruesa que la hidrodinámica (de velocidad). Para valores de $Pr$ extremadamente bajos (p. ej., $Pr=0.01$), el campo de temperatura penetra profundamente en el dominio del fluido, resultando en un decaimiento térmico muy gradual.
2. Perfiles de Velocidad y Temperatura: En fluidos de bajo $Pr$, el pico de velocidad se desplaza alejándose de la pared caliente y la capa límite de cantidad de movimiento se ensancha. Esto ocurre porque las fuerzas de flotabilidad se distribuyen sobre una región térmica más amplia, lo que conduce a una aceleración local más débil cerca de la pared. Por el contrario, los fluidos de mayor $Pr$ exhiben gradientes térmicos y picos de velocidad más agudos cerca de la superficie.
3. Tasas de Transferencia de Calor: El número de Nusselt local ($Nu_x$) disminuye a medida que el número de Prandtl cae. El predominio de la difusión térmica dispersa la energía térmica rápidamente hacia el fluido, reduciendo el gradiente de temperatura normal a la pared y, consecuentemente, disminuyendo las tasas de transferencia de calor por convección.
4. Evolución Transitoria: La convección de bajo $Pr$ exhibe una respuesta térmica retardada. Los fluidos con $Pr=0.01$ muestran una convergencia asintótica lenta hacia el equilibrio térmico, mientras que los fluidos de $Pr$ moderado y alto se estabilizan mucho más rápido.
5. Leyes de Escalamiento: Los resultados confirman una relación de escalamiento donde el espesor de la capa límite térmica ($\delta_t$) es inversamente proporcional a la raíz cuadrada del número de Prandtl ($\delta_t \sim Pr^{-1/2}$), lo cual es consistente con las predicciones teóricas para regímenes dominados por la difusión.

Significancia y Reivindicaciones
El artículo afirma proporcionar un marco numérico robusto capaz de capturar con precisión los mecanismos de transporte térmico fuertemente difusivos inherentes a la convección natural de bajo $Pr$. Al utilizar un MDF explícito con estrategias de estabilización específicas (upwinding y adherencia estricta a CFL), el estudio demuestra que es posible resolver el complejo acoplamiento entre los campos de velocidad y temperatura en estos regímenes.

Los autores sostienen que sus hallazgos ofrecen perspectivas físicas relevantes para aplicaciones de ingeniería que involucran sistemas de enfriamiento por metales líquidos, gestión térmica de reactores nucleares y tecnologías avanzadas de energía térmica. El trabajo enfatiza que comprender la evolución transitoria de los campos térmicos en fluidos de baja difusividad es crítico para el diseño de sistemas de enfriamiento pasivo y de intercambiadores de calor donde la difusión térmica dicta el comportamiento del transporte. El estudio concluye que el modelo numérico propuesto es confiable para investigar estos fenómenos, como lo evidencia su concordancia con los datos de referencia y su capacidad para reproducir las tendencias físicas esperadas, tales como el engrosamiento de las capas térmicas y la reducción de los coeficientes de transferencia de calor en condiciones de bajo $Pr$.