• Tunnel and Underground Space
• /
• v.29 no.5
• /
• pp.318-331
• /
• 2019

A three-dimensional(3D) equivalent continuum modeling was performed to analyze hydraulic behavior of rock mass considering discontinuities by using DFN model and smeared fracture model. DFN model was generated by FLAC3D and smeared fracture model was applied by using FISH functions, which is built-in functions in FLAC3D, for equivalent continuum model of fractured rock mass. Comparative analysis with 3DEC, which is for discontinuum analysis, was conducted to verify reliability of equivalent continuum analysis by using FLAC3D. Similar results of hydraulic analysis under the same conditions could be achieved. Equivalent continuum analysis of fractured rock mass by using DFN model was implemented to compare with existing analytical methods for inflow into the tunnel.

• 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.25 no.6
• /
• pp.556-567
• /
• 2015

The evaluation of THM coupling behavior in deep underground repository condition is essential for the long term safety and stability assessment of high-level radioactive waste repository. In order to develop reliable THM analysis techniques effectively, an international cooperation project, DECOVALEX, is carried out. In DECOVALEX-2015 Task B2, the in situ THM experiment planned to be conducted by JAEA was modeled by the research teams from the participating countries. In this study, a THM coupling technique combining TOUGH2 and FLAC was developed and applied to 1 dimensional THM modeling, in which rock, buffer, and heater are considered. The results were compared with those from other research teams.

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

• Tunnel and Underground Space
• /
• v.14 no.4
• /
• pp.295-303
• /
• 2004

For the deep underground excavation site with the geological complexity of soft seam and hard rock, the numerical analysis and in-situ measurement on the behaviors of roadway and surrounding rock according to stepwise excavation of the steep soft seam are carried out. The strata behavior is modeled using elasto-plastic FEM considering the empirical failure criteria of Hoek ＆ Brown and the strain-softening model. Hydraulic pressure capsule, MPBX and tape extensometer are installed around the roadway for the in-situ measurement of rock stress and deformation. Despite the complexity of geology and excavation procedure, the elasto-plastic analysis considering the empirical failure criteria of Hoek ＆ Brown and the strain-softening model shows good agreement with the in-situ measurement. Comparison of numerical modeling with in-situ measurement enables to predict the behaviors of the roadway and to obtain design parameters for the excavation and support at depth.

• Proceedings of the Korean Society for Rock Mechanics Conference
• /
• 2000.09a
• /
• pp.201-210
• /
• 2000

Rock bridge in rock masses can be considered as one of several types of opening-mode fractures, and also it has been known to have a great influence on the stability of structures in rock mats. In the beginning of researching a rock bridge it used to be studied only in characteristics of its behavior, as considering resistance of material itself. However the distribution pattern of rock bridges, which can affect the stability of rock structures, is currently researched with a fracture mechanical approach in numerical studies. For investigating the effect of rock bridges on the development pattern of hydraulic fractures, the author analyzed numerically the stress state transition in rock bridges and their phenomena with a different pattern of the rock bridge distributons. From the numerical studies, a two-crack configuration could be defined to be representative of the most critical conditions for rock bridges, only when cracks are systematic and same in their length and angle. Moreover, coalescence stresses and onset of propagation stresses could be known to increase with decreasing s/L ratio or increasing d/L ratio. The effect of pre-existing crack on hydraulic fracturing was studied also in numerical models. Different to the simple hydraulic fracturing modeling in which the fractures propagated exactly parallel to the maximum remote stress, the hydraulic fractures with pre-existing cracks dial not propagate parallel to the maximum remote stress direction. These are representative of the tendency to change the hydraulic fractures direction because of the existence of pre-existing crack. Therefore s/L, d/L ratios will be identical as a function effective on hydraulic fractures propagation, that is, the $K_{I}$ vague increase with decreasing s/L ratio or increasing d/L ratio and its magnification from onset to propagation increases with decreasing s/L ratio. The scanline is a commonly used method to estimate the fracture distribution on outcrops. The data obtained from the scanline method can be applied to the evaluation of stress field in rock mass.s.

• Tunnel and Underground Space
• /
• v.24 no.1
• /
• pp.46-54
• /
• 2014

Direct shear tests have been initiated to understand the characteristics of joints which crucially affect the stability of rock mass. In this research, numerical approach in direct shear tests has been initiated using 3DEC on the basis of 3D distinct element method. Normal loads were altered in four different levels on artificial joint tests depending on the sawtooth angle and strengths on constant normal stress conditions, measuring the peak shear strength according to the direct shear tests under laboratory condition. Also results obtained from mechanical properties through laboratory test were used to perform numerical modeling, and shear strength obtained from the modeling was used to compare with laboratory direct shear test. As a result numerical analysis from distinct element method can simulate well on the shear behavior of rockmass.

• 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.10 no.2
• /
• pp.165-172
• /
• 2000

When tunnels or underground structures are constructed in anisotropic rock mass, designers and constructors have to consider the anisotropic characteristics in rock mass because their physical and mechanical properties are depended on the anisotropic angles(${\beta}$). In this study, therefore, we have first investigated the mechanical behavior of the gneiss specimen from lab. tests, and then have analysed the behavior of specimens for to the transversely isotropic model in elastic medium using the FLAC program. The results of this study were summarized as follows; 1) In the result of the variation tests, in general, the properties of strength were depended on the angle of inclination in spite of the hard rock. And except for the shear strength test, the lowest and peak stress were appeared at 60$^{\circ}$ and 90$^{\circ}$respectively. 2) The results of specimen modeling analysis using FDM well indicated the mechanical behaviors of the specimen of transversely isotropic model.

• Tunnel and Underground Space
• /
• v.10 no.2
• /
• pp.237-248
• /
• 2000

Based on the preliminary results from the therm analysis, which is currently carrying, three-dimensional computer simulations using a finite element code, ABAQUS Ver. 5.8, were designed to determine the mechanically stable cavern and deposition hole spacing. Linear elastic modeling for the cases with different cavern and deposition hole spacing were carried out under three different in situ stress conditions. From the simulations, the response of the rock to the stress redistribution after the excavation of the openings could be investigated. Also the optimum cavern and deposition hole spacing could be estimated based on the factor of safety. When the in situ stress determined from the actual stress measurements in Korea were used, the case with cavern spacing of 40m and deposition hole spacing of 3m was in very stable condition, because the factor of safety was calculated as 3.42., When the in situ stress conditions for Sweden and Canada were used, the previous case, they seem to be in stable condition, since the factors of safety are still higher than 1.0. From these results, it was concluded that the rock will not fail even after the stress redistribution.