Conjugate heat transfer of transient thermal conduction,natural convection and thermal radiation with internal heat sources is simulated numerically using PHOENICS(Version 1.4) with the false density method.The flow field and temperature distribution are calculated using both the body fitted coordinates and Cartesian coordinates.The variation of the maximum temperature in the substance with internal heat sources is presented as function of the environmental temperature.