Numerical simulation of interaction of cellular detonation wave with systems of inert porous filters

Numerical modeling of the interaction of cellular detonation in a hydrogen-air mixture with several systems of porous filters covering part of the channel was carried out. The main regimes of detonation propagation and the critical conditions for detonation failure in the filter systems were obtained for each system. A map of detonation regimes was constructed, from which it follows that with an increase in the concentration of particles in the filters, it is possible to increase the distance between the filters or reduce the number of filters in the system necessary to successfully suppress detonation. A comparison of various filter systems in terms of blockage ratio and the total surface area of particles was made, from which it was concluded that the system of two filters was the most efficient to suppress detonation.     

Numerical study on the interaction characterization of rotating detonation wave and turbine rotor blades

Rotating detonation as a kind of pressure gain combustion is expected to greatly improve efficiency when applied to gas turbine engines. In this paper, the operation of rotating detonation combustor and turbine rotor blade was studied. Firstly, the analysis of the interaction between detonation wave and turbine blade shows that the compression of gas by detonation wave and reflected wave will lead to a sharp increase in the temperature at the wall of blade. When the detonation wave propagates, the oscillation amplitudes of pressure and temperature at the turbine inlet are 70% and 75% respectively, and the detonation oblique shock will change the flow trajectory of the air flow, resulting in the flow direction deviating from the incident angle. Then the comparison between detonation and deflagration shows that the total pressure of detonation is higher and will have greater work potential. The torque generated by the blades under detonation has the characteristics of high-frequency oscillation, which may be detrimental to the operation of the engine.     

Numerical research of rotating detonation initiation processes with different injection patterns

Three-dimensional numerical simulations are performed to study the initiation process of the rotating detonation combustor initiated by a pre-detonator. Navier-Stokes equations are solved with a 7-species and 8-steps H₂/Air chemical reaction mechanism. Three different injection patterns are considered, including the inner-slit injection, outer-slit injection and mid-slit injection. A single stable rotating detonation wave is established finally for all three injection patterns. Firstly, the initiator dynamics of the pre-detonator is discussed. The initial detonation wave spreads into the combustor from the pre-detonator and two oppositely propagating detonation waves are formed in the combustor in all cases, causing wave collisions. Similar fresh gas injection disturbance is seen, suggesting the little effect of different injection patterns at the beginning. After the recovery of the injection, ‘L-form’ fresh mixture layers are formed in the inner-slit/outer-slit injection patterns, leading to a rapid establishment of stable propagation mode. While in the mid-injection pattern, transition of fresh mixture layer structure is seen from ‘I-form’ to ‘L-form’. Severe lateral expansion exists during this period, resulting in a long-term wave initiation.     

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.     

Experimental study of a hydrogen-air rotating detonation engine with variable air-inlet slot

Rotating detonation engines have attracted considerable attentions in recent years. In this study, the experiments of initiating rotating detonation waves were performed on a H₂/air rotating detonation wave with the variable air-inlet slot. The results showed that the stability of detonation-wave pressure and velocity both initially increased and then decreased with the increase of slot width, and it could improve the stability of detonation-wave velocity via increasing the equivalence ratio. The intensity of reflected wave was strong for the tests of d = 0.5 mm, which leaded to the advance ignition of fresh mixture and a velocity deficit reaching up to 20%. The strong interaction between air plenum and combustor and bad mixing effect may be the reasons of forming unstable detonation wave for the tests of large-scale slots. The air-inlet slot of d = 1 mm, which got a best experiment results relative to other tests, had a wide equivalence-ratio scope to produce stable detonation wave.     

Numerical investigation on H2/Air non-premixed rotating detonation engine under different equivalence ratios

In this paper, three-dimensional numerical simulations are performed to investigate the formation and propagation characteristics of rotating detonation wave in a non-premixed engine. By changing the mass flow rate of H₂ and fixing air mass flow rate, the effects of equivalence ratio involving fuel lean and rich operating conditions are mainly discussed. Numerical results show that equivalence ratio plays a very critical role in the formation process and propagation mode, which further affects the propulsion performance of rotating detonation engine significantly. For current numerical geometry and operating conditions, the lean limit of equivalence ratio for formatting a stable RDW is about 0.4, dual-wave mode (at equivalence ratio of 0.6, 0.8, 1.0 and 1.4) and single-wave mode (at equivalence ratio of 1.2) are obtained, respectively. When equivalence ratio is 1.0, rotating detonation engine can exhibit excellent operating performance with the shortest formation time, best propagation stability, middling class thrust and specific impulse. Besides, the pressure contour analysis indicates that the effects of equivalence ratio and mass flow rate of H₂ on the collision strength and times during the re-initiation process are the main mechanisms for determining the formation possibility and propagation mode of rotating detonation wave. Besides, the intensity of accumulated pressure wave and distributions of equivalence ratio are two important factors for the generation of new detonation wave front. Furthermore, it is also detected from the comparisons of the propulsion performance that the effects of equivalence ratio on thrust and specific impulse under fuel lean conditions are more significant than those under fuel-rich conditions.     

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

Experimental investigation on propagation characteristics of liquid-fuel/preheated-air rotating detonation wave

The rotating detonation combustor can be applied to the turbine engine to develop into a new power device, and the liquid-fuel/air rotating denotation has important research significance for engine applications. In this research, the propagation characteristics of liquid-fuel/air rotating detonation wave were experimentally investigated. A hydrocarbon mixture—liquid gasoline was employed for the fuel, the oxidizer was high-temperature air preheated by a hydrogen-oxygen heater, and the rotating detonation wave was initiated via a hydrogen-oxygen pre-detonator. The effects of the equivalence ratio, ignition pressure, and air total temperature on the propagation characteristics of the liquid-fuel rotating detonation wave were analyzed. The liquid-fuel/air continuous rotating detonation wave can be successfully obtained with a single-wave mode, and the velocity and peak pressure of the rotating detonation waves increase as the equivalence ratio increases. As the detonation-wave pressures at the outlet of the pre-detonator increase, the establishment time of the rotating detonation wave gradually decreases, and the average establishment time is 4.01 ms. Stable rotating detonation waves are obtained with the air total temperature of 600–800 K, but the intensity of the detonation wave has a large deficit due to some instabilities.     

Numerical investigation of the effect of equivalence ratio on the propagation characteristics and performance of rotating detonation engine

The equivalence ratio is one of the significant factors affecting the propagation characteristics and performance of the rotating detonation engine (RDE). Using the compressible reacting flow solver based on the OpenFOAM open-source platform, the effect of different equivalence ratios of premixed H₂/Air gases on the propagation characteristics and performance of RDE under different total inlet pressures (P₀) is investigated. The reactants are injected through the discrete inlet to mimic the spatial inhomogeneity of the reactants in the actual RDE combustor. The results indicate a Y-shaped flow field structure is formed behind the rotating detonation wave (RDW) using the discrete inlet. There is only one RDW in the flow field with the change of the equivalence ratio when the P₀ is 0.5 MPa, and the primary factors affecting the RDW propagation velocity differ under fuel-lean and fuel-rich conditions. The RDW propagation mode switches from single-wave to co-directional double-wave and double-wave collision with the change of the equivalence ratio when the P₀ is 1.2 MPa. The velocity deficit of RDW in the double-wave mode is larger than that in the single-wave mode. Additionally, in the single-wave mode, the specific impulse decreases as the equivalence ratio increases, but the specific thrust increases as the equivalence ratio increases. When the propagation mode changes, the specific impulse and specific thrust show different trends with the change of the equivalence ratio. It demonstrates that multiple co-directional RDWs have the effect of stabilizing the thrust.     

Numerical investigation of the effect of inlet mass flow rates on H2/air non-premixed rotating detonation wave

In this paper, a three dimensional numerical investigation was carried out to study the formation and propagation characteristics of non-premixed rotating detonation wave using H₂/air as reactive mixtures. At a constant global equivalence ratio, the effects of inlet mass flow rates of H₂ and air on various performance parameters of rotating detonation wave and based on it combustor were analyzed in detail. On this basis, the mode switching process of rotating detonation wave caused by transiently changing the inlet mass flow rates was also discussed. The numerical results showed that inlet mass flow rates of H₂ and air played a very critical role in the formation, propagation and mode switching of rotating detonation wave. With the increase of inlet mass flow rates, rotating detonation wave could be switched from single wave to double waves. The propagation direction of double waves depended on the changing process of inlet mass flow rates. Meanwhile, compared to the single wave, double waves or its based combustor had the obvious advantages in formation time, stability and thrust, but had disadvantage in pressure ratio. In addition, both fill characteristics and mixing quality of fresh reactive mixtures are the underlying important mechanisms to explain the effects of inlet mass flow rates on rotating detonation waves.     

Numerical investigation of two-wave collision and wave structure evolution of rotating detonation engine with hollow combustor

A three-dimensional numerical simulation of rotating detonation engine (RDE) with hollow combustor is performed to analyze wave structure evolution systematically. Wave structure evolution is classified into four categories, namely two-wave collision (counter-rotating waves), abscission of detonation tail, and shock wave to detonation transition. Two-wave collision consists of symmetric detonation collision, asymmetric detonation collision, and detonation/shock collision. Two symmetric detonation waves turn into shock waves after collision. Collision of asymmetric detonation waves creates single detonation wave. The detonation/shock collision decreases the detonation wave intensity. Abscission of detonation tail and shock to detonation transition can both create single detonation wave or two opposite-direction detonation waves, depending on the wave hitting angle and the amount of fresh gas. All phenomena mentioned above affect the number of detonation waves in the combustion chamber.     

Effects of injection parameters on propagation patterns of hydrogen-fueled rotating detonation waves

Two-dimensional rotating detonation waves (RDWs) with separate injections of hydrogen and air are simulated using the Navier–Stokes equations together with a detailed chemical mechanism. The effects of injection stagnation temperature and slot width on the detonation propagation patterns are investigated. Results find that extremely high temperatures can lead to a chaotic mode in which detonation waves are generated and extinguished randomly. Increasing the slot width can reduce the number of detonation waves and finally trigger detonation quenching at a low injection stagnation temperature. But increasing the slot width can change the RDW propagation pattern from a chaotic to a stable mode under high injection temperature. Furthermore, the kinetic parameter τ (representing the chemical reactivity of the mixture) and the kinematic parameter α (representing the mixing efficiency of hydrogen and oxygen) are introduced to distinguish the RDW propagation patterns.     

Numerical simulation of radial-stratified rotating detonation flow field structures with different injection patterns

Rotating detonation engine has been widely studied in recent years because of its high theoretical efficiency and heat release rate. In many numerical simulations, the combustible mixture is injected and fully filled at the head of the combustor. In this paper, annular injection slits are proposed and three representative injection patterns are simulated by changing the injection directions. Stable single-wave modes are formed in all three patterns and two kinds of combustible mixture layer structures are found, namely “L-shape” and “T-shape” structures. Following the combustible mixture layer, the detonation wave is not fully filled in the radial direction, thus radial and circumferential shock waves are induced from the detonation wave, forming more complex wave structures. After the radial shock wave, velocity vortex and significant deflagration are found and propagate with the shock wave, thus maintaining a higher pressure and temperature there.     

Effects of boundary layer on wedge-induced oblique detonation structures in hydrogen-air mixtures

Oblique detonation wave (ODW) structures are studied widely in recent years, but most of them are solved by the Euler equations without considering viscosity and then effects of boundary layer. In this study, the Navier-Stokes Equations are used to simulate the wedge-induced ODWs in hydrogen-air mixtures, and the two types of ODW transition structures at different incident Mach number M_i are analyzed to clarify the effects of viscosity and hence the boundary layer. Results show that the effect of boundary layer on ODW structures should be classified by the types of ODW transition patterns. As for the smooth transition pattern of ODW at high Mach numbers, the effect of boundary layer can be neglected, but for the abrupt transition pattern of ODW at low Mach numbers, the effect of boundary layer is large and it changes the ODW structure greatly. Resulting from the interaction of shock and boundary layer, a recirculation zone is formed within the viscous ODW layer at M_i = 7, which leads to the phenomenon that the straight oblique shock wave evolves into two sections, with the downstream one having a larger shock angle. Additionally, the corresponding transition position moves upstream, and the initiation length becomes only one third of that in inviscid ODW. The great importance of considering viscosity in ODW simulations and future designs of combustor of oblique detonation engine has been addressed.     

Investigation of counter-rotating shock wave phenomenon and instability mechanisms of rotating detonation engine with hollow combustor and Laval nozzle

The counter-rotating shock wave phenomenon and instability mechanisms in the hollow rotating detonation combustor with a Laval nozzle are investigated. The in-house solver BYRFoam based on the OpenFOAM platform and the detailed chemical reaction kinetic mechanism are used. The coupling of the detonation wave and the shock wave is revealed, with the continuous presence of the counter-rotating shock wave in the combustor that propagates in the opposite direction to the detonation wave and oscillates in intensity. The intensities of both the detonation wave and the counter-rotating shock wave are found to be constantly oscillating, and this instability is referred to as the collision-oscillation instability. It is caused by the complex interaction between the detonation wave, the fresh gas and the counter-rotating shock wave. The velocity difference between the detonation wave and the counter-rotating shock wave is found to lead to the migration of the collision point, which in turn leads to the peak periodic undulation phenomenon of the pressure curve, that is, the wave packet. This instability is called the counter-rotating instability, which co-exists with the collision-oscillation instability for a long time. The results of this study are compared with the experimental data, and the reasons for the oscillation of the experimental pressure signal are explained.     

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.     

等离子体点火对燃烧转爆轰影响的数值计算

采用CE/SE方法对耦合的流体力学方程与麦克斯韦方程求解,对等离子体点火和气液两相爆轰过程进行数值模拟.讨论了3种不同点火位置对燃烧转爆轰(DDT)的影响.结果表明,等离子体点火可以在短时间内点燃爆轰管内汽油/空气混合物,经过一系列复杂的波系与波系、波系与管壁之间的相互作用最终形成稳定的爆轰波.壁面点火比中心点火具有优...     

热管蓄冰过程的数值模拟研究．

通过建立热管的简化热阻模型,采用数值计算的方法模拟了热管的动态蓄冰。模拟结果与实验结果吻合。同时,对热管冷凝段和蒸发段的长度比为1：2和1：4两种形式进行了数值模拟比较。结果表明,1：4的布置形式更优。