Internal erosion induces alterations in the initial microstructure of soils, simultaneously affecting physical, hydraulic, and mechanical properties. The initial soil composition plays a crucial role in governing the initiation and progression of seepage-induced suffusion. This study employs the controlled variable method to develop granular soil models with varying particle size ratios, initial fine particle contents, and coarse particle shapes. Seepage suffusion simulations coupled with microstructural analyses are conducted using the CFD-DEM approach. Results demonstrate that particle size ratio, fine particle content, and coarse particle shape exert distinct influences on cumulative erosion mass, fine particle distribution, contact fabric, and mechanical redundancy at both macroscopic and microscopic scales. This numerical investigation advances the fundamental understanding of internal erosion mechanisms and informs the development of micro-mechanical constitutive models. Furthermore, for binary granular media composed of coarse and fine particles, careful control of the particle size ratio and fine content is recommended when utilizing gap-graded soils in embankment and dam construction to improve structural resilience and resistance to internal erosion.

Hydraulic structures such as embankments and dams are essential for water storages, flood control, and transportation, but are vulnerable to suffusion under complex loading conditions. This study investigates the effect of suffusion on the cyclic shear behavior of gap-graded soils using the coupled computational fluid dynamics and discrete element method (CFD-DEM). A series of seepage infiltration and drained cyclic shear tests are conducted on specimens with varying mean stresses and initial stress anisotropy to systematically evaluate the mechanical consequences of suffusion. The findings reveal that the higher mean stress and initial stress anisotropy significantly exacerbate fines loss and deformation, particularly along principal seepage directions during suffusion. Furthermore, the eroded specimens exhibit substantial stiffness degradation and microstructural changes, including the deteriorated interparticle contacts and more pronounced fabric anisotropy. Notably, fines loss intensifies the load-bearing reliance on coarse particles during cyclic loading. These results provide new micromechanical insights into suffusion-induced degradation, offering valuable implications for developing advanced constitutive model of gap-graded soils accounting for suffusion-induced fines loss and cyclic loading conditions.

The complex phenomenon of suffusion is the selective erosion of the fine fraction under the effect of seepage flow within the matrix of coarser particles. Three processes are involved simultaneously: detachment, transport, and partial filtration of the fine particles. With the objective to characterize the influence of the stress state on suffusion-related parameters, downward seepage flow tests were conducted under hydraulic-gradient controlled conditions. Four stress states are investigated: triaxial isotropic, triaxial compression, triaxial extension and rigid vertical boundaries. Also, four different cohesionless gap-graded soils were tested, from underfilled to overfilled microstructures. The entire erosion process can be divided into four phases: onset, self-filtration, blow-out and steady state. The definitions of several suffusion-related parameters are given for each suffusion phase, in terms of hydraulic gradient, hydraulic conductivity variation, cumulative expended energy, erosion resistance index and Darcy velocity. The results demonstrate that the suffusion kinetics of soils in transition between underfilled and overfilled microstructures are more affected by the stress state than others.

Suffusion, a process whereby water gradually carries away fine particles from soil, is thought to be one of the possible reasons for the settlement or inclination of bridge piers after a major flood (delayed displacement). The aim of this study is to offer fresh insights into suffusion and its mechanical impact on the affected soil, with a specific focus on how it relates to bridge pier failures. Riverbed material replicated with relatively larger fine particles than those used in past studies which focused on soil in embankments or dikes. Through both monotonic and cyclic loading tests on soil samples with varying initial fines contents, while maintaining a constant relative density of 79%, several important discoveries are made. The small strain stiffness of suffused soil fluctuates as erosion occurs, along with a decrease in shear strength and an increase in soil contraction under monotonic stress. Furthermore, the research simulates the train loading exerted on the base soil of bridge piers susceptible to suffusion by subjecting the soil samples to cyclic loading both before and after erosion, mirroring practical conditions. The key findings of this study reveal that the stiffness of soil drops during erosion with no significant deformation of the soil. This leads to a large strain accumulation in the soil specimens under subsequent cyclic traffic loading. These findings highlight that the delayed settlement or inclination of bridge piers under cyclic or train loading after major flood is possibly due to suffusion in the base soil of the piers. (c) 2024 Production and hosting by Elsevier B.V. on behalf of The Japanese Geotechnical Society. This is an open access article under the CC BY NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).

In recent years, local scour has occurred on the pier foundations of river bridges during heavy rain and river flooding, often resulting in bridge collapse and outflow. This study focused on the characteristic displacement, called delayed displacement, of the river bridge pier. The critical displacement of the piers was first observed several days after the flood when the train passed and not immediately after the flood. The authors hypothesized that one of the possible reasons for the delayed displacement is the suffusion of the supporting ground beneath the pier foundation during the flood, followed by a compressive behavior due to the collapse of the soil skeleton under repeated traffic loads. Accordingly, this study performed erosion tests simulating flood and cyclic loading tests simulating train passage using a triaxial test apparatus to check the validity of this hypothesis. In some test cases, suffusion without any deformation occurred in the erosion test but deformed in the cyclic loading test just after the erosion test. This behavior matches the behavior of delayed displacement. It was also suggested that the risk of the delayed displacement becomes high when the soil skeleton was assumed to primarily comprise fine particles, and the void ratio and hydraulic gradient were high. By contrast, when the soil skeleton was assumed to primarily comprise coarse particles, suffusion occurred in the erosion test, but did not deform in the subsequent cyclic loading test. Thus, the risk of delayed displacement is considered to be low when coarse particles are dominant. Furthermore, clear relationship between suffusion and the consequent reduction in soil stiffness cannot be observed. This result indicates that no significant change in the stiffness occur in the supporting ground of the pier foundation at the stage immediately before the delayed displacement. Thus, identifying the deterioration in the stability of the piers through impact loading test, which is based on the concept that local scour reduces the natural frequency of the bridge pier, is difficult.

This study investigates, for the first time ever, the suffusion on gap-graded granular soils under torsional shear conditions from a microscopic perspective. A numerical model of the hollow cylinder torsional shear test (HCTST) using the discrete element method (DEM) is first developed, where an algorithm for simulating the real inner and outer rubber membranes of the hollow cylinder apparatus (HCA) is introduced. After the validation, the computational fluid dynamics (CFD) approach is introduced for the coupling between the particle and fluid phases. Then, a series of the coupled CFD-DEM suffusion simulations considering the rotation of the major principal stress axis (alpha) and intermediate principal stress ratio (b) are conducted. It is found that more fine particles are eroded in cases having smaller alpha and b, and the clogging phenomenon in the middle zones becomes more significant as both alpha and b increase. From the microscopic perspective, the specimens whose contact anisotropy principal direction is close to the fluid direction will lose more fines, and the anisotropy magnitude also plays an important role. In addition, the differences in structure and vertical connectivity of the pores in HCTST samples under various complex loading conditions cause fine particles to have different migration paths, further resulting in different fines mass loss.

Seepage -induced suffusion involves the migration of fine particles within a soil matrix. Seepage flow is affected by the soil permeability anisotropy of anisotropic soil fabric; however, suffusion anisotropy is unclear because of the limited function of existing permeameters. In recent studies, the effect of seepage direction has been investigated under only low hydraulic gradients because the control of seepage direction relies merely on gravity. In this study, a new, large -sized permeameter is developed with which suffusion tests can be conducted along horizontal or vertical seepage directions under high hydraulic gradients. Correspondingly, the permeameter can accommodate a specimen of 540 x 500 x 470 or 540 x 540 x 440 mm3 (length x width x height). The seepage direction is switched by changing the boundary conditions of the specimen with detachable perforated plates that allow pressurized water originating from different inlets to flow along horizontal or vertical directions. Two repeated pairs of tests were performed on a gap -graded clayey gravel to investigate the suffusion anisotropy of saturated clayey gravel. The results show that the maximum relative deviations of measurements for initial hydraulic conductivity, initiation, and failure hydraulic gradients are less than 3.5 %, demonstrating satisfactory reliability. The ratio of the initial horizontal hydraulic conductivity to vertical hydraulic conductivity for the test soil is 13.87, indicating a significantly anisotropic fabric induced by compaction. The ratios of horizontal initiation and failure hydraulic gradients to vertical initiation and failure hydraulic gradients are 0.52 and 0.59, respectively. This implies that suffusion anisotropy should not be neglected for evaluating the internal instability of anisotropic soils.

Bei dem Beitrag handelt es sich um die erweiterte Fassung der gleichnamigen Keynote Lecture auf der 4. Bodenmechanik Tagung im Rahmen der Fachsektionstage Geotechnik, die auf Anregung der Fachsektionsleitung auch in der Zeitschrift geotechnik veroffentlicht werden soll. Bei den Phanomenen der inneren Erosion in durchstromtem Boden und in Erdbauwerken geht es um das Losen, den Transport und die Ablagerung bevorzugter Fraktionen mit der Folge einer anderung der Bodeneigenschaften. Die Phanomene der inneren Erosion werden als Kontakterosion, Suffosion, Kolmation und ruckschreitende Erosion charakterisiert. Die Kinematik dieser physikalischen Prozesse ergibt sich mit der Energie einer Sickerstromung aus der Bewegung des Einzelkorns im Porenraum, den moglichen Freiheitsgraden beim Transport. Der Artikel gibt einen uberblick uber die Art und Bedingungen der verschiedenen Phanomene sowie uber deren spezifische Kinematik innerhalb der Bodenstruktur. Die relevanten international verwendeten Nachweismethoden und Kriterien werden aufgefuhrt und in ihrer Aussagekraft bewertet. Die kennzeichnenden Einflussparameter werden aufgezeigt. Fur die einzelnen Phanomene der inneren Erosion werden Strategien zur Bewertung und Beherrschung des Erosionsrisikos diskutiert. Phenomena, kinematics and risk assessment strategies of internal erosionInternal soil erosion due to seeping water in natural sediments as well as in earthworks can lead to a significant change in soil properties and could even destroy the structural integrity. The physical process of erosion always is induced by loosening, migration, and deposition of predominant fine particles within the soil structure. Depending on the kinematics, the phenomena are divided into contact erosion, suffusion, colmation and backward erosion piping. The kinematics is controlled by the energy of a seepage on the one hand, on the other by the degrees of freedom and the boundary conditions of an individual grain movement within the pore space. This article provides an overview of the characteristic, specifics, and conditions of the different phenomena considering their kinematics within the soil structure. Internationally used approaches and methods of assessment are listed, their significance and their limitations will be evaluated. The impact parameters that control the different processes are shown. Strategies for assessing and controlling the risk of structural damage are discussed for the different phenomena of internal soil erosion.

Water level rise and fall can lead to changes in seepage flow in the bank slopes of reservoirs. Cyclic seepage flow can result in the redistribution and even loss of fine particles in the slope, which eventually influences the slope stability. In this study, numerical simulations were conducted based on a multispecies transport finite element method (FEM) to investigate the influence of periodic water level fluctuations on the reservoir bank slope. The unsaturated soil was treated as a five-constituent mixture, by which erosion and deposition behaviours were captured by phase transition between deposited and fluidized fines. A constitutive model for suffusion considering the hydraulic fluctuation effect was employed to quantify the erosion process. The change in the seepage field and the corresponding evolution of the erosion rate and fines eroded ratio in the slope during the fluctuation period were investigated. The influence of the water level fluctuation frequency on the fines redistribution in the slope was analysed. Furthermore, the evolution patterns of the slope stability under different fluctuation frequencies were compared. The numerical results indicated that a higher fluctuation frequency could promote fine particle loss in the zones near the slope surface, consequently inducing greater deterioration in the slope stability.

Internal erosion can induce significant changes in the mechanical properties of soils, posing various hazards to dam and dike structures. Despite its importance, our current understanding of this phenomenon remains incomplete. The influence of pre-shearing stress conditions on the mechanical behaviours of soils during the internal erosion process is particularly challenging, as existing experiments have not been able to maintain the constant pre-shearing stress ratios. To bridge this gap in knowledge, this paper presents a series of discrete element method (DEM) simulations focused on gap-graded cohesionless soil. The primary objective of these simulations is to investigate two specific cases of internal erosion: suffusion and suffosion processes. Soil specimens are subjected to different pre-shearing stress ratios in the standard triaxial tests before being submitted to different levels of erosion to study their constitutive responses. The results show that erosion-induced deformation (i.e. suffosion) only starts after a specific amount of mass loss. This mass loss and the pre-shearing stress ratio form a well-defined criterion for triggering suffosion, which is named suffosion surface. The volumetric strain is shown to be a better indicator to describe the suffosion process than the commonly used void ratio. The pre-shearing stress ratio significantly influences the suffosion response of the soil sample, with a higher preshearing stress ratio facilitating soil failure. Furthermore, soil specimens undergo both deviatoric and volumetric responses during the suffosion process. To this end, new DEM-based statistical equations were proposed to describe the observed mechanisms, which are helpful for the future development of constitutive models to describe internal soil erosion.