The Direct Quadrature Method of Moments (DQMOM) has been presented for the solution of population balance equation in the wide range of the multi-phase flows. This method has the inherently interesting features which can be easily applied to the multi-inner variable equation. In addition, DQMOM is capable of easily coupling the gas phase with the discrete phases while it requires the relatively low computational cost. Soot inception, subsequent aggregation, surface growth and oxidation are described through a population balance model solved with the DQMOM for soot formation. This approach is also able to represent the evolution of the soot particle size distribution. The turbulence-chemistry interaction is represented by the laminar flamelet model together with the presumed PDF approach and the spherical harmonic P-1 approximation is adopted to account for the radiative heat transfer.