• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.151-154
• /
• 2002

The twin impulse wave leads to very complicated flow fields, such as Mach stem, spherical waves, and vortex ring. The twin impulse wave discharged from the exits of the two tubes placed in parallel is investigated to understand detailed flow physics associated with the twin impulse wave, compared with those in a single impulse wave. In the current study, the merging phenomena and propagation characteristics of the impulse waves are investigated using a shock tube experiment and by numerical computations. The Harten-Yee's total variation diminishing (TVD) scheme is used to solve the unsteady, two-dimensional, compressible, Euler equations. The Mach number $M_{s}$, of incident shock wave is changed below 1.5 and the distance between two-parallel tubes, L/d, is changed from 1.2 to 4.0. In the shock tube experiment, the twin impulse waves are visualized by a Schlieren optical system for the purpose of validation of computational work. The results obtained show that on the symmetric axis between two parallel tubes, the peak pressure produced by the twin-impulse waves and its location strongly depend upon the distance between two parallel tubes, L/d and the incident shock Mach number, $M_{s}$. The predicted Schlieren images represent the measured twin-impulse wave with a good accuracy.

• Coupled systems mechanics
• /
• v.12 no.1
• /
• pp.41-68
• /
• 2023

The present work is concerned with the study of the influence of inhomogeneous initial stresses in a hollow cylinder containing a compressible inviscid fluid on the propagation of axisymmetric longitudinal waves propagating in this cylinder. The study is carried out using the so-called three-dimensional linearized theory of elastic waves in bodies with initial stresses to describe the motion of the cylinder and using the linearized Euler equations to describe the flow of the compressible inviscid fluid. It is assumed that the inhomogeneous initial stresses in the cylinder are caused by the internal pressure of the fluid. To solve the corresponding eigenvalue problem, the discrete-analytic solution method is applied and the corresponding dispersion equation is obtained, which is solved numerically, after which the corresponding dispersion curves are constructed and analyzed. To obtain these dispersion curves, parameters characterizing the magnitude of the internal pressure, the ratio of the sound velocities in the cylinder material and in the fluid, and the ratio of the material densities of the fluid and the cylinder are introduced. Based on these parameters, the influence of the inhomogeneous initial stresses in the cylinder on the dispersion of the above-mentioned waves in the considered hydro-elastic system is investigated. Moreover, based on these results, appropriate conclusions about this influence are drawn. In particular, it is found that the character of the influence depends on the wavelength. Accordingly, the inhomogeneous initial stresses before (after) a certain value of the wavelength lead to a decrease (increase) of the wave propagation velocity in the zeroth and first modes.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03a
• /
• pp.367-373
• /
• 2008

Spinning detonations propagating in a circular tube were numerically investigated with a one-step irreversible reaction model governed by Arrhenius kinetics. Activation energy is used as parameter as 10, 20, 27 and 35, and the specific heat ratio and the heat release are fixed as 1.2 and 50. The time evolution of the simulation results was utilized to reveal the propagation mechanism of single-headed spinning detonation. The track angle of soot record on the tube wall was numerically reproduced with various levels of activation energy, and the simulated unique angle was the same as that of the previous reports. The maximum pressure histories of the shock front on the tube wall showed stable pitch at Ea=10, periodical unstable pitch at Ea=20 and 27 and unstable pitch consisting of stable, periodical unstable and weak modes at Ea=35, respectively. In the weak mode, there is no Mach leg on the shock front, where the pressure level is much lower than the other modes. The shock front shapes and the pressure profiles on the tube wall clarified the mechanisms of these stable and unstable modes. In the stable pitch at Ea=10, the maximum pressure history on the tube wall remained nearly constant, and the steady single Mach leg on the shock front rotated at a constant speed. The high and low frequency pressure oscillations appeared in the periodical unstable pitch at Ea=20 and 27 of the maximum pressure history. The high frequency was one cycle of a self-induced oscillation by generation and decay in complex Mach interaction due to the variation in intensity of the transverse wave behind the shock front. Eventually, sequential high frequency oscillations formed the low frequency behavior because the frequency behavior was not always the same for each cycle. In unstable pitch at Ea=35, there are stable, periodical unstable and weak modes in one cycle of the low frequency oscillation in the maximum pressure history, and the pressure amplitude of low frequency was much larger than the others. The pressure peak appeared after weak mode, and the stable, periodical unstable and weak modes were sequentially observed with pressure decay. A series of simulations of spinning detonations clarified that the unsteady mechanism behind the shock front depending on the activation energy.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.10a
• /
• pp.367-373
• /
• 2008

Spinning detonations propagating in a circular tube were numerically investigated with a one-step irreversible reaction model governed by Arrhenius kinetics. Activation energy is used as parameter as 10, 20, 27 and 35, and the specific heat ratio and the heat release are fixed as 1.2 and 50. The time evolution of the simulation results was utilized to reveal the propagation mechanism of single-headed spinning detonation. The track angle of soot record on the tube wall was numerically reproduced with various levels of activation energy, and the simulated unique angle was the same as that of the previous reports. The maximum pressure histories of the shock front on the tube wall showed stable pitch at Ea=10, periodical unstable pitch at Ea=20 and 27 and unstable pitch consisting of stable, periodical unstable and weak modes at Ea=35, respectively. In the weak mode, there is no Mach leg on the shock front, where the pressure level is much lower than the other modes. The shock front shapes and the pressure profiles on the tube wall clarified the mechanisms of these stable and unstable modes. In the stable pitch at Ea=10, the maximum pressure history on the tube wall remained nearly constant, and the steady single Mach leg on the shock front rotated at a constant speed. The high and low frequency pressure oscillations appeared in the periodical unstable pitch at Ea=20 and 27 of the maximum pressure history. The high frequency was one cycle of a self-induced oscillation by generation and decay in complex Mach interaction due to the variation in intensity of the transverse wave behind the shock front. Eventually, sequential high frequency oscillations formed the low frequency behavior because the frequency behavior was not always the same for each cycle. In unstable pitch at Ea=35, there are stable, periodical unstable and weak modes in one cycle of the low frequency oscillation in the maximum pressure history, and the pressure amplitude of low frequency was much larger than the others. The pressure peak appeared after weak mode, and the stable, periodical unstable and weak modes were sequentially observed with pressure decay. A series of simulations of spinning detonations clarified that the unsteady mechanism behind the shock front depending on the activation energy.

• Proceedings of the KSME Conference
• /
• 2003.04a
• /
• pp.1572-1577
• /
• 2003

High-power pulsed laser ablation under atmospheric pressure is studied utilizing numerical and experimental methods with emphasis on recondensation ratio, and the dynamics of the laser induced vapor flow. In the numerical calculation, the temperature pressure, density and vaporization flux on a solid substrate are first obtained by a heat-transfer computation code based on the enthalpy method, and then the plume dynamics is calculated by using a commercial CFD package. To confirm the computation results, the probe beam deflection technique was utilized for measuring the propagation of a laser induced shock wave. Discontinuities of properties and velocity over the Knudsen layer were investigated. Related with the analysis of the jump condition, the effect of the recondesation ratio on the plume dynamics was examined by comparing the pressure, density, and mass fraction of ablated aluminum vapor. To consider the effect of mass transfer between the ablation plume and air, unlike the most previous investigations, the equation of species conservation is simultaneously solved with the Euler equations. Therefore the numerical model computes not only the propagation of the shock front but also the distribution of the aluminum vapor. To our knowledge, this is the first work that employed a commercial CFD code in the calculation of pulsed ablation phenomena.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.15 no.8 s.101
• /
• pp.962-967
• /
• 2005

The twin-impulse wave discharged from two parallel tubes is investigated to see flow patterns, compared with the single impulse wave. In the present study, the merging phenomena and propagation characteristics of the impulse waves are investigated by experiment and numerical computation. The Harten-Yee's total variation diminishing scheme is used to solve the unsteady, two-dimensional, compressible, Euler equations. The Mach number Ms of incident shock wave is lower than 1.5 and the distance between the tubes is between 1.2 and 4.0. In the shock tube experiment, the twin impulse waves are visualized by a Schlieren optical system in order to validate the computational result. It is shown that on the symmetric axis between two parallel tubes, the peak pressure produced by the twin impulse waves and its location strongly depend upon the tube distance and the incident shock Mach number, Ms. The predicted Schlieren images show a good agreement with the measured twin-impulse wave.

• Transactions of the Korean Society of Machine Tool Engineers
• /
• v.10 no.1
• /
• pp.133-140
• /
• 2001

Recently, to assure the integrity of a structural components such as piping pressure vessels and thinning structure, Lamb wave inspection technique has been used in material evaluation. It is very important to select the optimal Lamb wave mode and to analyze the signal accurately because of its unique dispersion properties grnerating several modes within the speci-men. It this study, the feasibility of material evaluation applications using wavelet analysis of Lamb wave has been veir-fied experimentally. These results show as follows; 1)dispersion characteristic of each mode in dispersion curve is demon-strated that A0 mode propagating material surface is useful mode having the lest energy loss and not sensitive to surface condition. 2) it can be detected even the micro defect ($1\times2mm$) fabricated in ultrasonic probe flaw distance (290mm) to axis direction. 3) the wavelet transform which is called "time-frequency analysis" shows the Lamb wave propagation due to the change of materials characterization can be evaluated at each frequency and experimental group velocity of Lamb wave agrees quite well with that of simulated dispersion curve.ion curve.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.2
• /
• pp.61-71
• /
• 1993

The stress waves generated by the mechanical energies by impact or the chemical energies by the explosions are transmitted through medium. The wave propagation process through medium is a very complicated procedure due to the reflections and refractions of the waves at the free surfaces and interfaces. In this study the pressure independent Von-Mises model is employed for the wave propagation analysis in the layered systems. Governing equations of this study are conservation equations of momentum and mass in Lagrangian coordinate system which is fixed to the material. Due to the shock-front which violates the continuity assumptions inherent in the differential equations numerical artificial viscosity is used to spread the shock front over several computational zones. These equations are solved by Finite Difference Method with discretized time and space coordinates. The associate normality flow rule as a plastic theory is implemented to find the plastic strains.

• Journal of Korean Society of Environmental Engineers
• /
• v.27 no.10
• /
• pp.1065-1071
• /
• 2005

Many techniques of leak detection in pipeline systems have developed based on the propagation wave speeds and wave attenuation. In this paper, the transient analysis methodology is used for calculating the wave speed in the plastic pipe and a frequency analysis methodology is developed for leakage detection in water pipe networks. Data acquisition system for pressurized pipeline system were designed md fabricated to obtain high frequency pressure data. The methodology properly handles the unavoidable uncertainties in measurement and modeling error. Based on information from head pressure test data, it provides leak prediction capability from the transient events with leakage.

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

Two-dimensional hydrodynamic analysis was performed to analyze the attenuating characteristics of shock waves generated by the detonation of the bulkhead initiator. Through the interlocking analysis between HNS and HMX stacking initiator and STS bulkhead, we have precisely simulated detonation growth and pressure wave attenuation phenomena. The free surface velocity at the surface of the bulkhead was measured for quantitative comparison with the test data by VISAR. As a result, it was confirmed that the pressure attenuating pattern of the shock wave exponentially decreased according to the bulkhead thickness. The observed inflection point at the particle velocity measured over time is due to the subsequent propagation of the shock wave due to the rapid spallation of the interface between the detonator and the bulkhead.