• Tunnel and Underground Space
• /
• v.25 no.2
• /
• pp.155-167
• /
• 2015

The thermal-hydrological-mechanical (T-H-M) behavior of rock mass surrounding a high-temperature cavern thermal energy storage (CTES) operated for a period of 30 years has been investigated by TOUGH2-FLAC3D simulator. As a fundamental study for the development of prediction and control technologies for the environmental change and rock mass behavior associated with CTES, the key concerns were focused on the hydrological-thermal multiphase flow and the consequential mechanical behavior of the surrounding rock mass, where the insulator performance was not taken into account. In the present study, we considered a large-scale cylindrical cavern at shallow depth storing thermal energy of $350^{\circ}C$. The numerical results showed that the dominant heat transfer mechanism was the conduction in rock mass, and the mechanical behavior of rock mass was influenced by thermal factor (heat) more than hydrological factor (pressure). The effective stress redistribution, displacement and surface uplift caused by heating of rock and boiling of ground-water were discussed, and the potential of shear failure was quantitatively examined. Thermal expansion of rock mass led to the ground-surface uplift on the order of a few centimeters and the development of tensile stress above the storage cavern, increasing the potential of shear failure.

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

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

• 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 Korea Water Resources Association Conference
• /
• 2023.05a
• /
• pp.222-222
• /
• 2023

오픈 소스 유체역학 소프트웨어인 OpenFOAM은 다양한 유체 흐름에 적용 가능한 프로그램들로 구성되어 있다. 이 중 interFoam은 밀도가 다른 두 유체(i.e., 물, 공기) 간의 경계를 추적하는 기법을 기반으로 한 프로그램으로, 파랑의 거동 모의에 주로 쓰이고 있다. 파생형 프로그램으로는 동격자(dynamic mesh) 및 중첩 격자 기법(overset grid method)을 interFoam에 추가한 overInterDyMFoam이 있다. 두 기법을 사용하면 각각 여러 영역에서 유체흐름과 다중 물체 간의 상호작용을 효율적으로 모의할 수 있다. 본 연구에서는 overInterDyMFoam을 사용하여 두 개 수문의 개방 움직임을 구현하고 생성된 파랑이 포말대(swash zone)에 접근하였을 때의 흐름 특성을 조사하였다. 수치모형실험 결과 수문 개방 속도가 댐 붕괴 파랑 흐름 전파속도에 영향을 미치는 사실을 발견하였다. 또한, 처오름과 처내림의 상호작용에 의한 난류 운동 특성을 조사하기 위해 수문 개방시간 간격을 0초~3초로 설정하였다. 수치모형실험 결과는 수리모형실험의 수면 변동 시계열과 속도 시계열 결과와 비교하여 모형의 정확성이 검증되었다.

• Tunnel and Underground Space
• /
• v.16 no.1 s.60
• /
• pp.38-49
• /
• 2006

LNG storage in lined rock cavern demands various techniques concerned with rock mechanics, thermo-mechanics and hydrogeology in design, construction and maintenance stage. LNG pilot cavern was constructed in Daejeon in order to verify these techniques. In this paper, evaluation of drainage system and ice ring formation was studied by numerical simulation. By Modflow analysis in the viewpoint of aquifer and Seep/W analysis in the viewpoint of flow system, it was verified that the drainage system in the pilot cavern was efficiently operated. Since ice ring formation can be simulated by interactive relation between heat transfer and water flow, coupled analysis of those was performed. In this analysis, the position of ice ring was presumed and it was demonstrated that the formation is affected by velocity and direction of groundwater flow.

• Tunnel and Underground Space
• /
• v.11 no.3
• /
• pp.257-269
• /
• 2001

Numerical simulation is performed to understand the structural behavior of an underground radwaste repository, assumed to be located at the depth of 500 m, in a granitic rock mats, in which a fault intersects the roof of the repository cavern. Two dimensional universal distinct element code, UDEC is used in the analysis. The numerical model includes a granitic rock mass, a canister with PWR spent fuels surrounded by the compacted bentonite inside the deposition hole, and the mixed bentonite backfilled in the rest of the space within the repository cavern. The structural behavior of three different cases, each case with a fault of an angle of $33^{\circ},\;45^{\circ},\;and\;58^{\circ}$ passing through the cavern roof-wall intersection, has been compared. And then fro the case with the $45^{\circ}$ fault, the hydro-mechanical, thermo-mechanical, and thermo-hydro-mechanical interaction behavior have been studied. The effect of the time-dependent decaying heat, from the radioactive materials in PWR spent fuels, on the repository and its surroundings has been studied. The groundwater table is assumed to be located 10m below the ground surface, and a steady state flow algorithm is used.

• Tunnel and Underground Space
• /
• v.32 no.1
• /
• pp.30-58
• /
• 2022

The deep geological repository for high-level radioactive waste disposal is a multi barrier system comprised of engineered barriers and a natural barrier. The long-term integrity of the deep geological repository is affected by the coupled interactions between the individual barrier components. Erosion and piping phenomena in the compacted bentonite buffer due to buffer-rock interactions results in the removal of bentonite particles via groundwater flow and can negatively impact the integrity and performance of the buffer. Rapid groundwater inflow at the early stages of disposal can lead to piping in the bentonite buffer due to the buildup of pore water pressure. The physiochemical processes between the bentonite buffer and groundwater lead to bentonite swelling and gelation, resulting in bentonite erosion from the buffer surface. Hence, the evaluation of erosion and piping occurrence and its effects on the integrity of the bentonite buffer is crucial in determining the long-term integrity of the deep geological repository. Previous studies on bentonite erosion and piping failed to consider the complex coupled thermo-hydro-mechanical-chemical behavior of bentonite-groundwater interactions and lacked a comprehensive model that can consider the complex phenomena observed from the experimental tests. In this technical note, previous studies on the mechanisms, lab-scale experiments and numerical modeling of bentonite buffer erosion and piping are introduced, and the future expected challenges in the investigation of bentonite buffer erosion and piping are summarized.

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