A hybrid scheme coupling the discrete element method （DEM） with the computational fluid dynamics （CFD） is developed to model solid-liquid flows. Instead of solving the pressure Poisson equation, we use the compressible volume-averaged continuity and momentum equations with an isothermal stiff equation of state for the liquid phase in our CFD scheme. The motion of the solid phase is obtained by using the DEM, in which the particle-particle and particle-wall interactions are modelled by using the theoretical contact mechanics. The two phases are coupled through the Newton＇s third law of motion. To verify the proposed method, the sedi-mentation of a single spherical particle is simulated in water, and the results are compared with experimental results reported in the literature. In addition, the drafting, kissing, and tumbling （DKT） phenomenon between two particles in a liquid is modelled and rea-sonable results are obtained. Finally, the numerical simulation of the density-driven segregation of a binary particulate suspension in-volving 10 000 particles in a closed container is conducted to show that the presented method is potentially powerful to simulate real particulate flows with large number of moving particles.
