• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.9 s.240
• /
• pp.997-1005
• /
• 2005

The extinction characteristics of low strain rate normal gravity (1-g) nonpremixed methane-air flames were studied numerically and experimentally. A time-dependent axisymmetric two-dimensional (2D) model considering buoyancy effects and radiative heat transfer was developed to capture the structure and extinction limits of 1-g flames. One-dimensional (1D) computations were also conducted to provide information on 0-g flames. A 3-step global reaction mechanism was used in both the 1D and 2D computations to predict the measured extinction limit and flame temperature. A specific maximum heat release rate was introduced to quantify the local flame strength and to elucidate the extinction mechanism. Overall fractional contribution by each term in the energy equation to the heat release was evaluated to investigate the multi-dimensional structure and radiative extinction of 1-g flames. Images of flames were taken for comparison with the model calculation undergoing extinction. The two-dimensional numerical model was validated by comparing flame temperature profiles and extinction limits with experiments and ID computation results. The 2D computations yielded insight into the extinction mode and flame structure of 1-g flames. Two combustion regimes depending on the extinction mode were identified. Lateral heat loss effects and multi-dimensional flame structure were also found. At low strain rates of 1-g flame ('Regime A'), the flame is extinguished from the weak outer flame edge, which is attributed to multi-dimensional flame structure and flow field. At high strain rates, ('Regime B'), the flame extinction initiates near the flame centerline due to an increased diluent concentration in reaction zone, which is the same as the extinction mode of 1D flame. These two extinction modes could be clearly explained with the specific maximum heat release rate.

• The Journal of Engineering Geology
• /
• v.7 no.2
• /
• pp.91-99
• /
• 1997

In general, there are various approaches available in literature to model discontinuous rock masses and engineers are often confused which one to use for designing structures in rock masses. Modelling rock masses can be classified mainly into two approaches. One is discrete modelling of intact rock and discontinuities and the other is the equivalent continuum modelling. In this study five models are selected ;(1) Crack tensor model, (2) Equivalent volume defect model, (3) Damage model, (4) Micro - structure model (Parallel model and Series model), and (5) Homogenization model. Most of these models are mainly concerned with how to define additional strain due to discontinuities over the representative elementary volume (REV) and how to relate the stress field of discontinuities to that acting on the REV. The characteristics of these models are clarified by comparing with results of some laboratory tests.

• The Journal of Engineering Geology
• /
• v.19 no.1
• /
• pp.9-14
• /
• 2009

In this work applicability of a sharp-interface model to well-scale problems is examined. In coastal areas one of the main concerns in groundwater development is saltwater intrusion. To determine the saltwater intrusion at a well one needs to evaluate the saltwater content from a pumping well. We conducted laboratory sand-tank experiments to determine if a sharp-interface model can simulate the saltwater content of a pumping well. Lateral saltwater intrusion experiments and freshwater lens experiments are conducted. A partially penetrating well and a fully penetrating well are used. Comparison between numerical results and observation results indicates that the sharp-interface model can determine saltwater contents reasonably well. We conclude that a sharp-interface model can be applied to well-scale groundwater flow problems in coastal areas.

• Journal of the Korean Society of Industry Convergence
• /
• v.15 no.1
• /
• pp.21-27
• /
• 2012

Decomposed granite soils are in a wide range of conditions depending on the degrees of weathering. This paper is intended to examine laboratory tests such as consolidation tests and conventional triaxial compression tests conducted in order to find out the mechanical properties of Cheongju granite soil. Along with the foregoing, the results of basic physical tests conducted in order to grasp the physical properties of Cheongju granite soil were described and based on the results, methods to calculate the mechanical parameters of numerical approaches using Lade's double surface work-hardening constitutive model were examined. Finally, it is intended to explain the stress properties of Cheongju granite soil used as a geotechnical material based on its shear behavior and critical state concept using the results of isotropic consolidation tests and triaxial compression tests. As a conclusion, it can be seen that in the relationship between confining stress and maximum deviator stress, the slope is maintained at a constant value of 2.95. In the drained CTC test, maximum deviator stress generally existed in a range of axial strain of 6~8% and larger dilatancy phenomena appeared when confining stress was smaller. Finally, based on the results of the CTC tests on Cheongju granite soil, although axial strain, deviator stress and pore water pressure showed mechanical properties similar to those of overconsolidated soil, Cheongju granite soil showed behavior similar to that of normally consolidated soil in terms of volumetric strain.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.9 s.240
• /
• pp.988-996
• /
• 2005

Flame structure and extinction mechanism of counterflow methane/air non-premixed flame diluted with nitrogen are studied by NASA 2.2 s drop tower experiments and two-dimensional numerical simulations with finite rate chemistry and transport properties. Extinction mechanism at low strain rate is examined through the comparison among results of microgravity experiment, 1D and 2D simulations with a finite burner diameter. A two-dimensional simulation in counterflow flame especially with a finite burner diameter is shown to be very important in explaining the importance of multidimensional effects and lateral heat loss in flame extinction, effects that cannot be understood using a one-dimensional flamelet model. Extinction mechanism at low strain rate is quite different from that at high strain rate. Low strain rate flame is extinguished initially at the outer flame edge, the flame shrinks inward, and finally is extinguished at the center. It is clarified from the overall fractional contribution by each term in energy equation to heat release rate that the contribution of radiation fraction with 1D and 2D simulations does not change so much and the overall fractional contribution is decisively attributed to radial conduction ('lateral heat loss'). The experiments by Maruta et at. can be only completely understood if multi-dimensional heat loss effects are considered. It is, as a result, verified that the turning point, which is caused only by pure radiation heat loss, has to be shifted towards much lower global strain rate in microgravity flame.

• Fire Science and Engineering
• /
• v.13 no.4
• /
• pp.20-29
• /
• 1999

This paper describes the smoke filling process of a fire field model based on a self-developed SMEP(Smoke Movement Estimating Program) code to the simulation of fire induced flows in the two types of atrium space containing a ceiling heat flux. The SMEP using PISO algorithm solves conservation equations for mass, momentum, energy and species, together with those for the modified k- epsilon turbulence model with buoyancy term. Compressibility is assumed and the perfect gas law is used. The results of the calculated upper-layer average temperature and smoke layer interface height has shown reasonable agreement compared with the zone models. The zone models used are the CFAST developed at the Building and Fire Research Laboratory NIST U.S.A. and the NBTC one-room of FIRECALC developed at CSIRO, Australia. The smoke layer interface heights that are important in fire safety were not as sensitive as the smoke layer temperature to the nature of ceiling heat flux condition.

• Journal of Electrical Engineering and Technology
• /
• v.6 no.6
• /
• pp.853-866
• /
• 2011

We propose a multi-objective space search algorithm (MSSA) and introduce the identification of fuzzy inference systems based on the MSSA and information granulation (IG). The MSSA is a multi-objective optimization algorithm whose search method is associated with the analysis of the solution space. The multi-objective mechanism of MSSA is realized using a non-dominated sorting-based multi-objective strategy. In the identification of the fuzzy inference system, the MSSA is exploited to carry out parametric optimization of the fuzzy model and to achieve its structural optimization. The granulation of information is attained using the C-Means clustering algorithm. The overall optimization of fuzzy inference systems comes in the form of two identification mechanisms: structure identification (such as the number of input variables to be used, a specific subset of input variables, the number of membership functions, and the polynomial type) and parameter identification (viz. the apexes of membership function). The structure identification is developed by the MSSA and C-Means, whereas the parameter identification is realized via the MSSA and least squares method. The evaluation of the performance of the proposed model was conducted using three representative numerical examples such as gas furnace, NOx emission process data, and Mackey-Glass time series. The proposed model was also compared with the quality of some "conventional" fuzzy models encountered in the literature.

• Tunnel and Underground Space
• /
• v.16 no.4 s.63
• /
• pp.334-346
• /
• 2006

The discrete element code in 2-D, PFC2D, has been used as a tool to simulate various phenomena in rock mechanics and rock engineering. However, the code has an disadvantage that procedure to determine micro-parameters, namely properties of particles and contacts is repetitive and time-consuming. In this study, we analyzed the effect of micro-parameters(for generation of a contact-bonded model) on macro-properties(that were measured numerically by uniaxial compressive test). Based on the analysis, also, the time-saving and reliable method was suggested to determine a complete set of micro-parameters. In order to verify the suggested method, numerical specimens were generated in PFC2D for 10 different rock types at home and abroad. By the two trials for each specimen, in the result, the Young's modulus, Poisson's ratio and uniaxial compressive strength could be reproduced with being in relative error by about 5% to the values obtained by laboratory tests.

• Tunnel and Underground Space
• /
• v.21 no.4
• /
• pp.287-296
• /
• 2011

We performed a numerical modeling study of thermodynamic and multiphase fluid flow processes associated with underground compressed air energy storage (CAES) in a lined rock cavern (LRC). We investigated air tightness performance by calculating air leakage rate of the underground storage cavern with concrete linings at a comparatively shallow depth of 100 m. Our air-mass balance analysis showed that the key parameter to assure the long-term air tightness of such a system was the permeability of both concrete linings and surrounding rock mass. It was noted that concrete linings with a permeability of less than $1.0{\times}10^{-18}\;m^2$ would result in an acceptable air leakage rate of less than 1% with the operational pressure range between 5 and 8 MPa. We also found that air leakage could be effectively prevented and the air tightness performance of underground lined rock cavern is enhanced if the concrete lining is kept at a higher moisture content.

• Tunnel and Underground Space
• /
• v.26 no.5
• /
• pp.409-417
• /
• 2016

The shale gas is emerging as one of the oil and gas resources which can replace the traditional oil and gas resources. As the shale layer where the shale gas is deposited has low permeability, the hydrofracturing method is required to improve the productivity. This study is designed to conduct the laboratory hydrofracturing test on the samples which are modeled after the drilling hole having the general drilling hole and spiral groove. And compare the initial fracturing pressure and fluid contact between them in order to the result of the hydrofracturing depending on the shape of the drilling hole. In addition, the results were compared with the numerical modeling values from 3DEC and they were also compared with the data from the advance researches. It was found from the study that rather than the contact area of the high pressures water, the force concentration depending on the form of guide hole was more effective in the hydrofracturing.