• Tunnel and Underground Space
• /
• v.34 no.2
• /
• pp.127-142
• /
• 2024

The robustness of a numerical method means that its computational performance is maintained under various modeling conditions. New numerical methods or codes need to be assessed for robustness through benchmark testing. The TOUGH-FLAC modeling approach has been applied to various fields such as subsurface carbon dioxide storage, geological disposal of spent nuclear fuel, and geothermal development both domestically and internationally, and the modeling validity has been examined by comparing the results with experimental measurements and other numerical codes. In the present study, a benchmark test of the TOUGH-FLAC approach was performed based on a coupled thermal-hydro-mechanical behavior problem with an analytical solution. The analytical solution is related to the temperature, pore water pressure, and mechanical behavior of a fully saturated porous medium that is subjected to a point heat source. The robustness of the TOUGH-FLAC approach was evaluated by comparing the analytical solution with the results of numerical simulation. Additionally, the effects of thermal-hydro-mechanical coupling terms, fluid phase change, and timestep on the computation of coupled behavior were investigated.

• 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.31 no.3
• /
• pp.167-183
• /
• 2021

It is essential to understand the complex thermo-hydro-mechanical-chemical (THMC) coupled behavior in the engineered barrier system and natural barrier system to secure the high-level radioactive waste repository's long-term safety. The heat from the high-level radioactive waste induces thermal pressurization and vaporization of groundwater in the repository system. Groundwater inflow affects the saturation variation in the engineered barrier system, and the saturation change influences the heat transfer and multi-phase flow characteristics in the buffer. Due to the complexity of the coupled behavior, a numerical simulation is a valuable tool to predict and evaluate the THMC interaction effect on the disposal system and safety assessment. To enhance the knowledge of THMC coupled interaction and validate modeling techniques in geological systems. DECOVALEX, an international cooperation project, was initiated in 1992, and KAERI has participated in the projects since 2008 in Korea. In this study, we introduced the main contents of all tasks in the DECOVALEX-2023, the current DECOVALEX phase, to the rock mechanics and geotechnical researchers in Korea.

• Journal of Nuclear Fuel Cycle and Waste Technology(JNFCWT)
• /
• v.17 no.2
• /
• pp.183-202
• /
• 2019

For design and performance assessment of a high-level radioactive waste (HLW) disposal system, it is necessary to understand the characteristics of coupled thermo-hydro-mechanical (THM) behavior. However, in previous studies for the Korean Reference HLW Disposal System (KRS), thermal analysis was performed to determine the spacing of disposal tunnels and interval of disposition holes without consideration of the coupled THM behavior. Therefore, in this study, TOUGH2-MP/FLAC3D is used to conduct THM modeling for performance assessment of the Korean Reference HLW Disposal System (KRS). The peak temperature remains below the temperature limit of $100^{\circ}C$ for the whole period. A rapid rise of temperature caused by decay heat occurs in the early years, and then temperature begins to decrease as decay heat from the waste decreases. The peak temperature at the bentonite buffer is around $96.2^{\circ}C$ after about 3 years, and peak temperature at the rockmass is $68.2^{\circ}C$ after about 17 years. Saturation of the bentonite block near the canister decreases in the early stage, because water evaporation occurs owing to temperature increase. Then, saturation of the bentonite buffer and backfill increases because of water intake from the rockmass, and bentonite buffer and backfill are fully saturated after about 266 years. The stress is calculated to investigate the effect of thermal stress and swelling pressure on the mechanical behavior of the rockmass. The calculated stress is compared to a spalling criterion and the Mohr-Coulumb criterion for investigation of potential failure. The stress at the rockmass remains below the spalling strength and Mohr-Coulumb criterion for the whole period. The methodology of using the TOUGH2-MP/FLAC3D simulator can be applied to predict the long-term behavior of the KRS under various conditions; these methods will be useful for the design and performance assessment of alternative concepts such as multi-layer and multi-canister concepts for geological spent fuel repositories.

• Tunnel and Underground Space
• /
• v.32 no.6
• /
• pp.518-529
• /
• 2022

In Task C of the DECOVALEX-2023 project, nine institutes from six nations are developing their numerical codes to simulate thermo-hydro-mechanical coupled behavior for the FE experiment performed at Mont Terri underground rock laboratory, Switzerland. Currently, Step 1 for comparing the simulation results to field data is the ongoing stage, and we used the OGS-FLAC simulator for a series of numerical simulations. As a result, temperature increase depending on the heating hysteresis was well simulated, and saturation variation in the bentonite depending on phase change was observed. However, due to the suction overestimation, relative humidity and temperature change in the bentonite and the pressure variation in the Opalinus clay showed a difference compared to the field data. From the observation, it is confirmed that the effect of the bentonite capillary pressure is dominant to the flow analysis in the disposal system. We further plan to draw improved results considering tunnel support material and accurate initial water pressure distribution. Additionally, the thermal, hydrological, and mechanical anisotropy of the Opalinus clay was well simulated. From the simulation results, we confirmed the applicability of the OGS-FLAC simulator in the disposal system analysis.

• Tunnel and Underground Space
• /
• v.29 no.1
• /
• pp.38-51
• /
• 2019

Barcelona Basic Model (BBM) can describe not only swelling owing to decrease in effective stress, but also wetting-induced swelling due to decrease in suction. And the BBM can also consider increase in cohesion and apparent preconsolidation stress with suction, and decrease in the apparent preconsolidation stress with temperature. Therefore, the BBM is widely used all over the world to predict and to analyze coupled thermo-hydro-mechanical behavior of bentonite which is considered as buffer materials at the engineered barrier system in the high-level radioactive waste disposal system. However, the BBM is not well known in Korea, so this paper introduce the BBM to Korean rock engineers and geotechnical engineers. In this study, Modified Cam Clay (MCC) model is introduced before all, because the BBM was first developed as an extension of the MCC model to unsaturated soil conditions. Then, the thermo-elasto-plastic version of the BBM is described in detail.

• Tunnel and Underground Space
• /
• v.25 no.2
• /
• pp.168-185
• /
• 2015

The thermal-hydrological-mechanical (T-H-M) behavior of rock mass surrounding a large-scale high-temperature cavern thermal energy storage (CTES) at a shallow depth has been investigated, and the effects of hydrological conditions such as water table and rock permeability on the behavior have been examined. The liquid saturation of ground water around a storage cavern may have a small impact on the overall heat transfer and mechanical behavior of surrounding rock mass for a relatively low rock permeability of $10^{-17}m^2$. In terms of the distributions of temperature, stress and displacement of the surrounding rock mass, the results expected from the simulation with the cavern below the water table were almost identical to that obtained from the simulation with the cavern in the unsaturated zone. The heat transfer in the rock mass with reasonable permeability ${\leq}10^{-15}m^2$ was dominated by the conduction. In the simulation with rock permeability of $10^{-12}m^2$, however, the convective heat transfer by ground-water was dominant, accompanying the upward heat flow to near-ground surface. The temperature and pressure around a storage cavern showed different distributions according to the rock permeability, as a result of the complex coupled processes such as the heat transfer by multi-phase flow and the evaporation of ground-water.

• Tunnel and Underground Space
• /
• v.29 no.6
• /
• pp.468-479
• /
• 2019

It is essential to comprehend coupled hydro-mechanical behavior to utilize subsurface for the recent demand for underground space usage. In this study, we developed a new simulator for numerical simulation as a tool for researching to consider the various domestic field and subsurface conditions. To develop the new module, we combined OpenGeoSys, one of the scientific software package that handles fluid mechanics (H), thermodynamics (T), and rock and soil mechanics (M) in the subsurface with FLAC3D, one of the commercial software for geotechnical engineering problems reinforced. In this simulator development, we design OpenGeoSys as a master and FLAC3D as a slave via a file-based sequential coupling. We have chosen Terzaghi's consolidation problem related to single-phase fluid flow at a saturated condition as a benchmark model to verify the proposed module. The comparative results between the analytical solution and numerical analysis showed a good agreement.

• 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.

• Journal of the Korean Geotechnical Society
• /
• v.38 no.9
• /
• pp.45-52
• /
• 2022

Because discontinuity in the rock mass and contact of soil-structure interaction exhibits coupled thermal-hydromechanical (THM) behavior, it is necessary to develop an interface element based on the full governing equations. In this study, we derive force equilibrium, fluid continuity, and energy equilibrium equations for the interface element. Additionally, we present a stiffness matrix of the elastoplastic mechanical model for the interface element. The developed interface element uses six nodes for displacement and four nodes for water pressure and temperature in a two-dimensional analysis. The fully coupled THM analysis for fluid injection into a fault can model the complicated evolution of injection pressure due to decreasing effective stress in the fault and thermal contraction of the surrounding rock mass. However, the result of hydromechanical analysis ignoring thermal phenomena overestimates hydromechanical variables.