• Journal of computational fluids engineering
• /
• v.10 no.1
• /
• pp.24-30
• /
• 2005

Unstructured grid system is suitable for flows of complex geometries. For problems with moving boundary walls, the grid system must be time-dependently changing and deforming according to the movement of the boundaries when we use a body fitted grid system. In this paper, a new moving-grid finite-volume method on unstructured grid system is proposed and developed for unsteady compressible flows with shock waves. To assure geometric conservation laws on moving grid system, a control volume on the space-time unified domain is adopted for estimating numerical flux. The method is described and applied for two-dimensional flows.

• Journal of computational fluids engineering
• /
• v.4 no.1
• /
• pp.53-61
• /
• 1999

Two-dimensional, unsteady, incompressible and compressible Navier-Stokes codes are developed for the computation of the viscous turbulent flow over high-lift airfoils. The compressible code involves a conventional upwind-differenced scheme for the convective terms and LU-SGS scheme for temporal integration. The incompressible code with pseudo-compressibility method also adopts the same schemes as the compressible code. Three two-equation turbulence models are evaluated by computing the flow over single and multi-element airfoils. The compressible and incompressible codes are validated by predicting the flow around the RAE 2822 transonic airfoil and the NACA 4412 airfoil, respectively. In addition, both the incompressible and compressible code are used to compute the flow over the NLR 7301 airfoil with flap to study the compressible effect near the high-loaded leading edge. The grid systems are efficiently generated using Chimera overlapping grid scheme. Overall, the κ-ω SST model shows closer agreement with experiment results, especially in the prediction of adverse pressure gradient region on the suction surfaces of high-lift airfoils.

• Journal of KSNVE
• /
• v.11 no.2
• /
• pp.285-291
• /
• 2001

The impingement of a weak shock wane discharged from the open end of a shock tube upon a flat plate was investigated using shock tube experiments and numerical simulations. Harten-Yee Total Variation Diminishing method was used to solve axisymmetric, unsteady, compressible flow governing equations. Experiments were carried out to validate the present computations. The effects of the flat plate and baffle plate sizes on the impinging flow field over the flat plate were investigated. Shock Mach number was varied in the range from 1.05 to 1.20. The distance between the plate and shock tube was changed to investigate the effect on the peak pressure. From both the results of experiments and computations we obtained a good empirical equation to predict the peak pressure on the flat plate.

• Proceedings of the KSME Conference
• /
• 2000.11b
• /
• pp.562-567
• /
• 2000

The Impingement of a weak shock wave discharged from the open end of a shock tube upon a flat plate was investigated using shock tube experiments and numerical simulations. Harten-Yee Total Variation Diminishing method was used to solve axisymmetric, unsteady, compressible flow governing equations. Experiments were carried out to validate the present computations. The effects of the flat plate and baffle plate sizes on the impinging flow field over the flat plate were investigated. Shock Mach number was vaned in the range from 1.05 to 1.20. The distance between the plate and shock tube was changed to investigate the effect on the peak pressure. From both the results of experiments and computations we obtained a good empirical equation to predict the peak pressure on the flat plate.

• Proceedings of the KSME Conference
• /
• 2003.11a
• /
• pp.550-555
• /
• 2003

When a shock wave arrives at an open end of tube, an impulse wave is discharged from the tube exit and complicated flow is formed near tube exit. The flow field is influenced by the cross-sectional geometry of tube exit, such as circular, square, rectangular, trapezoid and etc. In the current study, three-dimensional propagation characteristics of impulse wave discharged from the tube exit with non-circular cross section are numerically investigated using a CFD method. Total variation diminishing (TVD) scheme is used to solve the three-dimensional, unsteady, compressible Euler equations. Computations are performed for the Mach numbers of the incident shock wave $M_{s}$ below 1.5. The results obtained show that the peak pressure of the impulse wave and propagation directivity depends on the cross-sectional geometry of tube exit and the Mach number of incident shock wave.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.11a
• /
• pp.191-196
• /
• 1998

A two-dimensional airfoil flying over a wavy wall is calculated by solving the unsteady Euler equation. Unsteady Transonic flow over an NACA00012 airfoil in pitching motion has been computed for code validation. Some numerical results for NACA6409 airfoil under different wave number, wave length, fly height are presented. The numerical results show the variation of lift and pitching moment coefficients are increased as wave length decrease.

• 한국전산유체공학회:학술대회논문집
• /
• 2001.05a
• /
• pp.121-126
• /
• 2001

An accelerating flow field through a compressor cascade is studied numerically by unsteady computational simulation. The two-dimensional Navier-Stokes equations for compressible flow is used for the study of unsteady high incidence angle flow, with preconditioning scheme to cover the wide range of Mach number and $\kappa-\omega$ model for the turbulent viscous flow analysis. A DCA(double circular arc) compressor blade is accelerated artificially in this study to understand the unsteady effect by comparing the present results with the existing steady-state experimental and computational results. Also, the accelerating flow field during the starting phase of gas turbine is studied with actual experimental data for the understanding of flow field and performance characteristics at off-design condition.

• 한국가시화정보학회:학술대회논문집
• /
• 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.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.31 no.3 s.258
• /
• pp.217-224
• /
• 2007

A numerical analysis for a silencer with three baffles of 120mm tank gun has been performed. The Reynolds-Averaged Wavier-Stokes equations with Baldwin-Lomax turbulence model were employed to compute unsteady, compressible flow inside the tank gun and the silencer. An axisymmetric computational domain was constructed by using 12 multi block chimera grids. The resolution of flow field is observed by depicting calculated pressure and muzzle brake force. The peak blast pressure and noise through the silencer reduced approximately 99% and 41dB in comparison to the tank gun without the silencer at near filed.

• 한국가시화정보학회:학술대회논문집
• /
• 2002.11a
• /
• pp.113-116
• /
• 2002

The shock wave discharged from an annular duct leads to very complicated flow features, such as Mach stem, spherical waves, and vortex rings. In the current study, the merging phenomenon and propagation characteristics of the shock wave are numerically investigated using a CFD method. The Harten-Yee's total variation diminishing (TVD) scheme is used to the unsteady, axisymmetric, two-dimensional, compressible Euler equations. The Mach number of incident shock wave $M_s$ is varied in the range below 2.0. The computational results are visualized to observe the major features of the annular shock waves discharged from the tube. On the symmetric axis, the peak pressure produced by the shock wave and its location depend upon strongly the radius of the annular tubes. A Mach stem is generated along the symmetric axis of the annular tubes.