CFD–DEM Modelling of Ground Collapse Induced by Underground Pipeline Leakage in Water-Rich Sand Layers

Urban underground pipeline aging and leakage can result in soil erosion and ground collapse, constituting a major threat to urban public safety. To investigate this disaster mechanism, this present study established a two-dimensional numerical model based on the computational fluid dynamics–discrete element method (CFD–DEM) two-way fluid–solid coupling approach, simulating and reproducing the entire process from soil erosion, soil arch evolution to ground collapse caused by underground pipeline leakage in water-rich sand layers. The simulation shows that under the action of seepage pressures, soil particles are eroded and lost, forming a cavity above the pipeline defect. As soil continues to be lost, the disturbed zone expands toward the ground surface, causing ground settlement, and in water-rich sand layers, a funnel-shaped sinkhole is eventually formed. The ground collapse process is closely related to the groundwater level and the thickness of the overlying soil layer above the pipeline. Rising groundwater levels reduce the effective stress and shear strength of the soil, significantly exacerbating seepage erosion. Increasing the thickness of the overlying soil layer can enhance the confining pressure, improve soil compactness, and promote the formation of soil stress arch, thereby effectively slowing down the rate of ground collapse. This study reproduces the process of ground collapse numerically and reveals the mechanism of ground collapse induced by underground pipeline leakage in water-rich sand layers.