Development of program for evaluating the temperature of Zr-U metallic fuel rod

A code for evaluating the temperature of Zr-U metallic rod has been developed. Finite element (FE) method is adopted for the developed code sharing the user subroutines which has been prepared for the ABAQUS commercial FE code. The developed program for the Zr-U metallic fuel rod corresponds to a nonlinear transient heat transfer problem, and uses a sparse matrix solver for FE equations during iterations at every time step. The verifications of the developed program were conducted using the ABAQUS code. Steady state and transient problems were analyzed for 1/8 rod model due to the symmetry of the fuel rod and full model. From the evaluation of temperature for the 1/8 rod model at steady state, maximal error of 0.18 % was present relative to the ABAQUS result. Analysis for the transient problem using the fuel rod model resulted in the same as the variation of centerline temperature from the ABAQUS code during a hypothetical power transient. The distribution of heat flux for the entire cross section and surface was almost identical for the two codes