Extension of discontinuous deformation analysis method for coupled reactive transport problems

Fluid flow and the accompanying reactant transport in fractured rock masses are key concerns for understanding the performance and safety assessments of subsurface engineering problems, particularly in waste disposal and tunneling project. Nevertheless, predicting solute transport behavior remains a great challenge due to computational sophistication and numerical instability if numerous discontinuities and chemical reaction is taken into consideration. In this study, the partial differential equation describing the reactive transport of solute is solved and coupled with the Discontinuous Deformation Analysis (DDA) method. First, the extended DDA method is validated by simulating an indoor dissolution test on the single rock fracture. Then, it is applied to establish a discrete fracture network (DFN), and be further employed to investigate the reactive transport behavior and geochemical evolution in a case of acidic mine drainage in a fractured stratum. The study analyzes the effects of stress field, reaction speed, and diffusion coefficient on the reactive transport. Numerical results indicate that the conclusions reached are consistent with current findings, combined with its prominent features in simulating the mechanical behavior of discontinuities, demonstrating that the extended method is a promising tool for the coupled hydro-mechanical-chemo analysis.