• Journal of Mechanical Science and Technology
• /
• v.17 no.12
• /
• pp.2066-2077
• /
• 2003

The time-dependent behavior of unsteady condensation of moist air through the Ludwieg tube is investigated by using a computational fluid dynamics (CFD) work. The two-dimensional, compressible, Navier-Stokes equations, fully coupled with the condensate droplet growth equations, are numerically solved by a third-order MUSCL type TVD finite-difference scheme, with a second-order fractional time step. Baldwin-Lomax turbulence model is employed to close the governing equations. The predicted results are compared with the previous experiments using the Ludwieg tube with a diaphragm downstream. The present computations represent the experimental flows well. The time-dependent unsteady condensation characteristics are discussed based upon the present predicted results. The results obtained clearly show that for an initial relative humidity below 30% there is no periodic oscillation of the condensation shock wave, but for an initial relative humidity over 40% the periodic excursions of the condensation shock occurs in the Ludwieg tube, and the frequency increases with the initial relative humidity. It is also found that total pressure loss due to unsteady condensation in the Ludwieg tube should not be ignored even for a very low initial relative humidity and it results from the periodic excursions of the condensation shock wave.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.12 no.1
• /
• pp.48-56
• /
• 2002

A computational work of the impulse wave which is discharged from the open end of a pipe is compared to the Lighthill\`s aeroacoustics theory. The second-order total variation diminishing(TVD) scheme is employed to solve the axisymmetric, compressible, unsteady Euler equations. The relationship between the initial compressure wave form and the resulting impulse wave is characterized in terms of the peak pressure. The overpressure, pressure gradient and wavelength of the initial compression wave are changed to investigate the influence of the initial compressure wave form on the peak pressure of impulse wave. The results obtained show that for the initial compression wave of a large wavelength and small pressure gradient the peak pressure of the impulse wave depends upon the wavelength and pressure gradient of compression wave, but for the initial compression wave of a short wavelength and large pressure gradient the peak pressure of the impulse wave is almost constant regardless of the wavelength and pressure gradient of compression wave. The peak pressure of the impulse wave is increased with an increase in the overpressure of the initial compression wave. The results from the numerical ana1ysis are well compared to the results from the aeroacoutics theory with a food agreement.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03b
• /
• pp.452-455
• /
• 2008

A new passive control technique of cavity-induced pressure oscillations has been investigated numerically for a supersonic two-dimensional flow over open rectangular cavities at Mach number 1.83 just upstream of a cavity, in which a sub-cavity system is installed on the backward-facing step of the main cavity. A third-order TVD (Total Variation Diminishing) finite difference scheme with MUSCL is used to discretize the spatial derivatives in the unsteady compressible Navier-Stokes equations. The results obtained show that the present sub-cavity system is very effective in reducing cavity-induced pressure oscillations. The results also showed that the resultant amount of attenuation of cavity-induced pressure oscillations was dependent on the length and thickness of the flat plate, and also on the depth of the sub-cavity used as an oscillation suppressor.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.26 no.7
• /
• pp.951-958
• /
• 2002

Rapid expansion of a moist air or a stream through a supersonic nozzle often leads to non-equilibrium condensation shock wave, causing a considerable energy loss in flow field. Depending on amount of latent heat released due to non-equilibrium condensation, the flow is highly unstable or a periodical oscillation accompanying the condensation shock wave in the nozzle. The unsteadiness of the condensation shock wave is always associated with several kinds of instabilities as well as noise and vibration of flow devices. In the current study, a passive control technique using a porous wall with a plenum cavity underneath is applied for the purpose of alleviation of the condensation shock oscillations in a transonic nozzle. A droplet growth equation is coupled with two-dimensional Navier-Stokes equation system. Computations are carried out using a third-order MUSCL type TVD finite-difference scheme with a second-order fractional time step. An experiment using an indraft wind tunnel is made to validate the present computational results. The results show that the oscillations of the condensation shock wave are completely suppressed by the current passive control method.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.26 no.7
• /
• pp.997-1004
• /
• 2002

The current study describes experimental and computational works on the passive control of the steady and unsteady condensation shock waves, which are generated in a transonic nozzle. The bleed slots are installed on the contoured wall of the transonic nozzle in order 10 control the magnitude of the condensation shock wave and its oscillations. For computations, a droplet growth equation is copuled with two-dimensional Navier-Stokes equation systems. Computations are carried out using a third-order MUSCL type TVD finite-difference scheme with a second-order fractional time step. Baldwin-Lomax turbulence model is employed to close the governing equations. An experiment using an indrafi transonic wind tunnel is made to validate the computational results. The current computations represented well the experimental flows. From both the experimental and computational results it is found that the magnitude of the condensation shock wave in the bleed slotted nozzle is signi ficantly reduced, compared with no passive control of solid wall. The oscillations of the condensation shock wave are successfully suppressed by a bleed slot system.

• Journal of Mechanical Science and Technology
• /
• v.20 no.11
• /
• pp.1961-1971
• /
• 2006

A soft recovery system (SRS) is a device that stops a high speed projectile without damaging the projectile. The SRS is necessary to verify the shock resistant requirements of microelectronics and electro-optic sensors in smart munitions, where the projectiles experience over 20,000 g acceleration inside the barrel. In this study, a computer code for the performance evaluation of a SRS based on ballistic compression decelerator concept has been developed. It consists of a time accurate compressible one-dimensional Euler code with use of deforming grid and a projectile motion analysis code. The Euler code employs Roe's approximate Riemann solver with a total variation diminishing (TVD) method. A fully implicit dual time stepping method is used to advance the solution in time. In addition, the geometric conservation law (GCL) is applied to predict the solutions accurately on the deforming mesh. The equation of motion for the projectile is solved with the four-stage Runge-Kutta time integration method. A small scale SRS to catch a 20 mm bullet fired at 500 m/s within 1,600 g-limit has been designed with the proposed method.

• Journal of Mechanical Science and Technology
• /
• v.18 no.3
• /
• pp.407-418
• /
• 2004

Prandtl-Meyer expansion flow with homogeneous condensation is investigated experimentally and by numerical computations. The steady and unsteady periodic behaviors of the diabatic shock wave due to the latent heat released by condensation are considered with a view of technical application to the condensing flow through steam turbine blade passages. A passive control method using a porous wall and cavity underneath is applied to control the diabatic shock wave. Two-dimensional, compressible Navier-Stokes with the nucleation rate equation are numerically solved using a third-order TVD (Total Variation Diminishing) finite difference scheme. The computational results reproduce the measured static pressure distributions in passive and no passive Prandtl-Meyer expansion flows with condensation. From both the experimental and computational results, it is found that the magnitude of steady diabatic shock wave can be considerably reduced by the present passive control method. For no passive control, it is found that the diabatic shock wave due to the heat released by condensation oscillates periodically with a frequency of 2.40㎑. This unsteady periodic motion of the diabatic shock wave can be completely suppressed using the present passive control method.

• Journal of The Korean Astronomical Society
• /
• v.34 no.4
• /
• pp.285-287
• /
• 2001

We have studied the nonlinear evolution of a magnetized disk of isothermal gas, which is sustained by its self-gravity. Our objective is to investigate how the Jeans, Parker, and convective instabilities compete with each other in structuring/de-structuring large scale condensations in such disk. The Poisson equation for the self-gravity has been solved with a fourth-order accurate Fourier method along with the Green function, and the MHD part has been handled by an isothermal TVD code. When large wavelength perturbations are applied, the combined action of the Jeans and Parker instabilities suppresses the development of the convection and forms a dense core of prolate shape in the mid-plane. Peripheral structures around it are filamentary. The low density filaments connect the dense core to the diffuse upper region. On the other hand, when small wavelength perturbations are applied, the disk develops into an equilibrium state which is reminiscent of the Mouschovias's 2-D non-linear equilibrium of the classical Parker instability under an externally given gravity.

• Proceedings of the KSME Conference
• /
• 2004.04a
• /
• pp.1643-1648
• /
• 2004

When a shock wave propagates into a Helmholtz resonator, very complicated wave phenomena are formed both inside and outside the resonator tube. Shock wave reflection, shock focusing phenomena and shock-vortex interactions cause strong pressure fluctuations inside the resonator, consequently leading to powerful sound emission. In the present study, the wave phenomena inside and outside the Helmholtz resonator are, in detail, investigated with a help of CFD. The Mach number of the incident shock wave is varied below 2.0 and several types of resonators are tested to investigate the influence of resonator geometry on the wave phenomena. A TVD scheme is employed to solve the axisymmetric, compressible, Euler equations. The results obtained show that the configuration of the Helmholtz resonator significantly affects the peak pressure of shock wave focusing, its location, the amplitude of the discharged wave and resonance frequency.

• Proceedings of the KSME Conference
• /
• 2004.04a
• /
• pp.1595-1600
• /
• 2004

When a shock wave arrives at a duct, an impulsive wave is discharged from the duct exit and causes serious noise and vibration problems. In the current study, the characteristics of the impulsive wave discharged from a partial closed duct exit is numerically investigated using a CFD method. The Yee-Roe- Davis's total variation diminishing(TVD) scheme is used to solve the axisymmetric, unsteady, compressible Euler equations. With several partial closed duct exits, the Mach number of the incident shock wave $M_s$ and the distance L/D between the duct exit and a flat plate are varied in the range of $M_s$ = 1.01 ${\sim}$ 1.50 and L/D = 1.0 ${\sim}$ 4.0, respectively. The results obtained show that the magnitude of the impulsive wave impinging upon the flat plate strongly depends upon $M_s$, L/D and the partial closure of duct exit. The impulsive wave on the flat plate can be considerably alleviated by the partial closure of duct exit and, thus, the present method can be a passive control for the impulsive wave.