Liquid metal turbulent heat transfer in a circular tube

It is well-known that the turbulent Prandtl number, Pr/sub t/, must be taken higher than unity in liquid metal turbulent heat transfer, but what is the correct relationship between the turbulent Prandtl number and the fluid Prandtl number, Pr, or the Peclet number, Pe, is still uncertain. A simple closure of v-barΘ-bar equation has been proposed. The closed form of the equation was tested with a local equilibrium assumption and by substituting experimental data into hydrodynamic turbulent quantities appearing in the equation presented. The tested equation is promising but cannot be counted as a prediction because the experimental data were used. In this article, the hydrodynamic turbulent quantities are calculated employing a full Reynolds stress model. Three types of full Reynolds stress models are tested: the model by Launder, Reece and Rodi, the model by Pope and Whitelaw, and one by Prud'homme and Elgobashi. All the governing equations are solved numerically with a finite difference scheme. The computation is executed for a fully developed state of liquid metal flow in a circular tube heated at uniform heat flux