Abstract A two fluid model is used to predict gas–solid particle turbulent dilute flows around an obstacle. The turbulence gas phase is treated with the Chen–Kim modified k– ε model. Two turbulence models are applied to the particulate phase: the A p model or the Hinze–Tchen's algebraic model, and the k p model. The first model assumes that a local equilibrium exists between the fluctuating motion of the particles and the gas turbulence, while the second model describes the particle turbulent properties considering the effect of convection, diffusion, production and destruction of the particle turbulence via a k p equation, where k p is the kinetic energy of the particles. These models are implanted in the CFD code PHOENICS. The 2D experiment of Fessler and Eaton (flow around a vertical backward-facing step), and the 3D experiment made at C.S.T.B. (flow around a cube) are used as references for validation.