• Journal of computational fluids engineering
• /
• v.2 no.2
• /
• pp.9-17
• /
• 1997

The performance of several turbulence models in computing an axisymmetric supersonic base flow is investigated. A compressible Navier-Stokes code, which incorporates k-ε, k-ω model and Reynolds stress closure with three kinds of pressure-strain correlation model, has been developed using implicit LU-SGS algorithm with second-order upwind TVD scheme. Numerical computations have been carried out for Herrin and Dutton's base flow. It is observed that the two-equation models give large backward axial velocity approaching to the base and somewhat larger variation of base pressure distribution than the Reynolds stress model. It is also found that the Reynolds stress model with third order pressure-strain model in the anisotropy tensor predicts most accurate mean flow field.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.5
• /
• pp.1280-1290
• /
• 1995

In order to investigate the compression wave propagating in a high-speed railway tunnel, a numerical calculation was applied to the wave phenomenon occurring in a model tunnel. Unsteady, one-dimensional inviscid or viscous flows were solved by an explicit TVD scheme, and the calculated flows were compared with the results of measurement in real tunnels. Tunnel noises caused by emission of the compression wave were characterized in terms of excess pressure of compression wave, pressure gradient in the wave front and width of the compression wave. Calculated attenuation, pressure gradient and width of compression wave with the propagating distance agreed with the results of measurement in the real tunnels. The results also show that tunnel noises are proportional to the train velocity entering the tunnel.

• Proceedings of the KSME Conference
• /
• 2000.11b
• /
• pp.544-549
• /
• 2000

The effects of nonequilibrium condensation on the shock boundary layer interaction over a transonic bump model were investigated experimentally and numerically. An experiment was conducted using a supersonic indraft wind tunnel. A droplet growth equation was incorporated into two-dimensional Navier-Stokes equation systems. Computations were carried out using a third-order MUSCL type TVD finite-difference scheme with a second-order fractional time step. Computations compared with the experimental results. Nonequilibirum condensation suppressed the boundary layer separation and the pressure fluctuations due to the shock boundary layer interaction. Especially the nonequilibrium condensation was helpful to suppress the high frequency components of the pressure fluctuations.

• 한국가시화정보학회:학술대회논문집
• /
• 2003.11a
• /
• pp.95-98
• /
• 2003

When a shock wave discharges from an open end of a duct, an impulse wave is generated outside the duct, causing serious noise and vibration problems. The magnitude of the impulse wave can be reduced by installing of a perforated duct. In the current study, the characteristics of the impulse wave discharged from the exit of a perforated duct are numerically investigated. A TVD (total variation diminishing) scheme is used to solve the unsteady, axisymmetric, compressible Euler equations. In computations, the porosity of a perforated pipe $(\sigma)$ and the Mach number of incident shock wave $(M_s)$ are varied in the range of $\sigma=0\~19\%\;and\;M_s=1.01\~1.50$, respectively. The results show that the directivity and magnitude of impulse wave strongly depend upon the Mach number of incident shock wave and the porosity of the perforated pipe. The present CFD results are in close agreement with experimental results.

• 한국가시화정보학회:학술대회논문집
• /
• 2004.11a
• /
• pp.72-75
• /
• 2004

The present study describes a computational work to investigate detailed behaviors of the twin shock waves discharged from the exits of two-parallel ducts. In computations, the Yee-Roe-Davis's TVD scheme was used to solve the unsteady, three-dimensional, inviscid, compressible, Euler equations. The distance between two ducts is varied and the Mach number of the incident shock wave is changed below 2.0. The results obtained show that on the symmetric axis between two-parallel ducts, the maximum pressure achieved by the merge of twin shock waves and its location strongly depend upon the distance between two-parallel ducts and the Mach number of the incident shock wave. It is also found that the twin shock waves discharged from the exits of two-parallel ducts leads to the complicated flow fields, such as Mach stem, spherical waves, and vertical structures.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.165-168
• /
• 2002

This paper describes the dynamics of the weak shock wave propagating inside some kinds of branched pipe bends. Computations are carried out by solving the two-dimensional, compressible, unsteady Euler Equations. The second-order TVD(Total Variation Diminishing) scheme is employed to discretize the governing equations. For computations, two types of branched pipe($90^{\circ}$ branch,$45^{\circ}$ branch) with a diameter of D are used. The incident normal shock wave is assumed at D upstream of the pipe bend entrance, and its Mach number is changed between 1.1 and 2.4. The flow fields are numerically visualized by using the pressure contours and computed schlieren images. The comparison with the experimental data performed for the purpose of validation of computational work. Reflection and diffraction of the propagating shock wave are clarified. The present computations predicted the experimented flow field with a good accuracy.

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

This study introduced a lattice Boltzmann computational scheme capable of modeling thermo hydrodynamic flows with simpler equilibrium particle distribution function compared with other models. The equilibrium particle distribution function is the local Maxwelian equilibrium function in this model, with all the constants uniquely determined. The characteristics of the proposed model is verified by calculation of the sound speeds, and the shock tube problem. In the lattice Boltzmann method, a thermal fluid or compressible fluid model simulates the reflection of a weak shock wave colliding with a sharp wedge having various angles $\theta$$\sub$w/. Theoretical results using LBM are satisfactory compared with the experimental result or the TVD.

• Proceedings of the KSME Conference
• /
• 2001.11b
• /
• pp.456-461
• /
• 2001

This paper depicts the weak shock wave propagating inside some kinds of pipe bends. Computational work is to solve the two-dimensional, compressible, unsteady Euler Equations. The second-order TVD scheme is employed to discretize the governing equations. For the computations, the incident normal shock wave is assumed at the entrance of the pipe bend, and its Mach number is changed between 1.1 and 1.7. The turning angle and radius of the curvature of the pipe bend are changed to investigate the effects on the shock wave structure. The present computational results clearly show the shock wave reflection and diffraction occurring in the pipe bend. In particular, the vortex generation, which occurs at the edge of the bend, and its shedding mechanism are discussed in details.

• Proceedings of the KSME Conference
• /
• 2001.11b
• /
• pp.462-467
• /
• 2001

The present study addresses a computational work of the weak shock wave propagatings inside a silencer system of automobile exhaust pipe. Four different types of the silencer systems and the initial shock wave Mach number $M_s$ of $1.01\sim1.30$ are applied to investigate their effects on the noise reduction and the flow field in a silencer system. The results obtained from the present computational work are compared with the experimental results. The second order total variation diminishing (TVD) scheme is employed to solve the two dimensional, compressible, unsteady, Euler equations. The present computational results predict the experimental results with a quite good accuracy. Of the four silencer systems applied, the most desirable silencer system to reduce the peak pressure at the exit of the exhaust pipe is discussed.

• The Bulletin of The Korean Astronomical Society
• /
• v.35 no.2
• /
• pp.46.2-46.2
• /
• 2010

A magnetosonic wave is a longitudinal wave propagating perpendicularly to the magnetic fields and involves compression and rarefaction of the plasma. Lee and Kim (2000) investigated the theoretical solution for the evolution of nonlinear magnetosonic waves in the homogeneous space which adopt the approach of simple waves. We confirm the solution using a one-dimensional MHD code with Total Variation Diminishing (TVD) scheme. Then we apply the solution for the solar wind profiles. We examined the properties of nonlinear waves for the various initial perturbations at near the Lagrangian (L1) point. Also we describe waves steepening process while the shock is being formed by assuming different timescales for a driving source.