This paper formulates a new numerical model for the dense gas-solid flows, combining Eulerian/Eulerian approach and Eulerian/Lagrangian approach together. The model takes account of inter-particle interactions in Eulerian coordinate using the kinetic theory of granular flow based on the Chapman-Enskog theory of dense gases. A stochastic particle dispersion model is incorporated in the model to calculate the turbulence diffusion of particles in Lagrangian coordinate. A comparison with published experimental results is made, and good agreement is shown.