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

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

• 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.33 no.2
• /
• pp.35-47
• /
• 2017

Thermal behavior of soils is mainly focused on thermal conduction, and the study of natural convection is very limited. Increase of soil temperature causes natural convection due to buoyancy from density change of pore water. The limitations of the analysis using fluid dynamics for natural convection in the porous media is discussed and a new numerical analysis is presented for natural convection in porous media using THM governing equations fully coupled in the macroscopic view. Numerical experiments for thermal probe show increase in the uncertainty of thermal conductivity estimated without considering natural convection, and suggest appropriate experimental procedures to minimize errors between analytical model and numerical results. Burial of submarine power cable should not exceed the temperature changes of $2^{\circ}C$ at the depth of 0.2 m under the seabed, but numerical analysis for high permeable ground exceeds this criterion. Temperature and THM properties of the seafloor are important design factors for the burial of power cable, and in this case effects of natural convection should be considered. Especially, in the presence of heat sources in soils with high permeability, natural convection due to the variation of density of pore water should be considered as an important heat transfer mechanism.

• Tunnel and Underground Space
• /
• v.33 no.5
• /
• pp.312-338
• /
• 2023

The underground disposal facility for spent nuclear fuel demands a specialized design, distinct from conventional practices, to ensure long-term thermal, mechanical, and hydraulic integrity, preventing the release of radioactive isotopes from high-temperature spent nuclear fuel. SKB has established design criteria for such facilities and executed practical design implementations for Forsmark. Moreover, in response to subsurface uncertainty, SKB has proposed an empirical approach involving monitoring and adaptive design modifications, alongside stepwise development. SKB has further introduced a unique support system, categorizing ground types and behaviors and aligning them with corresponding support types to confirm safety through comparative analyses against existing systems. POSIVA has pursued a comparable approach, developing a support system for Onkalo while accounting for distinct geological characteristics compared to Forsmark. This demonstrates the potential for domestic implementation of spent nuclear fuel disposal facility designs and the establishment of a support system adapted to national attributes.