Road Collapse Induced by Pipeline Leakage in Water-Rich Sand: Experiments and Computational Fluid Dynamics-Discrete Element Method Simulations

To investigate the mechanism of road collapse induced by structural defects in underground drainage/sewerage pipelines in water-rich sands, laboratory physical model tests were conducted to reproduce the macroscopic development of surface subsidence. A computational fluid dynamics-discrete element method (CFD-DEM) model was then established and validated against the tests to assess its reliability. Using the validated model, we examined the effects of defect size and groundwater level on the progression of groundwater-ingress-driven internal erosion and tracked the evolution of vertical stress and intergranular contacts around the pipe. Results show that internal erosion proceeds through three stages—initial erosion, slow settlement, and collapse—culminating in an inverted-cone collapse pit. After leakage onset, the vertical stress in the surrounding soil exhibits a short-lived surge followed by a decline on both sides above the pipe. The number of intergranular contacts decreases markedly; erosion propagates preferentially in the horizontal direction, where the reduction in contacts is most pronounced. Within the explored range, higher groundwater levels and larger defects accelerate surface settlement and yield deeper and wider collapse pits. Meanwhile, soil anisotropy strengthens with increasing groundwater level but peaks and then slightly relaxes as defect size grows. These qualitative findings improve understanding of the leakage-induced failure mechanism of buried pipelines and offer references for discussions on monitoring, early warning, and risk awareness of road collapses.