• 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.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.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.31 no.6
• /
• pp.561-577
• /
• 2021

The compacted bentonite buffer in a geological repository for high-level radioactive waste disposal is saturated due to groundwater inflow. Saturation of the bentonite buffer results in bentonite swelling and bentonite penetration into the rock discontinuities present around the disposal hole. The penetrated bentonite is exposed to groundwater flow and can be eroded out of the repository, resulting in bentonite mass loss which can affect the physical integrity of the engineered barrier system. Hence, the evaluation of buffer-rock interactions and coupled behavior due to groundwater inflow and bentonite penetration is necessary to ensure long-term disposal safety. In this study, the effects of the bentonite penetration and swelling on the physical properties of jointed rock mass were evaluated using the quasi-static resonant column test. Jointed rock specimens with bentonite penetration were manufactured using Gyeongju bentonite and hollow cylindrical granite rock discs obtained from the KAERI underground research tunnel. The effects of vertical stress and saturation were assessed using the P-wave and S-wave velocities for intact rock, jointed rock and jointed rock with bentonite penetration specimens. The joint normal and joint shear stiffnesses of each joint condition were inferred from the wave velocity results assuming an equivalent continuum. The joint normal and joint shear stiffnesses obtained from this study can be used as input factors for future numerical analysis on the performance evaluation of geological waste disposal considering rock discontinuities.

• 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.33 no.3
• /
• pp.189-207
• /
• 2023

In the present study, the thermoshearing experiment on a rough rock fracture were modeled using a three-dimensional grain-based distinct element model (GBDEM). The experiment was conducted by the Korea Institute of Construction Technology to investigate the progressive shear failure of fracture under the influence of thermal stress in a critical stress state. The numerical model employs an assembly of multiple polyhedral grains and their interfaces to represent the rock sample, and calculates the coupled thermo-mechanical behavior of the grains (blocks) and the interfaces (contacts) using 3DEC, a DEM code. The primary focus was on simulating the temperature evolution, generation of thermal stress, and shear and normal displacements of the fracture. Two fracture models, namely the mated fracture model and the unmated fracture model, were constructed based on the degree of surface matedness, and their respective behaviors were compared and analyzed. By leveraging the advantage of the DEM, the contact area between the fracture surfaces was continuously monitored during the simulation, enabling an examination of its influence on shear behavior. The numerical results demonstrated distinct differences depending on the degree of the surface matedness at the initial stage. In the mated fracture model, where the surfaces were in almost full contact, the characteristic stages of peak stress and residual stress commonly observed in shear behavior of natural rock joints were reasonably replicated, despite exhibiting discrepancies with the experimental results. The analysis of contact area variation over time confirmed that our numerical model effectively simulated the abrupt normal dilation and shear slip, stress softening phenomenon, and transition to the residual state that occur during the peak stress stage. The unmated fracture model, which closely resembled the experimental specimen, showed qualitative agreement with the experimental observations, including heat transfer characteristics, the progressive shear failure process induced by heating, and the increase in thermal stress. However, there were some mismatches between the numerical and experimental results regarding the onset of fracture slip and the magnitudes of fracture stress and displacement. This research was conducted as part of DECOVALEX-2023 Task G, and we expect the numerical model to be enhanced through continued collaboration with other research teams and validated in further studies.

• Coupled systems mechanics
• /
• v.13 no.3
• /
• pp.201-217
• /
• 2024

This study aims to develop a new model to obtain the minimum area in circular isolated footings with eccentric column taking into account that the surface in contact with the ground works partially in compression, i.e., a part of the contact area of the footing is subject to compression and the other there is no pressure (pressure zero). The new model is formulated from a mathematical approach based on a minimum area, and it is developed by integration to obtain the axial load "P", moment around the X axis "M_x" and moment around the Y axis "M_y" in function of σ_max (available allowable soil pressure) R (radius of the circular footing), α (angle of inclination where the resultant moment appears), y₀ (distance from the center of the footing to the neutral axis measured on the axis where the resultant moment appears). The normal practice in structural engineering is to use the trial and error procedure to obtain the radius and area of the circular footing, and other engineers determine the radius and area of circular footing under biaxial bending supported on elastic soils, but considering a concentric column and the contact area with the ground works completely in compression. Three numerical problems are given to determine the lowest area for circular footings under biaxial bending. Example 1: Column concentric. Example 2: Column eccentric in the direction of the X axis to 1.50 m. Example 3: Column eccentric in the direction of the X axis to 1.50 m and in the direction of the Y axis to 1.50 m. The new model shows a great saving compared to the current model of 44.27% in Example 1, 50.90% in Example 2, 65.04% in Example 3. In this way, the new minimum area model for circular footings will be of great help to engineers when the column is located on the center or edge of the footing.

• Tunnel and Underground Space
• /
• v.31 no.4
• /
• pp.289-308
• /
• 2021

This study conducts coupled thermo-hydro-mechanical numerical modeling to investigate the maximum temperature and conditions for securing mechanical stability of the high-level radioactive waste repository when temperature criteria of bentonite buffer are 100℃ and 125℃, respectively. In case of temperature criterion of buffer as 100℃, the maximum temperatures at the interface between canister and buffer are calculated to be 99.4℃ and 99.8℃, respectively for a case with disposal tunnel spacing of 40 m and deposition hole spacing of 5.5 m and for the other case with disposal tunnel spacing of 30 m and deposition hole spacing of 6.5 m. In case of temperature criterion of buffer as 125℃, spacings of disposal tunnel and deposition hole could be decreased to 30 m and 4.5 m, respectively, which reduces the disposal area up to 55% compared to the disposal area of KRS⁺. According to analysis of mechanical stability for various disposal spacings, RMR of rock mass for KRS⁺ should be larger than 72.4 which belongs to good rock in RMR classification to prevent failure of rock mass. As disposal spacing is decreased, required RMR of rock mass is increased. In order to prevent failure of rock mass for a case with disposal tunnel spacing of 30 m and deposition hole spacing of 4.5 m, RMR larger than 87.3 is needed. However, mechanical stability of the repository is secured for all cases with RMR over 75 considering the enhancement of rock strength due to confining stress induced by swelling of the bentonite buffer and backfill.

• Tunnel and Underground Space
• /
• v.32 no.1
• /
• pp.30-58
• /
• 2022

The deep geological repository for high-level radioactive waste disposal is a multi barrier system comprised of engineered barriers and a natural barrier. The long-term integrity of the deep geological repository is affected by the coupled interactions between the individual barrier components. Erosion and piping phenomena in the compacted bentonite buffer due to buffer-rock interactions results in the removal of bentonite particles via groundwater flow and can negatively impact the integrity and performance of the buffer. Rapid groundwater inflow at the early stages of disposal can lead to piping in the bentonite buffer due to the buildup of pore water pressure. The physiochemical processes between the bentonite buffer and groundwater lead to bentonite swelling and gelation, resulting in bentonite erosion from the buffer surface. Hence, the evaluation of erosion and piping occurrence and its effects on the integrity of the bentonite buffer is crucial in determining the long-term integrity of the deep geological repository. Previous studies on bentonite erosion and piping failed to consider the complex coupled thermo-hydro-mechanical-chemical behavior of bentonite-groundwater interactions and lacked a comprehensive model that can consider the complex phenomena observed from the experimental tests. In this technical note, previous studies on the mechanisms, lab-scale experiments and numerical modeling of bentonite buffer erosion and piping are introduced, and the future expected challenges in the investigation of bentonite buffer erosion and piping are summarized.

• Tunnel and Underground Space
• /
• v.30 no.4
• /
• pp.382-393
• /
• 2020

The engineered barrier system of high-level radioactive waste disposal must maintain its performance in the long term, because it must play a role in slowing the rate of leakage to the surrounding rock mass even if a radionuclide leak occurs from the canister. In particular, it is very important to clarify gas dilation flow phenomenon clearly, that occurs only in a medium containing a large amount of clay material such as a bentonite buffer, which can affect the long-term performance of the bentonite buffer. Accordingly, DECOVALEX-2019 Task A was conducted to identify the hydraulic-mechanical mechanism for the dilation flow, and to develop and verify a new numerical analysis technique for quantitative evaluation of gas migration phenomena. In this study, based on the conventional two-phase flow and mechanical behavior with effective stresses in the porous medium, the hydraulic-mechanical model was developed considering the concept of damage to simulate the formation of micro-cracks and expansion of the medium and the corresponding change in the hydraulic properties. Model verification and validation were conducted through comparison with the results of 1D and 3D gas injection tests. As a result of the numerical analysis, it was possible to model the sudden increase in pore water pressure, stress, gas inflow and outflow rate due to the dilation flow induced by gas pressure, however, the influence of the hydraulic-mechanical interaction was underestimated. Nevertheless, this study can provide a preliminary model for the dilation flow and a basis for developing an advanced model. It is believed that it can be used not only for analyzing data from laboratory and field tests, but also for long-term performance evaluation of the high-level radioactive waste disposal system.