A nonlinear oscillator concept for one-dimensional pulsating detonations

To interpret the results of direct numerical simulations for the one-dimensional pulsating detonation, a nonlinear oscillator model is proposed based on the integral conservation considerations. A gauging procedure is suggested in which the results from the direct numerical simulation are used to compute the coefficients of the oscillator equation. Various terms of the oscillator equation obtained are compared with those of mechanical nonlinear oscillators in order to illustrate their analogy. Among many important features captured by this nonlinear oscillator equation, the oscillatory behavior of the detonation wave front can be interpreted as resonant excitation of the chemical energy release. Received 24 March 1997 / Accepted 7 May 1998     

The influence of detonation cell size and regularity on the propagation of gaseous detonations in granular materials

This paper presents results from an experimental study of transmission of gaseous detonation waves through various granular filters. Spherical glass beads of 4 and 8 mm diameter and crushed rock of 7.5 mm volume averaged diameter were used as filter material. Varying the initial pressure of the detonating gas mixture controlled the detonation cell size. Dilution with argon was used to vary the detonation cell regularity. The complete range from almost no detonation velocity deficit to complete extinction of the combustion wave was observed. The existing correlation for gaseous detonation velocity deficit where is the critical diameter for the gaseous detonation and is the pore size, is found to be applicable for both smooth spherical particles and irregular crushed rock when considering irregular detonation structures. Soot films and pressure measurements show that as the detonation cell size is increased, reinitiation of a reanular filter until it finally no longer occurs at . Complete extinction of the combustion wave occurs at . These two limits appear to be about the same for irregular and regular detonation cell structures. For irregular structures without argon dilution, can be found for detonation wave failure, and can be found for complete extinction of the combustion wave. For argon dilution these limits are changed to and , respectively. The data are a bit scarce as a basis for proposing a new correlation for regular structures, but as a first approximation log is suggested for regular structures. The detonation or combustion wave is found to approach a constant velocity in the granular filter if not extinguished. Received 31 October 2001 / Accepted 15 July 2002 Published online 4 November 2002 Correspondence to: T. Slungaard (e-mail: slung@maskin.ntnu.no) 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     

Direct initiation of detonations in unconfined gasoline sprays

Large scale experiments on detonation initiation in gasoline-air by two different sources were carried out at stoichiometric conditions. Unconfined clouds of volume generated by a special facility had a shape of semicylinder of 15–17 m in length and 6–8 m in radius. Both the charge of condensed HE and the charge of stoichiometric propane-air were used to initiate detonation in the mixture. In case of initiation by a propane-air charge the critical initiation energy was up to 7 times as large as that for HE initiation. The detonation cell size for gasoline-air was determined as 0.04–0.05 m. It was shown, that the well-known correlation between the critical energy of point blast initiation and the cell size failed for this system. The cell size obtained is close to one of propane-air, but no direct transfer of detonation from one mixture to another was observed. Received 10 March 1995 / Accepted 12 March 1995     

Recent work on gaseous detonations

The paper reviews recent progress in the field of gaseous detonations, with sections on shock diffraction and reflection, the transition to detonation, hybrid, spherically-imploding, and galloping and stuttering fronts, their structure, their transmission and quenching by additives, the critical energy for initiation and detonation of more unusual fuels. The final section points out areas where our understanding is still far from being complete and contains some suggestions of ways in which progress might be made. Received 9 September 1999 / Accepted 10 May 2001     

Experimental and numerical studies on transmission of gaseous detonation to a less sensitive mixture

A phenomenon of detonation transmission from one gaseous mixture (donor) to another of lower sensitivity (acceptor) was studied experimentally and numerically. The objective was to study effects of a donor mixture length and acceptor mixture sensitivity on the possibility of detonation transmission. Experiments were carried out in detonation tube 9.5–12 m long and 174 mm id. Three types of donor mixtures were used in the driver: stoichiometric acetylene/air, stoichiometric hydrogen/air, and 20% of hydrogen/air. Air mixtures with 14–29.6% of hydrogen were used as acceptors. Driver length varied from 0.17 to 5.6 m. Detonation transmission was studied for an abrupt opening of interface between two mixtures. Series of 1D and 2D calculations were made to simulate the problem numerically. Both, results of experiments and calculations revealed a set of parameters that effect transmission process. Critical conditions were determined as minimum driver length expressed in terms of characteristic chemical reaction length scales of acceptor mixture. They were shown to depend on differences in reaction rates and energy contents of donor and acceptor mixture. Received 6 January 1997 / Accepted 20 March 1997     

Experimental study on detonation parameters and cellular structures of fuel cloud

In this paper,detonation parameters of fuel cloud,such as propylene oxide(PO),isopropyl nitrate(IPN),hexane,90 # oil and decane were measured in a self-designed and constructed vertical shock tube.Results show that the detonation pressure and velocity of PO increase to a peak value and then decrease smoothly with increasing equivalence ratio.Several nitrate sensitizers were added into PO to make fuel mixtures,and test results indicated that the additives can efficiently enhance detonation velocity and pressure of fuel cloud and one type of additive n-propyl nitrate(NPN) played the best in the improvement.The critical initiation energy that directly initiated detonation of all the test liquid fuel clouds showed a U-shape curve relationship with equivalence ratios.The optimum concentration lies on the rich-fuel side(φ > 1).The critical initiation energy is closely related to molecular structure and volatility of fuels.IPN and PO have similar critical values while that of alkanes are larger.Detonation cell sizes of PO were respectively investigated at 25 C,35 C and 50 C with smoked foil technique.The cell width shows a U-shape curve relationship with equivalence ratios at all temperatures.The minimal cell width also lies on the rich-fuel side(φ > 1).The cell width of PO vapor is slightly larger than that of PO cloud.Therefore,the detonation reaction of PO at normal temperature is controlled by gas phase reaction.     

On detonation behavior of mixed fuels

An experimental study of detonation was carried out for mixed fuels (, and ). The results obtained demonstrate that detonation velocities of mixed fuels agree very well with normal C-J velocity and that detonation limits were not extended by the addition of a small amount of . As for mixed fuel, detonation velocities and detonation limits as a function of the equivalence ratio of the whole mixture of mixed fuel coincided with those of single fuel. However, the detonation induction time of a mixture of with a lower sensitivity to detonation decreased considerably with the addition of a small amount of . In addition, the influence of mixed fuels on the quasi-detonation and the fast deflagration in the process of DDT was investigated, showing interesting and complicated features in mixed fuel. Received 20 December 1994 / Accepted 20 November 1995     

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     

Analytical study of idealized two-dimensional cellular detonations

In this study, the idealized two-dimensional detonation cells were decomposed into the primary units referred to as sub-cells. Based on the theory of oblique shock waves, an analytical formula was derived to describe the relation between the Mach number ratio through triple-shock collision and the geometric properties of the cell. By applying a modified blast wave theory, an analytical model was developed to predict the propagation of detonation waves along the cell. The calculated results show that detonation wave is, first, strengthened at the beginning of the cell after triple-shock collision, and then decays till reaching the cell end. The analytical results were compared with experimental data and previous numerical results; the agreement between them appears to be good, in general. Received 13 February 2001 / Accepted 2 August 2001     

Critical reactions for the hydrazine vapor detonations

A detailed kinetic model devoted to the hydrazine vapor detonation has been built. It consists of 33 reversible reactions and 13 species. A reduced kinetic model has been proposed by using the Principal Component Analysis of matrix F (PCAF) method as implemented in KINALC. It is constituted of 26 reactions and 11 species. This model has been shown to be valid over a pressure range of 0.1 to 10 atm. However, the predictions of the models are significantly affected by changes in the enthalpy of formation of NH. With the help of the full kinetic model, a value of A, the proportionality factor in the ZND model between the induction distance in the detonation wave and the detonation cell size, of has been derived if one considers that the collision efficiency of NH on the thermal decomposition of hydrazine is equal to the one of N. The value of A for pure hydrazine detonation is shown to be strongly dependent on the value of the collision efficiency of NH. Received 24 December 2000 / Accepted 11 April 2001     

Effect of reflection type on detonation initiation at shock-wave focusing

Abstract. From practical and theoretical standpoints, the initiation of combustion in gaseous media due to the shock waves focusing process at various reflectors is a subject of much current interest. The complex gas flowfield coupled with chemical kinetics provides a wide spectrum of possible regimes of combustion, such as fast flames, deflagration, detonation etc. Shock wave reflection at concave surfaces or wedges causes converging of the flow and produces local zones with extremely high pressures and temperatures. The present work deals with the initiation of detonation due to shock waves focusing at parabolic and wedge reflectors. Particular attention has been given to the determination of the critical values of the incident shock wave (ISW) Mach number, parameters of the combustible mixture, and geometrical sizes of reflector at which different combustion regimes could be generated. Received 30 August 1999 / Accepted 23 February 2000     

Pressure profiles in detonation cells with rectangular and diagonal structures

Experimental results presented in this work enable us to classify the three-dimensional structure of the detonation into two fundamental types: a rectangular structure and a diagonal structure. The rectangular structure is well documented in the literature and consists of orthogonal waves travelling independently from each another. The soot record in this case shows the classical diamond detonation cell exhibiting ‘slapping waves’. The experiments indicate that the diagonal structure is a structure with the triple point intersections moving along the diagonal line of the tube cross section. The axes of the transverse waves are canted at 45 degrees to the wall, accounting for the lack of slapping waves. It is possible to reproduce these diagonal structures by appropriately controlling the experimental ignition procedure. The characteristics of the diagonal structure show some similarities with detonation structure in round tube. Pressure measurements recorded along the central axis of the cellular structure show a series of pressure peaks, depending on the type of structure and the position inside the detonation cell. Pressure profiles measured for the whole length of the two types of detonation cells show that the intensity of the shock front is higher and the length of the detonation cell is shorter for the diagonal structures. Received 17 May 2000 / Accepted 29 November 2000     

Dynamics study of detonation-wave cellular structure 1. Statistical properties of detonation wave front

We have investigated the evolution of cellular detonation-wave structure as a gaseous detonation travels along a round tube and measured cell lengths as a function of the initial pressure of the gas. We have tested acetylene-containing combustible gas mixtures with different degrees of regularity. Along with the smoked-foil technique, an emission method has been used to the measure current and average values of the detonation cell length. The method is based on the detection of an emission spectrum behind the detonation front in the spectral range corresponding to local gas temperatures that are much higher than those for the Chapman-Jouguet equilibrium condition. This technique provides quasi-continuous cell-length measurements along the normal to the detonation front over the length of several factors of ten times the tube. Our study has experimentally identified the steady states of detonation structure in round tubes, referred to here as the single detonation modes. When the state of a single mode is fully established, then both the flow structure and the energy release at detonation front develop strictly periodically along the tube at a constant frequency inversely proportional to the cell length of the mixture. The mixture regularity has had no influence on the occurrence of the detonation mode, which is defined by the value of initial pressure or the total energy release of the mixture. Outside of the pressure range where a detonation mode was most likely to occur, the detonation front is unstable and may exhibit an irregular cellular pattern. Monitoring the evolution of cells over a long distance revealed that the local gas emissivity, which is time dependent and corresponds to axial pulsations of the detonation structure, has the appearance of a superposition of separate harmonics describing the states of emissivity oscillations and cell structure of single detonation modes. Received 18 October 1999 / Accepted 10 June 2001     

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     

Numerical modelling of detonation structure in two-phase flows

A one-dimensional physical model and a numerical method for the simulation of heterogeneous detonation were proposed based on an Eulerian approach for heterogeneous flows. The combination of modern shock-capturing schemes in combination with a dynamically moving, adaptive grid ensure the properresolution of both reaction zones and flow discontinuities. Numerical examples illustrate the effect of the heat release due to heterogeneous combustion. Received August 4, 1995 / Accepted December 12, 1995     

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     

Non-equilibrium model of steady detonations in aluminum particles - oxygen suspensions

Steady-state detonation regimes are studied on the basis of the mathematical model of detonation of aluminum particles in oxygen taking into account differences in velocities and temperatures of the mixture components. The final steady state is analyzed by determining the types of final singularities in the plane of relaxation parameters (the ratios of characteristic times of thermal and velocity relaxations and combustion). The regions of existence of steady-state regimes are found numerically, depending on the detonation wave velocity and relaxation parameters. Numerical illustrations of various flow types are presented, and the properties of the detonation wave structure caused by velocity nonequilibrium are examined. Qualitative agreement of data obtained with frozen relaxation parameters and their dependence on the flow parameters is shown. Received 5 July 1997 / Accepted 13 July 1998     

Numerical predictions of oblique detonation stability boundaries

Oblique detonation stability was studied by numerically integrating the two-dimensional, one-step reactive Euler equations in a generalized, curvilinear coordinate system. The integration was accomplished using the Roe scheme combined with fractional stepping; nonlinear flux limiting was used to prevent unphysical solution oscillations near discontinuities. The method was verified on one- and two-dimensional flows with exact solutions, and its ability to correctly predict one-dimensional detonation stability boundaries was demonstrated. The behavior of straight oblique detonations attached to curved walls was then considered. Using the exact, steady oblique detonation solution as an initial condition, the numerical simulation predicted both steady and unsteady oblique detonation solutions when a detonation parameter known as the normal overdrive, , was varied. For a standard test case with a specific heat ratio of , a dimensionless activation energy of , and dimensionless heat release of , an oblique detonation with a constant dimensionless component of velocity tangent to the lead shock, , underwent transition to unstable behavior at . This is slightly higher than the threshold of predicted by one-dimensional theory; thus, the two-dimensionality renders the flow slightly more susceptible to instability. Received 4 August 1996 / Accepted 5 March 1996     

Some observations of the controlled generation and onset of detonation

Detailed observation of deflagration to detonation transition (DDT) is inherently difficult. This is primarily due to the stochastic nature of flame acceleration and shock formation processes that in most practical situations give rise to the conditions required for detonation to emerge. The present paper describes how shock tube techniques have been used to control the conditions required for the onset of detonation. The paper first outlines some initial experiments involving turbulent flame acceleration before concentrating on experiments in a reflected shock mode. To aid interpretation of the observations the paper also presents a simple gasdynamic analysis of particle trajectories and considers the various physical and chemical processes that could lead to the onset of detonation. Received 27 November 2001 / Accepted 28 January 2002