Abstract
About 50% of the world's dam failures are triggered by piping, which is a primary cause of embankment breaks. The phenomenon of piping results from the erosion of soil particles and their transport within a soil mass. In this paper, a numerical method is proposed to analyze the erosion within soils and the transport of eroded soil particles by adopting the concept of the erosion rate of soils. In such an analysis, the saturated-unsaturated seepage flow of the pore liquid, the detachment of the soil particles from the soil fabric, and the migration of the eroded particles are taken into consideration, and the equations related to the conservation of the pore liquid and the eroded soil particles are numerically solved. This numerical simulation allows for the procurement of the temporal alteration and the spatial distribution of the porosity, the particle size distribution, and the concentration of the detached soil particles in the pore liquid, as well as the distribution of pore liquid pressure. The results have revealed that the method can reproduce the experimental data from previous studies on the internal erosion of soils and that it qualitatively predicts, from the numerical experiments, the typical development of piping within soils, such as soil blocks and embankments, as the solutions to the initial and the boundary value problems of the governing equations.
