• 한국연소학회:학술대회논문집
• /
• 2005.10a
• /
• pp.15-19
• /
• 2005

Three-dimensional structures of unsteady detonation wave propagating through a square-shaped tube is studied using computational method and parallel processing. Inviscid fluid dynamics equations coupled with variable-${\gamma}$ formulation and simplified one-step Arrhenius chemical reaction model were analysed by a MUSCL-type TVD scheme and four stage Runge-Kutta time integration. Results in three dimension show the two unsteady detonation wave propagating mode, the Rectangular and diagonal mode of detonation wave instabilities. Two different modes of instability showed the same cell length but different cell width and the geometric similarities in smoked-foil record.

• International Journal of Fluid Machinery and Systems
• /
• v.2 no.3
• /
• pp.248-253
• /
• 2009

This paper presents a numerical study for unsteady flows in a high-pressure steam turbine with a partial admission stage. Compressible Navier-Stokes equations are solved by the high-order high-resolution finite-difference method based on the fourth-order compact MUSCL TVD scheme, Roe's approximate Riemann solver, and the LU-SGS scheme. The SST-model is also solved for evaluating the eddy-viscosity. The unsteady two-dimensional flows through whole nozzle-rotor cascade channels considering a partial admission are numerically investigated. 108 nozzle passages with two blockages and 60 rotor passages are simultaneously calculated. The influence of the flange in the nozzle box to the lift of rotors is predicted. Also the efficiency of the partial admission stage changing the number of blockages and the number of nozzles is parametrically predicted.

• 유체기계공업학회:학술대회논문집
• /
• 2000.12a
• /
• pp.169-175
• /
• 2000

A design method for transonic turbine blades is developed based on Navier-Stokes equations. The present computing process is done on the four separate steps, 1.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. Applications are made to the VKI transonic rotor blades. Computed results are analyzed with respect to the aerodynamic performance and are compared with the experimental data.

• 한국전산유체공학회:학술대회논문집
• /
• 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 Mechanical Science and Technology
• /
• v.16 no.4
• /
• pp.557-563
• /
• 2002

In this paper, the characteristics of flow and spray motions affected by from piston head configurations were investigated numerically. Calculations were carried out from intake process to the end of compression. GTT (Generalized Tank and Tube method) code, which includes a third order upwind Chakravarthy-Osher TVD scheme and k-$\varepsilon$ turbulence model with fuel spray analysis was used for the calculations. As a results, piston heads with smaller radii of curvature were found to give stronger reverse tumble than those with larger radii of curvature. Similar results are shown in the convection and diffusion of fuel sprays.

• Proceedings of the KSME Conference
• /
• 2003.11a
• /
• pp.544-549
• /
• 2003

When a shock wave arrives at an open end of tube, an impulse wave is discharged from the tube exit and causes serious noise and vibration problems. In the current study, the effect of the cross-sectional area of tube on the impulse wave is numerically investigated using a CFD method. The Harten-Yee's total variation diminishing(TVD) scheme is used to solve the axisymmetric, two-dimensional, unsteady, compressible Euler equations. With three different cross-sectional areas of tube, the Mach number of the incident shock wave $M_{s}$ is varied between 1.01 and 1.5. The results obtained show that the directivity and magnitude of impulse wave strongly depend upon the Mach number of incident shock wave and are influenced by the tube area. It is also known that the tube cross-sectional area significantly affects the magnitude of impulse wave at or near the tube axis.

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

• Journal of The Korean Astronomical Society
• /
• v.28 no.2
• /
• pp.223-243
• /
• 1995

We describe the implementation of a multi-dimensional numerical code to solve the equations for idea! magnetohydrodynamics (MHD) in cylindrical geometry. It is based on an explicit finite difference scheme on an Eulerian grid, called the Total Variation Diminishing (TVD) scheme, which is a second-order-accurate extension of the Roe-type upwind scheme. Multiple spatial dimensions are treated through a Strang-type operator splitting. Curvature and source terms are included in a way to insure the formal accuracy of the code to be second order. The constraint of a divergence-free magnetic field is enforced exactly by adding a correction, which involves solving a Poisson equation. The Fourier Analysis and Cyclic Reduction (FACR) method is employed to solve it. Results from a set of tests show that the code handles flows in cylindrical geometry successfully and resolves strong shocks within two to four computational cells. The advantages and limitations of the code are discussed.

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

• Journal of computational fluids engineering
• /
• v.21 no.2
• /
• pp.1-11
• /
• 2016

Considerable efforts are required to develop a monotone, robust and stable high-order numerical scheme for solving the hyperbolic system. The discontinuous Galerkin(DG) method is a natural choice, but elimination of the spurious oscillations from the high-order solutions demands a new development of proper limiters for the DG method. There are several available limiters for controlling or removing unphysical oscillations from the high-order approximate solution; however, very few studies were directed to analyze the exact role of the limiters in the hyperbolic systems. In this study, the performance of the several well-known limiters is examined by comparing the high-order($p^1$, $p^2$, and $p^3$) approximate solutions with the exact solutions. It is shown that the accuracy of the limiter is in general problem-dependent, although the Hermite WENO limiter and maximum principle limiter perform better than the TVD and generalized moment limiters for most of the test cases. It is also shown that application of the troubled cell indicators may improve the accuracy of the limiters under some specific conditions.