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

• Tunnel and Underground Space
• /
• v.23 no.1
• /
• pp.13-30
• /
• 2013

Increase of pore fluid pressure resulting from injection of $CO_2$ may reactivate pre-existing faults, and the induced seismic activities can raise the safety issues such as seal integrity, restoration of storage capacity, and, in the worst case, removal of previously injected $CO_2$. Thus, fault stability and potential for $CO_2$ leakage need to be assessed at the stage of site selection and planning of injection pressure, based on the results of large-scale site investigations and numerical modeling for various scenarios. In this report, studies on the assessment of fault stability during injection of $CO_2$ were reviewed. The seismic activities associated with an artificial injection of fluids or a release of naturally trapped high-pressure fluids were first examined, and then site investigation methods for the magnitude and orientation of in situ stresses, the distribution and change of pore fluid pressure, and the location of faults were generally summarized. Recent research cases on possibility estimation of fault reactivation, prediction of seismic magnitude, and modeling of $CO_2$ leakage through a reactivated fault were presented.

• Tunnel and Underground Space
• /
• v.27 no.4
• /
• pp.217-229
• /
• 2017

This study focuses on the improvement of the working environment in domestic collieries where temperature is increasing due to heat of the earth that is caused by the long-term mining. In order to improve the working environment of the mine, a ventilation evaluation was carried out for Hwasoon Mining Industry. In order to increase the ventilation efficiency of the mine, numerical analysis of the effect on temperature was carried out by using climsim, a temperature prediction program. The analysis shows that A coal mine needs $6,152m^3/min$ for in-flow ventilation rate but the total input air flowrate is $4,710m^3/min$, $1,442m^3/min$ of in-flow ventilation rate shortage. The 93 m hypothetical ventilation shaft from -395 ML to -488 ML could result about $3^{\circ}C$ temperature drop in the coal mine of -488 ML far. As a result of predicting the $CO_2$ concentration at -523 ML development using artificial neural network, the emission of $CO_2$ increased as the amount of coal and coal bed thickness increased. The factors that have the greatest effect on the amount of $CO_2$ emissions were coal layer thickness and coal mining. And, as the air quantity increases, it has a great effect on the decrease of carbon dioxide concentration.

• Journal of Korea Water Resources Association
• /
• v.50 no.2
• /
• pp.89-97
• /
• 2017

Turbulent flow structure in the high amplitude meandering channel is complex due to secondary recirculation with helicoidal motions and shear layers formed by flow separation from the curved sidewall. In this work, the secondary flow and the superelevation of the water surface produced in the high-amplitude Kinoshita channel are reproduced by the unsteady Reynolds-averaged Navier-Stokes (RANS) computations using the VOF technique for resolving the variation of water surface elevation and three statistical turbulence models ($k-{\varepsilon}$, RNG $k-{\varepsilon}$, $k-{\omega}$ SST). The numerical results computed by a second-order accurate finite volume method are compared with an existing experimental measurement. Among applied turbulence models, $k-{\omega}$ SST model relatively well predicts overall distribution of the secondary recirculation in the Kinoshita channel, while all three models yield similar prediction of water superelevation transverse slope. The secondary recirculation driven by the radial acceleration in the upstream bend affects the flow structure in the downstream bend, which yields a pair of counter-rotating vortices at the bend apex. This complex flow pattern is reasonably well reproduced by the $k-{\omega}$ SST model. Both $k-{\varepsilon}$ based models fail to predict the clockwise-rotating vortex between a pair of counter-rotating vortices which was observed in the experiment. Regardless of applied turbulence models, the present computations using the VOF method appear to well reproduce the superelevation of water surface through the meandering channel.

• Magazine of the Korean Society of Agricultural Engineers
• /
• v.40 no.5
• /
• pp.91-99
• /
• 1998

The widespread use of thin shell structures has created a need for a systematic method of analysis which can adequately account for arbitrary geometric form and boundary conditions as well as arbitrary general type of loading. Therefore, the stress and analysis of thin shell has been one of the more challenging areas of structural mechanics. A wide variety of numerical methods have been applied to the governing differential equations for spherical and cylindrical structures with a few results applicable to practice. The analysis of axisymmetric spherical shell is almost an every day occurrence in many industrial applications. A reliable and accurate finite element analysis procedure for such structures was needed. Dynamic loading of structures often causes excursions of stresses well into the inelastic range and the influence of geometry changes on the response is also significant in many cases. Therefore both material and geometric nonlinear effects should be considered. In general, the shell structures designed according to quasi-static analysis may fail under conditions of dynamic loading. For a more realistic prediction on the load carrying capacity of these shell, in addition to the dynamic effect, consideration should also include other factors such as nonlinearities in both material and geometry since these factors, in different manner, may also affect the magnitude of this capacity. The objective of this paper is to demonstrate the dynamic characteristics of spherical shell. For these purposes, the spherical shell subjected to uniformly distributed step load was analyzed for its large displacements elasto-viscoplastic static and dynamic response. Geometrically nonlinear behaviour is taken into account using a Total Lagrangian formulation and the material behaviour is assumed to elasto-viscoplastic model highly corresponding to the real behaviour of the material. The results for the dynamic characteristics of spherical shell in the cases under various conditions of base-radius/central height(a/H) and thickness/shell radius(t/R) were summarized as follows : The dynamic characteristics with a/H. 1) AS the a/H increases, the amplitude of displacement in creased. 2) The values of displacement dynamic magnification factor (DMF) were ranges from 2.9 to 6.3 in the crown of shell and the values of factor in the mid-point of shell were ranged from 1.8 to 2.6. 3) As the a/H increases, the values of DMF in the crown of shell is decreased rapidly but the values of DMF in mid-point shell is increased gradually. 4) The values of DMF of hoop-stresses were range from 3.6 to 6.8 in the crown of shell and the values of factor in the mid-point of shell were ranged from 2.3 to 2.6, and the values of DMF of stress were larger than that of displacement. The dynamic characteristics with t/R. 5) With the thickness of shell decreases, the amplitude of the displacement and the period increased. 6) The values of DMF of the displacement were ranged from 2.8 to 3.6 in the crown of shell and the values of factor in the mid-point of shell were ranged from 2.1 to 2.2.

• Korean Journal of Environmental Agriculture
• /
• v.10 no.1
• /
• pp.67-75
• /
• 1991

The QUAL-II E Model was applied to predict the water quality of the Mankyong drainage System, and lead to following conclusion. 1. The difference between computed and measured BOD at the M-3 (Bakgugeong) station was within 10%, indicating that the application of the QUAL-IIE Model for the prediction of water quality was satisfactory thus far. 2. The application of the model states that the discharge of concentrated pollutants at the M-1 station on the Jeonju stream, located 41Km upstream from the estuary, causes the worst problems. The sluice which extends residence time and enlarges watery surface improves water quality by a Self-purification process at the M-3 station, 28km upstream from the estuary. 3. The accuracy of the model diminished when this model was applied on the estuary downstream of the sluice. Hence, the application of the model on the estuary needs to be used with caution. 4. Among the conputed water quality parameters, BOD is the worst problem. At the M-3 station, BOD is computed to be 26.6 mg/1 in 1996, 30.7 mg/1 in 2,001, 33.0 mg/l in 2006, and 37.5 mg/1 in 2011. When preventive measures against water pollution are not properly exercised, severe problems in irrigation and water resources are expected. This study will be of used in the selection of irrigation water intake points, the criteria of effluent treatment, the management of water resources, and the establishment of water quality managemont policy.

• Journal of the Korean GEO-environmental Society
• /
• v.9 no.5
• /
• pp.45-52
• /
• 2008

The present study has extended OpenSees, which is an open-source software framework DOS program for developing applications to idealize geotechnical and structural problems, for the static analysis of axial load capacity and settlement of single piles in MS Windows environment. The Windows version of OpenSees as improved by this study has enhanced the DOS version from a general purpose software program to a special purpose program for driven and bored pile analysis with additional features of pre-processing and post-processing and a user friendly graphical interface. The method used in the load capacity analysis is the numerical methods based on load transfer functions combined with finite elements. The use of empirical nonlinear T-z and Q-z load transfer curves to model soil-pile interaction in skin friction and end bearing, respectively, has been shown to capture the nonlinear soil-pile response under settlement due to load. Validation studies have shown the static load capacity and settlement predictions implemented in this study are in fair agreement with reference data from the static loading tests.

• Journal of the Society of Naval Architects of Korea
• /
• v.30 no.1
• /
• pp.125-133
• /
• 1993

Most studies on the vibration analysis of a beam-column with ends elastically restrained and various intermediate constraints have been based on the Euler beam theory, which is inadequate for beam-columns of low slenderness ratios. In this paper, analytical methods for vibration and dynamic sensitivity of a Timoshenko beam-column with ends elastically restrained and various intermediate constraints are presented. Firstly, an exact solution method is shown. Since the exact method requires considerable computational effort, a Rayleigh-Ritz analysis is also investigated. In the latter two kinds of trial functions are examined for comparisions : eigenfunctions of the base system(the system without intermediate constraints) and polynomials having properties corresponding to the eigenfunctions of the base system. The results of some numerical Investigations show that the Rayleigh-Ritz analysis using the characteristic polynomials is competitive with the exact solutions in accuracy, and that it is much more efficient in computations than using the eigenfunctions of the base system, especially in the dynamic sensitivity analysis. In addition, the prediction of the changes of natural frequencies due to the changes of design variables based on the first order sensitivity is in good agreements with that by the ordinary reanalysis as long as the changes of design variables are moderate.

• Fire Science and Engineering
• /
• v.32 no.1
• /
• pp.16-23
• /
• 2018

In a closed-compartment with various configurations, the correlation that can predict the maximum heat release rate (HRR) with the changes in internal volume and fire growth rate was investigated numerically. The volume of the compartment was controlled by varying the length ratio based on the bottom surface shape of the ISO 9705 fire room, where the ceiling height was fixed to 2.4 m. As a main result, the effect of a change in ceiling height on the maximum HRR was examined by a comparison with a previous study that considered the change in ceiling height. In addition, a more generalized correlation equation was proposed that could predict the maximum HRR in closed-compartments regardless of the changes in ceiling height. This correlation had an average error of 7% and a maximum error of 19% for various fire growth rates when compared with the numerical results. Finally, the applicability of the proposed correlation to representative fire compartments applied to the domestic performance-based design (PBD) was examined. These results are expected to provide useful information on predicting the maximum HRR caused by flashover in closed-compartments as well as the input information required in a fire simulation.

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

This study presents the prediction methodology of debris flow occurrence areas using the SINMAP model. Former studies used a single calibration region applying some of the soil test results to predict debris flow occurrence in SINMAP model, which couldn't subdivide the soil properties for the target areas. On the other hands, a multi-calibration region using a detailed soil map and soil strength parameters (c, ${\phi}$) for each soil series to make up for limitation of former studies is proposed. In this process, soils with soil erodibility factor (K) are classified into three types: 1) gravel and gravelly soil. 2) sand and sandy soil, and 3) silt and clay. In addition, T/R estimation method using mean elevation of target area instead of T/R method using actual occurrence time is suggested in this study. The suggested method is applied to Seobyeok-1 ri area, Bonghwa-gun where debris flow occurred. As a result of comparison between two T/R estimation method, both T/R estimations are almost equal. Therefore, the suggested methodologies in this study will contribute to set up the national-wide mitigation plan against debris flow occurrence.