• Journal of the Korean Professional Engineers Association
• /
• v.24 no.6
• /
• pp.97-105
• /
• 1991

The cautious blasting works had been used with emulsion explosion electric M/S delay caps. Drill depth was from 3m to 6m with Crawler Drill ø70mm on the calcalious sand stone (soft-moderate-semi hard Rock). The total numbers of fire blast were 88 round. Scale distance were induces 15.52-60.32. It was applied to propagation Law in blasting vibration as follows. Propagation Law in Blasting Vibration (Equation omitted) where V : Peak partical velocity(cm/sec) D : Distance between explosion and recording sites(m) W : Maximum Charge per delay-period of eighit milliseconds o. more(kg) K : Ground transmission constant, empirically determind on the Rocks, Explosive and drilling pattern ets. b : Charge exponents n : Reduced exponents Where the quantity D / W$^n$ is known as the Scale distance. Above equation is worked by the U.S Bureau of Mines to determine peak particle velocity. The propagation Law can be catagrorized in three graups. Cubic root Scaling charge per delay Square root Scaling of charge per delay Site-specific Scaling of charge per delay Charge and reduction exponents carried out by multiple regressional analysis. It's divided into under loom and over 100m distance because the frequency is verified by the distance from blast site. Empirical equation of cautious blasting vibration is as follows. Over 30 ‥‥‥under 100m ‥‥‥V=41(D/$^3$√W)$\^$-1.41/ ‥‥‥A Over 100 ‥‥‥‥under 100m ‥‥‥V=121(D/$^3$√W)$\^$-1.56/ ‥‥‥B K value on the above equation has to be more specified for furthur understang about the effect of explosives, Rock strength. And Drilling pattern on the vibration levels, it is necessary to carry out more tests.

• 한국연소학회:학술대회논문집
• /
• 1998.10a
• /
• pp.139-154
• /
• 1998

Under certain circumstance, premixed turbulent flame can be treated as wrinkled thin laminar flame and its motion in a hydrodynamic flow field has been investigated by employing G-equation. Past studies on G-equation successfully described certain aspects of laminar flame propagation such as effects of stretch on flame speed. In those studies, flames were regarded as a passive interface that does not influence the flow field. The experimental evidences, however, indicate that flow field can be significantly modified by the propagation of flames through the volume expansion of burned gas. In the present study, a new method to be used with G -equation is described to include the effect of volume expansion in the flame dynamics. The effect of volume expansion on the flow field is approximated by Biot-Savart law. The newly developed model is validated by comparison with existing analytical solutions of G -equation to predict flames propagating in hydrodynamic flow field without volume expansion. To further investigate the influence of volume expansion, present method was applied to initially wrinkled or planar flame propagating in an imposed velocity field and the average flame speed was evaluated from the ratio of flame surface area and projected area of unburned stream channel. It was observed that the initial wrinkling of flame cannot sustain itself without velocity disturbance and wrinkled structure decays into planar flame as the flame propagates. The rate of decay of the structure increased with volume expansion. The asymptotic change in the average burning speed occurs only with disturbed velocity field. Because volume expansion acts directly on the velocity field, the average burning speed is affected at all time when its effect is included. With relatively small temperature ratio of 3, the average flame speed increased 10%. The combined effect of volume expansion and flame stretch is also considered and the result implied that the effect of stretch is independent of volume release.

• The Journal of the Acoustical Society of Korea
• /
• v.15 no.2
• /
• pp.72-78
• /
• 1996

This thesis deals with a method to solve a two-and-one-half-dimensional ($2\frac12$ D) problem, which means that the ocean environment is two-dimensional whereas the source is fully three-dimensionally propagating, including three-dimensional refraction phenomena and three-dimensional back-scattering, using two-dimensional two-way parabolic equation method combined with Fourier synthesis. Two dimensional two-way parabolic equation method uses Galerkin's method for depth and Crank-Nicolson method and alternating direction for range and provides a solution available to range-dependent problem with wave-field back-scattered from discontinuous interface. Since wavenumber, k, is the function of depth and vertical or horizontal range, we can reduce a dimension of three-dimensional Helmholtz equation by Fourier transforming in the range direction. Thus transformed two-dimensional Helmholtz equation is solved through two-way parabolic equation method. Finally, we can have the $2\frac12$ D solution by inverse Fourier transformation of the spectral solution gained from in the last step. Numerical simulation has been carried out for a canonical ocean environment with stair-step bottom in order to test its accuracy using the present analysis. With this spectral parabolic equation method, we have examined three-dimensional acoustic propagation properties in a specified site in the Korean Straits.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.29 no.6
• /
• pp.583-590
• /
• 2016

This paper presents a dynamic crack propagation algorithm with Rayleigh damping effect based on the MLS(Moving Least Squares) Difference Method. Dynamic equilibrium equation and constitutive equation are derived by considering Rayliegh damping and governing equations are discretized by the MLS derivative approximation; the proportional damping, which has not been properly treated in the conventional strong formulations, was implemented in both the equilibrium equation and constitutive equation. Dynamic equilibrium equation including time relevant terms is integrated by the Central Difference Method and the discrete equations are simplified by lagging the velocity one step behind. A geometrical feature of crack is modeled by imposing the traction-free condition onto the nodes placed at crack surfaces and the effect of movement and addition of the nodes at every time step due to crack growth is appropriately reflected on the construction of total system. The robustness of the proposed numerical algorithm was proved by simulating single and multiple crack growth problems and the effect of proportional damping on the dynamic crack propagation analysis was effectively demonstrated.

• Journal of Ocean Engineering and Technology
• /
• v.9 no.1
• /
• pp.92-100
• /
• 1995

The effect of different anisotropy and stress ratio on fatigue crack propagation behavior was investigated under various stress ratio(R=-0.4, -0.2, 0.2, 0.2, 0.4) using pure titanium sheet used in aerospace, chemical and food industry. The rack closure behavior under constant load amplitude fatigue crack propagation test was examined. Fatigue crack propagation rate da/dN was estimated in terms of effective stress intensity factor range, $\Delta$K$_{eff}$, regardless of various stress ratio but was influenced by anisotropy. Also, it was found that the effect of anisotropy was considerably decreased but still not negligible when he da/dN was evaluated by a conventional parameter, $\Delta$$K_{eff}$/E and when the modified da/dN.$\sqrt{\varepsilon}_f$ was evaluated by $\Delta$$K_{eff}$/E. On the other hand, da/dN could be evaluated uniquely by effective new parameter, $\Delta$K$_{eff}$/$sigma_{ys}$, regardless of anisotropy, as int he following equation da/dN=C''[\frac{{\Delta}K_{eff}}{{\sigma}_{ys}}]^{n''}. And effective stress intensity factor range ratio, U was estimated by the following equation with respect to the ratio of reversed plastic zone size, $\Delta r_{p}$ to monotonic plastic zone size, $r_p$ regardless of stress ratio and anisotropy. U=-4.45$(\Delta r_{p}/r_{p})^{2}$+4.1$(\Delta r_{p}/r_{p})$+0.245_{p})$+0.245

• Nuclear Engineering and Technology
• /
• v.43 no.3
• /
• pp.271-278
• /
• 2011

When performing fire PSA in a nuclear power plant, an event mapping method, using an internal event PSA model, is widely used to reduce the resources used by fire PSA model development. Feasible initiating events and component failure events due to fire are identified to transform the fault tree (FT) for an internal event PSA into one for a fire PSA using the event mapping method. A surrogate event or damage term method is used to condition the FT of the internal PSA. The surrogate event or the damage term plays the role of flagging whether the system/component in a fire compartment is damaged or not, depending on the fire being initiated from a specified compartment. These methods usually require explicit states of all compartments to be modeled in a fire area. Fire event scenarios, when using explicit identification, such as surrogate or damage terms, have two problems: (1) there is no consideration of multiple fire propagation beyond a single propagation to an adjacent compartment, and (2) there is no consideration of simultaneous fire propagations in which an initiating fire event is propagated to multiple paths simultaneously. The present paper suggests a fire propagation equation to identify all possible fire event scenarios for an explicitly treated fire event scenario in the fire PSA. Also, a method for separating fire events was developed to make all fire events a set of mutually exclusive events, which can facilitate arithmetic summation in fire risk quantification. A simple example is given to confirm the applicability of the present method for a $2{\times}3$ rectangular fire area. Also, a feasible asymptotic approach is discussed to reduce the computational burden for fire risk quantification.

• Proceedings of the KSME Conference
• /
• 2004.11a
• /
• pp.1617-1622
• /
• 2004

Two-phase flow of liquid and gas through pipe lines are frequently encountered in nuclear power plant or industrial facility. Pressure waves which can be generated by a valve operation or any other cause in pipe lines propagate through the two-phase flow, often leading to severe noise and vibration problems or fatigue failure of pipe line system. It is of practical importance to predict the propagation characteristics of the pressure waves for the safety design for the pipe line. In the present study, a theoretical analysis is performed to understand the propagation characteristics of a weak shock wave in a bubbly flow. A wave equation is developed using a small perturbation method to analyze the weak shock wave through a bubbly flow with comparably low void fractions. It is known that the elasticity of pipe and void fraction significantly affect the propagation speed of shock wave, but the frequency of relaxation oscillation which is generated behind the shock wave is not strongly influenced by the elasticity of pipe. The present analytical results are in close agreement with existing experimental data.

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

A numerical study is carried out for the detonation wave propagation through a T-branch. The T-branch is a crucial part of the combustion wave igniter, a novel concept of rocket ignition system aimed for the simultaneous ignition of multiple combustion chambers by delivering detonation waves. Euler equation and induction parameter equation are used as governing equations with a reaction term modeled from the chemical kinetics database obtained from a detailed chemistry mechanism. Second-order accurate implicit time integration and third-order space accurate TVD algorithm were used for solution of the coupled equations. Over two-million grid points enabled the capture of the dynamics of the detonation wave propagation including the degeneration and re-initiation phenomena, and some of the design factors were be obtained for the CWI flame tubes.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.31 no.4
• /
• pp.432-440
• /
• 2007

In this study, we have investigated the closure behavior of fatigue cracks in SAPH440 steel plates under a mixed-mode I+II loading. A crack image capturing system as a direct measuring method was used to measure the closure levels at a crack tip. The crack closure levels in the fluctuation and stable sections were increased with the increase of the mode mixture ratio. The mode mixture ratio independent fatigue crack propagation rates equation was calculated by considering mixed-mode crack closure levels. The equation was examined according to the application method of crack opening ratio. The fracture surface analysis by C-scan method was also performed in order to investigate the closure mechanism and propagation mode of crack under the mixed-mode I+II loading. The crack closure under the mixed mode I+II is confirmed as a surface roughness closure by the quantitative analysis of fracture surface using the proposed surface roughness parameter.

• The Transactions of The Korean Institute of Electrical Engineers
• /
• v.58 no.10
• /
• pp.2011-2015
• /
• 2009

In an insulating dielectric liquid such as transformer oil, space charge injection and propagation were analyzed under the Fowler-Nordheim and Richardson-Dushman's thermal emission charge injection conditions for blade-plane electrodes stressed by a step voltage. The governing equations were composed of all five equations such as the Poisson's equation for electric fields, three continuity equations for electrons, negative, and positive ions, and energy balanced equation for temperature distributions. The governing equations for each carrier, the continuity equations, belong to the hyperbolic-type PDE of which the solution has a step change at the space charge front resulting in numerical instabilities. To decrease these instabilities, the governing equations were solved simultaneously by the Finite Element Method (FEM) employing the artificial diffusion scheme as a stabilization technique. Additionally, the terminal current was calculated by using the generalized energy method which is based on the Poynting's theorem, and represents more reliable and stable approach for evaluating discharge current. To verify the proposed method, the discharge phenomena were successfully applied to the blade~plane electrodes, where the radius of blade cap was $50{\mu}m$.