气相爆轰波传播过程中的自点火效应

基于基元反应模型和单步反应模型,对直管道中H₂-air混合气体中爆轰波的传播过程进行了数值模拟,揭示了气相爆轰波传播过程中的自点火效应。利用数值模拟方法计算了不同爆轰模型的点火延迟时间,并得到了爆轰波三波点的传播过程以及所形成胞格结构的尺寸。结果表明,胞格宽度与点火延迟时间成正比;爆轰波诱导区内气体的点火延迟时间与三波点的运动周期基本一致。进一步对结果分析可知,爆轰波的自维持传播取决于点火延迟时间(表征化学反应的特征时间)和三波点的运动周期(表征流动的特征时间)的匹配;当二者相匹配时,经过前导激波压缩后形成的高温高压爆轰气体,在短时间内实现了自点火,同时释放出大量的能量推动了爆轰波的前进,即爆轰波的稳定自维持传播依靠其自点火机制。     

STUDY ON THE RELATIONSHIP BETWEEN IGNITION DELAY TIME AND GASEOUS DETONATION CELL SIZE

讨论点火延迟时间和爆轰波胞格尺度的内在关系,将点火延迟时间作为特征参量来模拟胞格尺度. 分别对两个总包单步化学反应模型和一个基元反应模型的点火延迟时间进行了数值模拟研究. 对于满足当量比的氢气/空气混合气体,分析了不同初始压力下点火延迟时间随初始温度的变化关系. 研究表明:总包单步反应模型的点火延迟时间不随压力变化,且与初始温度呈线性关系. 基元反应模型的点火延迟时间随压力变化,而且存在理论上的S 型曲线,但是在拐点区域和低温区域与CHEMKIN 计算的结果相差1~3 个量级. 现有模型模拟的胞格尺度普遍偏小,其相应的点火延迟时间也偏小,胞格尺度与点火延迟时间具有正相关性. 入射激波后的气体的点火延迟时间与三波点的运动周期一致,是定量化模拟胞格的关键因素.     

气相爆轰波胞格尺度与点火延迟时间关系研究

讨论点火延迟时间和爆轰波胞格尺度的内在关系,将点火延迟时间作为特征参量来模拟胞格尺度. 分别对两个总包单步化学反应模型和一个基元反应模型的点火延迟时间进行了数值模拟研究. 对于满足当量比的氢气/空气混合气体,分析了不同初始压力下点火延迟时间随初始温度的变化关系. 研究表明:总包单步反应模型的点火延迟时间不随压力变化,且与初始温度呈线性关系. 基元反应模型的点火延迟时间随压力变化,而且存在理论上的S 型曲线,但是在拐点区域和低温区域与CHEMKIN 计算的结果相差1~3 个量级. 现有模型模拟的胞格尺度普遍偏小,其相应的点火延迟时间也偏小,胞格尺度与点火延迟时间具有正相关性. 入射激波后的气体的点火延迟时间与三波点的运动周期一致,是定量化模拟胞格的关键因素.      

用于爆轰驱动的射流起爆实验研究

爆轰驱动激波风洞的自由来流模拟范围与驱动气体的爆轰极限密切相关,爆轰极限越宽则模拟范围越大。驱动气体一般是通过点火管进行起爆的,提高点火管的起爆能力可以拓宽爆轰极限。为了提高点火管起爆能力,就点火管口径、点火气体爆轰敏感性和单/双点火管3种因素的影响进行了实验研究。在不同的点火管初始条件下,对驱动段波速进行了测量。结论如下:(1)提高点火管口径可以显著提升起爆能力;(2)点火气体爆轰敏感性对起爆能力有影响,点火管为缩径内型面时,低敏感性气体起爆能力更强,点火管为等径内型面时则低敏感性气体和高敏感性气体的起爆能力大体持平;(3)在保证射流同步的前提下,双点火管能够提高起爆能力,为保证射流同步性需使用化学恰当比的氢氧混气等爆轰敏感性强的点火气体。     

Theoretical investigation of ignition and detonation of coal-particle gas mixtures

A mathematical model of ignition of a coal-particle gas suspension is developed within the framework of mechanics of reacting heterogeneous media. Some qualitative features are studied, which reveal different typical variants of heat dynamics of the mixture: heterogeneous ignition caused by the coke residue ignition reaction, homogeneous ignition induced by volatile oxidation in the gas phase, and hybrid ignition due to the simultaneous effect of surface and volume reactions. Verification of the model is performed using available experimental data on ignition delay times for coal-particle suspensions in air and oxygen behind reflected shock waves. The model is extended to detonation processes in the suspension, verified by experimental data on the relation between the propagation velocity and initial particle concentration. Stationary detonation regimes are studied, and specific features of detonation structures are discussed.This paper was presented at the 18th International Colloquium on the Dynamics of Explosions and Reactive System, Seattle, Washington, USA, August, 2001     

Rapid detonation initiation by sparks in a short duct: a numerical study

Rapid onset of detonation can efficiently increase the working frequency of a pulse detonation engine (PDE). In the present study, computations of detonation initiation in a duct are conducted to investigate the mechanisms of detonation initiation. The governing equations are the Euler equations and the chemical kinetic model consists of 19 elementary reactions and nine species. Different techniques of initiation have been studied for the purpose of accelerating detonation onset with a relatively weak ignition energy. It is found that detonation ignition induced by means of multiple sparks is applicable to auto-ignition for a PDE. The interaction among shock waves, flame fronts and the strip of pre-compressed fresh (unburned) mixture plays an important role in rapid onset of detonation.     

Natural gas ignition delay times behind reflected shock waves: Application to modelling and safety

Ignition delay times have been measured behind reflected shock waves using infrared and ultraviolet emission spectroscopies in the wide ranges of temperature (1485-1900 K), pressure (0.3-1.3 MPa), equivalence ratio (0.5-2) and dilution with argon (95-99%). For a mixture composed of several alkanes, it has been shown that the ignition delay can be evaluated from a combination of individual expressions obtained for each alkane (methane, ethane and propane). This method is applied to ignition delay time calculations for mixtures of several natural gases and validated with a particular Algerian natural gas. The measured ignition delay time are compared with those computed using a detailed kinetic mechanism. Agreement is generally satisfactory but the deviation between measurements and calculations increases when the ethane relative fraction increases in natural gas. Dynamic parameters of detonation in air are evaluated from these correlations. The predicted kinetic parameters of detonation (cell width) with a natural gas from Abu Dhabi, especially rich in ethane, are about half as those estimated with a Russian natural gas.Received: 16 April 2002, Accepted: 4 October 2002, Published online: 8 May 2003     

爆轰波平掠惰性气体界面及其解耦现象的数值研究

采用基元化学反应模型和迎风TVD格式,数值研究了爆轰波平掠惰性气体界面时的物理现象及其作用机制,并用点隐方法克服化学反应源项引起的计算刚性。数值结果显示,当爆轰波平掠过惰性气体界面时,形成了爆轰波、界面、透射激波以及稀疏波相互作用的现象。在高N2比例稀释的可燃混合气体情况下,当爆轰波平掠过特定惰性气体界面时,它与惰性气体界面相互作用产生的稀疏波可以导致爆轰波的解耦。     

Mechanism of deflagration-to-detonation transitions above repeated obstacles

Experiments are carried out to investigate the mechanism of the deflagration-to-detonation transition (DDT). Because, this mechanism has relevance to safety issues in industries, where combustible premixed gases are in general use. A stoichiometric gas of oxygen and hydrogen (oxy-hydrogen) is ignited in a tube, repeated obstacles are installed, and the DDT behaviours are visualized using a high-speed video camera. The pitch and height of the repeated obstacles and the initial pressure of the oxy-hydrogen premixed gas are varied in an attempt to obtain the optimum conditions that cause DDT a short distance from the ignition source. The experiments identified DDT as being essentially caused by one of the following mechanisms: (1) A deflagration wave is accelerated in terms of a vortex, which is generated behind the obstacle, and the flame acceleration induces a secondary shock wave. Eventually, the shock–flame interaction ahead of the obstacle causes DDT via a very strong local explosion. (2) Each shock wave generated by relatively weak local explosions between the obstacles is not sufficient to cause DDT directly, but DDT results from an accumulation of shock waves. The detonation induction distance is also examined, taking into account the physical and chemical parameters of the obstacles and the oxy-hydrogen premixed gas.     

气体-燃料液滴两相系统爆轰的数值模拟

用两相流体力学模型对气体 燃料液滴系统进行了研究。数值模拟了点火后两相系统爆轰波的发展过程，得到爆轰波的结构和参数。数值模拟结果表明气体 燃料液滴系统爆轰波有较宽的反应区，因而两相爆轰波的曲率对爆速的影响效应十分明显。进行了燃料液滴尺寸对爆轰波的结构和参数的影响的数值模拟。除了很小的液滴外，燃料液滴在爆轰波前导激波面和ＣＪ面间不能完全气化。随着液滴尺寸的增加，燃料液滴在爆轰波前导激波面和ＣＪ面间释放出的能量随之减少，爆轰参数也随之下降。     

用环形激波聚焦实现爆轰波直接起爆的数值模拟

利用基元反应模型和有限体积法对环形激波在可燃气体中聚焦实现爆轰波直接起爆进行了数值模拟。研究结果表明，标准状态下的氢气-空气混合气体在马赫数为3.1以上的环形激波聚焦产生的高温高压区作用下会诱发可燃气体的直接起爆形成爆轰波，爆轰波与激波和接触间断相互作用产生了复杂的波系结构；爆轰波爆点位置在对称轴上并不是固定的点，而是随着初始激波马赫数的变化而发生移动；可燃气体初始温度和压力对起爆临界马赫数都有影响，但是初始温度的影响大得多。     

Detonation formation in H $_2$-O $_2$/He/Ar mixtures at elevated initial pressures

In this paper the formation of detonation in H-O/He/Ar mixtures at elevated initial pressures was investigated in an initiation tube for a detonation driver with an exploding wire as the ignition source. In most experiments the detonation wave was formed by a DDT process in which a reactive shock wave accelerates behind the leading shock wave and eventually leads to the onset of detonation. The onset position was found to be at the leading shock wave or behind it. Only in very sensitive mixtures at high initial pressure the direct initiation of detonation was observed. The influence of ignition energy, initial pressure and composition on the detonation induction distance was determined. The results show that the detonation induction distance increases with the decrease of ignition energy and initial pressure and with the increase of the mole fraction of helium or argon. With the same mole fraction, argon increases the induction distance more than helium. In the facility utilized the DDT upper and lower limits of hydrogen in H-O mixtures are in the ranges from 36 to 40 % and from 78 to 82 %, respectively, and the upper limits for helium and argon in stoichiometric H-O mixtures are 40 % and 36 %, respectively. High pressure peaks generated by the DDT process were measured, especially in mixtures near the DDT limits. Statistical results show that such peak pressures can be up to 6 times of the CJ-pressures. Received 1 March 2000 / Accepted 25 May 2000     

Behavior of detonation waves at low pressures

With respect to stability of gaseous detonations, unsteady behavior of galloping detonations and re-initiation process of hydrogen-oxygen mixtures are studied using a detonation tube of 14 m in length and 45 mm i.d. The arrival of the shock wave and the reaction front is detected individually by a double probe combining of a pressure and an ion probe. The experimental results show that there are two different types of the re-initiation mechanism. One is essentially the same as that of deflagration to detonation transition in the sense that a shock wave generated by flame acceleration causes a local explosion. From calculated values of ignition delay behind the shock wave decoupled from the reaction front, the other is found to be closely related with spontaneous ignition. In this case, the fundamental propagation mode shows a spinning detonation. Received 10 March 1997 / Accepted 8 June 1997     

Similarity motions of a gas in the case of local supply of mass and energy in a fuel mixture

One-dimensional nonstationary similarity motions of a gas with exothermic reactions behind shock waves are analyzed. The thickness of the region of chemical reactions is ignored. New solutions are obtained for the problems of flows of a chemically active gas with the formation of shock and detonation waves. In particular, it is shown that in the framework of the adopted schemes of the combustion process a solution with five strong-discontinuity surfaces can be constructed. The results are given of numerical solutions for supercompression detonation and Chapman-Jouguet detonation.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 70–77, July–August, 1979.     

Ignition and combustion study of JP-8 fuel in a supersonic flowfield

The ignition and combustion process of fuels in a supersonic combustion chamber plays an important role in the design of hypersonic propulsion system. However, it is a quite complicated process, due to the large variation of inlet air velocity, temperature, oxygen concentration, and shocks in the supersonic combustion chamber. The purpose of this paper is to observe the ignition delay and combustion phenomenon of the JP-8 fuel droplets in a supersonic flowfield experimentally. A shock tube is used as a basic test facility to create a high-speed and high-temperature flowfield as a supersonic combustor. In the experiments, several test parameters are controlled, such as shock velocity, gas temperature, fuel droplet size and distance, initial fuel temperature, and oxygen concentration, etc. The test results show the influence of these parameters on ignition delay, ignition limitation, and detonation. The most important factor in the experiment is the initial fuel temperature effect, which is influenced by the altitude variation during a flight. Received 4 August 1995 / Accepted 12 December 1995     

可燃气体中激波聚焦的点火特性

数值模拟了二维平面激波从抛物面上反射在可燃气体中聚焦的过程,研究了形 成爆轰波的点火特性. 对理想化学当量比氢气/空气混合气体,在初始压强20kPa的条件下, 马赫数2.6-2.8的激波聚焦能产生两个点火区：第1个点火区是反射激波会聚引起的,第 2个点火区是由入射激波在抛物面上发生马赫反射引起的. 这种条件下流场中会出现爆燃转 爆轰,起爆点分别分布在管道壁面、抛物反射面和第2点火区附近. 起爆机理分别为激波管 道壁面反射、点火诱导激波的抛物面反射和点火诱导的激波与第2点火区产生的爆燃波的相 互作用. 不同的点火和起爆过程导致了不同的流场波系结构,同时影响了爆轰波传播的波动 力学过程.     

Effects of nitrates on hydrocarbon-air flames and detonations

Abstract. The subject of hydrocarbon sensitization by nitrates under conditions of a heterogeneous spray has been of interest due to its applicability in promoting ignition. To gain insight into the mechanisms of the nitrate sensitization effect, the present work was limited to vapour phase studies at elevated temperatures in order to avoid the influence of heterogeneous factors. The experiments performed included studies of flammability, flame propagation, shock ignition and detonation. The mixtures used were composed of air, hexane, and isopropyl nitrate (IPN) with IPN concentrations ranging from 0 to 100 mole % in hydrocarbon-IPN. In addition, methane and propane were also included in the flame experiments. For the shock ignition and detonation experiments, the measured ignition delay and detonation cell size had minimum values for IPN-air and maximum values for hexane-air. With increases in the IPN concentration, the ignition delay and detonation cell size fell monotonically between the values for hexane and IPN. This monotonic behaviour was explained to be the result of mixing the hydrocarbon with the more sensitive nitrate whose energetics are larger than or comparable to the hydrocarbon when mixed with air. For the slow combustion mode, the results also confirmed the monotonic behavior and showed that the lean flammability limit and the flame velocity fell between those of the hydrocarbon and IPN. Received 10 September 1999 / Accepted 27 July 2000     

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