This study developed a novel, 3D coupled computational fluid dynamics (CFD)-discrete element method (DEM) model by coupling two software programs, OpenFOAM and PFC^3D, to solve problems related to fluid–particle interaction systems. The complete governing equations and the flow chart of the coupling calculations are clearly presented herein. The coupled CFD-DEM model was first benchmarked using two classic geo-mechanics problems, for which the analytical solutions are available. Then, the CFD-DEM model was employed to investigate the flow characteristics of a particle heap subjected to the effects of water inside a pipe under different conditions. The results showed that particle size and pipe inclination angle can significantly affect the particle flow morphology, total kinetic energy and erosion rate for mono-sized particles, whereas polydisperse particles had a slight effect. This model can accurately describe the flow characteristics of particles inside a pipe, and the results of this study were consistent with those of previous studies. The reliability of this model was further demonstrated, which showed that this model can provide valuable references for solving similar problems such as soil erosion and bridge scour problems.

The cement manufacturing process causes environmental problems such as carbon dioxide emissions. Bentonite-cement-water glass grouting (BCG) is a grout with low environmental impact, designed to reduce the use of cement. The significance of this study is to explore the fundamental properties of BCGs with the bentonite contents of up to 50%, and the antiwash-out and permeability characteristics of BCG in comparison with cement-water glass grout (CG). The effects of water-to-solid ratios (W/S), bentonite content, water glass content and the volume ratio of bentonite-cement slurry to water glass slurry (V_a: V_b) on the fundamental properties of BCGs were tested in the laboratory. The results show that the BCGs have excellent characteristics of high viscosity, low bleeding, short and controllable setting time, and high early strength. The BCGs have been applied to the pre-grouting of the water-rich sandy pebble stratum in the Wangfujing North Station section of Beijing Metro Line 8 and achieved better water sealing and reinforcement effects than CG.

Due to rainfall infiltration, slope instability becomes frequent, which is the main reason for landslide disasters. In this study, the stability of slope affected by rainfall was analyzed using an indoor model test and geo-studio simulation method, and the variation law of phreatic line, seepage field, the most dangerous sliding surface, and safety factor with time were studied under rainfall infiltration. Research results showed that under the effect of rainfall, the slope failure presented a typical traction development mode. With the increase of time, the phreatic line of the slope kept rising, the water head keeps increasing, the seepage depth in the slope became deeper, and the slope stability worsened until the slope was damaged. The water head height decreased gradually from the slope left boundary to the right, and the water head width decreased gradually. The soil at the slope back edge was damaged, and the sliding soil accumulated at the slope foot, forming a gentle slope, which increased the shear strength of the slope, making the slope finally reach a stable state. In this process, the overlying soil changed from an unsaturated state to a saturated state, the pore water pressure and soil pressure increased, and then the slope was damaged, both of which decreased. Under high rainfall intensity, the slope was damaged, the soil in the slope was rapidly saturated, and the time required to produce the sliding area was short. When the rainfall intensity was the same, the smaller the slope angle was, the smaller the safety factor was. When the slope angle was the same, the greater the rainfall intensity was, the smaller the safety factor was.

River floods, dammed lake flood discharge, reservoir discharge, seawater recession, etc. all cause the water level in front of a slope to drop, which changes the original steady-state seepage field in the soil, leading to harmful slope instability. To study this phenomenon, a numerical model was established through theoretical analysis combined with the coupling of the Seep/W and Slope/W modules of the GeoStudio finite element software, and the numerical model was verified by the model test results of indoor medium sand and silt. This paper focuses on the effects when the water level in front of a slope drops at different speeds, different drop ratios, different initial water levels, different filling materials, and matrix suction on the seepage field and slope stability. The conclusions are as follows: (1) the greater the speed at which the water level in front of a slope falls, the greater the downward seepage force formed by the seepage field of the slope to the slope body; (2) the change curve of the safety factor at a higher speed is steeper when the water level falls at different speeds, and the safety factor value when the water level in front of the slope is constant is smaller; (3) the safety factor of the slope decreases with an increase in the drop ratio; when the drop ratio is the same, the loss of stability is worse if the initial water level is lower; (4) when there is a drawdown of water levels in front of the slope, the non-cohesive medium sand slope is more prone to instability failure than the cohesive silt slope; and (5) when this modeling method is applied to matrix suction, the effect of matrix suction increases the safety factor of the slope.

In practical engineering applications, silt is prone to liquefaction and quicksand. This paper mainly studies the improvement effects of urease, lignin and their mixture on the strength and liquefaction resistance of silt. Based on the results and phenomena of an unconfined compressive strength and dynamic triaxial test, the improvement effects of the compressive strength, deformation resistance and liquefaction resistance of silt under different improvement schemes are analyzed, and the optimal values of the cement or lignin when enzyme-induced calcium carbonate precipitation (EICP) technology, lignin alone or EICP and lignin are obtained. The results show that the optimum concentration of the constant temperature and humidity sample (referred to as the constant humidity sample) and the constant temperature immersion sample (referred to as the soaking sample) of urease in the unconfined compressive strength test is 1.0 mol/L, and the compressive strength of the soaking sample is 4.9 MPa, which is 1.56 times that of the constant humidity sample; the optimum addition ratio of the lignin-improved constant humidity sample is 3%, and its compressive strength is 2.07 Mpa; the optimum addition ratio of the samples immersed at constant temperature is 4%, and the compressive strength is 3.05 MPa; when urease combines with lignin to improve silt, 4% is the best lignin addition ratio, the compressive strength of the constant humidity sample reaches 1.57 Mpa and the compressive strength of the soaking sample reaches 3.75 MPa; in the dynamic triaxial multi-stage cyclic load test, all samples were cured at constant humidity sample, and in the urease modified silt scheme, 1.0 mol/L was the optimal cement concentration; in the scheme of improving silt with lignin, 3% is the optimal addition ratio; when 1.25 mol/L cementation solution plus urease crude extract is combined with different ratios of lignin in the experimental scheme, 3% is the best lignin addition ratio.