To better understand the dynamic process of rock avalanches blocking rivers, a novel numerical approach based on the coupled Eulerian-finite-discrete element method (CEFDEM) is proposed. The Samaoding paleolandslide blocking river event, which occurred at the upstream of the Jinsha River was used as a case study to further validate the new numerical approach. Field investigations, thermoluminescence dating, and geomorphological analysis were conducted to determine basic geological conditions and provide data for the numerical simula-tions. Then a calibrated 3D landslide blocking river simulation based on the CEFDEM was conducted. The landslide blocking river lasted for 70 s and landslide scale from the numerical simulation agrees well with the field investigation. The maximum overall feature speed of the entire sliding mass is 35 m/s, while part of the sliding mass can reach 69 m/s. Dynamic fragmentation of the rock slide is stratified such that the bottom of the sliding body has higher fragmentation degree than the top. The variation of kinetic energy, accumulated friction dissipation, and fracture energy of the sliding mass are also shown. The impulsive water wave is triggered immediately after sliding mass runs into river, and its maximum height is 132 m, while part of the wave can reach a speed of 64 m/s. The river water will be pushed by the subsequent sliding mass movement. A comparison of the CEFDEM model and particle flow code (PFC) model in landslide blocking river was conducted, and the advantages and limitations of CEFDEM model were discussed in detail.