Abstract: By solving the Navier-Stokes equations and a Eularian-Lagrangian multiphase model, three-dimensional simulations of solid-liquid two-phases flows in centrifugal pumps were made by using two-way coupling of the two phases and modifying the momentum equations of discrete solid phase with empirical formulas. The Alhert erosion model combined with the Mclaury wear test results was adopted to calculate the erosion of flow components, and the flow features and variation in pump external characteristics were analyzed for working conditions of different flow rates and different volume concentrations of solid phase in the pump inflow. Calculations indicate that the volume concentration has strong influence on pressure distribution and streamlines. With the inflow volume concentration increasing, average pressure on both sides of the impeller blade is increased while pressure difference between blade inlet and outlet is gradually reduced. Existence of solid particles suppresses flow turbulences and consequently weakens the second flows in impeller passages, leading to even an increase in pumping efficiency in certain working conditions. Under small flow condition, a swirl is generated in the flow before its entering the impeller inlet, and solid particles under centrifugal force are impacting and cutting the solid surfaces, which causes wear erosion of impeller ring. An increase in the inflow volume concentration of solid phase will increase its volume fraction distributed over flow components and hence aggravate their wear erosion.