Effects of slip velocity and volume fraction of slip spheres on the momentum transfer characteristics of assemblages of slip spheres are numerically investigated. The fluid slip along the surface of the sphere is considered by Navier's linear slip model. The dimensionless governing continuity and momentum equations are solved using a semi-implicit marker and cell method implemented on a staggered grid arrangement in spherical coordinates. The convection and viscous terms of momentum equations are discretized by means of the QUICK scheme and a second-order central differencing scheme, respectively. The present numerical solver is benchmarked via grid independence and comparisons with the existing literature values. Results were obtained over a wide range of pertinent dimensionless numbers such as the Reynolds number, volume fraction of the dispersed phase, and dimensionless slip parameter.