Effect of Cladding Surface Roughness on Thermal-Hydraulic Response of Nuclear Fuel Rod of Advanced Gas-Cooled Reactor

Abstract

A thermal-hydraulic study of an isolated Advanced Gas-Cooled Reactor (AGR) nuclear fuel rod with smooth and rough cladding surfaces is carried out by computational fluid dynamics simulation and analytical calculation. Square transverse ribs of various pitch/height ratios (6:12) are considered for the rough surface. Parameters of the rough cladding surface show greater values than those for the smooth surface except for surface heat flux. It is found that only the average surface heat flux increases with an increasing pitch/height ratio. On the other hand, the average values of wall shear stress, Darcy friction factor, skin friction factor, convective heat transfer coefficient, Nusselt number, and thermal-hydraulic performance decrease with an increasing pitch/height ratio. The simulated results are found to be very close to the values obtained from an analytical calculation. Also, square and circular ribs are compared. The circular ribs show lower values of convective heat transfer coefficient and wall shear stress but permit high surface heat flux. The results of this study will help researchers comprehend the effect of cladding surface roughness on fuel rod thermal behavior.

Keywords:

• Thermal hydraulics
• rib
• pitch/height ratio
• computational fluid dynamics

Nomenclature

A =	=	cross-sectional area of the channel
$C_p$ =	=	specific heat
c =	=	skin friction factor
d =	=	diameter
dh =	=	hydraulic diameter, 4A/PT
f =	=	Darcy friction factor
h =	=	convective heat transfer coefficient
k =	=	height of rib
$\mathrm{k}$ =	=	thermal conductivity
LC =	=	characteristic length
Nu =	=	Nusselt number, $\frac{h{L}_C}{k}$
p =	=	pitch or distance between two ribs
Pr =	=	Prandtl number, $\frac{\mu C_p}{k}$
PT =	=	total wetted perimeter
Re =	=	Reynolds number, $\frac{\rho u_m{L}_C}{\mu }$
$T_b$ =	=	bulk fluid temperature
$T_s$ =	=	temperature of the surface
UT =	=	shear velocity, $\sqrt{{\tau }_w/\rho }$
$u_m$ =	=	mean velocity
y =	=	wall-normal distance
y+ =	=	normalized distance, $y\frac{\rho U_T}{\mu }$

Greek

$\alpha$ =	=	helix angle of rib
$\eta$ =	=	thermal-hydraulic performance
$\mu$ =	=	dynamic viscosity
$\rho$ =	=	density
${\tau }_w$ =	=	wall shear stress, $\frac{f}{4}\frac{\rho u_m^2}{2}$