Solid particle erosion is a serious issue in centrifugal pumps that may result in economic losses. Erosion prediction in centrifugal pump is complex because the flow field inside it is three-dimensional (3D) unsteady and erosion can be affected by numerous factors. In this study, solid particle erosion of the entire centrifugal pump for liquid–solid flow is investigated numerically. Two-way coupled Eulerian–Lagrangian approach is adopted to calculate the liquid–solid interaction. The reflection model proposed by Grant and Tabakoff and the erosion model proposed by the Erosion/Corrosion Research Center are combined to calculate the erosion rate and predict the erosion pattern. Results show that for the baseline case, the inlet pipe is the least eroded component, whereas the impeller is the most eroded component. The highest average and maximum erosion rates occur at the hub of impeller. The most severe erosion region of a blade is the leading edge with a curvature angle that varies from 55 deg to 60 deg. The most severe erosion region of a volute is in the vicinity of a curvature angle of 270 deg. The impeller erosion pattern, especially the middle part of the hub and the vicinity of the blade pressure side, can be greatly influenced by operation parameters, such as flow rate, particle concentration, and particle size.