• Journal of the Korean Geotechnical Society
• /
• v.36 no.8
• /
• pp.49-60
• /
• 2020

A fully coupled thermo-hydro-mechanical model is established to evaluate frost heave behaviour of saturated frost-susceptible soils. The method is based on mass conservation, energy conservation, and force equilibrium equations, which are fully coupled with each other. These equations consider various physical phenomena during one-dimensional soil freezing such as latent heat of phase change, thermal conductivity changes, pore water migration, and the accompanying mechanical deformation. Using the thermo-hydro-mechanical model, a sensitivity analysis study is conducted to examine the effects of the geotechnical parameters and external conditions on the amount of frost heave and frost heaving rate. According to the results of the sensitivity analysis, initial void ratio significantly affects each objective as an individual parameter, whereas soil particle thermal conductivity and temperature gradient affect frost heave behaviour to a greater degree when applied simultaneously. The factors considered in this study are the main factors affecting the frost heaving amount and rate, which may be used to determine the frostbite sensitivity of a new sample.

• Journal of the Korean Geotechnical Society
• /
• v.29 no.1
• /
• pp.81-92
• /
• 2013

This paper presents the results of a numerical investigation into the behavior of retaining wall subject to cycles of freezing and thawing due to seasonal temperature change. The thermo-hydro-mechanical coupled finite element modeling strategy was first established to simulate the wall behavior. A series of finite element analyses were then performed on a range of conditions representing seasonal temperature change characteristics. The results indicated that the average freezing temperature and the number of cycles of freezing and thawing were the primary influencing factors for the wall behavior. Also revealed was that the duration of freezing period does not significantly affect the wall displacement and the lateral earth pressure, and that the earth pressure on the wall does not significantly change due to the freezing and thawing action suggesting that the increase in the wall displacement during the freezing and thawing action may be attributed to degradation of backfill due to the freezing and thawing action.

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

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

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