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

Resumo Técnico: Características de Fluxo e Transferência de Calor de Convecção Forçada sobre um Cilindro Circular Isofluxo em Galinstan para Números de Reynolds até 600

Definição do Problema
Este estudo apresenta uma investigação numérica da transferência de calor por convecção forçada estacionária de um cilindro circular sujeito a uma condição de contorno de fluxo de calor constante (isofluxo). O meio fluido é o Galinstan, uma liga metálica líquida caracterizada por um baixo número de Prandtl ($Pr = 0,025$) e condutividade térmica excepcionalmente alta. A pesquisa foca no comportamento hidrodinâmico e térmico do escoamento para uma faixa de número de Reynolds ($Re$) de $1 \le Re \le 600$. Embora a literatura anterior tenha abrangido extensivamente cilindros isotérmicos e fluidos de número de Prandtl moderado a alto, este trabalho aborda os desafios específicos de simular metais líquidos de baixo número de Prandtl sob condições de isofluxo, um regime relevante para sistemas avançados de gestão térmica, como o resfriamento de eletrônicos e sistemas de energia baseados em metais líquidos.

Metodologia
As equações governantes — compreendendo as equações de função de corrente, vorticidade e energia em coordenadas cilíndricas — são resolvidas utilizando um esquema de diferenças finitas compactas de quarta ordem (FOCS–CC). Para garantir estabilidade e convergência, esta discretização espacial é acoplada a uma técnica de iteração de tempo pseudo (PTI).

Principais características metodológicas incluem:

• Transformação de Coordenadas: O domínio físico é mapeado para um domínio computacional usando a transformação $r = e^{\pi\lambda}$ e $\theta = \pi\eta$, permitindo uma resolução de malha eficiente próximo à superfície do cilindro.
• Precisão de Alta Ordem: O esquema utiliza formulações compactas de quarta ordem para derivadas espaciais. Uma característica notável destacada pelos autores é a facilidade de implementação; o método estende as estruturas existentes de segunda ordem ao incorporar coeficientes de correção específicos ($\alpha, b_1, \Gamma_1$, etc.) sem alterar a estrutura numérica subjacente.
• Algoritmo de Solução: O sistema de equações resultante é integrado usando uma abordagem fatorada para reduzir grandes matrizes de banda em sistemas tridiagonais, que são resolvidos eficientemente usando o Algoritmo de Matriz Tridiagonal (TDMA).
• Validação: O código numérico é validado através de estudos de independência de malha e pela comparação do número de Nusselt médio ($Nu$) calculado e do coeficiente de arrasto total ($C_D$) com resultados numéricos estabelecidos de Dennis e Chang (1970), Fornberg (1980), Ahmad (1996) e Morgan (1975).

Principais Resultados e Observações
O estudo analisa sistematicamente os campos de fluxo e térmico através da faixa especificada de número de Reynolds:

1. Estrutura do Fluxo: Em baixos números de Reynolds ($Re=1$), o escoamento permanece totalmente aderido à superfície do cilindro, sem separação, dominado por forças viscosas. À medida que o $Re$ aumenta para 200 e 600, ocorre a separação do escoamento, formando uma região de esteira atrás do cilindro. O tamanho desta esteira expande-se com o aumento do $Re$, indicando uma crescente dominância das forças de inércia.
2. Campo Térmico: Para $Re=1$, as isotermas formam anéis quase concêntricos, significando que a transferência de calor é governada primordialmente pela difusão térmica. À medida que o $Re$ aumenta, o transporte convectivo intensifica-se, fazendo com que as isotermas se agrupem mais densamente no ponto de estagnação a montante e se estendam a jusante, formando uma pluma térmica. Apesar deste alongamento convectivo, o campo térmico permanece relativamente difuso em comparação com fluidos convencionais devido à alta difusividade térmica do Galinstan.
3. Métricas de Transferência de Calor: O número de Nusselt local exibe variação angular mínima em baixos $Re$, mas mostra variação significativa em $Re$ mais altos devido à influência da esteira separada. O número de Nusselt médio aumenta monotonicamente com o $Re$ em todo o intervalo investigado, embora a taxa de aumento diminua em números de Reynolds mais elevados.
4. Correlação Empírica: Com base nos dados numéricos, uma nova correlação empírica para o número de Nusselt médio é proposta para o Galinstan entre $1 \le Re \le 600$:
   $$Nu = \alpha + \beta\sqrt{Re} + \gamma Re$$
   onde $\alpha = 0,26133$, $\beta = 0,13302$ e $\gamma = -5,42758 \times 10^{-4}$. Esta correlação demonstra excelente concordância com os dados numéricos, apresentando um coeficiente de determinação ($R^2$) de 0,99939.

Significância e Alegações
Os autores afirmam que o esquema FOCS–CC proposto, integrado com a técnica PTI, fornece um arcabouço numérico robusto e preciso para simular a convecção forçada ao redor de cilindros circulares em metais líquidos de baixo número de Prandtl. O estudo estabelece que:

• Soluções estacionárias para as equações de Navier–Stokes incompressíveis podem ser computadas com sucesso até $Re = 600$ para esta configuração, estendendo-se além do início da instabilidade de fluxo não estacionário.
• A interação entre o transporte convectivo intensificado e a difusão térmica dominante é uma característica definidora do Galinstan, distinguindo seu comportamento térmico de fluidos com números de Prandtl moderados ou altos.
• A alta precisão do esquema compacto de quarta ordem permite a previsão confiável de ambas as características hidrodinâmicas (arrasto) e térmicas (número de Nusselt), conforme confirmado pela validação contra a literatura existente.

O artigo conclui que o método numérico e a correlação desenvolvidos oferecem ferramentas valiosas para analisar a transferência de calor em aplicações de metais líquidos, destacando especificamente os retornos decrescentes na melhoria da transferência de calor por convecção em números de Reynolds mais elevados devido à persistente influência da alta difusividade térmica.