• Journal of Korean Tunnelling and Underground Space Association
• /
• v.19 no.6
• /
• pp.959-972
• /
• 2017

In evaluation of blast-induced vibration, peak particle velocity (PPV) is generally calculated by using attenuation relation curve. Calculated velocity is compared with the value in legal requirements or the standards to determine the stability. Attenuation relation curve varies depending on frequency of test blasting, geological structure of the site and blasting condition, so it is difficult to predict accurately using such an equation. Since PPV is response value from the ground, direct evaluation of the structure is impractical. Because of such a limit, engineers tend to use the commercial numerical analysis program in evaluating the stability of the structure more accurately. However, when simulate the explosion process using existing numerical analysis program, it's never easy to accurately simulate the complex conditions (fracture, crushing, cracks and plastic deformation) around blasting hole. For simulating such a process, the range for modelling will be limited due to the maximum node count and it requires extended calculation time as well. Thus, this study is intended to simulate the elastic energy after fractured zone only, instead of simulating the complex conditions of the rock that results from the blast, and the analysis of response characteristics of the velocity depending on shape and size of the fractured zone was conducted. As a result, difference in velocity and attenuation character were calculated depending on fractured zone around the blast source appeared. Propagation of vibration tended to spread spherically as it is distanced farther from the blast source.

• Tunnel and Underground Space
• /
• v.21 no.2
• /
• pp.128-137
• /
• 2011

The shut-in pressure calculated in common hydrofracturing test for vertical borehole equals generally to the minimum horizontal principal stress, so it should be considered as an essential parameter for determining the in-situ stress regime around the rock mass. It shows usually an ambiguous value in pressure-time history curves, however, because of the relationship between the behavior of hydraulic fractures and the condition of remote stress regime. In this study, a series of numerical analyses have been carried out to compare several methods for determining the shut-in pressure during hydrofracturing. The hydraulic-mechanical coupling has been applied to numerical analysis for simulating the fracture propagation by hydraulic pressure, and the different discontinuity geometry has been considered in numerical models to examine the effect of numerical element shape on fracture propagation pattern. From the numerical simulations with the four different discontinuity geometries, it was revealed that the shut-in pressure obtained from graphical methods rather than statistical method was relatively small. Consequently a care should be taken in selecting a method for determining the shut-in pressure when a stress anomaly around borehole and a fracture propagation with complicate mechanism are considered.

• Journal of the Korean Society of Visualization
• /
• v.15 no.3
• /
• pp.27-33
• /
• 2017

Direct numerical simulations of turbulent channel flows with moving wall conditions on the top wall are performed to examine the effects of the moving wall on the turbulent characteristics. The moving wall velocity only applied to the top wall with the opposite direction to the main flow is systematically varied to reveal the sustained-mechanism for turbulence. The turbulence statistics for the Couette-Poiseuille flow, such as mean velocity, root mean square of the velocity fluctuations, Reynolds shear stress and pre-multiplied energy spectra of the velocity fluctuations, are compared with those of canonical turbulent channel flows. The comparison suggests that although the turbulent activity on the top wall increases with increasing the Reynolds number, that on the bottom wall decreases, contrary to the previous finding for the canonical turbulent channel flows. The increase of the turbulent energy on the top wall is attributed to not only the increase of the Reynolds number but also elongation of the logarithmic layer due to increase of the wall layer on the top wall. However, because the logarithmic layer is shortened on the bottom wall due to the decrease of the wall layer, the turbulence energy on the bottom wall decreases despite of the increase of the Reynolds number.

• Journal of Energy Engineering
• /
• v.8 no.2
• /
• pp.333-340
• /
• 1999

A numerical study of turbulent condensing vapor jet submerged in subcooled liquids has been conducted. A physical model of the process is presented employing the locally homogeneous flow approximation of two phase flow in conjunction with a $\kappa$-$\varepsilon$-g model of turbulence properties. In this model the turbulence is represented by differential equations for its kinetic energy and dissipation. A differential equation for the concentration fluctuations is solved and a clipped normal probability distribution function is proposed for the mixture fraction. Effects of steam mass flux, pool temperature and nozzle internal diameter on the condensing vapor jet are also analyzed. The model is evaluated using existing data for turbulent condensing vapor jets. The agreement between the predictions and the available experimental data is good.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.39 no.9
• /
• pp.805-815
• /
• 2011

Nonlinear Parabolized Stability Equations(NSPE) can be effectively used to study more throughly the transition process. NPSE can efficiently analyze the stability of a nonlinear region in transition process with low computational cost compared to Direct Numerical Simulation(DNS). In this study, NPSE in general coordinate system is formulated and a computer code to solve numerically the equations is developed. Benchmark problems for incompressible and compressible boundary layers over a flat plate are analyzed to validate the present code. It is confirmed that the NPSE methodology constructed in this study is an efficient and effective tool for nonlinear stability analysis.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2017.05a
• /
• pp.451-455
• /
• 2017

The reactive Eulerian-Lagrangian code is utilized to simulate combustion of ZPP/THPP in a closed vessel. In the paper, ignition delay of THPP is mainly studied since ZPP and THPP are isolated by a boron nitride wall. Only a numerical case study is conducted as experimental observation is inaccessible. Results showed THPP ignition delay affects initial shock strength thus not only the first peak become weak, but also the frequency of a pressure oscillation is slowed.

• Proceedings of the KSME Conference
• /
• 2007.05b
• /
• pp.2378-2383
• /
• 2007

Direct Numerical Simulation was carried out to predict mass transfer in turbulent flow around a rotating stepped cylinder. This investigation is a follow-up study of Nesic et al. [Corrosion, Vol. 56, No. 10, pp. 1005 - 1014] The original motivation of this work stemmed from the efforts to design a simple device which can generate flows of high turbulence intensity at low cost for corrosion researchers. Two cases were considered; Sc=1 and 10 both at Re=335. Here, Sc and Re stand for Schmidt number and Reynolds number, respectively, based on the step height and the surface speed of the cylinder upstream the step. Main focus was placed on the correlation between turbulent fluctuation and concentration field. The spatio-temporal evolution of concentration field is discussed. The numerical results are qualitatively compared with those of the experiment conducted with the same flow configuration.

• 한국전산유체공학회:학술대회논문집
• /
• 2004.03a
• /
• pp.166-171
• /
• 2004

Direct numerical simulations are peformed to investigate the physics of a spatially developing turbulent boundary layer flow suddenly subjected to spanwise oscillating electro-magnetic forces in the near-wall region. The Reynolds number based on the inlet momentum thickness and free-stream velocity is $Re_\theta=300$. A fully-implicit fractional step method is employed to simulate the flow. The mean flow properties and the Reynolds stresses are obtained to analyze the near-wall turbulent structure. It is found that skin-friction and turbulent kinetic energy can be reduced by the electro-magnetic forces. Instantaneous flow visualization techniques are used to observe the response of streamwise vortices to spanwise oscillating forces. The near-wall vortical structures are clearly affected by spanwise oscillating electro-magnetic forces.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.10
• /
• pp.1317-1325
• /
• 2000

A new efficient numerical method for computing three-dimensional, unsteady, incompressible flows is presented. To eliminate the restriction of CFL condition, a fully-implicit time advancement in which the Crank-Nicolson method is used for both the diffusion and convection terms, is adopted. Based on an approximate block LU decomposition method, the velocity -pressure decoupling is achieved. The additional decoupling of the intermediate velocity components in the convection term is made for the fully -implicit time advancement scheme. Since the iterative procedures for the momentum equations are not required, the velocity components decouplings bring forth the reduction of computational cost. The second-order accuracy in time of the present numerical algorithm is ascertained by computing decaying vortices. The present decoupling method is applied to minimal channel flow unit with DNS (Direct Numerical Simulation).

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.31 no.9
• /
• pp.799-806
• /
• 2007

Direct Numerical Simulation was carried out to predict mass transfer in turbulent flow around a rotating stepped cylinder. This investigation is a follow-up study of DNS of turbulent flow in Nesic et al. [Corrosion, Vol. 56, No. 10, pp. 1005 - 1014] The original motivation of this work stemmed from the efforts to design a simple device which can generate flows of high turbulence intensity at low cost for corrosion researchers. Two cases were considered; Sc=1 and 10 both at Re=335. Here, Sc and Re stand for Schmidt number and Reynolds number, respectively, based on the step height and the surface speed of the cylinder upstream of the step. Main focus was placed on the correlation between turbulence and mass transfer. The spatio-temporal evolution of concentration field is discussed. The numerical results are qualitatively compared with those of the experiment conducted with a similar flow configuration.