Numerical analysis on turbulent flow and heat transfer in the rocket nozzle has been studied using the mass-weighted-averaged full Navier-Stokes equations, the Morkovin hypothesis on turbulent flow, the $ extsc{k}$-$varepsilon$ turbulence model with the wall function specially designed to be able to consider the effects of pressure gradients, heat transfer and compressibility, and the numerical scheme of Karki. The present results are in good agreement with the experiments of Back et al.