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

A new underground LNG storage concept in the rock mass has been developed by combining underground cavern construction and new ice-ring harrier technologies with the conventional cryogenic insulation system. Technical feasibility of the storage system has been verified through construction and operation of the pilot storage cavern and a full-scale project is expected to start in the near future. One of the most important issues in the LNG storage system is the operational efficiency of the storage to minimize heat loss during a long period of operation due to the cryogenic heat transfer. This paper presents several important results of heat transfer and coupled hydro-thermal analyses by a finite element code Temp/W and Seep/W. A series of heat transfer analyses for full-scale caverns were performed to determine design parameters such as boil-off gas ratio (BOR), insulation thickness and pillar width. The result of the coupled hydro-mechanical analysis showed that BOR for underground storage system remains at about 0.04 %/day during the early stage of the operation. This value could be even much lower when the discontinuities in the rock masses are taken into consideration.

• Journal of Nuclear Fuel Cycle and Waste Technology(JNFCWT)
• /
• v.17 no.2
• /
• pp.183-202
• /
• 2019

For design and performance assessment of a high-level radioactive waste (HLW) disposal system, it is necessary to understand the characteristics of coupled thermo-hydro-mechanical (THM) behavior. However, in previous studies for the Korean Reference HLW Disposal System (KRS), thermal analysis was performed to determine the spacing of disposal tunnels and interval of disposition holes without consideration of the coupled THM behavior. Therefore, in this study, TOUGH2-MP/FLAC3D is used to conduct THM modeling for performance assessment of the Korean Reference HLW Disposal System (KRS). The peak temperature remains below the temperature limit of $100^{\circ}C$ for the whole period. A rapid rise of temperature caused by decay heat occurs in the early years, and then temperature begins to decrease as decay heat from the waste decreases. The peak temperature at the bentonite buffer is around $96.2^{\circ}C$ after about 3 years, and peak temperature at the rockmass is $68.2^{\circ}C$ after about 17 years. Saturation of the bentonite block near the canister decreases in the early stage, because water evaporation occurs owing to temperature increase. Then, saturation of the bentonite buffer and backfill increases because of water intake from the rockmass, and bentonite buffer and backfill are fully saturated after about 266 years. The stress is calculated to investigate the effect of thermal stress and swelling pressure on the mechanical behavior of the rockmass. The calculated stress is compared to a spalling criterion and the Mohr-Coulumb criterion for investigation of potential failure. The stress at the rockmass remains below the spalling strength and Mohr-Coulumb criterion for the whole period. The methodology of using the TOUGH2-MP/FLAC3D simulator can be applied to predict the long-term behavior of the KRS under various conditions; these methods will be useful for the design and performance assessment of alternative concepts such as multi-layer and multi-canister concepts for geological spent fuel repositories.

• 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.29 no.6
• /
• pp.468-479
• /
• 2019

It is essential to comprehend coupled hydro-mechanical behavior to utilize subsurface for the recent demand for underground space usage. In this study, we developed a new simulator for numerical simulation as a tool for researching to consider the various domestic field and subsurface conditions. To develop the new module, we combined OpenGeoSys, one of the scientific software package that handles fluid mechanics (H), thermodynamics (T), and rock and soil mechanics (M) in the subsurface with FLAC3D, one of the commercial software for geotechnical engineering problems reinforced. In this simulator development, we design OpenGeoSys as a master and FLAC3D as a slave via a file-based sequential coupling. We have chosen Terzaghi's consolidation problem related to single-phase fluid flow at a saturated condition as a benchmark model to verify the proposed module. The comparative results between the analytical solution and numerical analysis showed a good agreement.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.8 no.4
• /
• pp.355-363
• /
• 2006

In this study, it was analyzed how the pore pressure behind a subsea tunnel influences on the stability of the tunnel. The tunnel is located in the soft rock layer, and a soft sandy layer and weathered soil layer are located on the top of it. Coupled numerical analyses are performed for both drained and undrained condition with varying coefficients of lateral earth pressure. In the case of undrained conditions, the stability of the tunnel was analyzed with different thicknesses of shotcrete. On the other hand, a sensitivity analysis was performed with different hydraulic conductivities and porosities of the shotcrete for the drained conditions. The stability of a subsea tunnel was evaluated in terms of safety factor suggested by You et al.(2000, 2001, 2005) based on the shear strength reduction technique. In this paper, the safety factor of a tunnel was calculated under steady state flow condition during hydro-mechanical coupled analysis. As a result, it was found that the stability of a subsea tunnel could be rather increased by allowing a proper amount of groundwater inflow into a subsea tunnel.

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

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.33 no.4
• /
• pp.1489-1498
• /
• 2013

The monolithically coupled finite element analysis for a deformable unsaturated soil slope is performed to investigate the effect of antecedent rainfall which is assumed by initial conditions varying degree of saturation (36, 51, 77%) in finite element analysis. The distributions of matric suction and deformation on slope surface obtained from numerical simulation show the instability of antecedent rainfall-induced unsaturated soil slope. Moreover, the numerical analysis using Drucker-Prager model can be checked if a soil slope has reached failure (trial failure criterion $f^{tr}$ >0, plastic behavior) or not (trial failure criterion $f^{tr}$ < 0, elastic behavior). It is found that displacement of slope surface layer increases and the matric suction on soil slope decreases with an increase of initial degree of saturation by antecedent rainfall. Especially, the matric suction of the soil slope in dry condition (S=36%) rapidly decreases rather than that in wet condition (S=51%) at the same rainfall duration. The results of the trial failure criterion ($f^{tr}$ > 0) show slope instability in the toe region and surface of the slopes.

• Tunnel and Underground Space
• /
• v.30 no.1
• /
• pp.39-62
• /
• 2020

In this study, Barcelona Basic Model (BBM) was implemented into TOUGH2-MP/FLAC3D for the numerical analysis of coupled thermo-hydro-mechanical (THM) behavior of unsaturated soils and the prediction of long-term behaviors. Similar to the methodology described in a previous study for the implementation of BBM into TOUGH-FLAC, the User Defined Model (UDM) of FLAC based on the Modified Cam Clay Model (MCCM) and the FISH function of FLAC3D were used to extend the existing MCCM module in FLAC3D for the implementation of BBM into TOUGH2-MP/FLAC3D. In the developed BBM module in TOUGH2-MP/FLAC3D, the plastic strains due to change in suction increase (SI) in addition to mean effective stress are calculated. In addition to loading-collapse (LC) yield surface, suction increase (SI) yield surface is changed by hardening rules in the developed BBM module. Several numerical simulations were conducted to verify and validate the implementation of BBM: using an example presented in the FLAC3D manual for the standard MCCM, simulation results using COMSOL, and experimental data presented in SKB Reports. In addition, the developed BBM analysis module was validated by simultaneously performing a series of modeling tests that were performed for the validation of the Quick tools developed for the purpose of effectively deriving BBM parameters, and by comparing the Quick tools and Code_Bright results reported in a previous study.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.12 no.2
• /
• pp.145-155
• /
• 2010

Recently interest in subsea tunnels is increasing nationwide and the construction of a subsea tunnel is taking place. For the stability of such a subsea tunnel, grouting is necessary for the water barrier and reinforcement of the tunnel. In this study, therefore, it was investigated how the grouting reinforcement had an effect on the stability of a subsea tunnel located in a great depth. To this end, Hydro-mechanical coupled analyses were performed for a sensitivity analysis in terms of different grouting range, rock class, shotcrete thickness, coefficient of lateral earth pressure, grouting thickness, and pumping existence for the rock classes I, III, and V. FLAC-2D ver. 5.0 was used for the numerical analyses. It was came to the conclusion that the effect of the increased water pressure due to the water barrier of the grouting should be considered as well as the strength improved effect in designing grouting reinforcement of subsea tunnels.

• Tunnel and Underground Space
• /
• v.29 no.4
• /
• pp.262-279
• /
• 2019

In the engineering barriers of high-level radioactive waste disposal, gases could be generated through a number of processes. If the gas production rate exceeds the gas diffusion rate, the pressure of the gas increases and gases could migrate through the bentonite buffer. Because people and the environment can be exposed to radioactivity, it is very important to clarify gas migration in terms of long-term integrity of the engineered barrier system. In particular, it is necessary to identify the hydro-mechanical mechanism for the dilation flow, which is a very important gas flow phenomenon only in medium containing large amounts of clay materials such as bentonite buffer, and to develop and validate new numerical approach for the quantitative evaluation of the gas migration phenomenon. Therefore, in this study, we developed a two-phase flow model considering the mechanical damage model in order to simulate the gas migration in the engineered barrier system, and validated with 1D gas flow modelling through saturated bentonite under constant volume boundary conditions. As a result of numerical analysis, the rapid increase in pore water pressure, stress, and gas outflow could be simulated when the dilation flow was occurred.