This research investigates a numerical simulation of swirling turbulent non-premixed combustion. The effects on the combustion characteristics are examined with three turbulence models:namely as the Reynolds stress model, spectral turbulence analysis and Re-Normalization Group. In addition, the P-1 and discrete ordinate (DO) models are used to simulate the radiative heat transfer in this model. The goving equations associated with the required boundary conditions are solved using the numerical model. The accuracy of this model is val-idated with the published experimental data and the comparison elucidates that there is a reasonable agreement between the obtained values from this model and the corresponding experimental quantities. Among different models proposed in this research, the Reynolds stress model with the Probability Density Function (PDF) ap-proach is more accurate (nearly up to 50％) than other turbulent models for a swirling flowfield. Regarding the effect of radiative heat transfer model, it is observed that the discrete ordinate model is more precise than the P-1 model in anticipating the experimental behavior. This model is able to simulate the subcritical nature of the isothermal flow as well as the size and shape of the intal recirculation induced by the swirl due to combustion.