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

Resumo Técnico: Solução Numérica de Convecção Natural em uma Placa Vertical Aquecida com Baixas Razões de Difusividade

Enunciado do Problema
Este estudo investiga o escoamento de convecção natural transiente ao longo de uma placa vertical isotérmica sob condições de baixas razões de difusividade, visando especificamente fluidos onde a difusividade térmica excede significativamente a difusividade de momento (baixos números de Prandtl, $Pr \ll 1$). Embora a convecção natural sobre placas verticais seja um parâmetro de referência bem estabelecido, a evolução transiente e as estruturas de camada limite em regimes de baixo $Pr$ — característicos de metais líquidos e sais fundidos — diferem substancialmente de fluidos de $Pr$ moderado ou alto (como ar ou óleos). O artigo aborda uma lacuna na literatura referente ao comportamento transiente específico das camadas limites térmicas e dos mecanismos de transferência de calor nessas condições de alta difusividade, onde a difusão térmica domina o transporte de momento impulsionado pela flutuabilidade.

Metodologia
A pesquisa emprega uma abordagem numérica baseada no Método das Diferenças Finitas (MDF) para resolver as equações adimensionais governantes derivadas da aproximação de Boussinesq. O sistema de equações inclui a equação de continuidade, as equações de momento e a equação de energia.

• Discretização: O domínio espacial utiliza uma grade cartesiana estruturada uniforme. Os termos difusivos são aproximados por diferenças finitas centrais de segunda ordem para preservar a precisão espacial, enquanto os termos convectivos não lineares são discretizados usando um esquema upwind de primeira ordem para suprimir oscilações numéricas espúrias em regiões com gradientes acentuados.
• Integração Temporal: Um esquema de Euler explícito é utilizado para o avanço temporal. Para garantir a estabilidade numérica, o passo de tempo ($\Delta\tau$) adere estritamente à condição de Courant–Friedrichs–Lewy (CFL), que é particularmente restritiva para fluidos de baixo $Pr$ devido à difusividade térmica aumentada.
• Validação: A metodologia numérica é validada através de um estudo de independência de malha (convergindo em uma malha de $100 \times 100$) e pela comparação dos resultados com correlações clássicas estabelecidas e dados recentes da literatura para convecção de baixo $Pr$. As simulações cobrem uma faixa de número de Prandtl de $0,01$a$0,71$, com números de Reynolds ($Re=10^2$) e Rayleigh ($Ra=10^5$) fixos.

Principais Resultados
As simulações numéricas revelam comportamentos físicos distintos em regimes de baixo $Pr$ em comparação com a convecção convencional:

1. Estrutura da Camada Limite: À medida que o número de Prandtl diminui, a camada limite térmica torna-se significativamente mais espessa que a camada limite hidrodinâmica (velocidade). Para $Pr$ extremamente baixos (ex: $Pr=0,01$), o campo de temperatura penetra profundamente no domínio do fluido, resultando em um decaimento térmico muito gradual.
2. Perfis de Velocidade e Temperatura: Em fluidos de baixo $Pr$, o pico de velocidade desloca-se para longe da parede aquecida, e a camada limite de momento alarga-se. Isso ocorre porque as forças de flutuabilidade são distribuídas sobre uma região térmica maior, levando a uma aceleração local mais fraca perto da parede. Inversamente, fluidos de $Pr$ mais altos exibem gradientes térmicos e picos de velocidade mais nítidos próximos à superfície.
3. Taxas de Transferência de Calor: O número de Nusselt local ($Nu_x$) diminui conforme o número de Prandtl cai. A dominância da difusão térmica espalha a energia térmica rapidamente para o fluido, reduzindo o gradiente de temperatura normal à parede e, consequentemente, diminuindo as taxas de transferência de calor por convecção.
4. Evolução Transiente: A convecção de baixo $Pr$ exibe uma resposta térmica retardada. Fluidos com $Pr=0,01$ apresentam uma convergência assintótica lenta em direção ao equilíbrio térmico, enquanto fluidos de $Pr$ moderado e alto estabilizam-se muito mais rapidamente.
5. Leis de Escalonamento: Os resultados confirmam uma relação de escalonamento onde a espessura da camada limite térmica é inversamente proporcional à raiz quadrada do número de Prandtl ($\delta_t \sim Pr^{-1/2}$), consistente com previsões teóricas para regimes dominados por difusão.

Significância e Alegações
O artigo alega fornecer um arcabouço numérico robusto capaz de capturar com precisão os mecanismos de transporte térmico fortemente difusivos inerentes à convecção natural de baixo $Pr$. Ao utilizar um MDF explícito com estratégias de estabilização específicas (upwinding e adesão estrita ao CFL), o estudo demonstra que é possível resolver o complexo acoplamento entre os campos de velocidade e temperatura nesses regimes.

Os autores afirmam que suas descobertas oferecem percepções físicas relevantes para aplicações de engenharia envolvendo sistemas de resfriamento por metal líquido, gestão térmica de reatores nucleares e tecnologias avançadas de energia térmica. O trabalho enfatiza que compreender a evolução transiente dos campos térmicos em fluidos de baixa difusividade é crítico para o projeto de sistemas de resfriamento passivo e trocadores de calor onde a difusão térmica dita o comportamento de transporte. O estudo conclui que o modelo numérico proposto é confiável para investigar esses fenômenos, conforme evidenciado pelo seu acordo com dados de referência e sua capacidade de reproduzir tendências físicas esperadas, como o espessamento das camadas térmicas e a redução dos coeficientes de transferência de calor em condições de baixo $Pr$.