An official website of the United States government
Abstract The susceptibility of a granular soil to suffusion is strongly dependent on its grain size distribution (GSD) and the mechanical and hydraulic conditions it is subjected to. This study investigates the onset of suffusion considering the effect of confining pressure and stress anisotropy using a fully resolved computational fluid dynamics and discrete element method (CFD–DEM). Three benchmarks, including the sedimentations of single and two adjacent spheres and the classic one‐dimensional (1D) consolidation are performed to demonstrate the capability of this method for high‐fidelity particle‐fluid simulations. A modified hydraulic criterion for the onset of suffusion considering stress anisotropy is presented. The microstructural changes of soil specimens before and during global suffusion are inspected, with emphasis on the evolutions of particle kinetic energy and displacements, force chain networks, and stress anisotropy. We found that the critical hydraulic gradient is negatively correlated with the confining pressure and the degree of stress anisotropy. Fine particles in the soil matrix are locally detached at small hydraulic gradients before the apparent global suffusion, as manifested by the variation of particle kinetic energy and coordination numbers. The roles of different contact types on force transmission and stress anisotropy in eroded specimens are also examined.
more »
« less
This content will become publicly available on July 1, 2026
CFD–DEM Analysis of Internal Soil Erosion Induced by Infiltration into Defective Buried Pipes
Internal soil erosion caused by water infiltration around defective buried pipes poses a significant threat to the long-term stability of underground infrastructures such as pipelines and highway culverts. This study employs a coupled computational fluid dynamics–discrete element method (CFD–DEM) framework to simulate the detachment, transport, and redistribution of soil particles under varying infiltration pressures and pipe defect geometries. Using ANSYS Fluent (CFD) and Rocky (DEM), the simulation resolves both the fluid flow field and granular particle dynamics, capturing erosion cavity formation, void evolution, and soil particle transport in three dimensions. The results reveal that increased infiltration pressure and defect size in the buried pipe significantly accelerate the process of erosion and sinkhole formation, leading to potentially unstable subsurface conditions. Visualization of particle migration, sinkhole development, and soil velocity distributions provides insight into the mechanisms driving localized failure. The findings highlight the importance of considering fluid–particle interactions and defect characteristics in the design and maintenance of buried structures, offering a predictive basis for assessing erosion risk and infrastructure vulnerability.
more »
« less
- Award ID(s):
- 2451951
- PAR ID:
- 10612302
- Publisher / Repository:
- MDPI
- Date Published:
- Journal Name:
- Geosciences
- Volume:
- 15
- Issue:
- 7
- ISSN:
- 2076-3263
- Page Range / eLocation ID:
- 253
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Rathje, Ellen; Montoya, Brina M.; Wayne, Mark H. (Ed.)Soil piping is the gradual and progressive erosion of soil grains, causing a void (open pipe) to form as water flows through the soil. In dam engineering, this type of internal erosion is often referred to as concentrated leak erosion and has been a cause of failure at multiple dams. Soil piping has also been observed in many landslides and contributes significantly to soil degradation in hillslopes and agricultural areas. Despite these many important impacts, there is still limited understanding of how soil pipes develop and progress and what factors control pipe stability. One of the significant challenges with analyzing soil piping, or concentrated leak erosion, is that it typically occurs in the vadose zone, where unsaturated conditions are present. However, most studies examining internal erosion have focused on saturated conditions, and few studies have examined the role unsaturated hydraulic properties (i.e., air entry value, matric suction, etc.) may play in the likelihood of internal erosion. Consequently, this study aims to explore the mechanisms controlling the erosion rate within soil pipes from the perspective of unsaturated soil mechanics. Bench-scale experiments were performed to examine the formation and progression of an eroded pipe in a small slope constructed at different water contents. Soil samples were also tested to measure its unsaturated hydraulic properties. The results show that the likelihood of pipe formation varies with the moisture content and, therefore, suction in the soil, as does the potential for pipe collapse. This demonstrates that unsaturated soil properties are key to understanding the formation and progression of piping in slopes.more » « less
-
Soil piping (concentrated leak erosion) is a major contributor to soil erosion in many parts of the world, and collapse of eroded pipes can result in the formation of gullies and sinkholes or trigger slope instability. Despite these significant impacts, there is little understanding of factors controlling pipe collapse, and how water within the pipe influences moisture levels within a slope. In this study, physical models were employed on unsaturated model slopes with pre-formed macropores to investigate how soil properties, pipe characteristics, and hydraulic conditions govern internal erosion processes and slope stability. Experiments simulated shallow field conditions (0.45 m overburden) using 4 mm and 12 mm pipes to establish preferential flow paths, while varying model parameters including initial compaction moisture content and density, pipe condition (absent, closed, or open), slope angle, and model width. Volumetric water content sensors monitored moisture evolution, while cameras captured slope responses to subsurface flow. Results demonstrate that initial compaction conditions (water content and density), pipe size, hydraulic connectivity, and pipe condition control internal erosion processes and slope stability.more » « less
-
Abstract Sinkholes are the most abundant surface features in karst areas worldwide. Understanding sinkhole occurrences and characteristics is critical for studying karst aquifers and mitigating sinkhole‐related hazards. Most sinkholes appear on the land surface as depressions or cover collapses and are commonly mapped from elevation data, such as digital elevation models (DEMs). Existing methods for identifying sinkholes from DEMs often require two steps: locating surface depressions and separating sinkholes from non‐sinkhole depressions. In this study, we explored deep learning to directly identify sinkholes from DEM data and aerial imagery. A key contribution of our study is an evaluation of various ways of integrating these two types of raster data. We used an image segmentation model, U‐Net, to locate sinkholes. We trained separate U‐Net models based on four input images of elevation data: a DEM image, a slope image, a DEM gradient image, and a DEM‐shaded relief image. Three normalization techniques (Global, Gaussian, and Instance) were applied to improve the model performance. Model results suggest that deep learning is a viable method to identify sinkholes directly from the images of elevation data. In particular, DEM gradient data provided the best input for U‐net image segmentation models to locate sinkholes. The model using the DEM gradient image with Gaussian normalization achieved the best performance with a sinkhole intersection‐over‐union (IoU) of 45.38% on the unseen test set. Aerial images, however, were not useful in training deep learning models for sinkholes as the models using an aerial image as input achieved sinkhole IoUs below 3%.more » « less
-
The addition of short carbon fibers to the feedstock of large-scale polymer extrusion/deposition additive manufacturing results in significant increases in mechanical properties dependent on the fiber distribution and orientation in the beads. In order to analyze those factors, a coupled computational fluid dynamics (CFD) and discrete element modeling (DEM) approach is developed to simulate the behavior of fibers in an extrusion/deposition nozzle flow after calibrations in simple shear flows. The DEM model uses bonded discrete particles to make up flexible and breakable fibers that are first calibrated to match Jeffery’s orbit and to produce interactions that are consistent with Advani-Tucker orientation tensor predictions. The DEM/CFD model is then used to simulate the processing of fiber suspensions in the variable flow and geometries present in extrusion/deposition nozzles. The computed results provide enhanced insight into the evolution of fiber orientation and distribution during extrusion/deposition as compared to existing models through individual fiber tracking over time and space on multiple parameters of interest such as orientation, flexure, and contact forces.more » « less
