An Eulerian multiphase approach based on kinetic theory of granular flow was used to simulate flow of coal-water slurries (CWS) in horizontal pipelines. The RNG k-ε turbulent model was incorporated in the governing equation to model turbulent two-phase flow with strong particle-particle interactions. In this model, the coal particles with bimodal distribution were considered as two solid-phase components, and the moment exchange between solid and liquid as well as that between solid and solid were accounted for. The model was firstly validated with pressure gradient and concentration profile data from the open literature, and then validated with pressure gradient data of the authors’ experiments. The effects of influx velocity, total influx concentration and grain composition were numerically investigated, and the results have displayed some important slurry flow characteristics, such as constituent particle concentration distribution and velocity distribution as well as pressure gradients, which are very difficult to display in the experiments. The results suggest that both gravity difference between large and small particles and strong particleparticle interaction had significant effects on concentration distribution as well as velocity distribution.