• 한국전산유체공학회:학술대회논문집
• /
• 1995.10a
• /
• pp.182-187
• /
• 1995

The present paper explains some advancement made by the authors for the compressible flow computation of the Euler equations based on the unstructured grid and vertex- centered finite volume method. Accurate solutions to the unsteady axisymmetric shock wave propagation problems and three-dimensional airplane flows have been obtained by a high-order upwind TVD and FCT schemes. Unstructured grid adaption is made for the unsteady shock wave problems by the dynamic h-refinement/unrefinement procedure and for the three-dimensional steady flows by the Delaunay point-insertion method to generate three-dimensional tetrahedral mesh enrichment. Some physics of the shock wave diffraction phenomena and three-dimensional airplane flow are discussed.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.10a
• /
• pp.141-142
• /
• 2003

When the under-expanded supersonic jet impinges on the perpendicular plate, it is well known that the self-induced flow oscillation occurs at the specific conditions. This phenomenon is related with the noise problems of aeronautical and other industrial engineering. But, the very complicated flow field is formed and it is difficult to clear the flow structure and the mechanism of oscillation. This paper aims to clear the characteristics of flow field and the wave pattern during the under-expanded supersonic jet impinges on the plate. The numerical calculation was carried out using the TVD numerical method. In this paper, the flow visualization, the pressure fluctuation on the surface of plate and the mechanism of oscillation are discussed.

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