The turbulent heat transfer to low Prandtl number fluid flow through rod bundles is analyzed using k-$varepsilon$ two-equation model. For the prediction of the turbulent flow field, an anisotropic eddy viscosity model is used. In the analysis of the temperature field, the effects of various parameters such as geometry, Reynolds and Prandtl numbers are considered. The calculation in made for Prandtl numbers from 0.001 to 0.1 in order to analyze the heat transfer to low Prandtl number fluid such as liquid metals. The numerical results show that for small P/D (Pitch/Diameter) geometries low Prandtl number makes severe changes of the rod surface temperature.