Assessment of RANS and improved near-wall modeling for forced convection at low Prandtl numbers based on LES ur to Rer=2000

Most promising Generation IV nuclear reactor concepts are based on a liquid metal coolant. However, at low Prandtl (Pr ) numbers such as those of liquid metal, classical approaches derived for unity, or close to unity, Pr fail to accurately predict the heat transfer. This paper assesses the RANS modeling of forced turbulent convection at low Pr and in channel flow. Reference results at high Reynolds (Re ) number are required to ensure that the Peclet number is sufficiently high. Therefore, new reference results were obtained by performing a wall-resolved Large-Eddy Simulation of turbulent channel flows at a friction Reynolds number Reτ=2000Reτ=2000 and at Pr=0.01Pr=0.01 and 0.0250.025 (this also corresponds to the highest Re Direct Numerical Simulation (DNS) available in the literature for the flow, but without heat transfer). The LES velocity statistics are in very good agreement with those of the DNS and, as validated by the authors in previous publications, the LES approach used here accurately predicts the temperature statistics at low Pr . The LES results are used to assess RANS heat transfer modeling based on the effective turbulent Prandtl number (PrtPrt) concept. Among existing PrtPrt correlations, the correlation by Kays (1994) [10] is shown to yield the best results. Since it is also shown that the near-wall temperature profile does not follow a log-law, a new “law of the wall for temperature” is here proposed, which does not use any blending function. Its use as a wall-function is also validated in actual RANS simulations.

[1]  A. J. Reynolds.  The prediction of turbulent Prandtl and Schmidt numbers , 1975 .

[2]  Byung Do Oh,et al.  Canonical fractional-step methods and consistent boundary conditions for the incompressible Navier—Stokes equations , 2001 .

[3]  Yuichi Matsuo,et al.  DNS of turbulent heat transfer in channel flow with respect to Reynolds and Prandtl number effects , 1999 .

[4]  K. Chien,et al.  Predictions of Channel and Boundary-Layer Flows with a Low-Reynolds-Number Turbulence Model , 1982 .

[5]  Konstantin Mikityuk,et al.  Heat transfer to liquid metal: Review of data and correlations for tube bundles , 2009 .

[6]  G. Grötzbach,et al.  Challenges in low-Prandtl number heat transfer simulation and modelling , 2013 .

[7]  G. Grötzbach.  Revisiting the resolution requirements for turbulence simulations in nuclear heat transfer , 2011 .

[8]  Bernhard Weigand,et al.  An extended Kays and Crawford turbulent Prandtl number model , 1997 .

[9]  Nam-il Tak,et al.  Investigation on turbulent heat transfer to lead–bismuth eutectic flows in circular tubes for nuclear applications , 2006 .

[10]  W. Kays,et al.  Turbulent Prandtl number : where are we ? , 1994 .

[11]  L. Chandra,et al.  CFD analyses of liquid metal flow in sub-channels for Gen IV reactors , 2011 .

[12]  Javier Jiménez,et al.  Scaling of the velocity fluctuations in turbulent channels up to Reτ=2003 , 2006 .

[13]  Laurent Bricteux,et al.  A multiscale subgrid model for both free vortex flows and wall-bounded flows , 2009 .

[14]  L. Bricteux,et al.  Numerical prediction of turbulent heat transfer at low Prandtl number , 2012 .

[15]  B. A. Kader.  Temperature and concentration profiles in fully turbulent boundary layers , 1981 .

[16]  Laurent Bricteux,et al.  Direct and large eddy simulation of turbulent heat transfer at very low Prandtl number: Application to lead–bismuth flows , 2012 .

[17]  F. Nicoud,et al.  Subgrid-Scale Stress Modelling Based on the Square of the Velocity Gradient Tensor , 1999 .

[18]  Yuichi Matsuo,et al.  Surface heat-flux fluctuations in a turbulent channel flow up to Reτ=1020 with Pr=0.025 and 0.71 , 2004 .

[19]  G. Grötzbach.  Numerical simulation of turbulent temperature fluctuations in liquid metals , 1981 .

[20]  H. Anglart,et al.  Numerical modeling of forced-convection heat transfer to lead–bismuth eutectic flowing in vertical annuli , 2013 .

[21]  Xu Cheng,et al.  CFD analysis of thermal–hydraulic behavior of heavy liquid metals in sub-channels , 2006 .

[22]  C. Jayatilleke,et al.  The influence of Prandtl number and surface roughness on the resistance of the laminar sub-layer to momentum and heat transfer , 1966 .

[23]  F. Menter.  Two-equation eddy-viscosity turbulence models for engineering applications , 1994 .

[24]  B. Launder,et al.  Application of the energy-dissipation model of turbulence to the calculation of flow near a spinning disc , 1974 .

[25]  Sandro Manservisi,et al.  A CFD four parameter heat transfer turbulence model for engineering applications in heavy liquid metals , 2014 .

[26]  Oleg V. Vasilyev.  High Order Finite Difference Schemes on Non-uniform Meshes with Good Conservation Properties , 2000 .