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

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

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.10 no.4
• /
• pp.313-327
• /
• 2008

Thermo-mechanical coupled behavior of an underground structure during a fire accident have not been fully understood yet. Moreover, when such a thermo-mechanical coupled behavior is not considered in numerical analyses based on conventional heat transfer theory, fire-induced damage zone in an underground structure can be considerably underestimated. This study aims to develop a FEM-based numerical technique to simulate the thermo-mechanical coupled behavior of an underground structure in a fire accident. Especially, an element elimination model is newly proposed to simulate fire-induced structural loss together with a convective boundary condition. In the proposed model, an element where the maximum temperature calculated from heat transfer analysis is over a prescribed critical temperature is eliminated. Then, the proposed numerical technique is verified by comparing numerical results with experimental results from real fire model tests. From a series of parametric studies, the key parameters such as critical temperature, element size and temperature-dependent convection coefficients are optimized for the RABT and the RWS fire scenarios.

• 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.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 Society for Railway
• /
• v.13 no.2
• /
• pp.194-200
• /
• 2010

Recently, a longitudinal tunnel construction has increased because of subway construction extension, geomorphological effect and the development of construction Technologies etc. When the fire occurs in a tunnel and an underground structure, the many damage of human life and the economic losses are caused. In Korea, fire resistance character study of a tunnel and an underground structure is proceeding. However, when a concrete is exposed to high temperature, study of load carrying capacity reduction and stability evaluation for spalling of a concrete is not enough. Therefore in this study, fire resistance character of a concrete evaluated according to time heating temperature curve(RABT and RWS) and a result compared on virtual fire accident in order to apply fire scenario. Also this study performed thermo-mechanical coupled analysis of a FEM-based numerical technique and estimated fire-induced damage of a tunnel and an underground structure.

• Composites Research
• /
• v.28 no.5
• /
• pp.277-284
• /
• 2015

This paper presents a method to detect the fiber property variation of laminated GFRP plates from natural frequency response data. The combined finite element analysis using ABAQUS and the inverse algorithm described in this paper may allow us not only to detect the deteriorated elements from the mirco-mechanical point of view but also to find their numbers, locations, and the extent of damage. To solve the inverse problem using the combined method, this study uses several natural frequencies instead of mode shapes in a structure as the measured data. Several numerical results show that the proposed system is computationally efficient in identifying fiber stiffness degradation for complex structures such as composites with various layup sequences.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2011.04a
• /
• pp.89-92
• /
• 2011

The modeling of thermodynamic non-idealities and transport anomalies is a crucial prerequisite to realistically simulate the mixing and combustion processes of liquid propellants injected above critical pressures. This study has developed a specific set of subroutines to calculate the thermodynamic and transport properties based on the generalized cubic equation of state (EoS) in a coupled manner with the standard chemical kinetics packages (Chemkin). The existing flamelet analysis code is extended with the real-fluid package and applied to numerical investigation of local flame structures of kerosene and liquid oxygen at high pressure conditions relevant to the actual rocket engines.

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

• Tunnel and Underground Space
• /
• v.25 no.2
• /
• pp.168-185
• /
• 2015

The thermal-hydrological-mechanical (T-H-M) behavior of rock mass surrounding a large-scale high-temperature cavern thermal energy storage (CTES) at a shallow depth has been investigated, and the effects of hydrological conditions such as water table and rock permeability on the behavior have been examined. The liquid saturation of ground water around a storage cavern may have a small impact on the overall heat transfer and mechanical behavior of surrounding rock mass for a relatively low rock permeability of $10^{-17}m^2$. In terms of the distributions of temperature, stress and displacement of the surrounding rock mass, the results expected from the simulation with the cavern below the water table were almost identical to that obtained from the simulation with the cavern in the unsaturated zone. The heat transfer in the rock mass with reasonable permeability ${\leq}10^{-15}m^2$ was dominated by the conduction. In the simulation with rock permeability of $10^{-12}m^2$, however, the convective heat transfer by ground-water was dominant, accompanying the upward heat flow to near-ground surface. The temperature and pressure around a storage cavern showed different distributions according to the rock permeability, as a result of the complex coupled processes such as the heat transfer by multi-phase flow and the evaporation of ground-water.