AUTHOR=Huang Xi , Pang Bo , Chai Xiang , Yin Yuan 

TITLE=Proposal of a turbulent Prandtl number model for Reynolds-averaged Navier–Stokes approach on the modeling of turbulent heat transfer of low-Prandtl number liquid metal

JOURNAL=Frontiers in Energy Research

VOLUME=10

YEAR=2022

URL=https://www.frontiersin.org/journals/energy-research/articles/10.3389/fenrg.2022.928693

DOI=10.3389/fenrg.2022.928693

ISSN=2296-598X

ABSTRACT=<p>Because of their high molecular heat conductivity, low-Prandtl number liquid metal is a promising candidate coolant for various designs of advanced nuclear systems such as liquid metal–cooled fast reactors and accelerator-driven sub-critical system (ADS). With the fast-growing computational capacity, more and more attention has been paid to applying computational fluid dynamics (CFD) methods in thermal design and safety assessment of such systems for a detailed analysis of three-dimensional thermal–hydraulic behaviors. However, numerical modeling of turbulent heat transfer for low-Prandtl number liquid metal remains a challenging task. Numerical approaches such as wall-resolved large eddy simulation (LES) or direct numerical simulation (DNS), which can provide detailed insight into the physics of the liquid metal flow and the associated heat transfer, were widely applied to investigate the turbulent heat transfer phenomenon. However, these approaches suffer from the enormous computational consumption and are hence limited only to simple geometrical configurations with low to moderate Reynolds numbers. The Reynolds-averaged Navier–Stokes (RANS) approach associated with a turbulent Prandtl number <inline-formula id="inf1"><mml:math id="m1" xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:mi>P</mml:mi><mml:msub><mml:mi>r</mml:mi><mml:mi>t</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> accounting for the turbulent heat flux based on Reynolds analogy is still, at least in the current state in most of the circumstances, the only feasible approach for practical engineering applications. However, the conventional choice of <inline-formula id="inf2"><mml:math id="m2" xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:mi>P</mml:mi><mml:msub><mml:mi>r</mml:mi><mml:mi>t</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> in the order of 0.9∼unity in many commercial computational fluid dynamics codes is not valid for the low-Prandtl number liquid metal. In this study, LES/DNS simulation results of a simple forced turbulent channel flow up to a friction Reynolds number <inline-formula id="inf3"><mml:math id="m3" xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:mi>R</mml:mi><mml:msub><mml:mi>e</mml:mi><mml:mi>τ</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> of 2000 at <inline-formula id="inf4"><mml:math id="m4" xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:mi>P</mml:mi><mml:mi>r</mml:mi></mml:mrow></mml:math></inline-formula> of 0.01 and 0.025 were used as references, to which the Reynolds-averaged Navier–Stokes approach with varying <inline-formula id="inf5"><mml:math id="m5" xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:mi>P</mml:mi><mml:msub><mml:mi>r</mml:mi><mml:mi>t</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> was compared. It was found that the appropriate <inline-formula id="inf6"><mml:math id="m6" xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:mi>P</mml:mi><mml:msub><mml:mi>r</mml:mi><mml:mi>t</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> for the RANS approach decreases with bulk Peclet number <inline-formula id="inf7"><mml:math id="m7" xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:mi>P</mml:mi><mml:msub><mml:mi>e</mml:mi><mml:mi>b</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> and approaches a constant value of 1.5 when <inline-formula id="inf8"><mml:math id="m8" xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:mi>P</mml:mi><mml:msub><mml:mi>e</mml:mi><mml:mi>b</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> becomes larger than 2000. Based on this calibrated relation with <inline-formula id="inf9"><mml:math id="m9" xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:mi>P</mml:mi><mml:msub><mml:mi>e</mml:mi><mml:mi>b</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, a new model for <inline-formula id="inf10"><mml:math id="m10" xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:mi>P</mml:mi><mml:msub><mml:mi>r</mml:mi><mml:mi>t</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> used in the RANS approach was proposed. Validation of the proposed model was carried out with available LES/DNS results on the local temperature profile in the concentric annulus and bare rod bundle, as well as with experimental correlations on the Nusselt number in a circular tube and bare rod bundle.</p>