CFD analysis of thermal-hydraulic behavior of heavy liquid metals in sub-channels

CFD analysis was carried out for thermal-hydraulic behavior of heavy liquid metal flows, especially lead-bismuth eutectic, in sub-channels of both triangular and square lattices. Effect of various parameters, e.g. turbulence models and pitch-to-diameter ratio, on the thermal-hydraulic behavior was investigated. Among the turbulence models selected, only the second order closure turbulence models reproduce the secondary flow. For the entire parameter range studied in this paper, the amplitude of the secondary flow is less than 1% of the mean flow. A strong anisotropic behavior of turbulence is observed. The turbulence behavior is similar in both triangular and square lattices. The average amplitude of the turbulent velocity fluctuation across the gap is about half of the shear velocity. It is only weakly dependent on Reynolds number and pitch-to-diameter ratio. A strong circumferential non-uniformity of heat transfer is observed in tight rod bundles, especially in square lattices. Related to the overall average Nusselt number, CFD codes give similar results for both triangular and square rod bundles. Comparison of the CFD results with bundle test data in mercury indicates that the turbulent Prandt1 number for HLM flows in rod bundles is close to 1.0 at high Peclet number conditions, and increases by decreasing Peclet number. Based on the present results, the SSG Reynolds stress model with semi-fine mesh structures (y(1)(+) < 15) is recommended for the application of HLM flows in rod bundle geometries. (c) 2006 Elsevier B.V. All rights reserved.