Boulders mixed with debris flows roll downstream under interactions with debris flow slurry, which poses a great threat to the people, houses, bridges, and other infrastructure encountered during their movement. The catastrophic debris flow in Zhouqu County, which occurred on 7 August 2010, was used as an example to study the motion and accumulation characteristics of boulders in debris flows. In this study, a fluid–solid coupling model utilizing the general moving objects collision model and the renormalization group turbulent model was used in the FLOW-3D software, treating boulders with different shapes in the Zhouqu debris flow as rigid bodies and the debris flow as a viscous flow. Numerical simulation results show that this method can be used to determine the motion parameters of boulders submerged in debris flows at different times, such as the centroid velocity, angular velocity, kinetic energy, and motion coordinates. The research method employed herein can provide a reference for studying debris flow movement mechanisms, impact force calculations, and aid in designing engineering control structures.