• 한국연소학회:학술대회논문집
• /
• 2013.06a
• /
• pp.95-96
• /
• 2013

Numerical simulations are conducted to investigate the mechanism of spontaneous ignition of hydrogen within a certain length of downstream pipe released by the failure of pressure boundaries of various geometric assumption. The results show that local ignition is developed in limited area such as boundary layer and the mixing of hydrogen and air is weak at the planar pressure boundary conditions, whereas the flame fronts at the contact region are developed at the pressure boundaries of the spherical shape.

• 한국전산유체공학회:학술대회논문집
• /
• 2001.05a
• /
• pp.189-194
• /
• 2001

In this paper, non-reacting and reacting flowfields were computed using a preconditioned Navier-Stokes solver. The preconditioning technique of Merkle et al. and TVD scheme or Chakravarthy and Osher was employed and the results obtained using developed code have a good agreement with the previous results and experimental data. The preconditioned Wavier-Stokes equation set with low Reynolds number $\kappa-\epsilon$ equation and species continuity equations, are discretized with strongly implicit manner and time integrated with LU-SSOR scheme. For the purpose of treating unsteady problem the duel-time stepping scheme was employed. For the validation of the code in incompressible flow regime, steady driven square cavity flow was considered and calculation result shows reasonably good agreement with the result of incompressible code. Shock wave/boundary layer interaction problem was considered to show the shock capturing performance of preconditioned-TVD scheme. To validate unsteady flow, acoustic oscillation problem was calculated, and supersonic premix flame of $H_2$-air reaction problem which is calculated with turbulence model, 9-species/18-reaction step reaction model, shows reasonable agreement with the previous results. As a result, the preconditioning method has an advantage to calculate incompressible and compressible flow through one code and preconditioned solver easily developed from standard compressible code with minor efforts. But additional computational time and computer memory is required due to preconditioning matrix.

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

Various jet engines (Turbine engine family and RAM Jet engine) have been developed for high speed aircrafts. but their application to hypersonic flight is restricted by principle problems such as increase of total pressure loss and thermal stress. Therefore, the development of next generation propulsion system for hypersonic aircraft is a very important subject in the aerospace engineering field, SCRAM Jet engine based on a key technology, Supersonic Combustion. is supposed as the best choice for the hypersonic flight. Since Supersonic Combustion requires both rapid ignition and stable flame holding within supersonic air stream, much attention have to be given on the mixing state between air stream and fuel flow. However. the wider diffusion of fuel is expected with less total pressure loss in the supersonic air stream. So. in this study the direction of fuel injection is inclined 30 degree to downstream and the total pressure of jet is controlled for lower penetration height than thickness of boundary layer. Under these flow configuration both streams, fuel and supersonic air stream, would not mix enough. To spread fuel wider into supersonic air an aerodynamic force, baroclinic torque, is adopted. Baroclinic torque is generated by a spatial misalignment between pressure gradient (shock wave plane) and density gradient (mixing layer). A wedge is installed in downstream of injector orifice to induce an oblique shock. The schlieren optical visualization from side transparent wall and the total pressure measurement at exit cross section of combustor estimate how mixing is enhanced by the incidence of shock wave into supersonic boundary layer composed by fuel and air. In this study non-combustionable helium gas is injected with total pressure 0.66㎫ instead of flammable fuel to clarify mixing process. Mach number 1.8. total pressure O.5㎫, total temperature 288K are set up for supersonic air stream.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.48 no.2
• /
• pp.135-146
• /
• 2020

To define an interface in a conventional level-set method, an analytical function must be revealed for an interfacial geometry. However, it is not always possible to define a functional form of level sets when interfaces become complex in a Cartesian coordinate system. To overcome this difficulty, we have developed a new level-set formalism that discriminates the interior from the exterior of a CAD modeled interface by parameterizing the stereolithography (STL) file format. The work outlined here confirms the accuracy and scalability of the hydrodynamic reactive solver that utilizes the new level set concept through a series of tests. In particular, the complex interaction between shock and geometrical confinements towards deflagration-to-detonation transition is numerically investigated. Also, a process of flame spreading and damages caused by point source detonation in a real-sized plant facility have been simulated to confirm the validity of the new method built for reactive hydrodynamic simulation in any complex three-dimensional geometries.

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