The method of two dimensional Reynolds Averaged Navier-Stokes (RANS) equations has been employed for the simulation of turbulent particulate flow. This approach was fitted with appropriate closure equations that take into account all the pertinent forces and effects on the solid particles, such as: particle-turbulence interactions: turbulence modulation; particle-particle interactions; particle-wall interactions: gravitation, viscous drag and lift forces. The flow domain in all cases was a cylindrical pipe and the computations were carried for upward pipe flow. The finite volume technique was used for the numerical solution of the governing and closure equations. The results show the effect of loading and particle density on the profiles of the velocity. the turbulence intensity and the solids concentration. (C) 2009 Elsevier B.V. All rights reserved.