kg-εg-kp-εp-kpg-θ turbulence model is proposed which considers particle-particle collision and gas-particle turbulence. This model includes turbulent kinetic energy equation, turbulent kinetic energy dissipation rate equation, particle pseudo-temperature transportation equation and the two-phase velocity correlation transport equation. To close the turbulence model, algebraic expressions of two-phase Reynolds stresses and two-phase velocity correlation variable are established which considered both gas-particle interaction and anisotropy. This model is used to simulate gas-particle in swirling sudden-expansion chamber. Comparing with kg-εg-kp-εp-θ model which is simply closed using a semi-empirical dimensional analysis, the present model has better predicted capability. It is shown that the present model gives simulation results in much better agreement with the experimental results than the kg-εg-kp-εp-θ model.