• Tunnel and Underground Space
• /
• v.7 no.1
• /
• pp.13-19
• /
• 1997

It is necessary to study on thermo-hydro-mechanical effect at rock mass performing project such as radiowaste disposal in deep rock mass. In this study, thermo-hydro-mechanical coupling analysis which is considered interaction and the variation of rock properties induced by temperature increase was performed for the circular shaft when appling temperature of 20$0^{\circ}C$ at the shaft wall. The shaft is diameter of 2 m and under hydrostatic stress of 5 MPa. In the cases, thermal expansion by temperature increase progress from the wall to outward and thermal expansion could induce tensile stress over the tensile strength of rock mass at the wall. When rock properties were given as a function of temperature, thermal expansion increased, tensile stress zone expanded. Lately, water flow is activated by increase of permeability and decrease of viscosity.

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

As large underground projects such as radioactive waste disposal, hot water and heat storage, and geothermal energy become influential, the study, which consider all aspects of thermics, hydraulics and mechanics would be needed. Thermo Hydro-Mechanical coupling analysis is one of the most complex numerical technique because it should be implemented with the combined three governing equations to analyze the behavior of rock mass. In this study, finite element code, which is based on Biot's consolidation theory, was developed to analyze the thermo-hydro-mechanical coupling in continuum rock mass. To verify the implemented program, one-dimensional consolidation model under the isothermal and non-isothermal conditions was analyzed and was compared with the analytic solution. The parametric study on two-dimensional consolidation was also performed and the effects of several factors such as poisson's ratio and hydraulic anisotropy on rock mass behavior were investigated. In the future, this program would be revised to be used for analysis of general discontinuous media with incorporating discrete joint model.

• Tunnel and Underground Space
• /
• v.10 no.3
• /
• pp.355-365
• /
• 2000

As large underground projects such as radioactive waste disposal, hot water and heat storage, and geothermal energy become influential, the study, which consider all aspects of thermics, hydraulics and mechanics would be needed. Thermo-Hydro-Mechanical coupling analysis is one of the most complex numerical technique because it should be implemented with the combined three governing equations to analyze the behavior of rock mass. In this study, finite element code, which is based on Biot's consolidation theory, was developed to analyze the thermo-hydro-mechanical coupling in continuum rock mass. To verify the implemented program, one-dimensional consolidation model under the isothermal and non-isothermal conditions was analyzed and was compared with the analytic solution. The parametric study on two-dimensional consolidation was also performed and the effects of several factors such as poisson's ratio and hydraulic anisotropy on rock mass behavior were investigated. In the future, this program would be revised to be used for analysis of general discontinuous media with incorporating discrete joint model.

• Tunnel and Underground Space
• /
• v.8 no.3
• /
• pp.184-193
• /
• 1998

In order to dispose radioactive wastes safely, it is needed to understand the mechanical, thermal, fluid behavior of rockmass and physico-chemical interactions between rockmass and water. Also, the knowledge about mechanical and hydraulic properties of rocks is required to predict and to model many conditions of geological structure, underground in-situ stress, folding, hot water interaction, intrusion of magma, plate tectonics etc. This study is based on researches about rock mechanics issues associated with a waste disposal in deep rockmass. This paper includes the mechanical and hydraulic behavior of rocks in varying temperature conditions, thermo-hydro-mechanical coupling analysis in rock mass and deformation behavior of discontinuous rocks. The mechanical properties were measured with Interaken rock mechanics testing systems and hydraulic properties were measured with transient pulse permeability measuring systems. In all results, rock properties were sensitive to temperature variation.

• Coupled systems mechanics
• /
• v.9 no.2
• /
• pp.125-145
• /
• 2020

In this paper, we present a 3D thermo-hydro-mechanical coupled discrete beam lattice model of structure built of the nonisothermal saturated poro-plastic medium subjected to mechanical loads and nonstationary heat transfer conditions. The proposed model is based on Voronoi cell representation of the domain with cohesive links represented as inelastic Timoshenko beam finite elements enhanced with additional kinematics in terms of embedded strong discontinuities in axial and both transverse directions. The enhanced Timoshenko beam finite element is capable of modeling crack formation in mode I, mode II and mode III. Mode I relates to crack opening, mode II relates to in-plane crack sliding, and mode III relates to the out-of-plane shear sliding. The pore fluid flow and heat flow in the proposed model are governed by Darcy's law and Fourier's law for heat conduction, respectively. The pore pressure field and temperature field are approximated with linear tetrahedral finite elements. By exploiting nodal point quadrature rule for numerical integration on tetrahedral finite elements and duality property between Voronoi diagram and Delaunay tetrahedralization, the numerical implementation of the coupling results with additional pore pressure and temperature degrees of freedom placed at each node of a Timoshenko beam finite element. The results of several numerical simulations are presented and discussed.

• Tunnel and Underground Space
• /
• v.9 no.1
• /
• pp.1-11
• /
• 1999

A finite element code was developed to analyze coupled thermo-hydro-mechanical phenomena. This code is based on the finite element formulation provided by Noorishad et al. (1984) and Joint behavior was simulated Goodman's joint constitutive model. The developed code was applied for T-H-M coupling analysis for two kinds of shaft models, with a joint or without a joint respectively. For a model without a joint, temperature increased from the shaft wall to outward evidently. The radial displacement showed opposite directions of outward and inward at some distance from shaft wall. For a model with a joint, closure of joint was found due to thermal expansion. The temperature distribution along a joint showed relatively lower than that of rock matrix because of low thermal conductivity and high specific heat of water. And it could be concluded that effects of thermal flow to joint were more than that of hydraulic flow in a rock mass.

• Tunnel and Underground Space
• /
• v.8 no.1
• /
• pp.17-25
• /
• 1998

This paper investigates coupled thermal, mechanical and hydraulic phenomena in deep rock mass especially for underground heat storage system. Firstly, concepts of underground heat storage were presented and coupling phenomena in this area were illustrated. In order to understand the basic mechanism of thermal, hydraulic and deformation behavior in rock cavern disturbed by thermal gradient about 10$0^{\circ}C$, various numerical experiments were conducted using several codes. The study involves the behavior of fractured rock mass including rock joint. In spite of the limitation of codes modelling fully coupled effects, these codes could be applied in analysis of underground heat storage. The heat loss in rock mass, which is a major factor in heat storage, is insignificant in all results.

• Proceedings of the Korean Geotechical Society Conference
• /
• 2010.09a
• /
• pp.808-812
• /
• 2010

A great deal of attention is focused on coupled Thermo-Hydro-Mechanical (THM) behavior of multiphase porous media 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 3 mass balances equations (solid, water, and air), energy balance equation, and force equilibrium equation. Finite element code is developed from the Galerkin formulation and time integration of these governing equations for 4 main variables (displacement $\underline{u}$, gas pressure $P_g$, liquid pressure $P_l$), and temperature T). The code is validated with theoretical solutions for linear material with simple boundary conditions.

• Tunnel and Underground Space
• /
• v.17 no.1 s.66
• /
• pp.32-42
• /
• 2007

Severe damage can occur around deposition holes due to complex interaction of thermo-hydro-mechanical (THM) loading during the long term operation of high level radioactive waste repository. Many candidate sites for repository are located in crystalline rock mass, therefore mechanism of damage follows the form of brittle fracture and failure. This paper briefly introduces major outcomes from 15 years international collaborative project, DECOVALEX, and presents major study results for current ongoing benchmark test study from DECOVALEX-THMC, to evaluate the effect of THM loading to rock mass in excavation damaged zone (EDZ) near deposition holes. Through benchmark test model by simplifying THM loading to boundary loading obtained numerical results are compared, and discrete fracture interaction after up to 1 million years operation is discussed.

• Journal of the Korean Geotechnical Society
• /
• v.29 no.8
• /
• pp.65-73
• /
• 2013

Ground-coupled heat pump system has attracted attention as a promising renewable energy technology due to its improving energy efficiency and eco-friendly mechanism for space cooling and heating. Pipes buried in the ground play a role of direct thermal interaction between circulating fluid inside the pipe and surrounding soils in the geothermal exchange system. However, both complexities of turbulent flow coupling thermal-hydraulic phenomena and very long aspect ratio of the pipe make it difficult to model the heat exchange system directly. Energy balance for fluid flow inside the pipe was derived to model thermal-hydraulic phenomena, and one-dimensional pipe element was proposed through Galerkin formation and time integration of the equation. Developed element is combined to pre-developed FEM code for THM phenomena in porous media. Numerical results of Thermal Response Test showed that line-source model overestimates equivalent thermal conductivity of surrounding soils due to thermal interaction between adjacent pipes and finite length of the pipe. Thus, inverse analysis for the TRT simulation was conducted to present optimal transformation matrix with utmost convergence.