This paper presents the vector strsam function-vorticity (Ψ-ω) methed for computing gas-particles 3-dimensional turbulent reacting flow: the gas-phase governing equations and the related turbulence models, the properties of vector stream function, the results obtained from computation of pulverized coal combustors and industrial boiler, and discussions on a lot of work for further research.