Numerical Study of Low Pr Flow in a Bare 19-Rod Bundle Based on an Advanced Turbulent Heat Transfer Model

Abstract

Compared to assuming a constant turbulent Prandtl number model, an advanced four-equation model has the potential to improve the numerical heat transfer calculation accuracy of low–Prandtl number (Pr) fluids. Generally, a four-equation model consists of a two-equation k−ε turbulence model and a two-equation kθ−εθ heat transfer model. It is essential to analyze the influence of dissimilar turbulence models on the overall calculation accuracy of the four-equation model. The present study aims to study the effect of using different turbulence models on the same kθ−εθ heat transfer model. First, based on the open-source computational fluid dynamics software OpenFOAM, an advanced two-equation kθ−εθ heat transfer model was introduced into the solver buoyant2eqnFoam, which was developed based on the self-solver buoyantSimpleFoam of OpenFOAM. In the solver buoyant2eqnFoam, various turbulence models built into OpenFOAM can be conveniently called to close the Reynolds stress and an advanced two-equation heat transfer model can be utilized to calculate the Reynolds heat flux of low-Pr fluids. Subsequently, the solver buoyant2eqnFoam was employed to study the fully developed flow heat transfer of low-Pr fluids in a bare 19-rod bundle. The numerical results were compared and analyzed with the experimental correlations and other simulation results to validate the effectiveness and feasibility of the solver buoyant2eqnFoam. Furthermore, the influence of combining different turbulence models with the same two-equation kθ−εθ heat transfer model was also presented in this study. The results show that the turbulence model has a considerable influence on the prediction of turbulent heat transfer in the high Peclet number range, suggesting that it should be prudent when picking a turbulence model in the simulations of low-Pr fluids.