Pulsating one-dimensional detonation in ammonia-hydrogen–air mixtures

In this study, one-dimensional detonations in ammonia/hydrogen-air mixtures are numerically investigated by solving the fully compressible Navier-Stokes equations with detailed chemistry. Pulsating instabilities with single-mode are observed during the detonation wave propagation, accompanied by periodic coupling and decoupling of the lead shock wave and the reactive front. The ratio between driver pressure and initial pressure determines the overdrive degree and thus the oscillatory mode of detonation for a premixture with certain composition. The effects of hydrogen dilution and mixture equivalence ratio on pulsating detonations are also examined under a constant driver pressure. The growing hydrogen fraction in fuel blends significantly increases the oscillation frequency. In addition, the pulsating detonation frequency rises with increasing equivalence ratio under fuel-lean conditions, peaks under stoichiometric conditions, and falls under fuel-rich conditions as the equivalence ratio increases further. Evolutions of reactants, main intermediate radicals, and products are analysed in both fuel-lean and fuel-rich conditions. A chemical explosive mode analysis further confirms the highly-autoignitive nature of the mixture in the induction zone between reaction front and shock front where thermal diffusion plays a negligible role.     

Numerical study on the characteristics of rotating detonation wave with multicomponent mixtures

In order to investigate the effects of gas mixture components on the combustion characteristics of rotating detonation wave, two-dimensional simulation is presented to simulate the propagation process of rotating detonation wave with different methane conversions. The results indicate that there are five propagation modes of rotating detonation wave with different components: single-wave mode, single wave with counter-rotating components mode, double-waves mode, triple-waves mode and quadruple-waves mode. The detonation wave propagates along the forward direction in all five modes. With the increase of methane conversion, multi-wave mode appears in the combustion chamber. The fuel component has a great influence on the heat release ratio of detonation combustion. The velocity of detonation wave decreases with the increase of methane conversion. With the increase of methane conversion, the chemical reaction rate gradually increases, which leads to the intensification of chemical reaction on the deflagration surface. The reaction on the deflagration surface develops to the unburned fuel zone, which eventually leads to the formation of compression waves and shock waves in the fuel refill zone. When the shock wave sweeps through the fresh premixed gas, the reactant is compressed to form a detonation point and then ignite the fuel. A new detonation wave is finally formed. The total pressure ratio decreases with the increasing methane conversion, and the uniformity of the total pressure of outlet decreases with increasing methane conversion.     

Effects of concentration gradient on pulsating and cellular detonations

In the field of explosion accidents and the combustion chamber of propulsion systems, gaseous mixtures are more likely to be highly non-uniform and a detonation usually propagates in a non-homogeneous medium. In this paper, one-dimensional (1D) pulsating detonations propagating in non-homogeneous medium with concentration gradient are studied by direct numerical simulation with a detailed chemistry. It is found that a periodically pulsating detonation can adjust to adapt to the environment with concentration gradient. The superficial concentration gradient can improve survivability of 1D detonation under highly unstable case, widening detonability limit corresponding to equivalence ratio; while the steep concentration gradient will quickly extinguish the detonation. The extend of detonability limit led by the inhomogeneity is reinforced by the presence of transverse wave in 2D detonations. The 2D detonation in the gradient has better viability compared with the 1D due to the coupling interaction of leading shock and transverse wave with the gradient and the transportation of H₂ and H along the gradient. This study also demonstrates that the simulations of detonations in the gradient requires detailed chemical mechanism.     

双点激光点火直接起爆的数值模拟

针对双点激光点火直接起爆过程中爆轰波的形成、发展和传播问题,采用高精度数值模拟方法求解带化学反应的二维欧拉方程组,研究了不同环境压力情况对流场结构与波系变化的影响.结果表明,环境压力会影响激波强度与爆轰波的传播速度,是决定双点激光点火形成的火核在碰撞过程中能否实现爆轰并维持爆轰波传播的重要因素,利用双激光点相互作用形成...     

Simulations of detonation wave propagation in rectangular ducts using a three-dimensional WENO scheme

This paper reports high resolution simulations using a fifth-order weighted essentially non-oscillatory (WENO) scheme with a third-order TVD Runge-Kutta time stepping method to examine the features of detonation front and physics in square ducts. The simulations suggest that two and three-dimensional detonation wave front formations are greatly enhanced by the presence of transverse waves. The motion of transverse waves generates triple points (zones of high pressure and large velocity coupled together), which cause the detonation front to become locally overdriven and thus form “hot spots.” The transversal motion of these hot spots maintains the detonation to continuously occur along the whole front in two and three dimensions. The present simulations indicate that the influence of the transverse waves on detonation is more profound in three dimensions and the pattern of quasi-steady detonation fronts also depends on the duct size. For a “narrow” duct (4L×4L where L is the half-reaction length), the detonation front displays a distinctive “spinning” motion about the axial direction with a well-defined period. For a wider duct (20L×20L), the detonation front exhibits a “rectangular mode” periodically, with the front displaying “convex” and “concave” shapes one following the other and the transverse waves on the four walls being partly out-of-phase with each other.     

Experimental and numerical studies on detonation failure and re-initiation behind a half-cylinder

Mach reflection causes the re-initiation of decoupled detonation owing to changes in the boundary. A complementary series of experiments and numerical simulations, illustrating detonation failure and subsequent reinitiation processes, have been presented. Immediately across the half-cylinder, the decoupled detonation owing to the diffraction effect wave is reflected on the bottom wall to form a regular reflection, and then changes into the Mach reflection, which further determines the detonation reinitiation. Two different reinitiation modes after detonation wave diffraction were observed for the stable mixtures: the direct Mach reflection re-initiation mode and Mach reflection combined with the transverse detonation. However, for unstable detonations, a different reinitiation mode was obtained, whereby the development of intrinsic instabilities resonates with the reflection on the bottom wall, rendering the Mach reflection randomly occurring or even absent. The critical limit of detonation failure is characterized by the radius of the half-cylinder and the cell size. In addition, the transition length from regular to Mach reflection was measured to reveal the length-scale effect on the process.     

The initiation and propagation of detonation in supersonic combustible flow with boundary layer

Adaptive simulations solving the Navier-Stokes equations have been conducted in order to get a better understanding on the detonation initiation and propagation in a stoichiometric H₂/O₂/Ar supersonic mixture with boundary layer. The detonation is initiated by a continuous hot jet. When reflecting on the wall, the jet induced bow shock interacts with the boundary layer and forms the shock boundary layer interaction phenomena, while in Euler result the bow shock forms Mach reflection. The investigation shows that the Navier-Stokes simulation result is structurally in better agreement with the experiment compared with that of the inviscid Euler simulation result. The bow shock interacts with the separation shock, forming the shock induced combustion behind the interaction zone. Then the combustion front couples with shock and forms Mach stem induced detonation. The Mach stem induced detonation continues to getting higher and propagating upstream, initiating the main flow. The initiated partial detonation exists with the separation shock induced combustion front, forming an “oblique shock induced combustion-partial detonation” structure in the main flow. The investigation on the influence of free stream Mach number further confirms that the boundary layer has an important influence on detonation initiation. The parametric studies also show that there exists a free stream Mach number range to initiate the partial detonation in supersonic combustible flow successfully.     

Direct observations of reaction zone structure in propagating detonations

We report experimental observations of the reaction zone structure of self-sustaining, cellular detonations propagating near the Chapman-Jouguet state in hydrogen-oxygen-argon/nitrogen mixtures. Two-dimensional cross sections perpendicular to the propagation direction were imaged using the technique of planar laser induced fluorescence (PLIF) and, in some cases, compared to simultaneously acquired schlieren images. Images are obtained which clearly show the nature of the disturbances in an intermediate chemical species (OH) created by the variations in the strength of the leading shock front associated with the transverse wave instability of a propagating detonation. The images are compared to 2-D, unsteady simulations with a reduced model of the chemical reaction processes in the hydrogen-oxygen-argon system. We interpret the experimental and numerical images using simple models of the detonation front structure based on the “weak” version of the flow near the triple point or intersection of three shock waves, two of which make up the shock front and the third corresponding to the wave propagating transversely to the front. Both the unsteady simulations and the triple point calculations are consistent with the creation of keystone-shaped regions of low reactivity behind the incident shock near the end of the oscillation cycle within the “cell.”     

Self-organized generation of transverse waves in diverging cylindrical detonations

Self-organized generation of transverse waves associated with the transverse wave instabilities at a diverging cylindrical detonation front was numerically studied by solving two-dimensional Euler equations implemented with an improved two-step chemical kinetic model. After solution validation, four mechanisms of the transverse wave generation were identified from numerical simulations, and referred to as the concave front focusing, the kinked front evolution, the wrinkled front evolution and the transverse wave merging, respectively. The propagation of the cylindrical detonation is maintained by the growth of the transverse waves that match the rate of increase in surface area of the detonation front to asymptotically approach a constant average number of transverse waves per unit length along the circumference of the detonation front. This cell bifurcation phenomenon of cellular detonations is discussed in detail to gain better understanding on detonation physics.     

Deflagration to detonation transition in aluminum dust-air mixture under weak ignition condition

Deflagration to detonation transition (DDT) in flake aluminum dust-air mixture was studied in a 199 mm inner diameter and 29.6 m long horizontal tube. 40 sets of dust dispersion system were used to disperse flake aluminum into the experimental tube. An electric spark of 40 J was used to ignite the aluminum-air mixture. Self-sustained detonation was observed and the characteristics of deflagration to detonation transition process were studied. Contributed to the transverse wave and cellular structure of detonation wave front in aluminum-air mixture, the propagation velocity of detonation wave ranges from 1480 m/s to 1820 m/s and the maximum overpressure oscillates between 40 and 102 bar. Single head spinning detonation wave in aluminum dust-air mixture was observed and the cell size was evaluated.     

Effects of confinement on detonation in H2–O2 mixture with transverse concentration gradient

Since hydrogen has wide flammability limit and low ignition energy, it could be easily ignited and be easy for the transition to a detonation, leading to extremely serious impacts in explosion accidents or extremely high combustion effeciency in the propulsion. In the field of explosion accidents and the combustion chamber of propulsion systems, hydrogen mixtures are more likely to be highly non-uniform and a detonation usually propagates in a non-homogeneous medium. The work studies behaviors of detonations in non-homogenous medium by a high-resolution simulations. We widen computational domain and steepen the gradient to weaken the role of transverse wave on cellular detonations propagating in the medium with transverse concentration gradient and to reveal the interaction of longitudinal shock with reaction wave. The results show that characteristics of detonations in nonuniform medium is controlled by coupling role of gradient and confinement. As a domain is sufficient wide, the reflection wave is rather weak so that the detonation takes on galloping propagation, with a single-head mode. As the width is mediate, detonation cell takes on highly irregularity, similar to that of highly unstable detonation. However, in the narrow domain, steepening gradient plays a key role while confinement becomes minor in detonation propagation.     

Numerical simulation on the operating characteristics of rotating detonation ramjet engines at high flight mach number

For high-Mach-number incoming flow circumstances, a rotating detonation ramjet engine configuration is proposed in this research. By installing supporting blocks at the rear of the combustor, this configuration achieves continuous rotating detonation operation. Based on the Comparison of the flow structures obtained from the engine configuration with and without the supporting block before and after detonation ignition respectively, we obtain the intrinsic mechanism of detonation wave's propagation and re-initiation under the action of the supporting block. The supporting block creates a deflagration wave that is almost stationary before detonation ignition. In the detonation-ignited state, the deflagration wave is continually formed and traveling upstream under the influence of the supporting block, which is analogous to the periodical before detonation ignition of a transverse wave structure. The dynamic deflagration wave will cause the incomplete reactants behind the detonation wave to burn as the intensity of the detonation wave decreases. As a result, the incident shock wave is transformed into a Mach stem to achieve the re-initiation of the detonation wave.     

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.     

Knocking combustion in spark-ignition engines

Knocking combustion research is crucially important because it determines engine durability, fuel consumption, and power density, as well as noise and emission performance. Current spark ignition (SI) engines suffer from both conventional knock and super-knock. Conventional knock limits raising the compression ratio to improve thermal efficiency due to end-gas auto-ignition, while super-knock limits the desired boost to improve the power density of modern gasoline engines due to detonation. Conventional combustion has been widely studied for many years. Although the basic characteristics are clear, the correlation between the knock index and fuel chemistry, pressure oscillations and heat transfer, and auto-ignition front propagation, are still in early stages of understanding. Super-knock combustion in highly boosted spark ignition engines with random pre-ignition events has been intensively studied in the past decade in both academia and industry. These works have mainly focused on the relationship between pre-ignition and super-knock, source analyses of pre-ignition, and the effects of oil/fuel properties on super-knock. The mechanism of super-knock has been recently revealed in rapid compression machines (RCM) under engine-like conditions. It was found that detonation can occur in modern internal combustion engines under high energy density conditions. Thermodynamic conditions and shock waves influence the combustion wave and detonation initiation modes. Three combustion wave modes in the end gas have been visualized as deflagration, sequential auto-ignition and detonation. The most frequently observed detonation initiation mode is shock wave reflection-induced detonation (SWRID). Compared to the effect of shock compression and negative temperature coefficient (NTC) combustion on ignition delay, shock wave reflection is the main cause of near-wall auto-ignition/detonation. Finally, suppression methods for conventional knock and super-knock in SI engines are reviewed, including use of exhaust gas recirculation (EGR), the injection strategy, and the integration of a high tumble - high EGR-Atkinson/Miller cycle. This paper provides deep insights into the processes occurring during knocking combustion in spark ignition engines. Furthermore, knock control strategies and combustion wave modes are summarized, and future research directions, such as turbulence-shock-reaction interaction theory, detonation suppression and utilization, and super-knock solutions, are also discussed.     

Effects of a turbine guide vane on hydrogen-air rotating detonation wave propagation characteristics

The rotating detonation engine is a new machine that can generate thrust via continuous rotating detonation waves (RDWs). In this study, experiments were performed on a structure combining a rotating detonation combustor (RDC) and a turbine guide vane to investigate the propagation characteristic of hydrogen-air RDW. The results showed that the velocity of detonation wave initially increased and then decreased with the increase of equivalence ratio, and it got a velocity of 84% Chapman-Jouguet value. The velocity of detonation wave generally rose by 4.31% comparing with the no guide vane tests, while the scope of steady-operation state became narrow. The oscillation pressure was reduced by 64% after passing through the guide vane, and the magnitude of pressure was only 0.4 bar at the guide vane exit. Meanwhile, part of the shock wave was reflected back to combustor resulting in some small pressure disturbances, and the propagation mode of reflected wave was related to the propagation direction of RDW.     

Hydrogen-fueled detonation ramjet model: Wind tunnel tests at approach air stream Mach number 5.7 and stagnation temperature 1500 K

The mode of continuous spinning detonation (CSD) combustion of hydrogen in the annular combustor of a model of a hydrogen-fueled detonation ramjet under conditions of approach air stream Mach number 5.7 and stagnation temperature 1500 K is registered experimentally in a short-duration (pulsed) wind tunnel at the overall air-to-hydrogen equivalence ratio (ER) ranging from 0.7 to 1.4. The maximum values of thrust and specific impulse of the ramjet model are attained at ER = 1.25 and are estimated as 1550 N and 3300 s, respectively. At 1.4 < ER < 1.6, the mode of longitudinally pulsating detonation (LPD) combustion is registered with somewhat lower values of thrust and specific impulse.     

Study of thrust vector control for the rotating detonation model engine

The detonation wave in a rotating detonation engine is highly adaptable to the incoming flow, making the wave easier to control. In this study, a numerical simulation method is used to analyze the working process and flow field structure of a rotating detonation model engine with dual cavity injection of an H₂/air mixture by controlling the injection pressure ratio of the dual cavity and the number of detonation wave heads. It is found that the rotating detonation engine offers the possibility to control the thrust vector with two different modes. The first is a one-cycle alternate control mode with a small injection pressure ratio. Here two deflections occur in different directions occur across one detonation wave propagation cycle, but the overall deflection direction is in the low-pressure region. The second is a one-way control mode, with a large injection pressure ratio, and the deflection direction towards the low-pressure region. For the multi wave-mode, it belongs to one-way control mode because of constant deflection direction in the low-pressure area. From the perspective of thrust distribution along the circumference, the one-way control strategy satisfies the ability of a rotating detonation thrust vector control.     

Stable detonation wave propagation in rectangular-cross-section curved channels

The detonation propagation phenomena in curved channels were experimentally studied in order to determine the stable propagation condition. A stoichiometric ethylene–oxygen mixture gas and five types of rectangular-cross-section curved channels with different inner radii of curvature were employed. The detonation waves propagating through the curved channels were visualized using a high-speed video camera. Multi-frame short-time open-shutter photography (MSOP) was developed in the present study to simultaneously observe the front shock shape of the detonation wave and the trajectories of triple points on the detonation wave. The detonation wave became more stable under the conditions of a higher filling pressure of the mixture gas and/or a larger inner radius of curvature of the curved channel. The critical condition under which the propagation mode of the detonation wave transitioned from unstable to stable was having an inner radius of curvature of the curved channel (r_i) equivalent to 21–32 times the normal detonation cell width (λ). In the stable propagation mode, the normal detonation velocity (D_n) increased with the distance from the inner wall of the curved channel and approached the velocity of the planar detonation propagating through the straight section of the curved channel (D_str). The smallest D_n was observed on the inner wall and decreased with decreasing r_i/λ. The distribution of D_n on the detonation wave in the stable mode was approximately formulated. The approximated D_n given by the formula agreed well with the experimental results. The front shock shape of the detonation wave could be reconstructed accurately using the formula. The local curvature of the detonation wave (κ) nondimensionalized by λ decreased with increasing distance from the inner wall. The largest λκ was observed on the inner wall and increased with increasing r_i/λ. D_n/D_str decreased with increasing λκ. This nondimensionalized D_n–κ relation was nearly independent of r_i/λ.     

Numerical study on transition of hydrogen/air flame triggered by auto-ignition under effect of pressure wave in an enclosed space

End gas auto-ignition and transition of flame front are considered as the main causes of severe pressure oscillation in spark-ignition engines, which is one of the major features of knock and super-knock. The knowledge of characteristics of auto-ignition, flame front development, propagation of pressure wave and relations between them, still needs to be maintained. In this study, flame front transition induced by pressure wave and auto-ignition are investigated using one-dimensional simulation with detailed chemistry in an enclosed space Calculation cases with different initial thermodynamic conditions are investigated. Mass fraction of OH is employed as indicator of auto-ignition progress under variable conditions caused by pressure wave. Different propagation modes of flame front, including subsonic deflagration, detonation and supersonic deflagration, are developed under the effects of both pressure wave and auto-ignition. Results show that mass fraction of OH could successfully reflect auto-ignition progress, thus indicating occurrence and sequence of auto-ignition at different locations. Transitions from deflagration to detonation and detonation to supersonic deflagration are found to be triggered by sequential auto-ignition with different gradient of auto-ignition progress ahead of flame front induced by pressure wave.     

Mitigation of H2/air gaseous detonation via utilization of PAN-based carbon fibre felt

This paper first studied the mitigation of the hydrogen/air mixture detonation wave by tuning the tube inner wall with the absorbing material of polyacrylonitrile (PAN)-based carbon fibre felt. The experimental tests were performed in a single-trial circular cross-section tube filled with premixed hydrogen and air detonative mixture. The pressure values and flame front propagation were measured by means of pressure transducers and photodiodes respectively. The attenuation regimes of detonation wave in walled tubes with different thicknesses and layouts of absorption material were compared. The PAN-based carbon fibre felt makes a significant attenuation on the detonation propagation. The decoupling of leading shock wave and flame front can be observed under the effect of this absorbing layer. The ultimate strength close to the tube end and propagation velocity of the combustion wave decrease with the increase of felt thickness. When the interval layout felt is adopted, the spacing distance has almost no impact on the attenuation effectiveness. When the sectional layout is adopted, the effectiveness of detonation mitigation is however improved for a higher proportion of the absorbing material.