Advances in Coupling Computational Fluid Dynamics and Discrete Element Method in Geotechnical Problems

In some cases, the water content in granular soil increases to the extent that it becomes saturated, which noticeably alters its responses. For example, the pore water pressure within saturated granular soil would increase rapidly under sudden external loading, which is equivalent to undrained or constant volume conditions. This reduces the effective stress in soil dramatically and may result in catastrophic failure. There have been different numerical approaches to analyse such a failure mechanism of soil to provide a deeper understanding of soil behaviour at the microscopic level. One of the most common numerical tools for such analysis is the discrete element method (DEM) due to its advantage in obtaining microscopic properties (e.g., statistics on particle contacts and fabric), reproducibility and simple feedback control. However, most DEM studies ignore the fluid phase and merely consider the solid particles while the fluid pressure is indirectly calculated by mimicking undrained condition to a constant volume condition. Note that fluid’s influence does not limit to the change of pore water pressure. For example, the external loading would induce the movement of fluid, and the fluid-solid interaction could subsequently drag the solid particles to shift within the system. In addition, the state of soil could change from solid to suspension under an excess hydraulic gradient. Therefore, the study of the fluid-solid mixture is essential as it is a typical scenario in geotechnical practice, and the simulations of saturated sand should be conducted in numerical forms in which both the solid and fluid phases can be modelled.