Fluid–structure interaction (FSI) research is crucial for applications in fields such as naval engineering, geological hazards, and biomechanics. Traditional grid-based methods (such as CFD) often face challenges in simulating large-deformation flow fields and complex boundary conditions, where mesh distortion can compromise simulation accuracy. Building upon the DualSPHysics5.2 framework, this study leverages the strengths of weakly compressible SPH (WCSPH) in modeling free surface flows and large-deformation fluids, as well as the discrete element method (DEM), for accurately describing particle collisions and fragmentation behaviors. We propose an improved MSPH-DEM coupling algorithm that incorporates moving least squares (MLS) correction for kernel function gradient optimization. This algorithm utilizes MLS-based gradient correction to achieve smoother fluid surfaces as well as bidirectional coupling between fluids and particles. Experimental validation demonstrates that in dam break simulations, this method reduces pressure errors. In the dam break impacting a cube experiment, it enhances accuracy, while in the dam break impacting a baffle experiment, the horizontal displacement of marker points closely aligns with the experimental values from Liao et al. This approach effectively improves the accuracy of the simulations of FSI problems, offering a more reliable numerical simulation methodology for engineering applications such as geological hazard prevention.
Zou et al. (2026) studied this question.