Numerical study of the detonation structure in rich H2?NO2/N2O4 and very lean H2?N2O mixtures

In some mixtures and under certain conditions, detonation soot records show substructures. In nitromethane and nitrogen tetroxide mixtures, particular cellular structures can be observed. This kind of structures has been reported as the so-called double cellular structure. One- and two-dimensional simulations of detonation have shown that the double cellular structure is related to a non-monotonous energy release. Two-step energy release is also observed in rich H₂−NO₂/N₂O₄ and in very lean H₂−N₂O mixtures. The present study aims at the investigation of the effect of the energy release profile on the detonation structure in these two mixtures through numerical simulations. The origin of the non-monotonous energy release is explained in both mixtures using one-dimensional simulations with detailed chemistry. Reduced kinetic schemes are obtained and used to perform two-dimensional simulations. It is shown that in rich H₂−NO₂/N₂O₄ mixtures, the double cellular structure appears, whereas in very lean H₂−N₂O mixtures, classical substructures are observed. Both behaviours are explained based on ZND calculations and previous stability results. Phenomenological considerations led the authors to link the formation of the double cellular structure with the appearance of a large scale instability mode (a super cellular structure).     

Interaction of detonation with inert gas zone

The results of experimental study on detonation interaction with the regions of low reactivity, generated by the injection of an inert gas into reactive mixture, are reported. A square cross-section 60×60 mm, 3.6 m long detonation channel was used. The experiments were done for stoichiometric H₂−O₂ mixture at 0.3 bar and 0.5 bar initial pressure and room temperature. The inert gas (Ar, He, N₂ or CO₂) was injected from 0.523 dm³ container into the main channel 1 s before ignition. The size of the inert zone was controlled by inert initial pressure. The behavior of detonation was studied using direct streak photography and pressure transducers. The study has shown that at low pressure of Ar, N₂ and CO₂ injection only a slight decrease of detonation velocity occurs. At higher injection pressures complete damping of detonation and flame extinguishment occur, followed by flame reigniton and DDT outside the inert zone. For low He injection pressures an increase in detonation velocity was recorded. For higher injection pressures, detonation damping occurred, followed by DDT process. The results have shown that CO₂ has the strongest effect on damping 2H₂+O₂ detonation, with N₂ and Ar in the next places, and He very far behind. The effectiveness of inert gas in detonation damping was found proportional to its molecular weight and reciprocal to its specific heat ratio. The numerical simulations of detonation propagation through inert gas zone were also performed using the one- dimensional Detonation Lagrangean code with simple energy release model. The results of simulations are in good qualitative agreement with experimental tendencies. In particular, the model has shown that the re-initiation of detonation is enhanced by smooth concentration gradients at inert/reactive interface. An abridged version of this paper was presented at the 15th Int. Colloquium on the Dynamics of Explosions and Reactive Systems at Boulder, Colorado, from July 30 to August 4, 1995     

Turbulent jet initiation of detonation in hydrogen-air mixtures

Large scale experiments (50 m³) have been carried out on the initiation of detonation by means of a jet of hot combustion products. The effects of hydrogen concentration (18–30% vol.), jet orifice diameter (100–400 mm), and the mixture composition in constant volume explosion chamber (25–50%) were investigated. Both high enough hydrogen concentration and large enough jet size are necessary for detonation initiation. The minimum values are within the ranges of 20 to 25% vol. H₂, and of 100 to 200 mm correspondingly. A minimum ratio of jet size and mixture cell width 12–13 is required for detonation initiation.     

Detonation wave initiated by explosive condensation of supersaturated carbon vapor

An experimental study of the influence of condensation of supersaturated carbon vapor formed behind reflected shock waves on the process of propagation of a shock wave and formation of a detonation wave of condensation is carried out. Highly supersaturated carbon vapor was formed from thermal decay of unstable carbon suboxide C₃O₂ → C + 2CO behind a shock wave in mixtures containing 10–30% C₃O₂ in Ar. This reaction was followed by fast growth of condensed carbon particles, accompanied by heat release. Experiments have shown a considerable temperature and pressure increase in the narrow zone behind the wave front, resulting in shock wave amplification and transition to a detonation-like regime. An analysis of the kinetics and heat release in the given conditions and calculations based upon one-dimensional detonation theory have shown that in a mixture of 10% C₃O₂ + Ar, insufficient heat release resulted in a regime of “overdriven detonation”. In a mixture of 20% C₃O₂ + Ar a very good coincidence of measured values of pressure and wave velocity with calculated Chapman–Jouguet parameters is observed. In a 30% C₃O₂ + Ar mixture, an excess heat release caused a slow down of the effective condensation rate and a regime of “underdriven detonation” is observed.     

A generalized Zel’dovich-Neumann-Döring model of steady gaseous detonation

A model of N-component gas mixture detonation whose combustion zone is formed by n (1 ≤ n ≤ N − 1) independent reactions satisfying the kinetic formulation of the mass action law is studied. In the framework of this model, it is found that there exist two possible types of self-sustaining detonation waves with asymptotically stable combustion zones. It is shown that it is possible to observe only a strong wave propagating in the Chapman-Jouguet regime relative to the flow on the critical surface. It is also shown that, for a perfect gas mixture, this wave appears only when a combustion zone is formed by one or two reactions.     

Gaseous nitromethane and nitromethane-oxygen mixtures: A new detonation structure

Detonation in gaseous nitromethane (NM) and mixed with O₂ has been studied. Experiments were performed in a preheated steel tube at an initial temperatureT ₀∼=390 K for different initial pressuresP ₀ (1.7≥P ₀≥5 10⁻² bar). Different selfsustained detonation regimes were obtained, from multiheaded mode to spinning and galloping mode in marginal conditions. These chemical systems were characterized by a specific detonation cellular structure very different from that currently observed with classical gaseous C _n H _m /O₂/N₂ mixtures. All modes of detonation propagation in rich NM/O₂ mixtures exhibit a double scale cellular structure. The pattern of this double scale structure is particularly clear in the case of spinning mode. An abridged version of this paper was presented at the 15th Int. Colloquium on the Dynamics of Explosions and Reactive Systems at Boulder, Colorado, from July 30 to August 4, 1995     

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.     

Two-cell detonation: losses effects on cellular structure and propagation in rich H2–NO2/N2O4–Ar mixtures

Detonation experiments in H₂–NO₂/N₂O₄–Ar mixtures (Equivalence ratio 1.2 and initial pressure lower than 0.1 MPa) confined in a tube of internal diameter 52 mm reveal two propagation regimes depending on initial pressure: (1) a quasi-CJ regime is observed along with a double cellular structure at high pressures; (2) at lower pressures, a low velocity detonation regime is observed with a single structure. Transition between this two regimes happens when the spinning detonation of the larger cell vanishes. Each detonation regime is characterized by velocity and pressure measurements and cellular structure records. Coherence between all experimental data for each experiment leads in assumption that losses are responsible for the transition between one regime to another. In a second part, we study such behaviour for a two-step mixture through numerical simulations using a global two-step chemical kinetics and a simple losses model. Numerical simulations qualitatively agree with experiments. Both detonation regimes with their own cellular structures are reproduced.     

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

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.     

Experimental observation of the onset of detonation downstream of a perforated plate

In this study, the onset of detonation downstream of a perforated plate subsequent to the reflection of a Chapman–Jouguet detonation upstream is investigated. The experiments were performed with C₃H₈ + 5O₂ and C₂H₂+2.5O₂+70%Ar. The former has a much more irregular transverse wave pattern whereas the latter is known to have a piecewise laminar structure with a regular cellular structure. The onset of detonation phenomenon was found to be significantly different for the two mixtures. For the high argon diluted mixtures, the onset of detonation occurs in the vicinity downstream of the perforated plate. However, if the onset of detonation does not occur close to the plate, the precursor shock decouples from the reaction zone and a deflagration results. For the propane–oxygen mixtures, the onset of detonation is found to occur relatively far from the perforated plate at critical conditions. The major difference between these two mixtures is that the metastable turbulent reaction front can be maintained for relatively long distances for the propane–oxygen mixture. This turbulent metastable regime is also observed to be able to maintain a relatively constant propagation velocity for many channel widths prior to the onset of detonation. For the propane–oxygen mixtures, the onset is caused by a strong local explosion within the turbulent reaction zone.     

Population inversion behind a detonation wave propagating in a variable-density medium

The creation of an active medium by means of detonation has been investigated on a number of occasions. It was suggested that one could use the expansion of the detonation products of an acetylene-air mixture in vacuum [1] or the cooling of the detonation products of a mixture of hydrocarbons and air through a nozzle [2, 3]. In [4], the detonation of a solid high explosive was used to produce population inversion in the gas mixture CO₂-N₂-He(H₂O). Stimulated emission from HF molecules was observed in [5] behind the front of an overdriven detonation wave propagating in an F₂-H₂-Ar mixture in a shock tube. Population inversion behind a detonation wave was studied in H₂-F₂-He mixtures in [6–8] and in H₂-Cl₂-He mixtures in [9] with energy release on a plane and on a straight line in a medium with constant density. Similar problems were solved for shock waves propagating in both a homogeneous gaseous medium [7, 10] and in the supersonic part of an expanding nozzle. In the present paper, we study theoretically population inversion behind an overdriven detonation wave propagating in a mixture (fine carbon particles + acetylene + air) which flows through a hypersonic nozzle. The propagation of detonation in media with variable density and initial velocity was considered, for example, in [11, 12]. Analysis of the gas parameters behind a detonation wave propagating in a medium with constant density (for a given fuel) showed that the temperature difference across the detonation front is insufficient to produce population inversion of the vibrational levels of the CO₂ molecule.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 65–71, January–February, 1980.I am grateful to V. P. Korobeinikov for a helpful discussion of the results.     

Maximum density effects on convective instability of horizontal plane poiseuille flows in the thermal entrance region

A linear stability analysis is used to study the conditions marking the onset of secondary flow in the form of longitudinal vortices for plane Poiseuille flow of water in the thermal entrance region of a horizontal parallel-plate channel by a numerical method. The water temperature range under consideration is 0∼30°C and the maximum density effect at 4°C is of primary interest. The basic flow solution for temperature includes axial heat conduction effect and the entrance temperature is taken to be uniform at far upstream location jackie=−∞ to allow for the upstream heat penetration through thermal entrance jackie=0. Numerical results for critical Rayleigh number are obtained for Peclet numbers 1, 10, 50 and thermal condition parameters (λ ₁, λ ₂) in the range of −2.0≤λ ₁≤−0.5 and −1.0≤λ ₂≤1.4. The analysis is motivated by a desire to determine the free convection effect on freezing or thawing in channel flow of water.     

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.     

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.     

ON THE CRITICAL POINTS OF THE MAP Fp:X→‖AXB-C‖pp

Suppose A,B and C are the bounded linear operators on a Hilbert space H, when A has a generalized inverse A^- such that (AA^-)^*=AA^- and B has a generalized inverse B^- such that (B^-B)^*=B^-B,the general characteristic forms for the critical points of the map F_p:X^→‖AXB-C‖_p^p(1p=2. Similarly, the same question has been discussed for several operators.     

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.     

The optimal dimensions of convective-radiating circular fins

The optimal dimensions of convective-radiating circular fins with variable profile, heat-transfer coefficient and thermal conductivity, as well as internal heat generation are obtained. A profile of the form y=(w/2) [1+(r _o/r) ⁿ ] is studied, while variation of thermal conductivity is of the form k=k _o[1+ɛ((T−T _∞)/ (T _b−T _∞)) ^m ]. The heat-transfer coefficient is assumed to vary according to a power law with distance from the bore, expressed as h=K[(r−r _o)/(r _e−r _o)]^λ. The results for λ=0 to λ=1.9, and −0.4≤ɛ≤0.4, have been expressed by suitable dimensionless parameters. A correlation for the optimal dimensions of a constant and variable profile fins is presented in terms of reduced heat-transfer rate. It is found that a (quadratic) hyperbolic circular fin with n=2 gives an optimum performance. The effect of radiation on the fin performance is found to be considerable for fins operating at higher base temperatures, whereas the effect of variable thermal conductivity on the optimal dimensions is negligible for the variable profile fin. It is also observed, in general, that the optimal fin length and the optimal fin base thickness are greater when compared to constant fin thickness. Received on 22 February 1999     

Bentonite slurries—zeta potential,yield stress,adsorbed additive and time-dependent behaviour

Bentonite clay is a vital ingredient of drilling mud. The time-dependent and high shear thinning yield stress behaviour of drilling mud is essential for maintaining wellbore stability and to remove cuttings, cool and clean the drill bit of debris. As-prepared 3, 5 and 7 wt.% bentonite clay slurries displayed time-dependent behaviour where the yield stress (measured after quick stirring) decreased with time of rest. An equilibrium value is reached after 24 h. Despite the low solids concentration, the yield stress is already relatively high and is displayed at all pH level. The yield stress is maximum at pH 2 and minimum at pH ∼ 7. This yield stress is due to the formation of gel structure by the swelling clay particles. However the addition of phosphate additives such as (PO₃)^19 −, (P₃O₁₀)^5 − and (P₂O₇)^4 − completely dispersed the clay slurries at pH above 6. At pH below 6, yield stress is still present but is 3-folds smaller than that of the pure bentonite slurry. With phosphate additives, the magnitude of the critical zeta potential at the complete dispersion pH is ca 48 mV. However for the pure bentonite, the slurry remained flocculated at zeta potential of >50 mV in magnitude. Interestingly, (P₂O₇)^4 − anions is more effective than the other two phosphate additives in reducing the yield stress at low pH, ∼ 2.0.