CFD–DEM Modeling of Stress–Damage–Seepage Coupling Mechanisms and Support Strategies in Subsea Tunnel Excavation

The stability of subsea tunnels is governed by the strong coupling among stress redistribution, damage evolution, and seepage flow (Stress–Damage–Seepage, SDS). The dynamic interplay, especially under high water pressure, often leads to catastrophic failures, yet its mechanisms, particularly the role of support timing, remain insufficiently understood due to limitations in conventional numerical methods. This study aims to unravel the SDS coupling mechanisms during tunnel excavation under high hydraulic head, and to quantitatively investigate how support timing influences the stability of the surrounding rock within this coupled system. A coupled Computational Fluid Dynamics and Discrete Element Method (CFD-DEM) framework was employed. In this approach, excavation-induced damage, crack propagation, and fluid–particle interactions are explicitly resolved at the particle scale, whereas the macroscopic permeability evolution is captured through an imposed empirical exponential relationship. Simulations were conducted under both steady-state and transient seepage conditions with varying stress ratios and water heads. High-head transient seepage intensifies SDS coupling, dynamically redistributing seepage forces to damage zone edges and amplifying damage. Support timing critically mediates this interaction: premature support risks tensile failure at the tunnel periphery, while delayed support allows a vicious cycle of shear failure and increased inflow. Optimal “timely” support, applied after initial deformation, diverts high seepage forces inward, minimizing final damage. The spatiotemporal synchronization of transient seepage forces with damage evolution is pivotal for stability. Support timing acts as a key control variable. The CFD-DEM framework effectively elucidates these micro-mechanisms, providing a scientific basis for the dynamic design of support in high-pressure subsea tunnels.