• Tunnel and Underground Space
• /
• v.30 no.4
• /
• pp.271-305
• /
• 2020

An effect of coupled thermo-hydro-mechanical and chemical (THMC) behavior is an essential part of the performance and safety assessment of geological disposal systems for high-level radioactive waste and spent nuclear fuel. Furthermore, numerical models and modeling techniques are necessary to analyze and predict the coupled THMC behavior in the disposal systems. However, phenomena associated with the coupled THMC behavior are nonlinear, and the constitutive relationships between them are not well known. Therefore, it is challenging to develop numerical models and modeling techniques to analyze and predict the coupled THMC behavior in the geological disposal systems. It is also difficult to verify and validate the development of the models and techniques because it requires expensive laboratory tests and in-situ experiments that need to be performed for a long time. DECOVALEX was initiated in 1992 to efficiently develop numerical models and modeling techniques and validate the developed models and techniques against the lab and in-situ experiments. In Korea, Korea Atomic Energy Research Institute has participated in DECOVALEX-2011, DECOVALEX-2015, and DECOVALEX-2019 since 2008. In this study, all tasks in the three DECOVALEX projects were introduced to the researcher in the field of rock mechanics and geotechnical engineering in Korea.

• Tunnel and Underground Space
• /
• v.31 no.4
• /
• pp.270-288
• /
• 2021

We proposed a numerical method to simulate the hydro-mechanical behavior of rock fracture using a grain-based distinct element model (GBDEM) in the paper. As a part of DECOVALEX-2023 Task G, we verified the method via benchmarks with analytical solutions. DECOVALEX-2023 Task G aims to develop a numerical method to estimate the coupled thermo-hydro-mechanical processes within the crystalline rock fracture network. We represented the rock sample as a group of tetrahedral grains and calculated the interaction of the grains and their interfaces using 3DEC. The micro-parameters of the grains and interfaces were determined by a new methodology based on an equivalent continuum approach. In benchmark modeling, a single fracture embedded in the rock was examined for the effects of fracture inclination and roughness, the boundary stress condition and the applied pressure. The simulation results showed that the developed numerical model reasonably reproduced the fracture slip induced by boundary stress condition, the fracture opening induced by fluid injection, the stress distribution variation with fracture inclination, and the fracture roughness effect. In addition, the fracture displacements associated with the opening and slip showed good agreement with the analytical solutions. We expect the numerical model to be enhanced by continuing collaboration and interaction with other research teams of DECOVALEX-2023 Task G and validated in further study experiments.

• Tunnel and Underground Space
• /
• v.31 no.6
• /
• pp.610-622
• /
• 2021

The DECOVALEX project is one of the representative international cooperative projects to enhance the understanding of the complex Thermo-Hydro-Mechanical-Chemical(THMC) coupled behavior in the high-level radioactive waste disposal system based on the numerical simulation. DECOVALEX-2023 is the current phase consisting of 7 tasks, and Task C aims to model the THM coupled behavior in the disposal system based on the Full-scale Emplacement (FE) experiment at the Mont-Terri underground rock laboratory. This study performs the numerical simulation based on the OGS-FLAC developed for the current study. In the numerical model, we emplaced the heater with constant power horizontally based on the FE experiment and monitored the pressure development, temperature increase, and mechanical deformation at the specific monitoring points. We monitored the capillary pressure as the primary effect inducing the flow in the buffer system, and thermal stress and pressurization were dominant in the surrounding rocks' area. The results will also be compared and validated with the other participating groups and the experimental data further.

• Tunnel and Underground Space
• /
• v.17 no.1 s.66
• /
• pp.32-42
• /
• 2007

Severe damage can occur around deposition holes due to complex interaction of thermo-hydro-mechanical (THM) loading during the long term operation of high level radioactive waste repository. Many candidate sites for repository are located in crystalline rock mass, therefore mechanism of damage follows the form of brittle fracture and failure. This paper briefly introduces major outcomes from 15 years international collaborative project, DECOVALEX, and presents major study results for current ongoing benchmark test study from DECOVALEX-THMC, to evaluate the effect of THM loading to rock mass in excavation damaged zone (EDZ) near deposition holes. Through benchmark test model by simplifying THM loading to boundary loading obtained numerical results are compared, and discrete fracture interaction after up to 1 million years operation is discussed.

• Tunnel and Underground Space
• /
• v.30 no.6
• /
• pp.573-590
• /
• 2020

This study presents the current status of DECOVALEX-2023 project Task G and our research results so far. Task G, named 'Safety ImplicAtions of Fluid Flow, Shear, Thermal and Reaction Processes within Crystalline Rock Fracture NETworks (SAFENET)' aims at developing a numerical method to simulate the fracture creation and propagation, and the coupled thermohydro-mechanical processes in fracture in crystalline rocks. The first research step of Task G is a benchmark simulation, which is designed for research teams to make their modelling codes more robust and verify whether the models can represent an analytical solution for displacements of a single rock fracture. We reproduced the mechanical behavior of rock and embedded single fracture using a three-dimensional grain-based distinct element model for the simulations. In this method, the structure of the rock was represented by an assembly of rigid tetrahedral grains moving independently of each other, and the mechanical interactions at the grains and their contacts were calculated using 3DEC. The simulation results revealed that the stresses induced along the embedded fracture in the model were relatively low compared to those calculated by stress analysis due to stress redistribution and constrained fracture displacements. The fracture normal and shear displacements of the numerical model showed good agreement with the analytical solutions. The numerical model will be enhanced by continuing collaboration and interaction with other research teams of DECOVALEX-2023 Task G and validated using various experiments in a further study.

• Tunnel and Underground Space
• /
• v.31 no.6
• /
• pp.593-609
• /
• 2021

We proposed a numerical method for the thermo-mechanical behavior of rock fracture using a grain-based distinct element model (GBDEM) and simulated thermally induced fracture slip. The present study is the benchmark simulation performed as part of DECOVALEX-2023 Task G, which aims to develop a numerical method to estimate the coupled thermo-hydro-mechanical processes within the crystalline rock fracture network. We represented the rock sample as an assembly of Voronoi grains and calculated the interaction of the grains (blocks) and their interfaces (contacts) using a distinct element code, 3DEC. Based on an equivalent continuum approach, the micro-parameters of grains and contacts were determined to reproduce rock as an elastic material. Then, the behavior of the fracture embedded in the rock was characterized by the contacts with Coulomb shear strength and tensile strength. In the benchmark simulation, we quantitatively examined the effects of the boundary stress and thermal stress due to heat conduction on fracture behavior, focusing on the mechanism of thermally induced fracture slip. The simulation results showed that the developed numerical model reasonably reproduced the thermal expansion and thermal stress increment, the fracture stress and displacement and the effect of boundary condition. We expect the numerical model to be enhanced by continuing collaboration and interaction with other research teams of DECOVALEX-2023 Task G and validated in further study experiments.

• Tunnel and Underground Space
• /
• v.32 no.6
• /
• pp.568-585
• /
• 2022

In the present study, we proposed a numerical method for simulating thermally induced fracture slip using a grain-based distinct element model (GBDEM). As a part of DECOVALEX-2023, the thermo-mechanical loading test on a saw-cut rock fracture conducted at the Korea Institute of Civil Engineering and Building Technology was simulated. In the numerical model, the rock sample including a saw-cut fracture was represented as a group of random Voronoi polyhedra. Then, the coupled thermo-mechanical behavior of grains and their interfaces was calculated using 3DEC. The key concerns focused on the temperature evolution, thermally induced principal stress increment, and fracture normal and shear displacements under thermo-mechanical loading. The comparisons between laboratory experimental results and the numerical results revealed that the numerical model reasonably captured the heat transfer and heat loss characteristics of the rock specimen, the horizontal stress increment due to constrained displacement, and the progressive shear failure of the fracture. However, the onset of the fracture slip and the magnitudes of stress increment and fracture displacement showed discrepancies between the numerical and experimental results. We expect the numerical model to be enhanced by continuing collaboration and interaction with other research teams of DECOVALEX-2023 Task G and validated in further study.

• Tunnel and Underground Space
• /
• v.30 no.4
• /
• pp.306-319
• /
• 2020

Characterizations of Excavation Damage Zone (EDZ), which is hydro-mechanical degrading the host rock, are the important issues on the geological repository for the spent nuclear fuel. In the DECOVALEX 2019 project, Task G aimed to model the fractured rock numerically, describe the hydro-mechanical behavior of EDZ, and predict the change of the hydraulic factor during the lifetime of the geological repository. Task G prepared two-dimensional fractured rock model to compare the characteristics of each simulation tools in Work Package 1, validated the extended three-dimensional model using the TAS04 in-situ interference tests from Äspö Hard Rock Laboratory in Work Package 2, and applied the thermal and glacial loads to monitor the long-term hydro-mechanical response on the fractured rock in Work Package 3. Each modelling team adopted both Finite Element Method (FEM) and Discrete Element Method (DEM) to simulate the hydro-mechanical behavior of the fracture rock, and added the various approaches to describe the EDZ and fracture geometry which are appropriate to each simulation method. Therefore, this research can introduce a variety of numerical approaches and considerations to model the geological repository for the spent nuclear fuel in the crystalline fractured rock.

• Tunnel and Underground Space
• /
• v.28 no.6
• /
• pp.670-691
• /
• 2018

This study presents the research results of the BMT(Benchmark Model Test) simulations of the DECOVALEX-2019 project Task B. Task B named 'Fault slip modelling' is aiming at developing a numerical method to predict fault reactivation and the coupled hydro-mechanical behavior of fault. BMT scenario simulations of Task B were conducted to improve each numerical model of participating group by demonstrating the feasibility of reproducing the fault behavior induced by water injection. The BMT simulations consist of seven different conditions depending on injection pressure, fault properties and the hydro-mechanical coupling relations. TOUGH-FLAC simulator was used to reproduce the coupled hydro-mechanical process of fault slip. A coupling module to update the changes in hydrological properties and geometric features of the numerical mesh in the present study. We made modifications to the numerical model developed in Task B Step 1 to consider the changes in compressibility, Permeability and geometric features with hydraulic aperture of fault due to mechanical deformation. The effects of the storativity and transmissivity of the fault on the hydro-mechanical behavior such as the pressure distribution, injection rate, displacement and stress of the fault were examined, and the results of the previous step 1 simulation were updated using the modified numerical model. The simulation results indicate that the developed model can provide a reasonable prediction of the hydro-mechanical behavior related to fault reactivation. The numerical model will be enhanced by continuing interaction and collaboration with other research teams of DECOVALEX-2019 Task B and validated using the field experiment data in a further study.

• Tunnel and Underground Space
• /
• v.28 no.5
• /
• pp.400-425
• /
• 2018

This study presents the research results and current status of the DECOVALEX-2019 project Task B. Task B named 'Fault slip modelling' is aiming at developing a numerical method to simulate the coupled hydro-mechanical behavior of fault, including slip or reactivation, induced by water injection. The first research step of Task B is a benchmark simulation which is designed for the modelling teams to familiarize themselves with the problem and to set up their own codes to reproduce the hydro-mechanical coupling between the fault hydraulic transmissivity and the mechanically-induced displacement. We reproduced the coupled hydro-mechanical process of fault slip using TOUGH-FLAC simulator. The fluid flow along a fault was modelled with solid elements and governed by Darcy's law with the cubic law in TOUGH2, whereas the mechanical behavior of a single fault was represented by creating interface elements between two separating rock blocks in FLAC3D. A methodology to formulate the hydro-mechanical coupling relations of two different hydraulic aperture models and link the solid element of TOUGH2 and the interface element of FLAC3D was suggested. In addition, we developed a coupling module to update the changes in geometric features (mesh) and hydrological properties of fault caused by water injection at every calculation step for TOUGH-FLAC simulator. Then, the transient responses of the fault, including elastic deformation, reactivation, progressive evolutions of pathway, pressure distribution and water injection rate, to stepwise pressurization were examined during the simulations. The results of the simulations suggest that the developed model can provide a reasonable prediction of the hydro-mechanical behavior related to fault reactivation. The numerical model will be enhanced by continuing collaboration and interaction with other research teams of DECOLVAEX-2019 Task B and validated using the field data from fault activation experiments in a further study.