The problem with simultaneous conductive, convective and radiative heat transfer is analyzed in two dimensional gradual expansion nozzle flow. Main objective is to generate a program to attack the radiation in curvilinear coordinates with momentum and energy equations involved. The governing equations are transformed into generalized coordinates. Whereas the SIMPLE algorithm developed in a body-fitted coordinate system is used to obtain flowfield solution, the finite volume method for radiation is applied to solve the radiative transfer equation. The accuracy and performance of the finite volume method for radiation are validated in a well-known combined heat transfer problem by using curvilinear grid system. Then, detailed thermal characteristics are investigated in a gradually expanding nozzle through such various parameters as Prandtl number, conduction to radiation parameter, wall emissivity and scattering albedo. Results show that the radiation plays a significant role in the gradually expanding nozzle flow.