• The Magazine of the Society of Air-Conditioning and Refrigerating Engineers of Korea
• /
• v.13 no.4
• /
• pp.215-222
• /
• 1984

An analytical study on the thermal instability of fluid in a vertical solt between two permeable walls has been carried out using fast converging power series solution method. For given values of prandtl number Pr and permeability paramter ${\sigma}$, the critical Grashof number $Gr_c$ and the critical wave number ac are found as eigenvalues of the problem formulated by the stability equations and the appropriate boundary conditions which are derived on the basis of linear stability theory. In the case of ${\sigma}\;>\;10^4$, the results approach those of solid boundary case, but in the case of ${\sigma}\;<\;10^3$, the decrease of $Gr_c$ and $a_c$become more prominent. In other words, the permeable walls cause the flow to be more unstable than the solid walls. This is considered to be due to the slip of the fluid on the wail, which decrease the friction force.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2003.05a
• /
• pp.91-93
• /
• 2003

A comprehensive numerical study is carried out to investigate for the understanding of the flow evolution and flame development in a supersonic combustor with normal injection of ncumally injecting hydrogen in airsupersonic flows. The formulation treats the complete conservation equations of mass, momentum, energy, and species concentration for a multi-component chemically reacting system. For the numerical simulation of supersonic combustion, multi-species Navier-Stokes equations and detailed chemistry of H2-Air is considered. It also accommodates a finite-rate chemical kinetics mechanism of hydrogen-air combustion GRI-Mech. 2.11[1], which consists of nine species and twenty-five reaction steps. Turbulence closure is achieved by means of a k-two-equation model (2). The governing equations are spatially discretized using a finite-volume approach, and temporally integrated by means of a second-order accurate implicit scheme (3-5).The supersonic combustor consists of a flat channel of 10 cm height and a fuel-injection slit of 0.1 cm width located at 10 cm downstream of the inlet. A cavity of 5 cm height and 20 cm width is installed at 15 cm downstream of the injection slit. A total of 936160 grids are used for the main-combustor flow passage, and 159161 grids for the cavity. The grids are clustered in the flow direction near the fuel injector and cavity, as well as in the vertical direction near the bottom wall. The no-slip and adiabatic conditions are assumed throughout the entire wall boundary. As a specific example, the inflow Mach number is assumed to be 3, and the temperature and pressure are 600 K and 0.1 MPa, respectively. Gaseous hydrogen at a temperature of 151.5 K is injected normal to the wall from a choked injector.A series of calculations were carried out by varying the fuel injection pressure from 0.5 to 1.5MPa. This amounts to changing the fuel mass flow rate or the overall equivalence ratio for different operating regimes. Figure 1 shows the instantaneous temperature fields in the supersonic combustor at four different conditions. The dark blue region represents the hot burned gases. At the fuel injection pressure of 0.5 MPa, the flame is stably anchored, but the flow field exhibits a high-amplitude oscillation. At the fuel injection pressure of 1.0 MPa, the Mach reflection occurs ahead of the injector. The interaction between the incoming air and the injection flow becomes much more complex, and the fuel/air mixing is strongly enhanced. The Mach reflection oscillates and results in a strong fluctuation in the combustor wall pressure. At the fuel injection pressure of 1.5MPa, the flow inside the combustor becomes nearly choked and the Mach reflection is displaced forward. The leading shock wave moves slowly toward the inlet, and eventually causes the combustor-upstart due to the thermal choking. The cavity appears to play a secondary role in driving the flow unsteadiness, in spite of its influence on the fuel/air mixing and flame evolution. Further investigation is necessary on this issue. The present study features detailed resolution of the flow and flame dynamics in the combustor, which was not typically available in most of the previous works. In particular, the oscillatory flow characteristics are captured at a scale sufficient to identify the underlying physical mechanisms. Much of the flow unsteadiness is not related to the cavity, but rather to the intrinsic unsteadiness in the flowfield, as also shown experimentally by Ben-Yakar et al. [6], The interactions between the unsteady flow and flame evolution may cause a large excursion of flow oscillation. The work appears to be the first of its kind in the numerical study of combustion oscillations in a supersonic combustor, although a similar phenomenon was previously reported experimentally. A more comprehensive discussion will be given in the final paper presented at the colloquium.

• The Journal of the Petrological Society of Korea
• /
• v.25 no.3
• /
• pp.211-230
• /
• 2016

Detailed mapping along the Keumwang fault reveals a complex history of multiple brittle reactivations following late Jurassic and early Cretaceous ductile shearing. The fault core consists of a 10~50 m thick fault gouge layer bounded by a 30~100 m thick damaged zone. The Pre-cambrian gneiss and Jurassic granite underwent at least six distinct stages of fault movements based on deformation environment, time and mechanism. Each stage characterized by fault kinematics and dynamics at different deformation environment. Stage 1 generated mylonite series along the Keumwang shear zone by sinistral ductile shearing during late Jurassic and early Cretaceous. Stage 2 was a mostly brittle event generating cataclasite series superimposed on the mylonite series of the Keumwang shear zone. The roundness of pophyroclastes and the amount of matrix increase from host rocks to ultracataclasite indicating stronger cataclastic flow toward the fault core. At stage 3, fault gouge layer superimposed on the cataclasite generated during stage 2 and the sedimentary basins (Umsung and Pungam) formed along the fault by sinistral strike-slip movement. Fragments of older cataclasite suspended in the fault gouge suggest extensive reworking of fault rocks at brittle deformation environments. At stage 4, systematic en-echelon folds, joints and faults were formed in the sedimentary basins by sinistral strike-slip reactivation of the Keumwang fault. Most of the shearing is accommodated by slip along foliations and on discrete shear surfaces, while shear deformation tends to be relatively uniformly distributed within the fault damage zone developed in the mudrocks in the sedimentary basins. Fine-grained andesitic rocks intruded during stage 4. Stage 5 dextral strike-slip activity produced shear planes and bands in the andesitic rocks. ESR(Electron Spin Resonance) dates of fault gouge show temporal clustering within active period and migrating along the strike of the Keumwang fault during the stage 6 at the Quaternary period.

• Journal of Advanced Marine Engineering and Technology
• /
• v.37 no.5
• /
• pp.500-509
• /
• 2013

This study reports a numerical analysis of the internal flow characteristics of the integrated urea-SCR muffler system with the various geometries of the multi-perforated tube which is set up between the muffler inlet and in front of SCR catalysts. The multi-perforated tube is generally used to disperse uniformly the urea-water solution spray and to make better use of the SCR catalyst, resulting in the increased $NO_x$ reduction and decreased ammonia slip. The effects of the multi-perforated tube orifice area ratios on the velocity distributions in front of the SCR catalyst, which is ultimately quantified as the uniformity index, were investigated for the optimal muffler system design. The steady flow model was applied by using a general-purpose commercial software package. The air at the room temperature was used as a working fluid, instead of the exhaust gas and urea-water solution spray mixture. From the analysis results, it was clarified that the multi-perforated tube geometry sensitively affected to the formation of the bulk swirling motion inside the plenum chamber set in front of the SCR catalyst and to the uniformity index of the velocity distribution produced at the inlet of the catalyst.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2012.05a
• /
• pp.34-34
• /
• 2012

Heavy storms rainfall has caused many landslides and slope failures especially in the mountainous area of the world. Landslides and slope failures are common geologic hazards and posed serious threats and globally cause billions in monetary losses and thousands of casualies each year so that studies on slope stability and its failure mechanism under rainfall are being increasing attention of these days. Rainfall-induced slope failures are generally caused by the rise in ground water level, and increase in pore water pressures and seepage forces during periods of intense rainfall. The effective stress in the soil will be decreased due to the increased pore pressure, which thus reduces the soil shear strength, eventually resulting in slope failure. During the rainfall, a wetting front goes downward into the slope, resulting in a gradual increase of the water content and a decrease of the negative pore-water pressure. This negative pore-water pressure is referred to as matric suction when referenced to the pore air pressure that contributes to the stability of unsaturated soil slopes. Therefore, the importance is the study of saturated unsaturated soil behaviors in evaluation of slope stability under heavy rainfall condition. In an actual field, a series of failures may occur in a slope due to a rainfall event. So, this study attempts to develop a numerical model to investigate this failure mechanism. A two-dimensional seepage flow model coupled with a one-dimensional surface flow and erosion/deposition model is used for seepage analysis. It is necessary to identify either there is surface runoff produced or not in a soil slope during a rainfall event, while analyzing the seepage and stability of such slopes. Runoff produced by rainfall may result erosion/deposition process on the surface of the slope. The depth of runoff has vital role in the seepage process within the soil domain so that surface flow and erosion/deposition model computes the surface water head of the runoff produced by the rainfall, and erosion/deposition on the surface of the model slope. Pore water pressure and moisture content data obtained by the seepage flow model are then used to analyze the stability of the slope. Spencer method of slope stability analysis is incorporated into dynamic programming to locate the critical slip surface of a general slope.

• The Journal of Engineering Geology
• /
• v.32 no.4
• /
• pp.643-659
• /
• 2022

Geologic sequestration technologies such as CCS (carbon capture and storage), EGS (enhanced geothermal systems), and EOR (enhanced oil recovery) have been widely implemented in recent years, prompting evaluation of the mechanical stability of storage sites. As fluid injection can stimulate mechanical instability in storage layers by perturbing the stress state and pore pressure, poroelastic models considering various injection scenarios are required. In this study, we calculate the pore pressure, stress distribution, and vertical displacement along a surface using commercial finite element software (COMSOL); fault slips are subsequently simulated using PyLith, an open-source finite element software. The displacement fields, are obtained from PyLith is transferred back to COMSOL to determine changes in coseismic stresses and surface displacements. Our sequential use of COMSOL-PyLith-COMSOL for poroelastic modeling of fluid-injection and induced-earthquakes reveals large variations of pore pressure, vertical displacement, and Coulomb failure stress change during injection periods. On the other hand, the residual stress diffuses into the remote field after injection stops. This flow pattern suggests the necessity of numerical modeling and long-term monitoring, even after injection has stopped. We found that the time at which the Coulomb failure stress reaches the critical point greatly varies with the hydraulic and poroelastic properties (e.g., permeability and Biot-Willis coefficient) of the fault and injection layer. We suggest that an understanding of the detailed physical properties of the surrounding layer is important in selecting the injection site. Our numerical results showing the surface displacement and deviatoric stress distribution with different amounts of fault slip highlight the need to test more variable fault slip scenarios.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.18 no.12
• /
• pp.679-684
• /
• 2017

Friction Stir Welding is a metal welding technique, in which friction heat between a welding tool and a welding material is used to weld parts at temperatures below the melting point of a material. In this study, the temperature and velocity changes in a magnesium alloy (AZ31) during the welding process were analyzed by computational flow dynamics technique while welding the material using a friction stir welding technique. For the analysis, the modeling and analysis were carried out using Fluent as a fluid analysis tool. First, the welding material was assumed to be a temperature-dependent Newtonian fluid with high viscosity, and the rotation region and the stationary region were simulated separately to consider the rotational flow generated by the rotation of the welding tool having a helical groove. The interface between the welding tool and welding material was given the friction and slip boundary conditions and the heat transfer effect to the welding tool was considered. Overall, the velocity and temperature characteristics of the welded material according to time can be understood from the results of transient analysis through the above flow analysis modeling.

• Tunnel and Underground Space
• /
• v.31 no.5
• /
• pp.309-332
• /
• 2021

Natural barrier systems surrounding the geological repository for the high-level radioactive waste should guarantee the hydraulic performance for preventing or delaying the leakage of radionuclide. In the case of the behavior of a crystalline rock, the hydraulic performance tends to be decided by the existence of discontinuities, so the coupled hydro-mechanical(HM) processes on the discontinuities should be characterized. The discontinuum modelling can describe the complicated behavior of discontinuities including creation, propagation, deformation and slip, so it is appropriate to model the behavior of a crystalline rock. This paper investigated the coupled HM processes in discontinuum modelling such as UDEC, 3DEC, PFC, DDA, FRACOD and TOUGH-UDEC. Block-based discontinuum methods tend to describe the HM processes based on the fluid flow through the discontinuities, and some methods are combined with another numerical tool specialized in hydraulic analysis. Particle-based discontinuum modelling describes the overall HM processes based on the fluid flow among the particles. The discontinuum methods that are currently available have limitations: exclusive simulations for two-dimension, low hydraulic simulation efficiency, fracture-dominated fluid flow and simplified hydraulic analysis, so it could be improper to the modelling the geological repository. Based on the concepts of various discontinuum modelling compiled in this paper, the advanced numerical tools for describing the accurate coupled HM processes of the deep geological repository should be developed.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.40 no.1
• /
• pp.105-110
• /
• 2020

This study provides a method and data for predicting the flight initiation wind speed of wind-borne debris. From the force equilibrium acting on debris including aerodynamic and inertia forces, the equation for predicting the flight initiation wind speeds are presented. Wind tunnel tests were carried out to provide necessary aerodynamic data in the equation for the debris with various aspect ratios. The proposed equation for flight initiation wind speeds was validated from free flying tests in the wind tunnel. The flights of debris were mostly initiated by slip when width to thickness was less than 10, otherwise overturning were dominant. The actual flight initiation speeds were lower than that of the computed ones. The surface boundary layer flow and the gap between the debris and surface might affect the prediction error.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.15 no.6
• /
• pp.1886-1896
• /
• 1991

This paper presents a method to predict the flying height of the head slider in a hard disk drive. Quantitative predicts of the flying height according to the variations of the external load and the disk velocity have been done by numerical computation. In addition, the magnitude of the external load to keep flying height constant were also suggested. The Modified Reynolds' equation driven from hydrodynamic lubrication theory under slip flow condition was used to describe air-bearing system under the slider. To solve the equation, a Finite Volume Method (FVM) has been applied. To determine the final minimum flying height and pitch angle of the head slider, the Secant iteration method is used which update initial guess of the minimum flying height and pitch angle of the slider. In this study, the model head slider has been selected from a real hard disk drive which is equipped in many commercial personal computers. As a result, as the disk velocity increases at constant external load, the minimum flying height and the pitch angle increase due to the in crease of the air-bearing force at the bottom of slider.