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

Résumé technique : Solution numérique de la convection naturelle sur une plaque verticale chauffée à faibles rapports de diffusivité

Énoncé du problème
Cette étude examine l'écoulement de la convection naturelle transitoire le long d'une plaque verticale isotherme dans des conditions de faibles rapports de diffusivité, ciblant spécifiquement les fluides où la diffusivité thermique excède significativement la diffusivité de quantité de mouvement (faibles nombres de Prandtl, $Pr \ll 1$). Bien que la convection naturelle sur des plaques verticales soit un cas de référence bien établi, l'évolution transitoire et les structures de la couche limite dans les régimes à faible $Pr$ — caractéristiques des métaux liquides et des sels fondus — diffèrent substantiellement de celles des fluides à $Pr$ modéré ou élevé (comme l'air ou les huiles). L'article comble une lacune dans la littérature concernant le comportement transitoire spécifique des couches limites thermiques et des mécanismes de transfert de chaleur dans ces conditions de haute diffusivité, où la diffusion thermique domine le transport de quantité de mouvement induit par la flottabilité.

Méthodologie
La recherche emploie une approche numérique basée sur la méthode des différences finies (FDM) pour résoudre les équations adimensionnelles dérivées de l'approximation de Boussinesq. Le système d'équations comprend l'équation de continuité, les équations de quantité de mouvement et l'équation de l'énergie.

• Discrétisation : Le domaine spatial utilise une grille de type de Cartesian uniforme structurée. Les termes diffusifs sont approximés par des différences finies centrales de second ordre pour préserver la précision spatiale, tandis que les termes convectifs non linéaires sont discrétisés à l'aide d'un schéma amont (upwind) de premier ordre afin de supprimer les oscillations numériques parasites dans les régions présentant des gradients abrupts.
• Intégration temporelle : Un schéma d'Euler explicite est utilisé pour l'avancement temporel. Pour garantir la stabilité numérique, le pas de temps ($\Delta\tau$) respecte strictement la condition de Courant–Friedrichs–Lewy (CFL), laquelle est particulièrement contraignante pour les fluides à faible $Pr$ en raison de la diffusivité thermique accrue.
• Validation : La méthodologie numérique est validée par une étude d'indépendance du maillage (convergeant vers une grille de $100 \times 100$) et par la comparaison des résultats avec des corrélations classiques établies et des données récentes de la littérature pour la convection à faible $Pr$. Les simulations couvrent une plage de nombres de Prandtl allant de $0,01$à$0,71$, avec des nombres de Reynolds ($Re=10^2$) et de Rayleigh ($Ra=10^5$) fixes.

Résultats clés
Les simulations numériques révèlent des comportements physiques distincts dans les régimes à faible $Pr$ par rapport à la convection conventionnelle :

1. Structure de la couche limite : À mesure que le nombre de Prandtl diminue, la couche limite thermique devient nettement plus épaisse que la couche limite hydrodynamique (vitesse). Pour des $Pr$ extrêmement bas (ex. $Pr=0,01$), le champ de température pénètre profondément dans le domaine fluide, entraînant une décroissance thermique très graduelle.
2. Profils de vitesse et de température : Dans les fluides à faible $Pr$, le pic de vitesse s'éloigne de la paroi chauffée et la couche limite de quantité de mouvement s'élargit. Cela se produit car les forces de flottabilité sont réparties sur une région thermique plus large, entraînant une accélération locale plus faible près de la paroi. Inversement, les fluides à $Pr$ élevé présentent des gradients thermiques et des pics de vitesse plus marqués près de la surface.
3. Taux de transfert de chaleur : Le nombre de Nusselt local ($Nu_x$) diminue lorsque le nombre de Prandtl chute. La dominance de la diffusion thermique propage rapidement l'énergie thermique dans le fluide, réduisant le gradient de température normal à la paroi et, par conséquent, abaissant les taux de transfert de chaleur par convection.
4. Évolution transitoire : La convection à faible $Pr$ présente une réponse thermique retardée. Les fluides avec $Pr=0,01$ montrent une convergence asymptotique lente vers l'équilibre thermique, tandis que les fluides à $Pr$ modéré ou élevé se stabilisent beaucoup plus rapidement.
5. Lois d'échelle : Les résultats confirment une relation d'échelle où l'épaisseur de la couche limite thermique ($\delta_t$) est inversement proportionnelle à la racine carrée du nombre de Prandtl ($\delta_t \sim Pr^{-1/2}$), ce qui est cohérent avec les prédictions théoriques pour les régimes dominés par la diffusion.

Signification et affirmations
L'article affirme fournir un cadre numérique robuste capable de capturer avec précision les mécanismes de transport thermique fortement diffusifs inhérents à la convection naturelle à faible $Pr$. En utilisant une méthode de différences finies explicite avec des stratégies de stabilisation spécifiques (amont et respect strict de la condition CFL), l'étude démontre qu'il est possible de résoudre le couplage complexe entre les champs de vitesse et de température dans ces régimes.

Les auteurs soutiennent que leurs conclusions offrent des perspectives physiques pertinentes pour les applications d'ingénierie impliquant les systèmes de refroidissement par métaux liquides, la gestion thermique des réacteurs nucléaires et les technologies avancées de l'énergie thermique. Le travail souligne que la compréhension de l'évolution transitoire des champs thermiques dans les fluides à faible diffusivité est cruciale pour la conception de systèmes de refroidissement passifs et d'échangeurs de chaleur où la diffusion thermique dicte le comportement du transport. L'étude conclut que le modèle numérique proposé est fiable pour l'investigation de ces phénomènes, comme en témoigne son accord avec les données de référence et sa capacité à reproduire les tendances physiques attendues, telles que l'épaississement des couches limites thermiques et la réduction des coefficients de transfert de chaleur dans les conditions de faible $Pr$.