• Tunnel and Underground Space
• /
• v.32 no.2
• /
• pp.144-159
• /
• 2022

The present study developed a sequential approach-based numerical simulator for modeling coupled thermal-hydrological-mechanical (THM) processes in the ground and investigated the computational performance of the coupling analysis algorithm. The present sequential approach linked the two different solvers: an open-source numerical code, OpenGeoSys for solving the thermal and hydrological processes in porous media and a commercial code, FLAC3D for solving the geomechanical response of the ground. A benchmark test of the developed simulator was carried out using a THM problem where an analytical solution is given. The benchmark problem involves the coupled behavior (variations in temperature, pore pressure, stress, and deformation with time) of a fully saturated porous medium which is subject to a point heat source. The results of the analytical solution and numerical simulation were compared and the validity of the numerical simulator was investigated.

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

• Tunnel and Underground Space
• /
• v.21 no.1
• /
• pp.66-81
• /
• 2011

In this study, it was aimed to develop a thermal-hydraulic-mechanical coupled fracture mechanics code that models a fracture initiation, propagation and failure of underground rock mass due to thermal and hydraulic loadings. The development was based on a 2D FRACOD (Shen & Stephasson, 1993), and newly developed T-M and H-M coupled analysis modules were implemented into it. T-M coupling in FRACOD employed a fictitious heat source and time-marching method, and explicit iteration method was used in H-M coupling. The validity of developed coupled modules was verified by the comparison with the analytical result, and its applicability to the fracture initiation and propagation behavior due to temperature changes and hydraulic fracturing was confirmed by test simulations.

• Advances in Energy Research
• /
• v.6 no.1
• /
• pp.75-90
• /
• 2019

Anthropogenic greenhouse gas emissions are rising rapidly despite efforts to curb release of such gases. One long term potential solution to offset these destructive emissions is the capture and storage of carbon dioxide. Partially depleted hydrocarbon reservoirs are attractive targets for permanent carbon dioxide disposal due to proven storage capacity and seal integrity, existing infrastructure. Optimum well completion design in depleted reservoirs requires understanding of prominent geomechanics issues with regard to rock-fluid interaction effects. Geomechanics plays a crucial role in the selection, design and operation of a storage facility and can improve the engineering performance, maintain safety and minimize environmental impact. In this paper, an integrated geomechanics workflow to evaluate reservoir caprock integrity is presented. This method integrates a reservoir simulation that typically computes variation in the reservoir pressure and temperature with geomechanical simulation which calculates variation in stresses. Coupling between these simulation modules is performed iteratively which in each simulation cycle, time dependent reservoir pressure and temperature obtained from three dimensional compositional reservoir models in ECLIPSE were transferred into finite element reservoir geomechanical models in ABAQUS and new porosity and permeability are obtained using volumetric strains for the next analysis step. Finally, efficiency of this approach is demonstrated through a case study of oil production and subsequent carbon storage in an oil reservoir. The methodology and overall workflow presented in this paper are expected to assist engineers with geomechanical assessments for reservoir optimum production and gas injection design for both natural gas and carbon dioxide storage in depleted reservoirs.

• Journal of Radiation Industry
• /
• v.17 no.4
• /
• pp.333-344
• /
• 2023

To validate the numerical model used in the study of deep disposal of spent nuclear fuel, we selected benchmark cases and performed numerical model validation. We selected the DECOVALEX-THMC Task D_THM1 FEBEX Type benchmark case, which was conducted from 2003 to 2007. We analyzed the thermal-hydraulic (TH) behavior using the finite element program CODE_BRIGHT and verified the results against previous studies. The temperature results were similar to the results of DECOVALEX-THMC Task D. The saturation results showed a similar trend to the results of DECOVALEX-THMC Task D, but the time to reach full saturation was different.

• Journal of the Korean Geotechnical Society
• /
• v.32 no.6
• /
• pp.17-26
• /
• 2016

Characteristic behaviors of geo-structure during freezing and thawing process have to be understood based on fundamental knowledge on phase change in porous soil and interaction between soil and structure. Inversion analysis using published one-dimensional soil freezing tests was conducted to suggest a mechanical model to consider an effect of the ice saturation on Young's modulus. Silty soil was more sensitive to temperature than weathered granite soil and sand, and weathered granite soil was more affected by initial water saturation in stiffness decrease than silty soil. Numerical simulations on chilled gas pipeline showed that shielding effect from surrounding frozen zone around the pipe decreases impact from external load onto the pipe. And a pipe installed in sand backfill showed more heaving due to relatively low stiffness of sand during freezing than that of surrounding in-situ weather granite soil. However, it had more stable stress condition due to effective stress redistribution from external load.

• 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 the Korean Geotechnical Society
• /
• v.19 no.5
• /
• pp.59-67
• /
• 2003

In this study, the change in hydraulic conductivity according to the changes in the contact area, aperture, confining pressure and temperature was observed to improve the reliability of the analysis of underground water flow. Also, the mechanical and thermal properties of domestic crystalline rocks in a great depth were obtained. It was found that the averaged intial aperture ranged from 544.33${\mu}{\textrm}{m}$ to 898.62${\mu}{\textrm}{m}$ and it followed a log-normal distribution. The hydraulic conductivity decreased with the increase of normal stress on the fracture surface and the hydraulic conductivity decreased as temperature increased. The change in hydraulic conductivity was strongly correlated with the change in contact area. It was verified by experiments that hydraulic conductivity was inversely proportional to the contact area. The measured mechanical and thermal properties were very close to the existing typical properties of domestic granites.

• Journal of the Korean Geotechnical Society
• /
• v.33 no.3
• /
• pp.49-62
• /
• 2017

In order to understand response of geo-structures to the freezing-thawing process in the ground, it is necessary to consider phase change of the pore water of the ground and also to understand soil interaction with structures. In this study, numerical analysis was carried out for freezing and thawing effect on the modular road system. Neumann's theoretical equation for freezing-thawing processes in porous media can be used to estimate frozen depth and heaving from basic soil properties and ground and surface temperature, but its application is limited to the case for the sediment with fully saturated condition and zero unfrozen water content. Numerical analysis of the modular road system was performed on various soil types and different ground water table as the varying freezing index. The amount of heaving in the silty soil was much larger than those in granite weathered soil or sandy soil, and lowering groundwater level reduced ground heaving induced by freezing. Numerical analysis for temperature history of the ground surface predicted residual heaving near the surface by the freeze-thaw process in silty soil. It ought to reduce stiffness and bearing capacity of the ground so that it will impair stability and serviceability of new road system. However, the amount of residual heaving was insignificant for the road system installed in weathered soil granite and sandy soil. Since modular road system is a pavement structure mounted on the supporting substructure unlike the prevalent road pavement system, strict criteria should be applied for uniform and differential settlement of the pavement system.