Axisymmetric numerical thermal analysis for a 3-dimensional regenerative cooling system in a rocket engine is carried out. To predict the accurate heat transfer with the stiff temperature distribution, several tests have been conducted for the grid size, the properties variation of the coolant and the combustion gas depending on temperature. The axisymmetric heat flux model is defined using fin efficiencies and is designed to be equivalent to the heat flux of the 3-dimensional coolant channel. For comparison purpose, the 1-dimensional analysis using Bartz equation is also conducted. The performance of the present model in predicting the cooling characteristics of a 3-dimensional regenerative cooling system is compared with the 3-dimensional results of RTE(Rocket Thermal Evaluation). It is found that the present method predicts much closer results to those of RTE code than 1-dimensional analysis.