• Journal of Korean Tunnelling and Underground Space Association
• /
• v.22 no.4
• /
• pp.419-434
• /
• 2020

In this study, the damage parameters of Mazars model for KURT (KAERI Underground Research Tunnel) granite are measured form uniaxial compressive and Brazilian tests under the simulated coupled condition of a deep geological disposal. The tests are conducted in three different temperatures (15℃, 45℃, and 75℃) and dry/saturated conditions. Major model parameters such as maximum effective tensile strain (𝜖_d0), A_t, B_t, A_c, and B_c differ from the typical reference values of concrete specimens. This is likely due to the difference in elastic modulus between rock and concrete. It is found that the saturation of specimens causes an increase in value of B_t and B_c while, the rise in temperature increases 𝜖_d0 and B_t and decreases B_c. The damage model obtained from this study will be used as the primary input parameters in the development of coupled Thermo-Hydro-Mechanical Damage numerical model in KAERI.

• 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.30 no.1
• /
• pp.39-62
• /
• 2020

In this study, Barcelona Basic Model (BBM) was implemented into TOUGH2-MP/FLAC3D for the numerical analysis of coupled thermo-hydro-mechanical (THM) behavior of unsaturated soils and the prediction of long-term behaviors. Similar to the methodology described in a previous study for the implementation of BBM into TOUGH-FLAC, the User Defined Model (UDM) of FLAC based on the Modified Cam Clay Model (MCCM) and the FISH function of FLAC3D were used to extend the existing MCCM module in FLAC3D for the implementation of BBM into TOUGH2-MP/FLAC3D. In the developed BBM module in TOUGH2-MP/FLAC3D, the plastic strains due to change in suction increase (SI) in addition to mean effective stress are calculated. In addition to loading-collapse (LC) yield surface, suction increase (SI) yield surface is changed by hardening rules in the developed BBM module. Several numerical simulations were conducted to verify and validate the implementation of BBM: using an example presented in the FLAC3D manual for the standard MCCM, simulation results using COMSOL, and experimental data presented in SKB Reports. In addition, the developed BBM analysis module was validated by simultaneously performing a series of modeling tests that were performed for the validation of the Quick tools developed for the purpose of effectively deriving BBM parameters, and by comparing the Quick tools and Code_Bright results reported in a previous 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.

• Journal of Nuclear Fuel Cycle and Waste Technology(JNFCWT)
• /
• v.4 no.2
• /
• pp.117-131
• /
• 2006

A coupled T-H-M(Thermo-Hydro-Mechanical) analysis was carried out for KENTEX (KAERI Engineering-scale T-H-M Experiment for Engineered Barrier System), which is a facility for validating the coupled T-H-M behavior in the engineered barrier system of the Korean reference HLW(high-level waste) disposal system. The changes of temperature, water saturation, and stress were estimated based on the coupled T-H-M analysis, and the influence of the types of mechanical constitutive material laws was investigated by using elastic model, poroelastic model, and poroelastic-plastic model. The analysis was done using ABAQUS, which is a commercial finite element code for general purposes. From the analysis, it was observed that the temperature in the bentonite increased sharply for a couple of days after heating the heater and then slowly increased to a constant value. The temperatures at all locations were nearly at a steady state after about 37.5 days. In the steady state, the temperature was maintained at $90^{\circ}C$ at the interface between the heater and the bentonite and at about $70^{\circ}C$ at the interface between the bentonite and the confining cylinder. The variation of the water saturation with time in bentonite was almost same independent of the material laws used in the coupled T-H-M processes. By comparing the saturation change of T-H-M and that of H-M(Hydro-Mechanical) processes using elastic and poroelastic material mod31 respectively, it was found that the degree of saturation near the heater from T-H-M calculation was higher than that from the coupled H-M calculation mainly because of the thermal flux, which seemed to speed up the saturation. The stresses in three cases with different material laws were increased with time. By comparing the stress change in H-M calculation using poroelasetic and poroelasetic-plastic model, it was possible to conclude that the influence of saturation on the stress change is higher than the influence of temperature. It is, therefore, recommended to use a material law, which can model the elastic-plastic behavior of buffer, since the coupled T-H-M processes in buffer is affected by the variation of void ratio, thermal expansion, as well as swelling pressure.

• International Journal of Concrete Structures and Materials
• /
• v.11 no.2
• /
• pp.199-213
• /
• 2017

Chloride penetration is among the main causes of corrosion initiation in reinforced concrete (RC) structures producing premature degradations. Weather and exposure conditions directly affect chloride ingress mechanisms and therefore the operational service life and safety of RC structures. Consequently, comprehensive chloride ingress models are useful tools to estimate corrosion initiation risks and minimize maintenance costs for RC structures placed under chloride-contaminated environments. This paper first presents a coupled thermo-hydro-chemical model for predicting chloride penetration into concrete that accounts for realistic weather conditions. This complete numerical model takes into account multiple factors affecting chloride ingress such as diffusion, convection, chloride binding, ionic interaction, and concrete aging. Since the complete model could be computationally expensive for long-term assessment, this study also proposes model simplifications in order to reduce the computational cost. Long-term chloride assessments of complete and reduced models are compared for three locations in France (Brest, Strasbourg and Nice) characterized by different weather and exposure conditions (tidal zone, de-icing salts and salt spray). The comparative study indicates that the reduced model is computationally efficient and accurate for long-term chloride ingress modeling in comparison to the complete one. Given that long-term assessment requires larger climate databases, this research also studies how climate models may affect chloride ingress assessment. The results indicate that the selection of climate models as well as the considered training periods introduce significant errors for mid- and long- term chloride ingress assessment.

• Tunnel and Underground Space
• /
• v.18 no.3
• /
• pp.165-174
• /
• 2008

Thermo-hydro-chemico-mechanical evolution of the radioactive waste repository and surrounding geological media is one of the key issues for the radioactive waste disposal. This article describes not only the basic context for the site selection but also a reasonable strategy for the repository related research based on the results of the French repository project carried out by ANDRA (National radioactive waste management agency). To have some alternatives for the determination of a preferable depth and geological media, it would be recommendable to establish a database system. The curing process of the fractures or microfissures in the EDZ (Excavation Disturbed Zone) during operation time has to be examined considering the evolution of the EDZ and the reversibility of the repository. It is prerequisite to carry out a feasibility study and to validate the design concept and design parameters in a properly constructed underground research laboratory (URL) in Korea.

• Journal of the Korean Geotechnical Society
• /
• v.27 no.11
• /
• pp.55-69
• /
• 2011

The importance of unsaturated state in various geo-engineering problems has led to the advance of mechanical constitutive model emulating behavior of unsaturated soils in response to thermo-hydro-mechanical loading. Elasto-plastic mechanical constitutive model for unsaturated soil is formulated based on Bishop's effective stress. Effective stress and temperature are main variables in constitutive equation, and incremental formulation of constitutive relationship is derived to compute stress update and stiffness tensor. Numerical simulations involving coupled THM processes are conducted to discuss numerical stability and applicability of developed constitutive model: one-dimensional test, tri-axial compression test, and clay-buffering at high level radioactive waste disposal. Numerical results demonstrated that developed model can predict very complex behavior of coupled THM phenomena and is applicable to geo-engineering problems under various environmental conditions, as well as interpret typical behavior of unsaturated soils.

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

• Geomechanics and Engineering
• /
• v.16 no.4
• /
• pp.363-373
• /
• 2018

The evaluation of Thermo-Hydro-Mechanical (THM) coupling behavior is important for the development of underground space for various purposes. For a high-level radioactive waste repository excavated in a deep underground rock mass, the accurate prediction of the complex THM behavior is essential for the long-term safety and stability assessment. In order to develop reliable THM analysis techniques effectively, an international cooperation project, Development of Coupled models and their Validation against Experiments (DECOVALEX), was carried out. In DECOVALEX-2015 Task B2, the in situ THM experiment that was conducted at Horonobe Underground Research Laboratory(URL) by Japan Atomic Energy Agency (JAEA), was modeled by the research teams from the participating countries. In this study, a THM coupling technique that combined TOUGH2 and FLAC3D was developed and applied to the THM analysis for the in situ experiment, in which rock, buffer, backfill, sand, and heater were installed. With the assistance of an artificial neural network, the boundary conditions for the experiment could be adequately implemented in the modeling. The thermal, hydraulic, and mechanical results from the modeling were compared with the measurements from the in situ THM experiment. The predicted buffer temperature from the THM modelling was about $10^{\circ}C$ higher than measurement near by the overpack. At the other locations far from the overpack, modelling predicted slightly lower temperature than measurement. Even though the magnitude of pressure from the modeling was different from the measurements, the general trends of the variation with time were found to be similar.