Sub-critical stable hydrogen-air premixed laminar flames in micro gaps

Sub-critical burning of lean hydrogen-air mixtures in micro gaps between two quartz disks was investigated both experimentally and numerically. Stationary regimes for different compositions and gap sizes were found when sub-critical flames remained in a stable position relative to the disk surfaces. The burning velocity in the micro gaps was observed to reach values much larger than the laminar burning velocity. A reaction-diffusion numerical model was proposed to corroborate experimental results. Different factors, such as boundary conditions for velocity, irradiation of the disk surfaces contacting the gas, and an increase in the chemical reaction rate near disk surfaces were modeled numerically in order to explain the increase in burning velocities. The best correlation between the numerical results and experimental data was observed in the scenario proposing as increased chemical reaction rate near the disk surfaces. Numerical simulations also showed that for large flame front velocities and wider sub-critical gaps, the flame front becomes unstable. The reason for this instability is the asynchronization of the combustion near the disk surfaces and the subsequent turbulization of the flame.     

Assessment of closure schemes in second-order conditional moment closure against DNS with extinction and ignition

The performance of second-order conditional moment closure (CMC) depends on models to evaluate conditional variances and covariances of temperature and species mass fractions. In this paper the closure schemes based on the steady laminar flamelet model (SLFM) are validated against direct numerical simulation (DNS) involving extinction and ignition. Scaling is performed to reproduce proper absolute magnitudes, irrespective of the origin of mismatch between local flamelet structures and scalar dissipation rates. DNS based on the pseudospectral method is carried out to study hydrogen-air combustion with a detailed kinetic mechanism, in homogeneous, isotropic, and decaying turbulent media. Lewis numbers are set equal to unity to avoid complication of differential diffusion. The SLFM-based closures for correlations among fluctuations of reaction rate, scalar dissipation rate, and species mass fractions show good comparison with DNS. The variance parameter in lognormal PDF and the constants in the dissipation term have been estimated from DNS results. Comparison is made for the resulting conditional profiles from DNS, first-order CMC, and second-order CMC with correction to the most critical reaction step according to sensitivity analysis. Overall good agreement ensures validity of the SLFM-based closures for modeling conditional variances and covariances in second-order CMC.     

Numerical investigation of mass flow rate effects on multiplicity of detonation waves within a H2/Air rotating detonation combustor

This paper presents results from numerical simulations of a non-premixed hydrogen-air rotating detonation combustor with radial injection. The fuel and air mass flow rate are varied in order to hold a unity global equivalence ratio. The calculations show that multiple detonation waves co-exist when the mass flow rate is increased. Conditional statistics of the detonation structure and combustion processes suggest similarities across co-existing waves. Quantification of the injection response to the rotation of a detonation indicates that at higher flow rate the refill time is short enough to allow for a quick and well mixed composition prior to the new front passage. Details of the combustion characteristics are analyzed. The results elucidate the correlation between initial injection conditions and detonation multiplicity on the overall physics within the combustor.     

Numerical simulation of detonation wave propagation through a rigid permeable barrier

The results of calculation of the detonation propagation in a porous medium for hydrogen-air mixture are presented. The porous medium was specified explicitly and consisted of sets of individual obstacles in the form of solid walls or the sets of finite-size plates. Various modes of detonation propagation depending on obstacle parameters are obtained: propagation in a cellular mode, stationary propagation with destruction of the cellular structure of the detonation front, propagation of a monotonically attenuating detonation wave with destruction of the cellular structure of the front. The possibility of reducing the detonation propagation velocity by replacing solid plates with finite-size ones was shown. The effect of the geometrical parameters of the plates and the step of it installation on the degree of detonation attenuation was estimated. It was determined that an increase in the number of plates leads to a stronger attenuation of the detonation.     

Premixed hydrogen-air flame front dynamics in channels with central and peripheral ignition

Experimental and analytical study of burning hydrogen-air mixtures with 12, 13, and 15 vol% hydrogen concentrations in channels with central and peripheral ignition was performed. Flame propagation speeds were determined by shadow and infrared high-speed imaging in the transverse and longitudinal directions, respectively. It was found that the increase in the flame front speed during the peripheral ignition reaches up to 1.7 times compared to the central ignition depending on mixture content. The pressure growth rate was examined in a closed channel. It was estimated that the time to reach a maximum pressure is 1.1 times less in the case or peripheral ignition than the central one. An analytical model was formed to describe the dynamics of the flame front in both cases. The model of a “reversed finger-flame” generated by a peripheral ignition was presented. The obtained results could be used in designing hydrogen-fueled combustible engines with the reduced knock-effect.     

Busemann diffuser for supersonic ramjet engine with detonation combustion of hydrogen-air mixture

The expediency of using a Busemann diffuser in a scramjet engine with detonation combustion of hydrogen-air mixtures is studied. A method is developed for calculating an air intake that provides spontaneous initiation of oblique detonation in a hydrogen-air mixture under conditions of rarefied atmosphere at hypersonic freestream Mach numbers. The geometry of the Laval nozzle (convergent-expanding) with the Busemann diffuser for the free flow Mach number 9 is determined. It is shown that at an altitude of about 40 km, the detonation combustion of the hydrogen-air mixture in this nozzle can provide more than 0.4 tons of thrust and more than 35% efficiency. The mathematical model is based on unsteady two-dimensional Euler equations for axisymmetric multicomponent reacting gas flow. To simulate chemical transformations, the detailed kinetic scheme is used that takes into account 33 non-equilibrium reactions for nine components of the mixture. The heat capacity, enthalpy, and entropy of a mixture are calculated using the reduced Gibbs energy of the gas components. Numerical modeling is carried out according to the modified Godunov scheme of the second order of approximation in spatial variables.     

Numerical study on the spark ignition characteristics of hydrogen-air mixture using detailed chemical kinetics

Hydrogen is a promising fuel and is expected to replace hydrocarbon fuels for its significant potentials to reduce the pollutants and greenhouse gases. It is very important to investigate Minimum ignition energy (MIE) on safety standards and ignition process of hydrogen-air mixtures. Even though the formation of flame kernels in quiescent hydrogen-air mixtures has been researched numerically and experimentally, the details of ignition mechanism have never been satisfactorily explained. In this study, the spark ignition of hydrogen-air mixture is investigated by using detailed chemical kinetics and considering the heat loss to the electrode. The purpose of this study is emphasized in the effects of the energy supply procedure, the radius of the spark channel, electrode size and electrode gap distance on the MIE. In addition, the effects of mixture temperature, electrode gap distance and electrode size on relationship between the equivalence ratio and the MIE are examined.     

Experimental study on induced accelerated combustion of premixed hydrogen-air in a confined environment

The study on induced accelerated combustion of premixed hydrogen-air in a confined environment is of great significance for the efficient utilization of hydrogen energy in internal combustion engines. The accelerated flame induced by the orifice plate is more stable and easy to control, which is beneficial to achieve controlled and rapid turbulent combustion. In this work, the accelerated combustion process induced by the orifice plate, and the influence of the orifice structure and initial conditions on the flame propagation and combustion characteristics were investigated by constant volume combustion bomb and schlieren method. The results show that the combustion process induced by the orifice plate consists of three stages: the initial stage of propagation, the accelerated stage of the orifice plate, and the end combustion stage. The reduction in aperture induces greater turbulence intensity and increases the perturbation of the orifice plate to the flame, resulting in a substantial increase in flame propagation speed through the orifice plate. As the initial pressure and the equivalence ratio increase, the velocity of turbulent flame induced by the orifice plate and the change rate of the velocity before and after the orifice plate increase. As the initial temperature increases, the turbulent flame propagation velocity does not change much, and the velocity change rate before and after the orifice plate decreases. The effect of the initial conditions on flame acceleration induced by the orifice plate is essentially the influence of flame propagation speed and instability. The greater the flame propagation speed and the stronger the flame instability, the stronger the induced turbulence and the greater the influence of the turbulent flow disturbance, and the greater the velocity of the turbulent flame induced by the orifice plate. There exists an optimum aperture for the shortest combustion duration at any initial conditions, but the optimal diameter is not sensitive to changes in initial conditions. The effect of orifice-induced combustion acceleration is remarkable, and the combustion durations induced by each orifice plate are shortened by more than 50%.     

Propagation process of H2/air rotating detonation wave and influence factors in plane-radial structure

The rotating detonation wave (RDW) propagation processes and influence factors are simulated in the plane-radial structure. The effects of inner radii of curvature, domain widths and stagnation pressures on propagation mode are studied. The RDW is initiated, and two kinds of propagation mode are obtained and analyzed. The flow field structure, parameters variation and influence factors on unstable propagation mode are explored in depth, and the geometrical and injection conditions of the unstable propagation are obtained. Results indicate that the decoupling and re-initiation occur repeatedly during the unstable propagation mode of the RDW, and the angular velocities of leading shock wave vary accordingly. When the domain width remains constant, the range of stagnation-pressure under unstable propagation mode increases as the inner radius increases. But the RDW propagates steadily when the inner radius increases to a certain value (Larger than 40 mm in this study). The effect of curvature radius and initial pressure ahead of detonation wave on the unstable propagation mode in this calculation model is similar to that in a curved channel. When r_i +0.464p_a > 80.932 or r_i ≥ 40 mm, the detonation wave can propagate steadily in the annular domain. When the curvature radius remains constant, the stagnation-pressure range of the unstable propagation mode decreases as the domain width increases.