Seepage failure prediction of breakwater using an unresolved ISPH-DEM coupling method enriched with Terzaghi’s critical hydraulic gradient

Advanced Modeling and Simulation in Engineering Sciences

Research Article

Mitsuteru Asai¹  Kiyonobu Kasama¹  Kumpei Tsuji¹ 

[1] Department of Civil Engineering, Kyushu University, 744, Motooka, Nishiku, 819-0395, Fukuoka, Japan;

关键词: ISPH;    DEM;    Unresolved coupling;    Seepage failure;    Terzaghi’s critical hydraulic gradient;   

DOI  :  10.1186/s40323-022-00239-3

 received in 2022-08-10, accepted in 2022-12-25,  发布年份 2022

来源: Springer

【 摘 要 】

This study develops a new numerical simulation model for rubble mound failure prediction caused by piping destruction under seepage flows. The piping has been pointed out as a significant cause of breakwater failure during tsunamis. Once boiling and heaving occur on the mound surface, the piping suddenly propagates in the opposite direction of seepage flow. For the seepage failure prediction, a coupled fluid-soil-structure simulator is developed by combining the ISPH for fluid and the DEM for rubble mounds and caisson blocks. The ISPH, a Lagrangian particle method for incompressible fluids, can simulate seepage and violent flows such as tsunamis. The DEM has been applied for discrete particle and rigid body simulations that include discontinuous deformation, as in the rubble mounds failure and large displacement of the caisson block. ISPH-DEM coupling simulations have already been proposed as a technique for multi-phase flows. Still, the technique cannot reproduce the sudden onset of piping from a stable mound. Two simple assumptions are applied to reduce the numerical cost for the fluid-soil-structure simulators of a breakwater structure composed of a rubble mound and the caisson block. Firstly, each rubble is modeled as an idealized spherical DEM particle with the mean diameter of the rubble. The ISPH particle size is assumed to be the same size as the DEM particle. Under these assumptions, the unresolved coupling model between rubble mound particles and fluid, which obtains the interaction through empirical drag force, should be applied. At the same time, the interaction between the fluid and the caisson block is fully resolved with the spatial resolution with the ISPH and DEM particle size. Our new contribution in this paper is how to model the interaction as an unresolved coupling between seepage flow simulated by ISPH and rubble mound particle modeled with DEM. Our original seepage failure experiment is simulated using the proposed ISPH-DEM coupling simulator. We identified the conventional drag force models as the unresolved coupling model are insufficient to initiate the boiling and piping observed in the experiment. It may be due in one part to excessive averaging of flow velocities caused by unresolved coupling. Therefore, Terzaghi’s critical hydraulic gradient is introduced to initiate the boiling and heaving. Unstable DEM particles, judged by Terzaghi’s critical hydraulic gradient, gradually lose their mass to represent unresolved suspended fine rubble mound particles. Our models qualitatively reproduce the sand boiling and backward erosion in the opposite direction of the seepage flow, as shown in the experiment.

【 授权许可】

CC BY   
© The Author(s) 2023

【 图 表 】

Fig. 1

Fig. 3

Fig. 3

42004_2022_800_Article_IEq108.gif

41116_2022_35_Article_IEq1.gif

Fig. 2

Fig. 6

Fig. 1

42004_2022_800_Article_IEq97.gif

Fig. 2

Fig. 3

Fig. 1

Scheme 1

Fig. 8

40249_2022_1049_Article_IEq7.gif

42004_2022_800_Article_IEq92.gif

Fig. 1

42004_2022_803_Article_IEq36.gif

42004_2022_800_Article_IEq91.gif

Fig. 6

Fig. 6

Fig. 1

42004_2022_800_Article_IEq86.gif

Fig. 4

42004_2022_800_Article_IEq84.gif

42004_2022_800_Article_IEq83.gif

42004_2022_800_Article_IEq82.gif

Fig. 5

Fig. 65

Fig. 2

Fig. 4

Fig. 64

42004_2022_800_Article_IEq75.gif

Fig. 2

42004_2022_800_Article_IEq73.gif

Fig. 3

Fig. 1

Fig. 63

42004_2022_800_Article_IEq69.gif

Fig. 2

Fig. 2

Fig. 1

Fig. 3

Fig. 3

Fig. 5

Fig. 62

40798_2022_490_Article_IEq91.gif

Fig. 1

40798_2022_490_Article_IEq87.gif

13690_2022_1010_Article_IEq14.gif

13690_2022_1010_Article_IEq12.gif

13690_2022_1010_Article_IEq10.gif

Fig. 2

13690_2022_1010_Article_IEq8.gif

13690_2022_1010_Article_IEq7.gif

13690_2022_1010_Article_IEq5.gif

13690_2022_1010_Article_IEq4.gif

Fig. 1

13690_2022_1010_Article_IEq1.gif

40798_2022_490_Article_IEq73.gif

Fig. 3

40798_2022_490_Article_IEq71.gif

40798_2022_490_Article_IEq70.gif

Fig. 2

Fig. 2

Fig. 1

Fig. 60

Fig. 1

Fig. 1

40798_2022_490_Article_IEq64.gif

40798_2022_490_Article_IEq63.gif

Fig. 3

Fig. 4

40798_2022_490_Article_IEq60.gif

Fig. 4

Fig. 59

Fig. 3

42004_2022_800_Article_IEq39.gif

Fig. 3

Fig. 1

Fig. 5

Fig. 1

Fig. 1

Fig. 5

42004_2022_800_Article_IEq32.gif

Fig. 4

Fig. 4

Fig. 2

40798_2022_490_Article_IEq58.gif

40798_2022_490_Article_IEq57.gif

40798_2022_490_Article_IEq56.gif

40798_2022_490_Article_IEq55.gif

40798_2022_490_Article_IEq54.gif

40798_2022_490_Article_IEq53.gif

40798_2022_490_Article_IEq52.gif

40798_2022_490_Article_IEq510.gif

Fig. 1

Fig. 2

40798_2022_490_Article_IEq48.gif

40798_2022_490_Article_IEq47.gif

40798_2022_490_Article_IEq46.gif

40798_2022_490_Article_IEq45.gif

40798_2022_490_Article_IEq44.gif

Fig. 4

40798_2022_490_Article_IEq42.gif

Fig. 1

40798_2022_490_Article_IEq41.gif

40798_2022_490_Article_IEq12.gif

40798_2022_490_Article_IEq11.gif

Fig. 58

40798_2022_490_Article_IEq9.gif

Fig. 5

Fig. 3

Fig. 5

40798_2022_490_Article_IEq5.gif

40798_2022_490_Article_IEq3.gif

40798_2022_490_Article_IEq2.gif

Fig. 1

40854_2022_439_Article_IEq93.gif

Fig. 3

Fig. 57

Fig. 4

Fig. 55

Fig. 2

Fig. 6

Fig. 5

Fig. 54

Fig. 2

Fig. 1

12864_2023_9114_Article_IEq10.gif

12888_2022_4500_Article_IEq3.gif

12888_2022_4500_Article_IEq2.gif

12888_2022_4500_Article_IEq1.gif

Fig. 2

Fig. 1

Fig. 7

Fig. 4

Fig. 11

Fig. 3

Fig. 10

Fig. 51

Fig. 9

Fig. 6

Fig. 1

Fig. 8

Fig. 3

Fig. 49

Fig. 4

Fig. 7

Fig. 4

Fig. 6

Fig. 2

Fig. 5

Fig. 48

Fig. 1

Fig. 4

Fig. 2

Fig. 5

Fig. 3

Fig. 5

41116_2022_35_Article_IEq659.gif

41116_2022_35_Article_IEq655.gif

Fig. 47

Fig. 2

41116_2022_35_Article_IEq649.gif

Fig. 1

41116_2022_35_Article_IEq645.gif

41116_2022_35_Article_IEq643.gif

41116_2022_35_Article_IEq641.gif

41116_2022_35_Article_IEq639.gif

41116_2022_35_Article_IEq635.gif

Fig. 46

41116_2022_35_Article_IEq632.gif

41116_2022_35_Article_IEq627.gif

41116_2022_35_Article_IEq626.gif

41116_2022_35_Article_IEq624.gif

41116_2022_35_Article_IEq569.gif

41116_2022_35_Article_IEq538.gif

Fig. 2

41116_2022_35_Article_IEq515.gif

41116_2022_35_Article_IEq504.gif

41116_2022_35_Article_IEq494.gif

Fig. 5

41116_2022_35_Article_IEq472.gif

41116_2022_35_Article_IEq459.gif

41116_2022_35_Article_IEq458.gif

Fig. 17

41116_2022_35_Article_IEq450.gif

41116_2022_35_Article_IEq446.gif