On the influence of low initial pressure and detonation stochastic nature on Mach reflection of gaseous detonation waves

The two-dimensional, time-dependent and reactive Navier–Stokes equations were solved to obtain an insight into Mach reflection of gaseous detonation in a stoichiometric hydrogen-oxygen mixture diluted by 25 % argon. This mixture generates a mode-7 detonation wave under an initial pressure of 8.00 kPa. Chemical kinetics was simulated by an eight-species, forty-eight-reaction mechanism. It was found that a Mach reflection mode always occurs for a planar detonation wave or planar air shock wave sweeping over wedges with apex angles ranging from \(5^\circ \) to \(50^\circ \) . However, for cellular detonation waves, regular reflection always occurs first, which then transforms into Mach reflection. This phenomenon is more evident for detonations ignited under low initial pressure. Low initial pressure may lead to a curved wave front, that determines the reflection mode. The stochastic nature of boundary shape and transition distance, during deflagration-to-detonation transition, leads to relative disorder of detonation cell location and cell shape. Consequently, when a detonation wave hits the wedge apex, there appears a stochastic variation of triple point origin and variation of the angle between the triple point trajectory and the wedge surface. As the wedge apex angle increases, the distance between the triple point trajectory origin and the wedge apex increases, and the angle between the triple point trajectory and the wedge surface decreases exponentially.     

Numerical simulation of Mach reflection of cellular detonations

The Mach reflection of cellular detonation waves on a wedge is investigated numerically in an attempt to elucidate the effect of cellular instabilities on Mach reflection, the dependence of self-similarity on the thickness of a detonation wave, and the initial development of the Mach stem near the wedge apex. A two-step chain-branching reaction model is used to give a thermally neutral induction zone followed by a chemical reaction zone for the detonation wave. A sufficiently large distance of travel of the Mach stem is computed to observe the asymptotic behavior in the far field. Depending on the scale at which the Mach reflection process occurs, it is found that the Mach reflection of a cellular detonation behaves essentially in the same way as a planar ZND detonation wave. The cellular instabilities, however, cause the triple-point trajectory to fluctuate. The fluctuations are due to interactions of the triple point of the Mach stem with the transverse waves of cellular instabilities. In the vicinity of the wedge apex, the Mach reflection is found to be self-similar and corresponds to that of a shock wave of the same strength, since the Mach stem is highly overdriven initially. In the far field, the triple-point trajectory approaches a straight line, indicating that the Mach reflection becomes self-similar asymptotically. The distance of the approach to self-similarity is found to decrease rapidly with decreasing thickness of the detonation front.     

气相爆轰波在障碍物上Mach反射的实验验证

本文公布了气相爆轰波沿收缩管道传播时发生Mach反射的实验证据。在爆轰波通过的管道中安装不同楔角的楔块,形成管道的收缩。爆轰波在通过楔块时会发生Mach反射。利用烟熏玻璃片记录到了爆轰波Mach反射时形成的三波点迹线及其两侧胞格尺寸和密度的变化。据我们掌握的资料,这是首次用胞格结构变化的记录证实,气相爆轰波与无化学反应的空气中的冲击波一样,在一定的入射条件下会发生Mach反射。这一实验结果可使我们更深入了解爆轰波的本质,也为数值模拟气相爆轰波在障碍物上Mach反射现象提供了可对比的依据。     

激波在不同三角楔面内传播特性的数值研究

激波在收缩管内的反射与聚焦会形成高温高压区,点燃可燃混合气并诱导爆轰,因此对爆轰发动机的点火具有重要意义。本文基于二维N-S方程,结合五阶WENO格式,对马赫数为6的正激波在三角形楔面内的反射与聚焦现象进行了数值研究。结果表明,楔面顶角的变化对激波的反射类型以及聚焦均有明显的影响:随着顶角的增加,激波的反射类型从马赫反射向过渡马赫反射和双马赫反射转变,且壁面上的前向射流更加明显;三波点第一次碰撞产生的高温高压区足够满足可燃混合气体的点火条件,且其温度与压力值随顶角的增加而增大;当激波在楔面上发生临界双马赫反射时,温度与压力达到最大;当顶角增加到一定值时,激波在楔面反射转变为常规反射,不会产生激波对碰,因而没有高温高压区。     

Numerical simulation of detonation reignition in H $_2$-O $_2$ mixtures in area expansions

Time-dependent, two-dimensional, numerical simulations of a transmitted detonation show reignition occuring by one of two mechanisms. The first mechanism involves the collision of triple points as they expand along a decaying shock front. In the second mechanism ignition results from the coalescence of a number of small, relatively high pressure regions left over from the decay of weakened transverse waves. The simulations were performed using an improved chemical kinetic model for stoichiometric H-O mixtures. The initial conditions were a propagating, two-dimensional detonation resolved enough to show transverse wave structure. The calculations provide clarification of the reignition mechanism seen in previous H-O-Ar simulations, and again demonstrate that the transverse wave structure of the detonation front is critical to the reignition process. Received 16 July 1998/ Accepted 17 November 1998     

Supersonic flow of a combustible gas mixture past a sphere

In recent years considerable interest has developed in the problems of steady-state supersonic flow of a mixture of gases about bodies with the formation of detonation waves and slow combustion fronts. This is due in particular to the problem of fuel combustion in a supersonic air stream.In [1] the problem of supersonic flow past a wedge with a detonation wave attached to the wedge apex is solved. This solution is based on using the equation of the detonation polar obtained in [2]-the analog of the shock polar for the case of an exothermic discontinuity. In [3] a solution is given of the problem of cone flow with an attached detonation wave, and [4] presents solutions of the problems of supersonic flow past the wedge and cone with the formation of attached adiabatic shocks with subsequent combustion of the mixture in slow combustion fronts. In the two latter studies two different solutions were also found for the problem of flow past a point ignition source, one solution with gas combustion in the detonation wave, the other with gas combustion in the slow combustion front following the adiabatic shock. These solutions describe two different asymptotic pictures of flow of a combustible gas mixture past bodies.In an experimental study of the motion of a sphere in a combustible gas mixture [5] it was found that the detonation wave formed ahead of the sphere splits at some distance from the body into an ordinary (adiabatic) shock and a slow combustion front. Arguments are presented in [6] which make it possible to explain this phenomenon and in certain cases to predict its occurrence.The present paper presents examples of the calculation of flow of a combustible gas mixture past a sphere with a detonation wave in the case when the wave does not split. In addition, the flow near the point at which the detonation wave splits is analyzed for the case when splitting occurs where the gas velocity behind the wave is greater than the speed of sound. This analysis shows that in the given case the flow calculation may be carried out without any particular difficulties. On the other hand, the calculation of the flow for the case when the point of splitting is located in the subsonic portion of the flow behind the wave (or in the region of influence of the subsonic portion of the flow) presents difficulties. This flow case is similar to the problem of the supersonic jet of finite width impacting on an obstacle.     

Wave dynamic processes in cellular detonation reflection from wedges

When the cell width of the incident detonation wave (IDW) is comparable to or larger than the Mach stem height, self-similarity will fail during IDW reflection from a wedge surface. In this paper, the detonation reflection from wedges is investigated for the wave dynamic processes occurring in the wave front, including transverse shock motion and detonation cell variations behind the Mach stem. A detailed reaction model is implemented to simulate two-dimensional cellular detonations in stoichiometric mixtures of H ₂/O ₂ diluted by Argon. The numerical results show that the transverse waves, which cross the triple point trajectory of Mach reflection, travel along the Mach stem and reflect back from the wedge surface, control the size of the cells in the region swept by the Mach stem. It is the energy carried by these transverse waves that sustains the triple-wave-collision with a higher frequency within the over-driven Mach stem. In some cases, local wave dynamic processes and wave structures play a dominant role in determining the pattern of cellular record, leading to the fact that the cellular patterns after the Mach stem exhibit some peculiar modes. The English text was polished by Yumming Chen.     

Numerical study of three-dimensional detonation structure transformations in a narrow square tube: from rectangular and diagonal modes into spinning modes

Three-dimensional (3-D) detonation structure transformations from rectangular and diagonal modes into spinning modes in a narrow square tube are investigated by high-resolution simulation. Numerical simulations are performed with a Riemann solver of the HLLC-type, new cell-based structured adaptive mesh refinement data structure, high-order, parallel adaptive mesh refinement reactive flow code. A simplified one-step kinetic reaction model is used to reveal the 3-D detonation structure. The four different types of initial disturbances applied in the ZND profiles lead to the structures of rectangular in phase, rectangular out of phase, rectangular partial out of phase and diagonal, respectively, during the initial stages of detonation propagation. Eventually, all these detonation structures evolve into the self-sustained spinning detonations. The asymmetric disturbance leads to a stable spinning detonation much faster than the rest. The important features in the formation of spinning detonation are revealed using a 3-D visualization, and a remarkable qualitative agreement with experimental and numerical results is obtained with respect to the transverse wave dynamics and detonation front structures. The transverse wave collisions produce the unburnt gas pockets and the energy to sustain the detonation front propagation and distortion. The periodic pressure oscillation of front plays a complex role as it shifts the reaction zone structure with an accompanying change in the driving energy of transition and the detonation parameters which result in the more distorted front and the unstable detonation. Eventually, the unstable distorted detonation evolves into a spinning detonation.     

Behavior of detonation propagation in mixtures with concentration gradients

Behavior of detonation waves in mixtures with concentration gradients normal to the propagation direction was studied experimentally. Mixtures with various concentration gradients were formed by sliding the separation plate which divides a detonation chamber from a diffusion chamber in which a diffusion gas was initially introduced. A stoichiometric hydrogen–oxygen mixture was charged in the detonation chamber, while oxygen or nitrogen was filled in the diffusion gas chamber. Temporal concentration measurement was conducted by the infrared absorption method using ethane as alternate of oxygen. Smoked foil records show a deformation of regular diamond cells to parallelogram ones, which well corresponds to local mixture concentration. Schlieren photographs reveal the tilted wave front whose angle is consistent with the deflection angle of the detonation front obtained from trajectories of the triple point. The local deflection angle increases with increase in local concentration gradient. Calculation of wave trajectory based on the ray tracing theory predicts formation of the tilted wave front from an initial planar front.      

爆轰波与激波对撞爆轰状态演化机理研究

对当量比氢氧混合气体中爆轰波与激波的正面对撞过程进行了二维数值研究. 采用了二阶精度NND差分格式与改进的两阶段化学反应模型,并以数值x-t纹影图以及烟迹图记录了对撞过程. 数值研究表明,透射爆轰波受到膨胀影响首先会衰减,甚至产生局部解耦现象;然后由于三波点的碰撞又能再次重新耦合. 在爆轰波对撞过程中,由于燃烧不均匀性而产生的弱横波对爆轰胞格的形成起着重要作用.      

The nature of the triple point singularity in the case of stationary reflection of weak shock waves

The structure of flow in the vicinity of a triple point in the problem of stationary irregular reflection of weak shock waves is numerically investigated within the framework of the Euler model, including the von Neumann paradox range. To improve the accuracy of the solution near singular points a new technology including a grid contracted toward the triple point and the discontinuity fitting is applied. It is shown that in the four-wave flow pattern the curvatures of the tangential discontinuity and the Mach front at the triple point are finite. The singularity is concentrated only in a sector between the reflected wave front and the expansion fan. When the three-wave flow pattern is realized, the curvatures of the tangential discontinuity and both wave fronts at the triple point are infinite. On the range of weak and moderate shock waves the logarithmic singularity in subsonic sectors near the triple point conserves up to transition to the regular reflection.     

Analysis of the shock structures in a regular detonation

Time-dependent two-dimensional numerical simulations have been used to investigate the detailed shock structures and patterns of energy release in the regions of the triple points and transverse waves in a planar detonation. As the system of shock triple points evolves between collisions, they trace a well shaped cellular pattern characteristic of detonations in argon-diluted, low-pressure mixtures of hydrogen and oxygen. In the region of the triple points, the shock structure evolves continuously from a single Mach structure to a double Mach structure and finally to a complex Mach structure characteristic of spinning detonations. Most of the energy released in the region of the triple points. The amount of energy release increases as the triple point comes closer to a collision with a wall or another triple point. Just before the collision, there is a large region of energy release that covers the length of the interacting transverse waves. The result is a rectangular high-energy region which boosts the propagation of the new detonation cell.     

Two-dimensional reactive flow dynamics in cellular detonation waves

This investigation deals with the two-dimensional unsteady detonation characterized by the cellular structure resulting from trajectories of triple-shock configurations formed by the transverse waves and the leading shock front. The time-dependent reactive shock problem considered here is governed by a system of nonlinear hyperbolic conservation laws coupled to a polytropic equation of state and a one-step Arrhenius chemical reaction rate with heat release. The numerical solution obtained allowed us to follow the dynamics of the cellular detonation front involving the triple points, transverse waves and unreacted pockets. The calculations show that the weak tracks observed inside the detonation cells around the points of collision of the triple-shock configurations arise from interactions between the transverse shocks and compression waves generated by the collision. The unreacted pockets of gas formed during the collisions of triple points change form when the activation energy increases. For the self-sustained detonation considered here, the unreacted pockets burn inside the region independent of the downstream rarefaction, and thus the energy released supports the detonation propagation. The length of the region independent of the downstream is approximately the size of one or two detonation cell. Received 13 February 1998 / Accepted 13 August 1998     

周期性非均匀介质中气相爆轰波演变模式研究

气相爆轰波在周期性非均匀介质中的起爆, 稳态传播和失效机制都极为复杂, 很多物理机制尚不明确, 是当前爆轰物理领域研究的热点和难点. 本文使用反应欧拉方程和两步化学反应模型对爆轰波在非均匀介质中的传播机理进行了数值模拟研究, 非均匀性由横向周期性分布的温度扰动体现, 重点分析不同波长、不同幅度的温度扰动对波阵面波系结构的影响. 计算结果表明, ZND爆轰波在温度扰动下向胞格爆轰波的转变主要受制于两种竞争性因素: 一是爆轰波内在的不稳定性; 二是温度扰动的波长和幅度, 前者是内因, 后者是外因. 温度扰动的存在抑制横波的发展, 延迟了ZND爆轰波向胞格爆轰波的演化, 并且内在不稳定性的增加可以减慢这种延迟现象. 这说明, 温度扰动可以在一定的范围内抑制胞格不稳定性的发展, 但是不能够终止这一过程. 温度的不连续性使得爆轰波阵面更为扭曲, 并在横波附近存在较弱的三波点结构, 即温度扰动可增加爆轰波固有的不稳定性, 改变爆轰波阵面的传播机理. 幅值较大的人工温度扰动可抑制爆轰波的传播和爆轰波自身的不稳定性. 爆轰波阵面胞格结构的形成取决于温度扰动与其自身的不稳定性.      

一维与二维爆轰传播的时空关联特性数值研究

一维爆轰传播的理论完备、计算准确, 二维斜爆轰传播由于壁面与黏性效应, 大尺度、高精度预测还有一定难度. 利用Euler方程和H₂-Air基元反应模型, 对二维有限长楔面诱导的斜爆轰和活塞驱动一维非定常正爆轰进行计算比较研究, 从时空两个维度方面, 分析了两者在起爆过程、稀疏波传播、爆轰波面演化中的关联特性. 研究结果表明: 在过驱动度相同的条件下, 经过时空变换的活塞驱动一维爆轰传播与二维驻定斜爆轰在起爆区波系结构、波面演化特征和主要参数分布规律方面无论定性或者定量对比均符合较好, 所以, 一维非定常爆轰和二维驻定斜爆轰具有时空相关性. 两者的差异主要体现在过驱动斜爆轰受稀疏波影响过渡到近Chapman-Jouguet (C-J)爆轰状态所需的弛豫时间不同, 原因可能是起源于活塞和壁面稀疏波强度的差异. 本文提出的一维与二维爆轰传播的时空关联方法不仅有助于认知斜爆轰起爆、过驱爆轰产生、胞格爆轰演化的三阶段规律, 还可以对比揭示壁面、边界层和黏性效应的影响, 应用在斜爆轰发动机燃烧室设计中能够有效节约计算时间和成本, 并降低复杂度.      

Investigation of the detonation regimes in gaseous mixtures with suspended starch particles

The existence of a secondary discontinuity at the rear of a detonation front shown in experiments by Peraldi and Veyssiere (1986) in stoichiometric hydrogen-oxygen mixtures with suspended 20-m starch particles has not been explained satisfactorily. Recently Veyssiere et al. (1997) analyzed these results using a one-dimensional (1-D) numerical model, and concluded that the heat release rate provided by the burning of starch particles in gaseous detonation products is too weak to support a double-front detonation (DFD), in contrast to the case of hybrid mixtures of hydrogen-air with suspended aluminium particles in which a double-front detonation structure was observed by Veyssiere (1986). A two-dimensional (2-D) numerical model was used in the present work to investigate abovementioned experimental results for hybrid mixtures with starch particles. The formation and propagation of the detonation has been examined in the geometry similar to the experimental tube of Peraldi and Veyssiere (1986), which has an area change after 2 m of propagation from the ignition point from a 69 mm dia. section to a 53 mm 53 mm square cross section corresponding to a 33% area contraction. It is shown that the detonation propagation regime in these experiments has a different nature from the double-front detonation observed in hybrid mixtures with aluminium particles. The detonation propagates as a pseudo-gas detonation (PGD) because starch particles release their heat downstream of the CJ plane giving rise to a non-stationary compression wave. The discontinuity wave at the rear of the detonation front is due to the interaction of the leading detonation front with the tube contraction, and is detected at the farthest pressure gauge location because the tube length is insufficient for the perturbation generated by the tube contraction to decay. Thus, numerical simulations explain experimental observations made by Peraldi and Veyssiere (1986). Received 5 July 1997 / Accepted 13 July 1998     

Structure of the detonation wave front in a mixture of nitromethane with acetone

It is shown that the leading front of an inhomogeneous detonation wave is a shock wave in which wave structures of the type of triple shock configurations are moving. It was experimentally found that the reaction in these inhomogeneities occurs in oblique shock waves. The reaction sites at the wave front are ring-shaped. In a 75: 25 mixture of nitromethane with acetone, up to 70% of the front surface is occupied by the reaction at the sites in the wave front. Measurements of the mass velocity profile indicate that afterburning takes place in the unloading area behind the Jouguet plane. Calculations of the heat release in the reaction mixture with a decrease in the mass velocity indicate that the material that have not reacted in the inhomogeneities can be ignited in the induction zone. It is suggested that the adiabatic flashes are a mechanism that generates inhomogeneities in the detonation wave front.     

The criterion of the existence or inexistence of transverse shock wave at wedge supported oblique detonation wave

A simplified theoretic method and numerical simulations were carried out to investigate the characterization of propagation of transverse shock wave at wedge supported oblique detonation wave.After solution validation,a criterion which is associated with the ratio Φ (u 2 /u CJ) of existence or inexistence of the transverse shock wave at the region of the primary triple was deduced systematically by 38 cases.It is observed that for abrupt oblique shock wave (OSW)/oblique detonation wave (ODW) transition,a transverse shock wave is generated at the region of the primary triple when Φ < 1,however,such a transverse shock wave does not take place for the smooth OSW/ODW transition when Φ > 1.The parameter Φ can be expressed as the Mach number behind the ODW front for stable CJ detonation.When 0.9 < Φ < 1.0,the reflected shock wave can pass across the contact discontinuity and interact with transverse waves which are originating from the ODW front.When 0.8 < Φ < 0.9,the reflected shock wave can not pass across the contact discontinuity and only reflects at the contact discontinuity.The condition (0.8 < Φ < 0.9) agrees well with the ratio (D ave /D CJ) in the critical detonation.     

A numerical simulation of reflection processes of a detonation wave on a wedge

Abstract. A two dimensional numerical simulation has been performed to study reflection processes of detonation waves on a wedge. The numerical scheme adopted is the flux corrected transport scheme and a two-step chemical reaction is assumed for a stoichiometric oxyhydrogen mixture diluted with argon. Transverse wave structures of the detonation are produced by artificial disturbances situated in front of a one-dimensional Chapman-Jouguet detonation wave. Numerical grids are generated by solving a Laplace equation. Results show that in the case where Mach reflection occurs, the cells in the Mach stem are smaller than those in the incident wave and are distorted in shape. There is also an initiating stage during which the cells in the Mach stem are created. The critical angle beyond which Mach reflection cannot occur is discussed. Received 15 October 1999 / Accepted 27 March 2000