Flow and Heat Transfer Characteristics of Forced Convection Past an Isoflux Circular Cylinder in Galinstan for Reynolds Numbers up to 600

Resumen Técnico: Características de Flujo y Transferencia de Calor de la Convección Forzada sobre un Cilindro Circular de Isoflujo en Galinstan para Números de Reynolds hasta 600

Planteamiento del Problema
Este estudio presenta una investigación numérica sobre la transferencia de calor por convección forzada en estado estacionario de un cilindro circular sometido a una condición de contorno de flujo de calor constante (isofluxo). El medio fluido es Galinstan, una aleación de metal líquido caracterizada por un número de Prandtl bajo ($Pr = 0.025$) y una conductividad térmica excepcionalmente alta. La investigación se centra en el comportamiento hidrodinámico y térmico del flujo sobre un rango de número de Reynolds ($Re$) de $1 \le Re \le 600$. Si bien la literatura previa ha cubierto extensamente cilindros isotérmicos y fluidos de número de Prandtl moderado a alto, este trabajo aborda los desafíos específicos de simular metales líquidos de bajo número de Prandtl bajo condiciones de isofluxo, un régimen relevante para sistemas avanzados de gestión térmica, tales como la refrigeración electrónica y los sistemas de energía basados en metales líquidos.

Metodología
Las ecuaciones gobernantes —que comprenden las ecuaciones de la función de corriente, vorticidad y energía en coordenadas cilíndricas— se resuelven utilizando un esquema de diferencias finitas compactas de cuarto orden (FOCS–CC). Para asegurar la estabilidad y convergencia, esta discretización espacial se acopla con una técnica de iteración de pseudo-tiempo (PTI).

Características metodológicas clave:

• Transformación de Coordenadas: El dominio físico se mapea a un dominio computacional utilizando la transformación $r = e^{\pi\lambda}$ y $\theta = \pi\eta$, lo que permite una resolución de malla eficiente cerca de la superficie del cilindro.
• Precisión de Alto Orden: El esquema utiliza formulaciones compactas de cuarto orden para las derivadas espaciales. Una característica notable destacada por los autores es la facilidad de implementación; el método extiende los marcos existentes de segundo orden mediante la incorporación de coeficientes de corrección específicos ($\alpha, b_1, \Gamma_1$, etc.) sin alterar la estructura numérica subyacente.
• Algoritmo de Solución: El sistema de ecuaciones resultante se integra mediante un enfoque factorizado para reducir las matrices de gran banda en sistemas tridiagonales, los cuales se resuelven eficientemente utilizando el Algoritmo de Matriz Tridiagonal (TDMA).
• Validación: El código numérico se valida mediante estudios de independencia de malla y comparando el número de Nusselt promedio ($Nu$) y el coeficiente de arrastre total ($C_D$) computados contra resultados numéricos establecidos de Dennis y Chang (1970), Fornberg (1980), Ahmad (1996) y Morgan (1975).

Resultados Clave y Observaciones
El estudio analiza sistemáticamente el campo de flujo y térmico a través del rango de número de Reynolds especificado:

1. Estructura del Flujo: Para números de Reynolds bajos ($Re=1$), el flujo permanece completamente adherido a la superficie del cilindro sin separación, dominado por fuerzas viscosas. A medida que $Re$ aumenta a 200 y 600, ocurre la separación del flujo, formando una región de estela detrás del cilindro. El tamaño de esta estela se expande con el aumento de $Re$, indicando un crecimiento en la dominancia de las fuerzas de inercia.
2. Campo Térmico: Para $Re=1$, las isotermas forman anillos casi concéntricos, lo que significa que la transferencia de calor está gobernada principalmente por la difusión térmica. A medida que $Re$ aumenta, el transporte convectivo se intensifica, causando que las isotermas se agrupen más estrechamente en el punto de estancamiento aguas arriba y se estiren hacia aguas abajo, formando una pluma térmica. A pesar de este estiramiento por convección, el campo térmico permanece relativamente difuso en comparación con fluidos convencionales debido a la alta difusividad térmica del Galinstan.
3. Métricas de Transferencia de Calor: El número de Nusselt local exhibe una variación angular mínima a bajo $Re$, pero muestra una variación significativa a $Re$ más altos debido a la influencia de la estela separada. El número de Nusselt promedio aumenta de forma monótona con $Re$ en todo el rango investigado, aunque la tasa de mejora disminuye en números de Reynolds más altos.
4. Correlación Empírica: Basándose en los datos numéricos, se propone una nueva correlación empírica para el número de Nusselt promedio para el Galinstan sobre $1 \le Re \le 600$:
   $$Nu = \alpha + \beta\sqrt{Re} + \gamma Re$$
   donde $\alpha = 0.26133$, $\beta = 0.13302$, y $\gamma = -5.42758 \times 10^{-4}$. Esta correlación demuestra una excelente concordancia con los datos numéricos, obteniendo un coeficiente de determinación ($R^2$) de 0.99939.

Significancia y Reivindicaciones
Los autores afirman que el esquema FOCS–CC propuesto, integrado con la técnica PTI, proporciona un marco numérico robusto y preciso para simular la convección forzada alrededor de cilindros circulares en metales líquidos de bajo número de Prandtl. El estudio establece que:

• Se pueden computar con éxito soluciones estacionarias para las ecuaciones de Navier–Stokes incompresibles hasta $Re = 600$ para esta configuración, extendiéndose más allá del inicio de la inestabilidad de flujo no estacionario.
• La interacción entre la transferencia de calor convectiva mejorada y la difusión térmica dominante es una característica definitoria del Galinstan, lo que distingue su comportamiento térmico del de los fluidos con números de Prandtl moderados o altos.
• La alta precisión del esquema compacto de cuarto orden permite la predicción confiable tanto de las características hidrodinámicas (arrastre) como térmicas (número de Nusselt), según lo confirmado por la validación contra la literatura existente.

El artículo concluye que la correlación y el método numérico desarrollados ofrecen herramientas valiosas para analizar la transferencia de calor en aplicaciones de metales líquidos, destacando específicamente los rendimientos decrecientes en la mejora de la transferencia de calor por convección a números de Reynolds más altos debido a la persistente influencia de la alta difusividad térmica.