In this paper a calculation procedure for simulating the combustion flow in the afterburner with the heat shield, flame stabilizer and the contracting nozzle is described and evaluated by comparison with experimental data. The modified two-equationk – ɛ model is employed to consider the turbulence effects, and thek – ɛ –g turbulent combustion model is used to determine the reaction rate. To take into account the influence of heat radiation on gas temperature distribution, heat flux model is applied to predictions of heat flux distributions. The solution domain spanned the entire region between centerline and afterburner wall, with the heat shield represented as a blockage to the mesh. The enthalpy equation and wall boundary of the heat shield require special handling for two passages in the afterburner. In order to make the computer program suitable to engineering applications, a subregional scheme is developed for calculating flow fields of complex geometries. The computational grids employed are 100×100 and 333×100 (non-uniformly distributed). The numerical results are compared with experimental data. Agreement between predictions and measurements shows that the numerical method and the computational program used in the study are fairly reasonable and appropriate for primary design of the afterburner.