Numerical Simulation of Sedimentation Behavior of Densely Arranged Particles in a Vertical Pipe Using Coupled SPH-DEM

This study develops a coupled Smoothed Particle Hydrodynamics (SPH) and the Discrete Element Method (DEM) framework to explore the sedimentation behavior of densely arranged particles in vertical pipes. An unresolved SPH-DEM model is proposed, which integrates porosity-dependent fluid governing equations through local averaging techniques to connect pore-scale interactions with macroscopic flow characteristics. Validated against single-particle settling experiments, the model accurately captures transient acceleration, drag equilibrium, and rebound dynamics. Systematic simulations reveal that particle number, arrangement patterns, and fluid domain geometry play critical roles in regulating collective settling: Increasing particle count induces nonlinear terminal velocity reduction. Systems of 16 particles show 50% lower velocity than single-particle cases due to enhanced shielding and energy dissipation. Particle configuration (compact layouts 4 × 8 vs. elongated arrangements 8 × 4) dictates hydrodynamic resistance, compact layouts facilitate faster settling by reducing cross-sectional blockage, while elongated arrangements amplify lateral resistance. The width of the fluid domain exerts threshold effects: narrow boundaries (0.03 m) intensify wall-induced drag and suppress vortices, whereas wider domains promote symmetric vortices that enhance stability. Additionally, critical transitions in multi-row/column systems are identified, where stress-chain redistribution and fluid-permeation thresholds govern particle detachment and velocity stratification. These findings deepen the understanding of granular–fluid interactions in confined spaces and provide a predictive tool for optimizing particle management in industrial processes such as wellbore cleaning and hydraulic fracturing.