• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.7
• /
• pp.2315-2325
• /
• 1996

A compression wave is attenuated or distorted as it propagates in a tube. The present study investigated the propagation characteristics of the compression waves which are generated by a train in a high-speed railway tunnel. A Total Variation Diminishing (TVD) difference scheme was applied to one-dimensional, unsteady viscous compressible flow. The numerical calculation involved the effects of wall friction, heat transfer and energy loss due to the friction heat in the boundary layer behind the propagating compression wave, and compared with the measurement results of a shock tube and a real tunnel. The present results show that attenuation of the compression wave in turbulent boundary layer is stronger than in laminar boundary layer, but nonlinear effect of the compression wave is greater in the laminar boundary layer. The energy loss due to the frictional heat had not influence on attenuation and distortion of the propagating compression waves.

• Journal of computational fluids engineering
• /
• v.8 no.4
• /
• pp.16-25
• /
• 2003

The design method for transonic turbine blades has been developed based on Wavier-Stokes equations. The present computing process is done on the four separate steps, i.e., determination of the blade profile, generation of the computational grids, cascade flow simulation and analysis of the computed results in the sense of the aerodynamic performance. The blade shapes are designed using the cubic polynomials under the control of the design parameters. Numerical methods for the flow equations are based on Van-Leer's FVS with an upwind TVD scheme on the finite volume. In the present study, numerical simulation has been done to investigate the effects of the design parameters on the aerodynamic peformance of the axial-flow turbine blades. Applications are made to the VKI transonic rotor blades. Computed results are analyzed with respect to four parameters and compared with the experimental data.

• Journal of computational fluids engineering
• /
• v.3 no.1
• /
• pp.11-21
• /
• 1998

A vorticity-based method for the numerical solution of the two-dimensional incompressible Navier-Stokes equations is presented. The governing equations for vorticity, velocity and pressure variables are expressed in an integro-differential form. The global coupling between the vorticity and the pressure boundary conditions is fully considered in an iterative procedure when numerical schemes are employed. The finite volume method of the second order TVD scheme is implemented to integrate the vorticity transport equation with the dynamic vorticity boundary condition. The velocity field is obtained by using the Biot-Savart integral. The Green's scalar identity is used to solve the total pressure in an integral approach similar to the surface panel methods which have been well established for potential flow analysis. The present formulation is validated by comparison with data from the literature for the two-dimensional cavity flow driven by shear in a square cavity. We take two types of the cavity now: (ⅰ) driven by non-uniform shear on top lid and body forces for which the exact solution exists, and (ⅱ) driven only by uniform shear (of the classical type).

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2008.05a
• /
• pp.275-280
• /
• 2008

Shock tube flow measurement has been often troubled with a finite opening time of diaphragm, but there is no systematic work to investigate its effect on the shock tube flow. In the present study, both the experimental and computational works have been performed on the shock tube flows at low pressure ratios. The computational analysis has been performed using the two-dimensional, unsteady, compressible Navier-Stokes equations, based upon a TVD MUSCL finite difference scheme. It is known that the present computational results reproduce the experimental data with good accuracy and simulate successfully the process of diaphragm opening as a function of time. The concept of an imaginary center is introduced to specify the non-centered expansion wave due to a finite opening time of diaphragm. The results obtained show that the diaphragm opening time is reduced as the initial pressure ratio of shock tube increases, leading to the effect of a finite opening time of diaphragm to be more remarkable at low pressure ratios.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.161-164
• /
• 2002

In the present study, a passive control method, using the porous wall and cavity system, is applied to the shock wave/boundary layer interactions in transonic moist air flow. The two-dimensional, unsteady, compressible Navier-Stokes equations, which are fully coupled with a droplet growth equation, are solved by the third-order MUSCL type TVD finite difference scheme. Baldwind-Lomax turbulence model is employed to close the governing equations. In order to investigate the effectiveness of the present control method, the total pressure losses of the flow and the time-dependent behaviour of shock motions are analyzed in detail. The computed results show that the present passive control method considerably reduces the total pressure losses due to the shock/boundary layer interaction in transonic moist air flow and suppresses the unsteady shock wave motions over the airfoil, as well. It is also found that the location of the porous ventilation significantly influences the control effectiveness.

• 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.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.26 no.5
• /
• pp.666-674
• /
• 2002

A rapid expansion of the moist air or stream through transonic nozzle often leads to not-equilibrium condensation shock, causing a considerable amount of energy loss to the entire flow field. Depending on amount of heat released, condensation shock wave occurs in the nozzle and interacts with the boundary layer flow. In the current study, a passive control technique using a porous wall with a plenum cavity underneath is applied for purpose of alleviation the condensation shock wave in a transonic nozzle. A droplet growth equation is incorporated into two-dimensional wavier-Stokes equation systems. 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 transonic wind tunnel is made to validate the present computational results. The results obtained show that the magnitude of condensation shock wave is reduced by the current passive control method.

• Proceedings of the KSME Conference
• /
• 2001.11b
• /
• pp.450-455
• /
• 2001

The time-dependent behavior of nonequilibrium condensation of moist air through the Ludwieg tube is investigated with a computational fluid dynamics(CFD) method. 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 computational results are compared with the previous experiments using the Ludwieg tube with a downstream diaphragm. The results 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 nonequilibrium condensation in the Ludwieg tube should not be ignored even for a very low initial relative humidity, and the periodic excursions of the condensation shock wave are responsible for the total pressure loss.

• Proceedings of the KSME Conference
• /
• 2001.11b
• /
• pp.432-437
• /
• 2001

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 analysis are well compared to the results from the aeroacoutics theory with a good agreement.

• Proceedings of the KSME Conference
• /
• 2001.11b
• /
• pp.640-645
• /
• 2001

The present study is to investigate the characteristics of the impulse waves discharged from the exit of the convergent and divergent pipes. An experiment is carried out using a shock tube with an open end and is compared to the computation of the axisymmetric, compressible, unsteady Euler equations, which are solved by the second-order total variation diminishing(TVD) scheme. For the computational work, some initial compression waves are assumed inside the pipe so that those are identical to the experimental ones of the shock tube. The results show that the peak pressures of the impulse waves discharged from the exit of convergent and divergent pipes decrease with an increase in the wavelength of the initial compression wave. All of the impulse waves have a strong directivity toward the pipe axis, regardless of the exit type of the pipe employed. The impulse waves discharged from the divergent pipe are stronger than those from the straight pipe, while the impulse waves of the convergent pipe are weaker than those from the straight pipe. It is believed that the convergent pipe can playa role of a passive control to reduce the peak pressure of the impulse wave. The present computations represent the experimented impulse waves with a good accuracy.