• Tunnel and Underground Space
• /
• v.31 no.4
• /
• pp.221-242
• /
• 2021

A numerical study of the performance assesment of coupled thermo-hydro-mechanical (THM) processes in improved Korean reference disposal system (KRS⁺) for high-level radioactive waste is conducted using TOUGH2-MP/FLAC3D simulator. Decay heat from high-level radioactive waste increases the temperature of the repository, and it decreases as decay heat is reduced. The maximum temperature of the repository is below a maximum temperature criterion of 100℃. Saturation of bentonite buffer adjacent to the canister is initially reduced due to pore water evaporation induced by temperature increase. Bentonite buffer is saturated 250 years after the disposal of high-level radioactive waste by inflow of groundwater from the surrounding rock mass. Initial saturation of rock mass decreases as groundwater in rock mass is moved to bentnonite buffer by suction, but rock mass is saturated after inflow of groundwater from the far-field area. Stress changes at rock mass are compared to the Mohr-Coulomb failure criterion and the spalling strength in order to investigate the potential rock failure by thermal stress and swelling pressure. Additional simulations are conducted with the reduced spacing of deposition holes. The maximum temperature of bentonite buffer exceeds 100℃ as deposition hole spacing is smaller than 5.5 m. However, temperature of about 56.1% volume of bentonite buffer is below 90℃. The methodology of numerical modeling used in this study can be applied to the performance assessment of coupled THM processes for high-level radioactive waste repositories with various input parameters and geological conditions such as site-specific stress models and geothermal gradients.

• Journal of the Korean Geotechnical Society
• /
• v.40 no.1
• /
• pp.55-63
• /
• 2024

In this study, a thermomechanical numerical model was developed to evaluate the stability of energy slabs. First, a wall-type energy slab was installed in a residential underground parking lot, and thermal performance tests were conducted. Based on the tests, a numerical thermohydraulics model of the energy slab was developed to accurately simulate the thermal behavior in thermal performance tests. Finally, utilizing the temperature data acquired using the developed model, a thermomechanical numerical model of the energy slab was established. The thermomechanical model was then used to simulate the thermal stresses induced by operating the energy slab. The results demonstrated a maximum thermal stress of 5,300 kPa, which highlights the need to utilize cement mortar with sufficient tensile strength to realize stable operation of the energy slab.

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

• Journal of the Korean Geotechnical Society
• /
• v.27 no.3
• /
• pp.75-83
• /
• 2011

Emerging issues related with fully coupled Thermo-Hydro-Mechanical (THM) behavior of unsaturated soil demand the development of a numerical tool in diverse geo-mechanical and geo-environmental areas. This paper presents general governing equations for coupled THM processes in unsaturated porous media. Coupled partial differential equations are derived from three mass balances equations (solid, water, and air), energy balance equation, and force equilibrium equation. With Galerkin formulation and time integration of these governing equations, finite element code is developed to find nonlinear solution of four main variables (displacement-u, gas pressure-$P_g$), liquid pressure-$P_1$), and temperature-T) using Newton's iterative scheme. Three cases of numerical simulations are conducted and discussed: one-dimensional drainage experiments (u-$P_g-P_1$), thermal consolidation (u-$P_1$-T), and effect of pile on surrounding soil due to surface temperature variation (u-$P_1$-T).

• 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.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.411-434
• /
• 2022

In deep geological disposal concept for spent nuclear fuel, it is well-known that rock mass at near-field experiences the thermal-hydraulic-mechanical (THM) coupled behavior. The mechanical properties of rock changes during the coupled process, and it is important to consider the changes into the analysis of numerical simulation and in-situ tests for long-term stability evaluation of nuclear waste disposal repository. This report collected the previous studies on indirect coupled behaviors of rock. The effects of water saturation and temperature on some mechanical properties of rock was considered, while the change in hydraulic conductivity of rock due to stress was included in the indirect coupled behavior.

• Tunnel and Underground Space
• /
• v.33 no.1
• /
• pp.10-28
• /
• 2023

Long-term shear behavior of the rock discontinuities should be analyzed and its stability should be evaluated to ensure the long-term stability of a high-level radioactive waste disposal repository. The long-term shear behavior of the discontinuities can be modeled with creep and RSF models. The shear creep test, velocity step test, and slide-hold-slide test can be performed to determine their model parameters or analyze the shear behavior by experiments under various conditions. Testing apparatuses for direct shear, triaxial compression, and biaxial shear were mainly used and improved to reproduce the thermo-hydro-mechanical conditions of local bedrock, and it was confirmed that the shear behavior could vary. In order to design a high-level radioactive waste disposal site in Korea, the long-term behavior of rock discontinuities should be investigated in consideration of rock types, thermo-hydro-mechanical conditions, metamorphism, and restoration of shear resistance.

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

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