• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2003.05a
• /
• pp.91-93
• /
• 2003

A comprehensive numerical study is carried out to investigate for the understanding of the flow evolution and flame development in a supersonic combustor with normal injection of ncumally injecting hydrogen in airsupersonic flows. The formulation treats the complete conservation equations of mass, momentum, energy, and species concentration for a multi-component chemically reacting system. For the numerical simulation of supersonic combustion, multi-species Navier-Stokes equations and detailed chemistry of H2-Air is considered. It also accommodates a finite-rate chemical kinetics mechanism of hydrogen-air combustion GRI-Mech. 2.11[1], which consists of nine species and twenty-five reaction steps. Turbulence closure is achieved by means of a k-two-equation model (2). The governing equations are spatially discretized using a finite-volume approach, and temporally integrated by means of a second-order accurate implicit scheme (3-5).The supersonic combustor consists of a flat channel of 10 cm height and a fuel-injection slit of 0.1 cm width located at 10 cm downstream of the inlet. A cavity of 5 cm height and 20 cm width is installed at 15 cm downstream of the injection slit. A total of 936160 grids are used for the main-combustor flow passage, and 159161 grids for the cavity. The grids are clustered in the flow direction near the fuel injector and cavity, as well as in the vertical direction near the bottom wall. The no-slip and adiabatic conditions are assumed throughout the entire wall boundary. As a specific example, the inflow Mach number is assumed to be 3, and the temperature and pressure are 600 K and 0.1 MPa, respectively. Gaseous hydrogen at a temperature of 151.5 K is injected normal to the wall from a choked injector.A series of calculations were carried out by varying the fuel injection pressure from 0.5 to 1.5MPa. This amounts to changing the fuel mass flow rate or the overall equivalence ratio for different operating regimes. Figure 1 shows the instantaneous temperature fields in the supersonic combustor at four different conditions. The dark blue region represents the hot burned gases. At the fuel injection pressure of 0.5 MPa, the flame is stably anchored, but the flow field exhibits a high-amplitude oscillation. At the fuel injection pressure of 1.0 MPa, the Mach reflection occurs ahead of the injector. The interaction between the incoming air and the injection flow becomes much more complex, and the fuel/air mixing is strongly enhanced. The Mach reflection oscillates and results in a strong fluctuation in the combustor wall pressure. At the fuel injection pressure of 1.5MPa, the flow inside the combustor becomes nearly choked and the Mach reflection is displaced forward. The leading shock wave moves slowly toward the inlet, and eventually causes the combustor-upstart due to the thermal choking. The cavity appears to play a secondary role in driving the flow unsteadiness, in spite of its influence on the fuel/air mixing and flame evolution. Further investigation is necessary on this issue. The present study features detailed resolution of the flow and flame dynamics in the combustor, which was not typically available in most of the previous works. In particular, the oscillatory flow characteristics are captured at a scale sufficient to identify the underlying physical mechanisms. Much of the flow unsteadiness is not related to the cavity, but rather to the intrinsic unsteadiness in the flowfield, as also shown experimentally by Ben-Yakar et al. [6], The interactions between the unsteady flow and flame evolution may cause a large excursion of flow oscillation. The work appears to be the first of its kind in the numerical study of combustion oscillations in a supersonic combustor, although a similar phenomenon was previously reported experimentally. A more comprehensive discussion will be given in the final paper presented at the colloquium.

• Aerospace Engineering and Technology
• /
• v.4 no.2
• /
• pp.27-35
• /
• 2005

In this study, a performance analysis code was developed for the off-design performance analysis of air turbo ramjet(ATR) engine, and the analyses were performed for the pre-designed ATR engine at several operating points in the envelope. Variable intake and thrust nozzle were assumed to cover the wide envelope. Mathematical models for each components were developed to calculate their off-design performance. Simple design formulas were introduced for some components to explore the performance variation versus the design parameters. As a result, the pre-defined engine couldn't cover the entire mission profile. And it was also found that the effect of the pre-cooler was not very great, especially in the region of low Mach number.

• The KSFM Journal of Fluid Machinery
• /
• v.11 no.5
• /
• pp.22-29
• /
• 2008

This study is performed to understand the effect of the variation in the passage area of a two-dimensional impeller on its performance characteristics. We observe the results with changing the area ratio of inlet to outlet about $1{\sim}2.8$. A comparison between the experimental and numerical results was performed for the same configuration in order to verify the reliability of the CFD code. Overall characteristics in the passages of impeller were analyzed in detail including streamline, Mach number, pressure and polytropic efficiency distribution. When the passage area ratio exceeds 2, the pressure ratio is high. An area ratio of 2.3 showed the highest efficiency. The results will be used as useful reference data to establish the design concept of two-dimensional impeller and to improve its performance.

• International Journal of Aeronautical and Space Sciences
• /
• v.16 no.2
• /
• pp.123-136
• /
• 2015

The present study numerically investigates the effect of shield on the flow characteristics of Hartmann whistle. The flow characteristics of un-shielded Hartmann whistle are compared with whistles of different shield heights 15 mm, 17 mm, 20 mm, 25 mm and 30 mm. The comparison of Mach number contours and transient velocity vectors of shielded Hartmann whistles with un-shielded ones for the same conditions reveal that the presence of shield causes the exiting jet to stick to the wall of the shield without causing spill-over around the cavity inlet, thus sustaining the shock oscillation as seen in the unshielded Hartmann whistle, which has intense flow/shock oscillation and spill-over around the cavity mouth. The velocity vectors indicate jet regurgitance in shielded whistles showing inflow and outflow phases like un-shielded ones with different regurgitant phases. The sinusoidal variation of mass flow rate at the cavity inlet in un-shielded Hartmann whistle indicates jet regurgitance as the primary operating mode with large flow diversion around the cavity mouth whereas the non-sinusoidal behavior in shielded ones represent that the jet regurgitance is not the dominant operating mode. Thus, this paper sufficiently demonstrates the effect of shield in modifying the flow/shock oscillations in the vicinity of the cavity mouth.

• International Journal of Aeronautical and Space Sciences
• /
• v.16 no.2
• /
• pp.295-310
• /
• 2015

In the present study, an unstructured mixed mesh flow solver was used to conduct a numerical prediction of the aerodynamic performance of the S-76 rotor in hover. For the present mixed mesh methodology, the near-body flow domain was modeled by using body-fitted prismatic/tetrahedral cells while Cartesian mesh cells were filled in the off-body region. A high-order accurate weighted essentially non-oscillatory (WENO) scheme was employed to better resolve the flow characteristics in the off-body flow region. An overset mesh technique was adopted to transfer the flow variables between the two different mesh regions, and computations were carried out for three different blade configurations including swept-taper, rectangular, and swept-taper-anhedral tip shapes. The results of the simulation were compared against experimental data, and the computations were also made to investigate the effect of the blade tip Mach number. The detailed flow characteristics were also examined, including the tip-vortex trajectory, vortex core size, and first-passing tip vortex position that depended on the tip shape.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.20 no.3
• /
• pp.103-110
• /
• 2016

A performance analysis program of a disk type fluidic valve was developed to predict a chamber pressure and a response time. A parametric study of this device was performed by using scattering plot method. A sensitivity of Mach number at a nozzle outlet showed the highest value about a outlet diameter of nozzle. An inlet flow rate is the most important parameter to design the fluidic valve because it has high sensitivity value both a outlet velocity and a response time.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.19 no.6
• /
• pp.81-90
• /
• 2015

In order to know the characteristics of the hypersonic air-breahting engine, high altitude and Mach number ground test is necessary. Therefore, high pressure and high temperature condition should be simulated to do ground test of the hypersonic air-breathing engine. In this paper, the hypersonic air-breathing engine ground test facility of the Korea Aerospace Research Institute was introduced and the composition and characteristics were described.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.8
• /
• pp.1122-1133
• /
• 2003

Numerical characteristics of implicit upwind schemes, such as upwind ADI, line Gauss-Seidel (LGS) and point Gauss-Seidel (LU) algorithms, for Navier-Stokes equations have been investigated. Time-derivative preconditioning method was applied for efficient convergence at low Mach/Reynolds number regime as well as at large grid aspect ratios. All the algorithms were expressed in approximate factorization form and von Neumann stability analysis was performed to identify stability characteristics of the above algorithms in the presence of high grid aspect ratios. Stability analysis showed that for high aspect ratio computations, the ADI and LGS algorithms showed efficient damping effect up to moderate aspect ratio if we adopt viscous preconditioning based on min-CFL/max-VNN time-step definition. The LU algorithm, on the other hand, showed serious deterioration in stability characteristics as the grid aspect ratio increases. Computations for several practical applications also verified these results.

• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.271-275
• /
• 2001

Three dimensional flow characteristics in a liquid fuel ramjet combustor are investigated using PIV method. The combustors have two rectangular inlets that form 90 degree each other. Three guide vane is installed in each rectangular inlet to improve the flow stability. We made three cases of test combustors in which those inlet angles are 30 degree, 45 degree and 60 degree. Each combustor easily changes the size of combustor's recirculation zone with the replacement of combustors dome. The experiments are performed in the water tunnel test with the same Reynolds number in the case of Mach 0.3 at inlet. PIV software is developed to measure the flow field in the combustor and the accuracy of developed PIV program is verified with rotating disk experiment and standard data. The experimental results show that the two main streams from rectangular inlet collide near the plane of symmetry and generate two large longitudinal vortex, A large and complex three-dimensional recirculating flow is measured in the recirculation zone.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2000.04a
• /
• pp.8-8
• /
• 2000

In this paper are presented main power and gasdynamic characteristics of C-l6VK hypersonic test cell of Research Test Center of CIAM. Gasdynamic adjustment of the C-l6VK test cell was carried out with the working section constructed on scheme of Ramjet/scramjet test in free stream. Gasdynamic adjustment was conducted stage by stage in tile following sequence. First, check and preparation of all technical systems and checking measuring system. Than determination of the characteristics of test cell on cold (without the heating of air at entrance) regime and determination of the characteristics of test cell on regimes with the heating of air. Finally determination of tile characteristics of test cell with the loading of the working part by object. On tile final stage of gasdynamic adjustment two experiments with tile axisymmetric Scramjet model loaded into the working part of test cell were conducted. The first experiment was conducted with the purpose of determination of flow parameters with the object leaded into the working part and verification of experiment cyclogram. The second experiment was conducted with injection of hydrogen into the combustion chamber of object, that is tile conditions on test cell simulated Scramjet flight Mach number M = 6. Such methodology of gasdynamic adjustment allows to determine influence of experimental object on flow parameters in the working part at different conditions of experiment (with the burning in combustion chamber of object and without the homing), and also to compare flow characteristics in the object duct.