• Tunnel and Underground Space
• /
• v.30 no.6
• /
• pp.573-590
• /
• 2020

This study presents the current status of DECOVALEX-2023 project Task G and our research results so far. Task G, named 'Safety ImplicAtions of Fluid Flow, Shear, Thermal and Reaction Processes within Crystalline Rock Fracture NETworks (SAFENET)' aims at developing a numerical method to simulate the fracture creation and propagation, and the coupled thermohydro-mechanical processes in fracture in crystalline rocks. The first research step of Task G is a benchmark simulation, which is designed for research teams to make their modelling codes more robust and verify whether the models can represent an analytical solution for displacements of a single rock fracture. We reproduced the mechanical behavior of rock and embedded single fracture using a three-dimensional grain-based distinct element model for the simulations. In this method, the structure of the rock was represented by an assembly of rigid tetrahedral grains moving independently of each other, and the mechanical interactions at the grains and their contacts were calculated using 3DEC. The simulation results revealed that the stresses induced along the embedded fracture in the model were relatively low compared to those calculated by stress analysis due to stress redistribution and constrained fracture displacements. The fracture normal and shear displacements of the numerical model showed good agreement with the analytical solutions. The numerical model will be enhanced by continuing collaboration and interaction with other research teams of DECOVALEX-2023 Task G and validated using various experiments in a further study.

• Tunnel and Underground Space
• /
• v.21 no.5
• /
• pp.394-405
• /
• 2011

In this paper, we applied coupled non-isothermal, multiphase fluid flow and geomechanical numerical modeling using TOUGH-FLAC coupled analysis to study the complex thermodynamic and geomechanical performance of underground lined rock caverns (LRC) for compressed air energy storage (CAES). Mechanical stress in concrete linings as well as pressure and temperature within a storage cavern were examined during initial and long-term operation of the storage cavern for CAES. Our geomechanical analysis showed that effective stresses could decrease due to air penetration pressure, and tangential tensile stress could develop in the linings as a result of the air pressure exerted on the inner surface of the lining, which would result in tensile fracturing. According to the simulation in which the tensile tangential stresses resulted in radial cracks, increment of linings' permeability and air leakage though the linings, tensile fracturing occurred at the top and at the side wall of the cavern, and the permeability could increase to $5.0{\times}10^{-13}m^2$ from initially prescribed $10{\times}10^{-20}m^2$. However, this air leakage was minor (about 0.02% of the daily air injection rate) and did not significantly impact the overall storage pressure that was kept constant thanks to sufficiently air tight surrounding rocks, which supports the validity of the concrete-lined underground caverns for CAES.

• Tunnel and Underground Space
• /
• v.28 no.1
• /
• pp.111-124
• /
• 2018

The purpose of this study is to analyze ventilation design factor for network-type double-deck road tunnel that have been developed actively around the world. A numerical analysis was carried out through computational fluid dynamics (CFD) to derive shock loss coefficient that occurs due to the change in cross sectional area at both merging section and diverging section. The model used for the numerical analysis is real-scale model and the reliability of the result is secured by comparing with the coefficient of the previous studies. As a result of this study, shock loss coefficient was calculated depending on the change in cross-sectional area ratio and was higher than the result of previous studies in case of both merging section and diverging section. It is considered that the characteristics of the geometrical structure of network-type double-deck road tunnel have a great impact on shock loss coefficient. Therefore, the result of this study is expected to be helpful for more accurate ventilation design of network-type double-deck road tunnel.

• Tunnel and Underground Space
• /
• v.22 no.6
• /
• pp.429-437
• /
• 2012

This study predicted temperature in the disposal tunnels using computational fluid dynamics based on natural ventilation quantity that comes from high altitude and temperature differences that are the characteristics of high level waste repository. The result of the previous study that evaluated quantitatively natural ventilation quantity using a hydrostatic method and CFD shows that significant natural ventilation quantity is generated. From the result, this study performed the prediction of temperature in disposal tunnels by natural ventilation quantity by the caloric values of the wastes, at both deep geological repository and surface repository. The result of analysis shows that deep geological repository is effective for thermal control in the disposal tunnels due to heat transfer to rock and the generation of sufficient natural ventilation quantity, while surface repository was detrimental to thermal control, because surface repository was strongly affected by external temperature, and could not generate sufficient natural ventilation quantity. Moreover, this study found that in the case of deep geological repository with a depth of 200 m, the heatof about $10^{\circ}C$ was transferred to the depth of 500 m. Thus, it is considered that if the high level waste repository scheduled to be built in the country is designed placing an emphasis on thermal control, deep geological repository rather than surface repository is more appropriate.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.37 no.1
• /
• pp.1-7
• /
• 2017

In general, the large loads which are applied from explosion, impact, earthquake and wind at a short time caused the materials of structures to large deformations, rotations and strains locally. If such phenomena will be analyzed, hydrocodes which can be considered fluid-structure interaction under computational continuum mechanics are inevitably needed. Also, the explosion mechanism is so complicated, it is reasonable that the behaviors of structure are predicted through explosion analyses and experiment at the same time. But, unfortunately, it is true that explosion experiments are limited to huge cost, large experiment facilities and safety problems. Therefore, in this study, it is shown that the results of explosion analyses using the AUTODYN are agreed with those of existing explosion experiments for reinforced concrete slabs within reasonable error limits. And the explosion damage of the same reinforced concrete slab are assessed for quite different reinforcement arrangement spacings, concrete cover depths, and vertical reinforcements. From the explosion analyses, it is known that the more the ratio of slab thickness to reinforcement arrangement spacing is increased, and small-diameter reinforcements are used than large-diameter reinforcements on the same reinforcement ratio, and vertical reinforcements are used, the more the anti-knock capacities are improved.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.19 no.1
• /
• pp.81-96
• /
• 2007

Nonlinear shoaling characteristics over surf zone are numerically investigated based on spatially averaged NavierStokes equation. We also test the validity of gradient model for turbulent stresses due to wave breaking using the data acquainted during SUPERTANK LABORATORY DATA COLLECTION PROJECT(Krauss et al., 1992). It turns out that the characteristics length scale of breaking induced current is not negligible, which firmly stands against ever popular gradient model, ${\kappa}-{\varepsilon}$ model, but favors Large Eddy Simulation with finer grid. Based on these observations, we model the residual stress of spatially averaged NavierStokes equation after Lagrangian Dynamic Smagorinsky(Meneveau et al., 1996). We numerically integrate newly proposed wave equations using SPH with Gaussian kernel function. Severely deformed water surface profile, free falling water particle, queuing splash after landing of water particle on the free surface and wave finger due to structured vortex on rear side of wave crest(Narayanaswamy and Dalrymple, 2002) are successfully duplicated in the numerical simulation of wave propagation over uniform slope beach, which so far have been regarded very difficult features to mimic in the computational fluid mechanics.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.17 no.5
• /
• pp.346-353
• /
• 2016

In the design conditions of some research reactors, the siphon phenomenon can cause continuous efflux of water during pipe rupture. A siphon breaker is a safety device that can prevent water efflux effectively. However, the analysis of the siphon breaking is complicated because many variables must be included in the calculation process. For this reason, a simulation program was developed with a user-friendly GUI to analyze the siphon breaking easily. The program was developed by MFC programming using Visual Studio 2012 in Windows 8. After saving the input parameters from a user, the program proceeds with three steps of calculation using fluid mechanics formulas. Bernoulli's equation is used to calculate the velocity, quantity, water level, undershooting, pressure, loss coefficient, and factors related to the two-phase flow. The Chisholm model is used to predict the results from a real-scale experiment. The simulation results are shown in a graph, through which a user can examine the total breaking situation. It is also possible to save all of the resulting data. The program allows a user to easily confirm the status of the siphon breaking and would be helpful in the design of siphon breakers.

• Journal of the Korean Society of Safety
• /
• v.11 no.4
• /
• pp.143-150
• /
• 1996

This is an explicit-Implicit, finite element analysis for linear as well as nonlinear hygrothermal stress problems. Additional features, such as moisture diffusion equation, crack element and virtual crack extension(VCE ) method for evaluating J-integral are implemented in this program. The Linear Elastic Fracture Mechanics(LEFM) Theory is employed to estimate the crack driving force under the transient condition for and existing crack. Pores in materials are assumed to be saturated with moisture in the liquid form at the room temperature, which may vaporize as the temperature increases. The vaporization effects on the crack driving force are also studied. The Ideal gas equation is employed to estimate the thermodynamic pressure due to vaporization at each time step after solving basic nodal values. A set of field equations governing the time dependent response of porous media are derived from balance laws based on the mixture theory Darcy's law Is assumed for the fluid flow through the porous media. Perzyna's viscoplastic model incorporating the Von-Mises yield criterion are implemented. The Green-Naghdi stress rate is used for the invariant of stress tensor under superposed rigid body motion. Isotropic elements are used for the spatial discretization and an iterative scheme based on the full newton-Raphson method is used for solving the nonlinear governing equations.

• Tunnel and Underground Space
• /
• v.23 no.3
• /
• pp.241-259
• /
• 2013

A packed bed thermal energy storage (TES) consisting of solid storage medium of rock or concrete through which the heat transfer fluid is circulated is considered as an attractive alternative for high temperature sensible heat storage, because of the economical viability and chemical stability of storage medium and the simplicity of operation. This study introduces the technologies of packed bed thermal energy storage, and presents a numerical model to analyze the thermal energy balance and the performance efficiency of the storage system. In this model, one dimensional transient heat transfer problem in the storage tank is solved using finite difference method, and temperature distribution in a storage tank and thermal energy loss from the tank wall can be calculated during the repeated thermal charging and discharging modes. In this study, a high temperature thermal energy storage connected with AA-CAES (advanced adiabatic compressed air energy storage) was modeled and analyzed for the temperature and the energy balance in the storage tank. Rock cavern type TES and above-ground type TES were both simulated and their results were compared in terms of the discharging efficiency and heat loss ratio.

• Tunnel and Underground Space
• /
• v.29 no.3
• /
• pp.197-213
• /
• 2019

We simulated the fault reactivation experiment conducted at 'Main Fault' intersecting the low permeability clay formations of Mont Terri Underground Research Laboratory in Switzerland using TOUGH-FLAC simulator. The fluid flow along a fault was modelled with solid elements and governed by Darcy's law with the cubic law in TOUGH2, whereas the mechanical behavior of a single fault was represented by creating interface elements between two separating rock blocks in FLAC3D. We formulate the hydro-mechanical coupling relation of hydraulic aperture to consider the elastic fracture opening and failure-induced dilation for reproducing the abrupt changes in injection flow rate and monitoring pressure at fracture opening pressure. A parametric study was conducted to examine the effects of in-situ stress condition and fault deformation and strength parameters and to find the optimal parameter set to reproduce the field observations. In the best matching simulation, the fracture opening pressure and variations of injection flow rate and monitoring pressure showed good agreement with field experiment results, which suggests the capability of the numerical model to reasonably capture the fracture opening and propagation process. The model overestimated the fault displacement in shear direction and the range of reactivated zone, which was attributed to the progressive shear failures along the fault at high injection pressure. In the field experiment results, however, fracture tensile opening seems the dominant mechanism affecting the hydraulic aperture increase.