Natural convection heat transfer from vertical rod bundles in liquid sodium

Natural convection heat transfer from vertical rod bundles in liquid sodium was numerically analyzed for three types of the bundle geometry (two parallel, equilateral triangle and equilateral square arrays). The unsteady laminar three dimensional basic equations for natural convection heat transfer caused by a step heat flux were numerically solved until the solution reaches a steady-state. The PHOENICS code was used for the calculation considering the temperature dependence of thermo-physical properties concerned. The 2 to 4 test rods for diameter (D=7.6 mm), heated length (L=200 mm) and L/d (=26.32) were used in this work. The surface heat fluxes for each cylinder were equally given for a modified Rayleigh number, R-f,R-L, ranging from 3.06x 10(4) to 3.14x 10(7) (q=1x10(4) to 7x 10(6) W/m(2)) in liquid temperature (T-L=673.15 K). The values of S/D, which are ratio of the diameter of flow channel for bundle geometry to the rod diameter, for the rod bundle were ranged from 1.4 to 3 on each bundle geometry. The spatial distributions of local and average Nusselt numbers, Nu(theta),(z) and (Nu(av,B))(N), on vertical rods of a bundle were clarified. The values of average Nusselt number, (Nu(av,B))(N,S/D), for three types of the bundle geometry with various values of S/D were calculated to examine the effect of the bundle geometry, S/D and R-f,R-L on heat transfer. The bundle geometry for the higher (Nu(av,B))(N) value under the condition of S/D=constant was examined. The correlation for (Nu(av,B))(N,S/)(D) for three types of bundle geometry above mentioned including the effects of R-f,R-L and S/D were developed. The correlations can describe the theoretical values of (Nu(av,B))(N,S/D) for three types of the bundle geometry for S/D ranging from 1.4 to 3 within -7.44 to 10.73 % difference.