Gases such as hydrogen can generate from the disposal canister in high-level radioactive waste disposal systems owing to the corrosion of cooper container in anoxic conditions. The gas can be accumulated in the voids of bentonite buffer around the disposal canister if gas generation rates become larger than the gas diffusion rate of bentonite buffer with the low-permeability. Continuous gas accumulations result in the increase in gas pressure, causing sudden dilation flow of gases with the gas pressure exceeding the gas breakthrough pressure. Given that the gas dilation flow can cause radionuclide leakage out of the engineered barrier system, it is necessary to consider possible damages affected by the radionuclide leakage and to properly understand the complicated behaviors of gas flow in the bentonite buffer with low permeability. In this study, the coupled hydro-mechanical model combined with the damage model that considers two-phase fluid flow and changes in hydraulic properties affected by mechanical deformations is applied to numerical simulations of 1-D gas injection test on saturated bentonite samples (refer to DECOVALEX-2019 Task A Stage 1A). To simulate the mechanical behavior of microcracks which occur due to the dilation flow caused by increase in gas pressure, a concept of elastic damage constitutive law is considered in the coupled hydro-mechanical model. When the TOUGH-FLAC coupling-based model proposed in this study is applied, changes in hydraulic properties affected by mechanical deformations combined with the mechanical damage are appropriately considered, and changes in gas injection pressure, pore pressures at radial filters and outlet, and stress recorded during the gas injection test are accurately simulated.

Discontinuum-based numerical methods can contain the multiple discontinuities in a model and reflect the thermal, hydraulic and mechanical characteristics of discontinuities. Therefore, discontinuum methods can be appropriate to simulate the model which require the detailed analysis of the coupled thermo-hydro-mechanical processes in fractured rock such as geothermal energy, CO2 geo-sequestration, and geological repository of the high-level radioactive waste. TOUGH-3DEC, the three-dimensional discontinuum simulators for the coupled thermo-hydro-mechanical analysis, was developed by linking the integral finite difference method TOUGH2 and the explicit distinct element method 3DEC to describe the coupled thermo-hydro-mechanical processes in both porous media and discontinuity. TOUGH2 handles thermo-hydraulic analysis by the internal simulation module, and 3DEC performs mechanical study based on the constitutive models of porous media and discontinuity with coupling the thermal and hydraulic response from TOUGH2. The thermal and hydraulic couplings are the key processes and should be carefully verified by sufficient cases, so this study performed the thermomechanical and hydro-mechanical simulations which are modelling the analytic solutions including the uniaxial consolidation, fracture static opening, and the heating of a hollow cylinder problems. Each thermo-mechanical and hydro-mechanical verification case is also validated by comparing with the results of the other continuum and discontinuum-based numerical methods. TOUGH-3DEC results follow the analytic solutions and show better accuracy than the continuum-based numerical methods in the static fracture opening problem. The developed TOUGH-3DEC simulator can be expanded to coupled thermo-hydro-mechanical-chemical analysis in fractured rock mass, and the simulator needs to be verified by more complicated coupled processes problems which require in the chemical coupling.