Abstract: To solve the strong coupling between landslide body and water, and the complex fluid-solid coupling problem caused by the evolution of nonlinear water waves, the smoothed particle flow method (SPH) was used to simulate the hydrodynamic characteristics of water, and the discrete element method (DEM) was used to simulate the large deformation and discontinuity of landslides.Through constructing SPH-DEM fluid-solid coupling algorithm, and developing a parallel computing program based on distributed MPI, the large-scale high-performance calculation of the SPH-DEM coupling process was realized.This program was applied to the simulation of Bairiya landslide surge in Jiangpinghe Reservoir of Hubei Province.The disaster chain process of landslide development, surge formation, propagation along river channel and impact on the dam was reproduced.The results showed that the maximum surge height in front of the dam was 474.00 m, which was lower than the dam elevation of 476.00 m.This program achieved more than 70% parallel acceleration efficiency when the processes number was 128.The research results can effectively analyze the dynamic process of landslide surge disaster, and the developed parallel program can provide reference for other domestic similar software.