The motivation of this study is to investigate the turbulence-chemistry interactions by using probability density function （PDF） method. A consistent hybrid Reynolds Averaged Navier-Stokes （RANS）/PDF method is used to simulate the turbulent non-reacting and reacting flows. The joint fluctuating velocity-frequency-composition PDF equation coupled with the Reynolds averaged density, momentum and energy equations are solved on unstructured meshes by the Lagrangian Monte Carlo （MC） method combined with the finite volume （FV） method. The simulation of the axisymmetric bluff body stabilized non-reacting flow fields is presented in this paper. The calculated length of the recirculation zone is in good agreement with the experimental data. Moreover, the significant change of the flow pattern with the increase of the jet-to-coflow momentum flux ratio is well predicted. In addition, comparisons are made between the joint PDF model and two different Reynolds stress models.
