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

• Tunnel and Underground Space
• /
• v.29 no.6
• /
• pp.439-455
• /
• 2019

A buffer is the major component of a high level radioactive waste repository. Due to their thermal conductivity and low permeability, bentonites have been considered as a key component of a buffer system in most countries. The deep geological condition generates ground water inflow and results in swelling pressure in the buffer and backfill. Investigation of swelling pressure of bentonite buffer is an important task for the safe disposal system. The swelling pressure that can be critical is affected by mechanical and hydro properties of the system. Therefore, in this study, a sensitivity analysis was conducted to examine the effect of hydro-mechanical (HM) behaviors in the MX-80 bentonite. Based on the results of the swelling pressure generation with HM model parameters, a coupled HM analysis of an unsaturated buffer and backfill in a deep geological repository was also carried out to investigate the major factor of the swelling pressure generation.

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

• 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.17 no.1 s.66
• /
• pp.56-65
• /
• 2007

A new underground LNG storage concept in the rock mass has been developed by combining underground cavern construction and new ice-ring harrier technologies with the conventional cryogenic insulation system. Technical feasibility of the storage system has been verified through construction and operation of the pilot storage cavern and a full-scale project is expected to start in the near future. One of the most important issues in the LNG storage system is the operational efficiency of the storage to minimize heat loss during a long period of operation due to the cryogenic heat transfer. This paper presents several important results of heat transfer and coupled hydro-thermal analyses by a finite element code Temp/W and Seep/W. A series of heat transfer analyses for full-scale caverns were performed to determine design parameters such as boil-off gas ratio (BOR), insulation thickness and pillar width. The result of the coupled hydro-mechanical analysis showed that BOR for underground storage system remains at about 0.04 %/day during the early stage of the operation. This value could be even much lower when the discontinuities in the rock masses are taken into consideration.

• Journal of the Korean Geotechnical Society
• /
• v.36 no.12
• /
• pp.17-25
• /
• 2020

The artificial ground freezing method is a ground amelioration technology that does not have a permanent effect on the ground. One of the key factors that determine the efficiency and design criteria of the artificial ground freezing is the groundwater flow. Therefore, in order to accurately evaluate the behavior of the artificial ground freezing, studies on the effect of water flow on the formation of ice walls must be preceded. In this paper, experimental and numerical analyses were conducted using only pure water to maximize the effect of water flow on the formation of ice walls. A hydro-thermal coupled model for freezing behavior was proposed and the accuracy of the model was verified. Through the numerical and experimental studies, the flow rate dominates not only the formation time but also the shape of the ice wall. In addition, this study proposes a method to indirectly predict the ice wall formation time, which is expected to be highly useful for a practical application where it is difficult to visually identify ice walls.

• 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.30 no.4
• /
• pp.359-381
• /
• 2020

Within the framework of DECOVALEX-2019 Task D, full-scale engineered barriers experiment (FEBEX) at Grimsel Test Site was numerically simulated to investigate an applicability of implemented Barcelona basic model (BBM) into TOUGH2-MP/FLAC3D simulator, which was developed for the prediction of the coupled thermo-hydro-mechanical behavior of bentonite buffer. And the calculated heater power, temperature, relative humidity, total stress, saturation, water content and dry density were compared with in situ data monitored in the various sections. In general, the calculated heater power and temperature provided a fairly good agreement with experimental observations, however, the difference between power of heater #1 and that of heater #2 could not captured in the numerical analysis. It is necessary to consider lamprophyre with low thermal conductivity around heater #1 and non-simplified installation progresses of bentonite blocks in the tunnel for better modeling results. The evolutions and distributions of relative humidity were well reproduced, but hydraulic model needs to be modified because the re-saturation process was relatively fast near the heaters. In case of stress evolutions due to the thermal and hydraulic expansions, the computed stress was in good agreement with the data. But, the stress is slightly higher than the measured in situ data at the early stage of the operation, because gap between rock mass and bentonite blocks have not been considered in the numerical simulations. The calculated distribution of saturation, water content, and dry density along the radial distance showed good agreement with the observations after the first and final dismantling. The calculated dry density near the center of the FEBEX tunnel and heaters were overestimated compared with the observations. As a result, the saturation and water content were underestimated with the measurements. Therefore, numerical model of permeability is needed to modify for the production of better numerical results. It will be possible to produce the better analysis results and more realistically predict the coupled THM behavior in the bentonite blocks by performing the additional studies and modifying the numerical model based on the results of this study.

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

• Journal of the Korean GEO-environmental Society
• /
• v.25 no.7
• /
• pp.29-34
• /
• 2024

Frost heave is one of the significant engineering characteristics of frozen ground and causes severe damages on geo-structures. Although thermo-hydro coupled analyses have been developed to predict frost heave behavior, these analyses involve excessive input parameters and have primarily been validated for frost heave in clayey soils. Frost heave mainly occurs in silty soils, which have relatively higher permeability compared to clayey soils, necessitating careful attention. This study introduces empirical models and verifies their reliability for silty soils. By using the validated model, the correlation of key input parameters is derived, which is expected to enhance the applicability of thermal-mechanical analysis for geo-structures on frozen ground in the future.