• 한국전산유체공학회:학술대회논문집
• /
• 1999.05a
• /
• pp.85-89
• /
• 1999

Numerical Analysis has been done for the supersonic off-design jet flow due to the pressure difference between the jet and the ambient fluid. The difference of pressure generates an oblique shock or an expansion wave at the nozzle exit, The waves reflect repeatedly at the center axis and on the sonic surface in the shear layer, and the pressure difference is resolved across these waves interacted with the turbulence mixing layer. In this paper, the axi-symmetric Navier-Stokes equation has been used with two equation $k-{\varepsilon}$ turbulence closure model. The second order TVD scheme with flux limiters, based on the flux vector split by the smooth eigenvalue split, has been used to capture internal shocks and other discontinuities. The correction term for the compressible flow and the damping function are used in the turbulence model. Numerical calculations have been done to analyze the off-design jet flow due to the pressure difference. The variation of pressure along the flow axis is compared with an experimental result and other numerical result. The characteristics of the interaction between the shock cell and the turbulence mixing layer have been analyzed.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.31 no.2 s.257
• /
• pp.109-115
• /
• 2007

Airframe-integrated Scramjet engines of NASA Langley type consist of a compressor, a combustion chamber and a nozzle. When some disturbances occur in one module of the engine, its influences are propagated to other modules. In this study, it is investigated numerically how shock waves were caused by thermal choking in one module propagate upstream and how they influence adjacent modules. The calculations are carried out in 2-dimensional supersonic viscous flow model using explicit TVD scheme in generalized coordinates. The adverse pressure gradient caused by heat addition brings about separation of the wall boundary layers and formation of the oblique shock wave that proceed to upstream. This moving shock wave formed one module blocks the flow coming into the adjacent modules, which makes the modules unstarted.

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

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2008.03a
• /
• pp.597-601
• /
• 2008

Three-dimensional structures of detonation wave propagating in circular tube were investigated. Inviscid fluid dynamics equations coupled with a conservation equation of reaction progress variable were analyzed by a MUSCL-type TVD scheme and four stage Runge-Kutta time integration. Variable-$\gamma$ formulation was used to account for the variable properties between unburned and burned states and the chemical reaction was modeled by using a simplified one-step irreversible kinetics model. The computational code was parallelized based on domain decomposition technique using MPI-II message passing library. The computations were carried out using a home made Windows based PC cluster having 160 AMD AthloxXP and Athlon64 processor. The computational domain consisted of through a roundshaped tube with wall conditions. As an initial condition, analytical ZND solution was distributed over the computational domain with disturbances. The disturbances has circumferential large gradient. The unsteady computational results in three-dimension show the detailed mechanisms of multi-cell mode of detonation wave instabilities resulting diamond shape in smoked-foil record.

• Journal of Mechanical Science and Technology
• /
• v.16 no.11
• /
• pp.1511-1521
• /
• 2002

High-speed moist air or steam flow has long been of important subject in engineering and industrial applications. Of many complicated gas dynamics problems involved in moist air flows, the most challenging task is to understand the nonequilibrium condensation phenomenon when the moist air rapidly expands through a flow device. Many theoretical and experimental studies using supersonic wind tunnels have devoted to the understanding of the nonequilibrium condensation flow physics so far. However, the nonequilibrium condensation can be also generated in the subsonic flows induced by the unsteady expansion waves in shock tube. The major flow physics of the nonequilibrium condensation in this application may be different from those obtained in the supersonic wind tunnels. In the current study, the nonequilibrium condensation phenomenon caused by the unsteady expansion waves in a shock tube is analyzed by using the two-dimensional, unsteady, Navier-Stokes equations, which are fully coupled with a droplet growth equation. The third-order TVD MUSCL scheme is applied to solve the governing equation systems. The computational results are compared with the previous experimental data. The time-dependent behavior of nonequilibrium condensation of moist air in shock tube is investigated in details. The results show that the major characteristics of the nonequilibrium condensation phenomenon in shock tube are very different from those in the supersonic wind tunnels.

• Journal of Mechanical Science and Technology
• /
• v.16 no.11
• /
• pp.1448-1456
• /
• 2002

The present study addresses the flow characteristics involved in the self-induced oscillations of the underexpanded jet impinging upon a cylindrical body. Both experiment and computational analysis are carried out to elucidate the shock motions of the self-induced oscillations and to find the associated major flow factors. The underexpanded sonic jet is made from a nozzle and a cylindrical body is placed downstream to simulate the impinging jet upon an obstacle. The computational analysis using TVD scheme is applied to solve the axisymmetric, unsteady, inviscid governing equations. A Schlieren system is employed to visualize the self-induced oscillations generated in flow field. The data of the shock motions are obtained from a high-speed video system. The detailed characteristics of the Mach disk oscillations and the resulting pressure variations are expatiated using the time dependent data of the Mach disk positions. The mechanisms of the self-induced oscillations are discussed in details based upon the experimental and computational results.

• The Bulletin of The Korean Astronomical Society
• /
• v.40 no.1
• /
• pp.65.2-65.2
• /
• 2015

We developed a three-dimensional magnetohydrodynamic (MHD) code to reproduce the structure of a solar wind, the properties of a coronal mass ejection (CME) and the interaction between them. This MHD code is based on the finite volume method incorporating total variation diminishing (TVD) scheme with an unstructured grid system. In particular, this grid system can avoid the singularity at the north and south poles and relax tight CFL conditions around the poles, both of which would arise in a spherical coordinate system (Tanaka 1994). In this model, we first apply an MHD tomographic method (Hayashi et al. 2003) to interplanetary scintillation (IPS) observational data and derive a solar wind from the physical values obtained at 50 solar radii away from the Sun. By comparing the properties of this solar wind to observational data obtained near the Earth orbit, we confirmed that our model captures the velocity, temperature and density profiles of a solar wind near the Earth orbit. We then insert a spheromak-type CME (Kataoka et al. 2009) into the solar wind to reproduce an actual CME event. This has been done by introducing a time-dependent boundary condition to the inner boundary of our simulation domain. On the basis of a comparison between a simulated CME and observations near the Earth, we discuss the physics involved in an ICME interacting with a solar wind.

• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.355-360
• /
• 2001

This study addresses a computational work of the impulsive wave which is discharged from the open end of a pipe. An initial compression wave inside the pipe is assumed to propagate toward atmosphere. The over pressure and wave-length of the initial compression wave are changed to investigate the characteristic values of the impulsive wave. The second order total variation diminishing (TVD) scheme is employed to solve the axisymmetric, compressible, unsteady Euler equations. The relationship between the initial compression wave form and impulsive wave is characterized in terms of the peak pressure of the impulsive wave and its directivity. The results obtained show that for the initial compression wave of a large wave-length the peak pressure of the impulsive wave does not depend on the over pressure of the initial compression wave, but for the initial compression wave of a very short wave-length, like a shock wave, the peak pressure of the impulsive wave is increased with an increase in the over pressure of the initial compression wave. The directivity of the impulsive wave to the pipe axis becomes significant with a decrease in the wave-length of the initial compression wave.

• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.299-304
• /
• 2001

The current study describes experimental and computational work on the passive control of the steady and unsteady condensation shock waves, which are generated in a transonic nozzle. The bleed slots are installed on the contoured wall of the transonic nozzle in order to control the magnitude of the condensation shock wave and its oscillations. For computations, a droplet growth equation is incorporated into the two-dimensional Navier-Stokes equation systems. Computations are carried out using a third-order MUSCL type TVD finite-difference scheme with a second-order tractional time step. Baldwin-Lomax turbulence model is employed to close the governing equations. An experiment using an indraft transonic wind tunnel is made to validate the computational results. The current computations represented well the experimental flows. From both the experimental and computational results it is found that the magnitude of the condensation shock wave in the bleed slotted nozzle is significantly reduced, compared with no passive control of solid wall. The oscillations of the condensation shock wave are successfully suppressed by a bleed slot system.

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