Submarine debris flows and their generated waves are common disasters in Nature that may destroy offshore infrastructure and cause fatalities. As the propagation of submarine debris flows is complex, involving granular material sliding and wave generation, it is difficult to simulate the process using conventional numerical models. In this study, a numerical model based on the smoothed particle hydrodynamics (SPH) algorithm is proposed to simulate the propagation of submarine debris flow and predict its generated waves. This model contains the Bingham fluid model for granular material, the Newtonian fluid model for the ambient water, and a multiphase granular flow algorithm. Moreover, a boundary treatment technique is applied to consider the repulsive force from the solid boundary. Underwater rigid block slide and underwater sand flow were simulated as numerical examples to verify the proposed SPH model. The computed wave profiles were compared with the observed results recorded in references. The good agreement between the numerical results and experimental data indicates the stability and accuracy of the proposed SPH model. © 2021 by the authors. Licensee MDPI, Basel, Switzerland.