The propagation mechanism of high speed turbulent deflagrations

The propagation regimes of combustion waves in a 30 cm by 30 cm square cross–sectioned tube with an obstacle array of staggered vertical cylindrical rods (with BR=0.41 and BR=0.19) are investigated. Mixtures of hydrogen, ethylene, propane, and methane with air at ambient conditions over a range of equivalence ratios are used. In contrast to the previous results obtained in circular cross–sectioned tubes, it is found that only the quasi–detonation regime and the slow turbulent deflagration regimes are observed for ethylene–air and for propane–air. The transition from the quasi–detonation regime to the slow turbulent deflagration regime occurs at (where D is the tube “diameter” and is the detonation cell size). When , the quasi–detonation velocities that are observed are similar to those in unobstructed smooth tubes. For hydrogen–air mixtures, it is found that there is a gradual transition from the quasi–detonation regime to the high speed turbulent deflagration regime. The high speed turbulent deflagration regime is also observed for methane–air mixtures near stoichiometric composition. This regime was previously interpreted as the “choking” regime in circular tubes with orifice plate obstacles. Presently, it is proposed that the propagation mechanism of these high speed turbulent deflagrations is similar to that of Chapman–Jouguet detonations and quasi-detonations. As well, it is observed that there exists unstable flame propagation at the lean limit where . The local velocity fluctuates significantly about an averaged velocity for hydrogen–air, ethylene–air, and propane–air mixtures. Unstable flame propagation is also observed for the entire range of high speed turbulent deflagrations in methane–air mixtures. It is proposed that these fluctuations are due to quenching of the combustion front due to turbulent mixing. Quenched pockets of unburned reactants are swept downstream, and the subsequent explosion serves to overdrive the combustion front. The present study indicates that the dependence on the propagation mechanisms on obstacle geometry can be exploited to elucidate the different complex mechanisms of supersonic combustion waves. Received 5 November 2001 / Accepted 12 June 2002 / Published online 4 November 2002 Correspondence to: J. Chao (e-mail: jenny.chao@mail.mcgill.ca) An abridged version of this paper was presented at the 18th Int. Colloquium on the Dynamics of Explosions and Reactive Systems at Seattle, USA, from July 29 to August 3, 2001.     

On some conditions for detonation initiation downstream of a perforated plate

An experimental study is presented of the influence of detonation wave parameters and detonation product composition upstream of a perforated plate on the onset of detonation downstream. Experiments were performed in detonation tube 106 mm in diameter, separated into two sections by a perforated plate combined with a diaphragm. The tube was equipped with pressure sensors and a semi-cylindrical smoked foil. Hydrogen–air mixtures with different hydrogen concentrations were used upstream and downstream of the perforated plate. It is shown for mixtures containing 25 and 34 % of hydrogen in air that the onset of detonation downstream depends on detonation parameters upstream of the perforated plate. An increase in the initial pressure upstream of the plate leads to detonation initiation immediately downstream. The variation of mixture composition upstream of a perforated plate does not affect on detonation initiation downstream under the present experimental conditions.     

Hydrogen detonation and fast deflagration triggered by a turbulent jet of combustion products

Initiation of detonation by a turbulent jet of combustion products has been studied in a detonation tube of 141 mm inner diameter. Jet formation techniques based on either a perforated plate or bursting membrane subjected to the impact of a stable detonation wave were utilized. Critical conditions of detonation initiation in hydrogen–air and hydrogen–oxygen–nitrogen mixtures have been found to depend on both the mixture sensitivity and the geometrical parameters of the arrangement. PACS 47.70.Fw; 82.33.Vx; 82.40.Fp This paper was based on work that was presented at the 19th Inter-national Colloquium on the Dynamics of Explosions and Reactive Systems, Hakone, Japan, July 27 - August 1, 2003     

A study on jet initiation of detonation using multiple tubes

A detonator consisting of a dense bundle of small-diameter tubes (4.4–19 mm) is tested experimentally using stoichiometric mixtures of hydrogen–oxygen and hydrogen–air. Tests are conducted in a 5,200-mm long detonation tube fitted with a schlieren photograph section and smoked foil to record the deflagration to detonation (DDT) transition. It is confirmed that the flame jet emanating from the tube assembly causes detonation initiation immediately downstream of the detonator, with little dependence on the size of the detonation tube. For the fuel–air mixture, the insertion of Shchelkin spirals into each of the smaller tubes enhances the development of the turbulent flame jet, leading to a shorter DDT distance. Multi-point spark ignition is also shown to provide a further reduction in the DDT distance compared to single-point ignition. PACS 47.40.-x; 47.40.Nm; 47.70.Fw; 82.40.-g; 82.40.Fp     

DDT in a smooth tube filled with a hydrogen–oxygen mixture

Results of experimental study on DDT in a smooth tube filled with sensitive mixtures having detonation cell size from 1 to 3 orders of magnitude smaller than the tube diameter are presented. Stoichiometric hydrogen–oxygen mixtures were used in the tests with initial pressure ranging from 0.2 to 8 bar. A dependence of the run-up distance to DDT on the initial pressure is studied. This dependence is found to be close to the inverse proportionality. It is suggested that the flow ahead of the flame results in formation of the turbulent boundary layer. This boundary layer controls the scale of turbulent motions in the flow. A simple model to estimate the maximum scale of the turbulent pulsations (boundary layer thickness) at flame positions along the tube is presented. The largest scale of the turbulent motions at the location of the onset of detonation is shown to be 1 order of magnitude greater than the detonation cell widths, λ, in all the tests. It is suggested that the onset of detonation is triggered during flame acceleration as soon as the maximum scale of the turbulent pulsations increases up to about 10 λ. The model to estimate the maximum size of turbulent motions, δ, and the correlation δ≈ 10λ, give a basis for estimations of the run-up distances to DDT in tubes with internal diameter D > 20λ. PACS 47.40.-x; 47.27.Nz This paper was based on work that was presented at the 19th Inter-national Colloquium on the Dynamics of Explosions and Reactive Systems, Hakone, Japan, July 27 - August 1, 2003     

Detonation diffraction from an annular channel

In this study, gaseous detonation diffraction from an annular channel was investigated with a streak camera and the critical pressure for transmission of the detonation wave was obtained. The annular channel was used to approximate an infinite slot resulting in cylindrically expanding detonation waves. Two mixtures, stoichiometric acetylene–oxygen and stoichiometric acetylene–oxygen with 70% Ar dilution, were tested in a 4.3 and 14.3 mm channel width (W). The undiluted and diluted mixtures were found to have values of the critical channel width over the cell size around 3 and 12 respectively. Comparing these results to values of the critical diameter (d _c ), in which a spherical detonation occurs, a value of critical d _c /W _c near 2 is observed for the highly diluted mixture. This value corresponds to the geometrical factor of the curvature term between a spherical and cylindrical diverging wave. Hence, the result is in support of Lee’s proposed mechanism [Lee in Dynamics of Exothermicity, pp. 321, Gordon and Breach, Amsterdam, 1996] for failure due to diffraction based on curvature in stable mixtures such as those highly argon diluted with very regular detonation cellular patterns.     

Effect of scale on the onset of detonations

Critical conditions for onset of detonations are compared at (1) two significantly different scales, (2) for a range of -air mixtures diluted with C, O, and (3) for two types of geometry – one a long obstructed channel and the other a room with a relatively small aspect ratios. For the range of scales, mixtures, and initial conditions tested, the detonation cell size was shown to be a reliable scaling parameter for characterization of detonation onset conditions. An experimental correlation is suggested for the critical detonation onset conditions. This correlation is based on a wide variety of available experimental data on DDT in mixtures of hydrogen and hydrocarbon fuels with air and on the use of detonation cell size as a scaling parameter characterizing the mixture. Received 14 November 1999 / Accepted 16 February 2000     

Measurement and chemical kinetic model predictions of detonation cell size in methanol–oxygen mixtures

In this study, detonation cell sizes of methanol–oxygen mixtures are experimentally measured at different initial pressures and compositions. Good agreement is found between the experiment data and predictions based on the chemical length scales obtained from a detailed chemical kinetic model. To assess the detonation sensitivity in methanol–oxygen mixtures, the results are compared with those of hydrogen–oxygen and methane–oxygen mixtures. Based on the cell size comparison, it is shown that methanol–oxygen is more detonation sensitive than methane–oxygen but less sensitive than hydrogen–oxygen.     

Measurement and scaling analysis of critical energy for direct initiation of gaseous detonations

In this paper, the critical energies required for direct initiation of spherical detonations in four gaseous fuels (C₂H₂, C₂H₄, C₃H₈ and H₂)–oxygen mixtures at different initial pressures, equivalence ratios and with different amounts of argon dilution are reported. Using these data, a scaling analysis is performed based on two main parameters of the problem: the explosion length R _o that characterizes the blast wave and a characteristic chemical length that characterizes the detonation. For all the undiluted mixtures considered in this study, it is found that the relationship is closely given by R_o ? 26 l{R_{\rm o} \approx 26 \lambda} , where λ is the characteristic detonation cell size of the explosive mixture. While for C₂H₂–2.5O₂ mixtures highly diluted with argon, in which cellular instabilities are shown to play a minor role on the detonation propagation, the proportionality factor increases to 37.3, 47 and 54.8 for 50, 65 and 70% argon dilution, respectively. Using the ZND induction length Δ _I as the characteristic chemical length scale for argon diluted or ‘stable’ mixtures, the explosion length is also found to scale adequately with R_o ? 2320 D_I{R_{\rm o} \approx 2320 \Delta_I} .     

爆轰波透射孔栅形成的高速爆燃波的结构和行为

以实验为主,研究光滑直管中乙炔氧气爆轰波透射孔栅形成的高速爆燃波结构和行为。实验以基于纹影平台的高速转鼓摄影记录孔栅下游近场内初始爆燃波的结构和发展,并以压力传感器跟踪其后继走向。研究发现,高速爆燃波是前驱激波和火焰的组合结构,随着初始压力的提升,它分别表现为紧随于前驱激波的层流和湍流燃烧火焰。入射爆轰胞格尺度小于或与孔栅扰动尺度相当时,湍流燃烧在下游近场迅速形成;层流结构的爆燃波通常无法抵制背景稀疏波而走向衰弱,而湍流燃烧结构可发生加速和向爆轰的转捩;他们之间存在一个不稳定的临界状态,高速爆燃波得以以0.5~0.6倍CJ爆轰速度传播较长距离,这一状态对应于双间断Rankine-Hugoniot关系的等容燃烧解。     

A Numerical Study Promoting Algebraic Models for the Lewis Number Effect in Atmospheric Turbulent Premixed Bunsen Flames

This numerical investigation carried out on turbulent lean premixed flames accounts for two algebraic – the Lindstedt–Vaos (LV) and the classic Bray–Moss–Libby (BML) – reaction rate models. Computed data from these two models is compared with the experimental data of Kobayashi et al. on 40 different methane, ethylene and propane Bunsen flames at 1 bar, where the mean flame cone angle is used for comparison. Both models gave reasonable qualitative trend for the whole set of data, in overall. In order to characterize quantitatively, firstly, corrections are made by tuning the model parameters fitting to the experimental methane–air (of Le = 1.0) flame data. In case of the LV model, results obtained by adjusting the pre-constant, i.e., reaction rate parameter, C_R, from its original value 2.6 to 4.0, has proven to be in good agreement with the experiments. Similarly, for the BML model, with the tuning of the exponent n, in the wrinkling length scale, L_y = C_l⋅ l_x(s_L/u′)ⁿ from value unity to 1.2, the outcome is in accordance with the measured data. The deviation between the measured and calculated data sharply rises from methane to propane, i.e., with increasing Lewis number. It is deduced from the trends that the effect of Lewis number (for ethylene–air mixtures of Le = 1.2 and propane–air mixtures of Le = 1.62) is missing in both the models. The Lewis number of the fuel–air mixture is related to the laminar flame instabilities. Second, in order to quantify for its influence, the Lewis number effect is induced into both the models. It is found that by setting global reaction rate inversely proportional to the Lewis number in both the cases leads to a much better numerical prediction to this set of experimental flame data. Thus, by imparting an important phenomenon (the Lewis number effect) into the reaction rates, the generality of the two models is enhanced. However, functionality of the two models differs in predicting flame brush thickness, giving scope for further analysis.     

Aluminum Particles–air Detonation at Elevated Pressures

The effect of initial pressure on aluminum particles–air detonation was experimentally investigated in a 13 m long, 80 mm diameter tube for 100 nm and 2 μm spherical particles. While the 100 nm Al–air detonation propagates at 1 atm initial pressure in the tube, transition to the 2 μm aluminum–air detonation occurs only when the initial pressure is increased to 2.5 atm. The detonation wave manifests itself in a spinning wave structure. An increase in initial pressure increases the detonation sensitivity and reduces the detonation transition distance. Global analysis suggests that the tube diameter for single-head spinning detonation or characteristic detonation cell size would be proportional to (d ₀: aluminum particle size, p ₀: initial pressure). Its application to the experimental data results in m ~ O(1) and n ~ O(1) for 1 to 2 μm aluminum–air detonation, thus indicating a strong dependence on initial pressure and gas-phase kinetics for the aluminum reaction mechanism in detonation. Hence, combustion models based on the fuel droplet diffusion theory may not be adequate in describing micrometric aluminum–air detonation initiation, transition and propagation. For 2 μm aluminum–air mixtures at 2 atm initial pressure and below, experiments show a transition to a “dust quasi-detonation” that propagates quasi-steadily with a shock velocity deficit nearly 40% with respect to the theoretical C–J detonation value. The dust quasi- detonation wave can propagate in a tube with a diameter less than 0.4–0.5 times the diameter required for a spinning detonation wave.     

Investigation of detonation initiation in aluminium suspensions

Detonation initiation is investigated in aluminium/oxygen and aluminium/air mixtures. Critical conditions for initiation of spherical detonations are examined in analogy with the criteria defined for gaseous mixtures, which correlate critical parameters of detonation initiation to the characteristic size of the cellular structure. However, experimental data on the detonation cell size in these two-phase mixtures are very scarce, on account of the difficulty to perform large-scale experiments. Therefore, 2D numerical simulations of the detonation cellular structure have been undertaken, with the same combustion model for Al/air and Al/O₂ mixtures. The cell size is found to be λ = 37.5 cm for a rich (r = 1.61) aluminium–air mixture, and λ = 7.5 cm for a stoichiometric aluminium-oxygen mixture, which is in reasonable agreement with available experimental data. Calculations performed in large-scale configurations (up to 25 m in length and 1.5 m in lateral direction) suggest that the critical initiation energy and predetonation radius for direct initiation of the unconfined detonation in the aluminium–air mixture are, respectively, 10 kg of TNT and 8 m. Moreover, numerical simulations reveal that the structure of the detonation wave behind the leading front is even more complicated than in pure gaseous mixtures, due to two-phase flow effects. This paper is based on work that was presented at the 21st International Colloquium on the Dynamics of Explosions and Reactive Systems, Poitiers, France, July 23–27, 2007.     

Optimization of the deflagration to detonation transition: reduction of length and time of transition

The aim of this experimental investigation is the study of Deflagration to Detonation Transition (DDT) in tubes in order to (i) reduce both run-up distance and time of transition (L _DDT and t _DDT) in connection with Pulsed Detonation Engine applications and to (ii) attempt to scale L _DDT with λ_CJ (the detonation cellular structure width). In DDT, the production of turbulence during the long flame run-up can lead to L _DDT values of several meters. To shorten L _DDT, an experimental set-up is designed to quickly induce highly turbulent initial flow. It consists of a double chamber terminated with a perforated plate of high Blockage Ratio (BR) positioned at the beginning of a 26 mm inner diameter tube containing a “Shchelkin spiral” of BR ≈ 0.5. The study involves stoichiometric reactive mixtures of H₂, CH₄, C₃H₈, and C₂H₄ with oxygen and diluted with N₂ in order to obtain the same cell width λ_CJ≈10 mm at standard conditions. The results show that a shock-flame system propagating with nearly the isobaric speed of sound of combustion products, called the choking regime, is rapidly obtained. This experimental set-up allows a L _DDT below 40 cm for the mixtures used and a ratio L _DDT/λ_CJ ranging from 23 to 37. The transition distance seems to depend on the reduced activation energy (E _a/RT _c) and on the normalized heat of reaction (Q/a ₀ ²). The higher these quantities are, the shorter the ratio L _DDT/λ_CJ is. PACS 47.40.Rs · 47.60.+i · 47.70.Pq · 47.80.CbThis paper was based on the work that was presented at the 19th International Colloquium on the Dynamics of Explosions and Reactive Systems, Hakone, Japan, July 27–August 1, 2003.     

Tube diameter effect on deflagration-to-detonation transition of propane–oxygen mixtures

An experimental study was conducted to study the tube diameter effect on deflagration- to-detonation run-up distance. The tube diameter effect is associated with the amplification factor, flame acceleration and heat loss. A simplified correlation of the run-up distance and tube diameter is proposed for the fuel-lean, stoichiometric and slight fuel-rich mixtures. The amplification factor, which is evaluated from the initial conditions of the propane–oxygen mixtures, might also be used to get a quick estimation of the run-up distance in tubes of larger diameter.

---

     

Influence of mixture sensitivity and pore size on detonation velocities in porous media

Experiments have been carried out to determine the dependence of the detonation velocity in porous media, on mixture sensitivity and pore size. A detonation is established at the top end of a vertical tube and allowed to propagate to the bottom section housing the porous bed, comprised of alumina spheres of equal diameter (1–32 mm). Several of the common detonable fuels were tested at atmospheric initial pressure. Results indicate the existence of a continuous range of velocities with change in Φ, spanning the lean and the rich propagation limits. For all fuels in a given porous bed, the velocity decreases from a maximum value at the most sensitive mixture near Φ≈1 (minimum induction length), toV/V _CJ≈0.3 at the limits. A decrease in pore size brings about a reduction inV/V _CJ and a narrowing of the detonability range for each fuel. For porous media comprised of spherical particles, it was possible to correlate the velocity data corresponding to a variety of different mixtures and for a broad range of particle sizes, using the following empirical expression:V/V _CJ=[1–0.35 log(d _c /d _p)]±0.1. The critical tube diameterd _c is used as a measure of mixture sensitivity andd _p denotes the pore diameter. An examination of the phenomenon at the composition limits, suggests that wave failure is controlled by a turbulent quenching mechanism.     

Detonability of THDCPD-exo–air mixtures

In the frame of industrial risk and propulsive application, the detonability study of JP10–air mixtures was performed. The simulation and measurements of detonation parameters were performed for THDCPD-exo/air mixtures at various initial pressure (1 bar < P ₀ < 3 bar) and equivalence ratio (0.8 < Φ < 1.6) in a heated tube (T ₀ ~ 375 K). Numerical simulations of the detonation were performed with the STANJAN code and a detailed kinetic scheme of the combustion of THDCPD. The experimental study deals with the measurements of detonation velocity and cell size λ. The measured velocity is in a good agreement with the calculated theoretical values. The cell size measurements show a minimum value for Φ ~ 1.2 at every level of initial pressure studied and the calculated induction length L _i corresponds to cell size value with a coefficient k = λ/L _i = 24 at P ₀ = 1 bar. Based on the comparison between the results obtained during this study and those available in the literature on the critical initiation energy E _c, critical tube diameter d _c and deflagration to detonation transition length L _DDT, we can conclude that the detonability of THDCPD–air mixtures corresponds to that of hydrocarbon–air mixtures.

---

This paper is based on the work presented at the 33rd International Pyrotechnics Seminar, IPS 2006, Fort Collins, July 16–21, 2006.     

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.      

两种不同气体中的高速爆燃波及其向爆轰的转变

对2C2H2 5O2及2C2H2 5O2 80%Ar两种可燃混合气体中的高速爆燃波及其向爆轰的转变过程进行实验研究.高速爆燃波由孔栅干涉爆轰波的方法直接生成,观测手段则以高速转鼓摄影获取孔栅近场流场x-t纹影图,以传感器追踪波面的后继发展.研究发现,两种气体中的爆燃波具有迥异的特性.前者燃烧波面在较低初压条件下为层流结构,而较高初压下为湍流结构,向爆轰转变点可以延伸至下游较长距离;后者在不同初压条件下燃烧波面无明显差异,爆轰的再次形成只能在孔栅下游近场内建立.两种气体中高速爆燃波的维持和爆轰转变过程均非纯粹激波压缩所致,湍流输运在其中起着必不可少的作用.分析显示,激波压缩效应对纯氧炔气体的高速爆燃和DDT贡献较小,湍流输运占主导地位;而氩气稀释气体较为稳定,缺乏自行衍生剧烈湍流燃烧的能力,因而激波压缩和外界扰动对其高速爆燃传播和爆轰转变起十分重要的作用.     

Comparison of critical conditions for DDT in regular and irregular cellular detonation systems

Abstract. The results of an experimental study of DDT in mixtures with regular and irregular detonation cellular structures are presented. Experiments were carried out in a tube 174 mm i. d. with obstacles (blockage ratios were 0.1, 0.3, and 0.6). Mixtures used were hydrogen–air and stoichiometric hydrogen–oxygen diluted with , Ar, and He. The critical conditions for DDT are shown to depend on the regularity of the cellular structure of test mixtures. The critical values of the cell sizes in Ar- and He-diluted mixtures are shown to be significantly smaller than those in -diluted mixtures. This means that systems with a highly regular detonation cellular structure have far less capacity for undergoing DDT compared to irregular ones with the same values of detonation cell sizes. Received 18 November 1999 / Accepted 15 May 2000