• 한국유체기계학회 논문집
• /
• 제15권6호
• /
• pp.85-91
• /
• 2012

Numerical simulations on a single stage transonic compressor which is developed by Korea Aerospace Research Institute are carried out and their results are compared with experimental data for cross validations. Comparisons between experimental data and numerical simulation results show good agreements on a performance curve, static pressure and total pressure distributions. CFD results show that there is a clear interaction between tip leakage flow and normal shock in the rotor passage. Tip leakage flows are almost dissipated after the strong normal shock and it forms a strong recirculation near the blade tip. Also a large separation region grows on the suction surface just after the normal shock. As the pressure ratio and blade loading increase, the normal shock line moves upstream and it starts to deviate from the blade leading edge. Then the tip leakage flow does not overcome the strong adverse pressure gradient and flow blockage originated from the tip recirculation region. As a result, the tip leakage flow heads for the neighboring blade leading edge, which results in a compressor stall.

• 한국전산유체공학회지
• /
• 제2권2호
• /
• pp.59-72
• /
• 1997

A numerical model is constructed to simulate the interactions of oblique shock wave / turbulent boundary layer on a strongly heated wall. The heated wall temperature is two times higher than the adiabatic wall temperature and the shock wave is strong enough to induce boundary layer separation. The numerical diffusion in the finite volume method is reduced by the use of a higher order convection scheme(UMIST scheme) which is a TVD version of QUICK scheme. The turbulence model is Chen-Kim two time scale model. The comparison of the wall pressure distribution with the experimental data ensures the validity of this numerical model. The effect of strong wall heating enlarges the separation region upstream and downstream. In order to eliminate the separation, wall suction is applied at the shock foot position. The bleeding slot width is about same as the upstream boundary layer thickness and suction mass flow is 10% of the flow rate in the upstream boundary layer. The final configuration of the shock reflection pattern and the wall pressure distribution approach to the non-viscous value when wall suction is applied.

• 대한수학회논문집
• /
• 제10권3호
• /
• pp.767-781
• /
• 1995

In this paper we present computation of a reflected shock in the hypersonic flow of air with chemical reactions. We consider two dimensional steady inviscid hypersonic flow of air around bodies including chemical reaction effects. At a high Mach number, a strong shock is formed in front of the body when a wedge is placed against the flow. In front of the shock, temperature and pressure increase greatly and the flow is in nonequilibrium state. If the shock hits a wall, then a reflected shock is formed in the nonequilibrium flow region. Behind this reflected shock, the temperature and pressure are very high. We carry out the computation of the reflected shock and the flow behind it. The jump conditions at the reflected shock are presented. A technique combining smooth transforms of domain and implicit difference methods is used to overcome numerical difficulties associated with the lack of resolution behind the shock and near the body.

• 한국가시화정보학회지
• /
• 제10권1호
• /
• pp.40-46
• /
• 2012

Recently micro shock tube is extensively being used in many diverse fields of engineering applications but the detailed flow physics involved in it is hardly known due to high Knudsen number and strong compressibility effects. Unlike the macro shock tube, the surface area to volume ratio for a micro shock tube is very large. This unique effect brings many complexities into the flow physics that makes the micro shock tube different compared with the macro shock tube. In micro shock tube, the inter- molecular forces of working gas can play an important role in specifying the flow characteristics of the unsteady shock wave flow which is essentially generated in all kinds of shock tubes. In the present study, a CFD method was used to predict and visualize the unsteady shock wave flows using the unsteady compressible Navier-Stokes equations, furnished with the no-slip and slip wall boundary conditions. Maxwell's slip equations were used to mathematically model the shock movement at high Knudsen number. The present CFD results show that the propagation speed of the shock wave is directly proportional to the initial pressure and diameter of micro shock tube.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2008년도 학술대회
• /
• pp.337-339
• /
• 2008

As a shock impinges into a vortex of variable strength, complex shock diffraction can occur. Since a vortex has a fixed rotating direction, the shock wave travelling in one direction creates strong asymmetry in the vortex flow field. The process is that first the shock is divided into two parts by the vortex. One part is moving in the adverse direction opposite to the vortex flow which is captured by the vortex center. The other part is moving in the favorable direction, namely, in the direction same as the vortex flow; it is swung around the vortex, accelerating the vortex flow. In this paper we have investigated numerically using ENO scheme how and why the shock-vortex interaction patterns appear so different for different parametric values. Conclusion is that there are three different types of shock-vortex interaction depending on two related parameters: shock Mach number and vortex Mach number. We present a parameter map by which we can discern what type of interaction pattern appears as a shock impinges into a vortex.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2008년 추계학술대회논문집
• /
• pp.337-339
• /
• 2008

Abstract: As a shock impinges into a vortex of variable strength, complex shock diffraction can occur. Since a vortex has a fixed rotating direction, the shock wave travelling in one direction creates strong asymmetry in the vortex flow field. The process is that first the shock is divided into two parts by the vortex. One part is moving in the adverse direction opposite to the vortex flow which is captured by the vortex center. The other part is moving in the favorable direction, namely, in the direction same as the vortex flow; it is swung around the vortex, accelerating the vortex flow. In this paper we have investigated numerically using ENO scheme how and why the shock-vortex interaction patterns appear so different for different parametric values. Conclusion is that there are three different types of shock-vortex interaction depending on two related parameters: shock Mach number and vortex Mach number. We present a parameter map by which we can discern what type of interaction pattern appears as a shock impinges into a vortex.

• 한국전산유체공학회지
• /
• 제14권3호
• /
• pp.76-85
• /
• 2009

The shock wave is deformed and the vortex is elongated simultaneously during the shock-vortex interaction. More precisely, the shock wave is deformed to a S-shape, consisting of a leading shock and a lagging shock by which the corresponding local vortex flows are accelerated and decelerated, respectively: the vortex flow swept by the leading shock is locally expanded and the one behind the lagging shock is locally compressed. As the leading shock escapes the vortex in the order of microseconds, the expanded flow region is quickly changed to a compression region due to the implosion effect. An induced shock is developed here and propagated against the vortex flow. This happens for a strong vortex because the tangential flow velocity of the vortex core is high enough to make the induced-shock wave speed supersonic relative to the vortex flow. For a weak shock, the vortex is basically subsonic and the induced shock wave is absent. For a vortex of intermediate strength, an induced shock wave is developed in the supersonic region but dissipated prematurely in the subsonic region. We have expounded these three shock-vortex interaction patterns that depend on the vortex flow regime using a third-order ENO method and numerical shadowgraphs.

• 한국군사과학기술학회지
• /
• 제14권1호
• /
• pp.147-153
• /
• 2011

A series of nonequilibrium molecular dynamics(NEMD) simulations are performed to investigate the kinetic energy and velocity distributions of molecules in shock waves. In the simulations, argon molecules are used as a medium gas through which shock waves are propagating. The kinetic energy distribution profiles reveals that as a strong shock wave whose Mach number is 27.1 is applied, 39.6% of argon molecules inside the shock wave have larger kinetic energy than molecular ionization energy. This indicates that an application of a strong shock wave to argon gas can give rise to an intense light. The velocity distribution profiles in z direction along which shock waves propagate clearly represent two Maxwell-Boltzmann distributions of molecular velocities in two equilibrium regions and one bimodal velocity distribution profile that is attributed to a nonequilibrium region. The peak appearing in the directional temperature in z direction is discussed on a basis of the bimodal velocity distribution in the nonequilibrium region.

• 한국초전도저온공학회:학술대회논문집
• /
• 한국초전도저온공학회 2000년도 KIASC Conference 2000 / 2000년도 학술대회 논문집
• /
• pp.43-47
• /
• 2000

Two modes of shock waves, a compression shock wave and a thermal shock wave, propagating in He II have been investigated. The shock waves are at a time generated by the impingement of a gasdynamic shock wave onto a He II free surface in the newly developed superfluid shock tube facility. Superconductive temperature sensors, piezo-type pressure transducers and visualization photograph were used for the measurement of them and the phenomena induced by them were investigated in detail. It is found that the compression by a compression shock wave in He II causes temperature drop because He II has negative thermal expansion coefficient. the thermal shock wave is found to be of a single triangular waveform with a limited shock strength. The waveform is similar to that generated by stepwise strong heating from an electrical heater for relatively long heating time. In the experiments at the temperatures near the lambda temperature, no thermal shock wave is sometimes detected in shock compressed He II. It can be understood that shock compression makes He Ii convert to He I in which no thermal shock wave is excited.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2003년도 춘계학술대회
• /
• pp.1742-1747
• /
• 2003

When a shock wave is discharged from the exit of a duct, complicated flow is formed near the duct exit. The flow field is much more complicated under the ground effects or any other objects near the exit of a duct, such as the circumstance near the exit of the high-speed railway tunnel. The resulting flow is essentially three-dimensional unsteady with the effects of strong compressibility. In the current study, three-dimensional flow fields of the weak shock wave which is discharged from the exit of a duct are numerically investigated using a CFD method. Computations are performed for the weak shock wave in the range below 1.5. The results obtained show that the directivity and magnitude of the weak shock discharged strongly depend upon the Mach number of initial shock wave and are significantly influenced by the ground effects.