In this paper, results from a study testing a DPD based simulation technique for modeling suspensions were presented. It is concluded that DPD can be used as an alternative computational tool for modeling a fairly wide variety of suspensions. Without significant modification, the method recovers well established predictions concerning the flow of a suspension for volume fractions in the dilute to semidilute systems regime. At higher volume fractions and Pe, modifications, such as a variable time step and explicit inclusion of lubrication forces are necessary to account for important phenomena that must be resolved at small time and length scales. This should also be true for other approaches like numerical solution of the Navier–Stokes equations at a similar resolution. Indeed, similar modifications have been employed for a lattice Boltzmann [Nguyen and Ladd (2002)] based model of suspensions. Further validation would be useful, e.g., testing how well the model describes sedimentation.