• 한국전산유체공학회:학술대회논문집
• /
• 1997.10a
• /
• pp.99-105
• /
• 1997

A finite-difference method based on conservative supra characteristic method type upwind flux difference splitting has been developed to study the nonequilibrium chemically reacting inviscid flow. For nonequilibrium air, NS-1 species equations were strongly coupled with flowfield equations through convection and species production terms. Inviscid nonequilibrium chemically reacting air mixture flows over Blunt-body were solved to demonstrate the capability of the current method. At low altitude flight conditions the nonequilibrium air models predicted almost the same temperature, density and pressure behind the shock as equilibrium flow: however, at high altitudes they showed substantial differences due to nonequilibrium chemistry effect. The new nonequilibrium chemically reacting upwind flux difference splitting mettled can be extended to viscous flow and multi-dimensional flow conditions.

• Journal of the Korea Institute of Military Science and Technology
• /
• v.4 no.2
• /
• pp.202-214
• /
• 2001

In this paper the upwind flux difference splitting Navier-Stokes method has been applied to study quasi one-dimensional nozzle flow and axisymmetric sphere-cone($5^{\circ}$) flow for the perfect gas, the equilibrium and the nonequilibrium chemically reacting hypersonic flow. The effective gamma(${ \tilde{\gamma}}$), enthalpy to internal energy ratio was used to couple chemistry with the fluid mechanics for equilibrium chemically reacting air. The influences of the various altitude(30km, 50km) at Mach number(15.0, 20.0) were investigated. The equilibrium shock position was located farthest downstream when compared with those of perfect gas in a quasi one-dimensional nozzle. The equilibrium shock thickness over the blunt body region was much thinner than that of perfect gas shock.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.30 no.4 s.247
• /
• pp.320-327
• /
• 2006

The computations of chemically reacting laminar and turbulent flows are performed using the preconditioned Navier-Stokes solver coupled with turbulent transport and multi-species equations. A low-Reynolds number $k-\varepsilon$ turbulence model proposed by Chien is used. The presence of the turbulent kinetic energy tenn in the momentum equation can materially affect the overall stability of the fluids-turbulence system. Because of this coupling effect, a fully coupled formulation is desirable and this approach is taken in the present study. Choi and Merkle's preconditioning technique is used to overcome the convergence difficulties occurred at low speed flows. The numerical scheme used for the present study is based on the implicit upwind ADI algorithm and is validated through the comparisons of computational and experimental results for laminar methane-air diffusion flame and $ H_2/O_2$ reacting turbulent shear flow. Preconditioning formulation shows better convergence characteristics than that of non-preconditioned system by approximately five times as much.

• Journal of computational fluids engineering
• /
• v.6 no.2
• /
• pp.32-39
• /
• 2001

Hypersonic analysis on the KSR-Ⅲ payload configuration has been performed using an axisymmetric Navier-Stokes code. A numerical code based on the Harten and Yee's upwind TVD scheme with simplified curve fits in the chemically reacting equilibrium air was developed. The carbuncle phenomenon on detached shock in front of the payload is controlled by using pressure gradients to tune the dissipation. Chemically reacting equilibrium computations for the reentry flight conditions of Mach No. 10.2, 8, 4.9 are presented and compared with the results of calorically perfect gas.

• 한국전산유체공학회:학술대회논문집
• /
• 2000.05a
• /
• pp.166-171
• /
• 2000

A two-dimensional/axisymmetric CSCM upwind flux difference splitting Wavier-Stokes method has been developed to study the finite rate chemically react-ing invisicd and viscous hypersonic flows over blunt-body. A upwind method was chosen due to its robustness in capturing the strong bow shock waves. For the nonequilibrium chemically reacting air, NS-I species conservation equations were strongly coupled with flowfield equations through convection and species production terms. The nonequilibrium wall pressure and heat transfer rate distributions along the vehicle were compared with those from equilibrium and perfect gas calculations. The nonequilibrium species distribution shows the reduced concentrations of O and N species when compared with equilibrium species distribution. The solutions resolved strong bow shock waves md heat transfer rate very accurately when compared with central difference schemes.

• Journal of the Korean Society of Combustion
• /
• v.1 no.1
• /
• pp.83-91
• /
• 1996

A numerical study is carried out to examine the ignition and propagation process of detonation wave in SCRam-accelerator operating in superdetonative mode. The time accurate solution of Reynolds averaged Navier-Stokes equations for chemically reacting flow is obtained by using the fully implicit numerical method and the higher order upwind scheme. As a result, it is clarified that the ignition process has its origin to the hot temperature region caused by shock-boundary layer interaction at the shoulder of projectile. After the ignition, the oblique detonation wave is generated and propagates toward the inlet while constructing complex shock-shock interaction and shock-boundary layer interaction. Finally, a standing oblique detonation wave is formed at the conical ramp.

• Journal of computational fluids engineering
• /
• v.1 no.1
• /
• pp.123-141
• /
• 1996

In this paper the upwind flux difference splitting Navier-Stokes method has been applied to study the perfect gas and the equilibrium chemically reacting hypersonic flow over an axisymmetric sphere-cone(5°) geometry. The effective gamma(γ), enthalpy to internal energy ratio was used to couple chemistry with the fluid mechanics for equilibrium chemically reacting air. The test case condition was at altitude(30km) and Mach number(15). The equilibrium shock thickness over the blunt body region was much thinner than that of perfect gas shock. The pressure difference between perfect gas and equilibrium gas was about 3 ∼ 5 percent. The heat transfer coefficient were also calculated. The results were compared with VSL results in order to validate the current numerical analysis. The results from current method were compared well VSL results ; however, not well at near nose. The proper boundary condition and grid system will be studied to improve the solution quality.

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

In this paper the upwind flux difference splitting Navier-Stokes method has been applied to study the perfect gas and the equilibrium chemically reacting hypersonic flow over an axisymmetric sphere-cone($5^{\circ}$) geometry. The effective gamma($\bar{r}$), enthalpy to internal energy ratio was used to couple chemistry with the fluid mechanics for equilibrium chemically reacting air. The test case condition was at altitude(30Km) and Mach number(15). The equilibrium shock thickness over the blunt body region was much thinner than that of perfect gas shock. The pressure difference between perfect gas and equilibrium gas was about $3\sim5$ percent. The skin friction coefficient and heat transfer coefficient were also calculated.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.23 no.3
• /
• pp.157-178
• /
• 2019

The present correspondence is conveyed on to consider the fascinating and novel characteristics of radiative hydromagnetic convective flow of a chemically reacting fluid over an exponentially accelerated inclined porous plate. Exact solutions for the fluid velocity, temperature and species concentration, under Boussinesq approximation, are obtained in closed form by the two term perturbation technique. The interesting parts of thermal dispersing outcomes are accounted in this correspondence. Graphical evaluation is appeared to depict the trademark direct of introduced parameters on non dimensional velocity, temperature and concentration profiles. Also, the numerical assortment for skin friction coefficient, Nusselt number and Sherwood number is examined through tables. The certification of current examination is confirmed by making an examination with past revelations available in composing, which sets a benchmark for utilization of computational approach.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2004.03a
• /
• pp.62-68
• /
• 2004

A comprehensive numerical analysis has been carried out for both non-reacting and reacting flows in a scramjet engine combustor with and without a cavity. The theoretical formulation treats the complete conservation equations of chemically reacting flows with finite-rate chemistry of hydrogen-air. Turbulence closure is achieved by means of a k-$\omega$ two-equation model. The governing equations are discretized using a MUSCL-type TVD scheme, and temporally integrated by a second-order accurate implicit scheme. Transverse injection of hydrogen is considered over a broad range of injection pressure. The corresponding equivalence ratio of the overall fuel/air mixture ranges from 0.167 to 0.50. The work features detailed resolution of the flow and flame dynamics in the combustor, which was not typically available in most of the previous studies. In particular, the oscillatory flow characteristics are captured at a scale sufficient to identify the .underlying physical mechanisms. Much of the flow unsteadiness is related not only to the cavity, but also to the intrinsic unsteadiness in the flow-field. The interactions between the unsteady flow and flame evolution may cause a large excursion of flow oscillation. The roles of the cavity, injection pressure, and heat release in determining the flow dynamics are examined systematically.