The present work is numerical simulat ion results from the modeling of a non-premixed natural gas flame performed in a co mbustor model. CFD studies using FLUENT code were carried out for d ifferent air swirl numbers and inlet thermal load with constant A/F. The isothermal flo w field was simulated using (SST) k-ω turbulence model and the reacting flow was simulated by the non-premixed combustion model with the P-1 radiation model available in the computational fluid dynamics package Fluent 6.3. The model geometry was created and meshing arrangement was generated using Gamb it pre-processing software. The domain of the model was based on the dimension of the combustor and burners. The case studied is a cylindrical enclosure of 0.1 m radius and 1.0 m length. Two reactant streams emerge fro m two separate coaxial jets producing a swirling diffusion flame. The natural gas is issued axially into the combustor fro m the annulus area between the swirler outer diameter and the burner hub diameter whereas the combustion air is introduced through an annular swirler having uniform swirl vanes. The results show a reasonable agreement of the measured and the calculated reverse flow zone sizes using the shear stress transport (SST) k-ω model. The boundary of the reverse flo w zone for weak air swirl nu mber of 0.5 is formed co mpletely inside the co mbustor with closed loop, while for air swirl numbers of 0.87 and 1.5 the boundaries fill the co mbustor and its size increased as the air swirl nu mber increased. Increasing the air swirl number, the high temperature regions size, the flame length, and the CO 2 and CO concentrations decreased while the O 2 concentration increased. Increasing the thermal load, the high temperature reg ions size, the flame length, CO 2 and CO concentrations increased, while the O 2 concentration decreased.