A direct modeling of colloidal suspensions consists of calculating trajectories of all suspended objects. Due to the large time computing and the large cost involved in such calculations, we consider in this paper another route. Colloidal suspensions are described on a mesoscopic level by a distribution function whose time evolution is governed by a Fokker-Planck-like equation. The difficulty encountered on this route is the high dimensionality of the space in which the distribution function is defined. A novel strategy is used to solve numerically the Fokker-Planck equation circumventing the curse of dimensionality issue. Rheological and morphological predictions of the model that includes both direct and hydrodynamic interactions are presented in different flows.