On the oscillating flame characteristics in nonpremixed laminar coflow-jets: An experimental and numerical study

This study investigates the characteristics of oscillating lifted flames in laminar coflow-jets experimentally and numerically by varying both fuel density (by varying propane and n-butane mixtures) and coflow density (by diluting air with N₂/He mixtures). Two different lifted flame oscillation behaviors are observed depending on these parameters: oscillating tribrachial lifted flame (OTLF) and oscillating mode-change lifted flame (OMLF), where a rapid increase in flame radius is observed. The regimes of the two flames are identified from experiments, which shows that OMLF occurs only when the effect of the negative buoyancy on the flow field by the fuel heavier than air becomes significant at low fuel jet velocity. OMLFs are also identified to distinguish OTLF regime from flame extinction, which implies that an OMLF can be extinguished when the positive buoyancy becomes weak, losing its stabilizing effect, or when the negative buoyancy becomes strong, further enhancing its destabilizing effect. Transient numerical simulations of both OTLF and OMLF reveal that the OMLF occurs by a strong toroidal vortex and a subsequent counterflow-like structure induced by relatively-strong negative buoyancy. Such a drastic flow redirection significantly changes the fuel concentration gradient such that the OMLF changes its mode from a tribrachial flame mode (decreasing edge speed with fuel concentration gradient) to the premixed flame-like transition mode when the fuel concentration gradient becomes very small (increasing edge speed with fuel concentration gradient). Again, a tribrachial flame mode is recovered during a rising period of flame edge and repeats an oscillation cycle.     

Stabilization mechanism of lifted flame edge in the near field of coflow jets for diluted methane

The stabilization mechanism of lifted flames in the near field of coflow jets has been investigated experimentally and numerically for methane fuel diluted with nitrogen. The lifted flames were observed only in the near field of coflow jets until blowout occurred in the normal gravity condition. To elucidate the stabilization mechanism for the stationary lifted flames of methane having the Schmidt number smaller than unity, the behavior of the flame in the buoyancy-free condition, and unsteady propagation characteristics after ignition were investigated numerically at various conditions of jet velocity. It has been found that buoyancy plays an important role for flame stabilization of lifted flames under normal gravity, such that the flame becomes attached to the nozzle in microgravity. The stabilization mechanism is found to be due to the variation of the propagation speed of the lifted flame edge with axial distance from the nozzle in the near field of the coflow as compared to the local flow velocity variation at the edge.     

Effects of fuel Lewis number on flame spread over solids

Using a detailed two-dimensional numerical model, a systematic investigation has been made to study the effect of fuel Lewis number (Le_F = α/D_F) and mass transfer on flame spread over thin solids. The fuel Lewis number affects the flame spread rates for both concurrent and opposed flames over thin fuels. The dependence of the flame spread rate on Le_F for these two spreading modes is, however, not the same. In opposed flame spreads (zero-gravity, self-propagation, and normal gravity downward propagation), the flame spread rate vs. Le_F curve is non-monotonic with a maximum value occurring at an intermediate value of Le_F = 0.5. In steady, concurrent spread in zero-gravity with low-speed flow and a constant flame length, the flame spread rate decreases with Le_F in a monotonic manner. By using the computational model as a tool, the effects of fuel mass diffusion perpendicular to and parallel with the solid surface are isolated to obtain more physical insight on the two-dimensional aspect of fuel mass transfer on flame spread. In addition, the model has also been used to decouple the solid evaporation process so that the fuel diffusion effect in the gas-phase can be isolated. Both of these theoretical exercises contribute to the understanding of mass transfer effects on the flame spreading phenomena over solids.     

Air/fuel mixing in jet flames

The paper examines eight diverse regimes in which fuels can mix and react with air. These comprise: (i) Lifted subsonic; and (ii) supersonic jet flames, with (iii) and without (iv) cross flows; (v) Rim-attached flames; (vi) Early Downwash flames; (vii) Downwash-attached jet flames; and (viii) Fire Whirls.Correlations of characteristics within these regimes are principally in terms of a dimensionless Flow Number, U*, Cross Flow Reynolds number, Re_c, and, for Fire Whirls, a dimensionless Critical Velocity, CV. Boundaries of seven of the eight regimes are identified, through plots of U*, against Re_c, and of the eighth through a plot of CV against U*. The circumstances of transitions between regimes are identified. The study involves a variety of CH₄ cross flow flame measurements, in a wind tunnel. Cross flows can initially create a small lee-side flame downwash, due to the depression in pressure. With increasing fuel flow this might extend 1.3 m downwards from the horizontal tip of the vertical burner. Jet flames can attach to the downwash, which can become significant above Re_c ≈ 2000. More extensive downwash might further delay blow-off. Regime boundaries are constructed on the U*/Re_c diagram covering lifted flames, early downwash, and downwash-attached flames. The most powerful flames tend to be lifted, choked, flames, with cross flow, and fire whirls. Combustion becomes less efficient at high Re_c and low U*, although CH₄ was efficiently reacted.Experimental values of the ratio of fuel to air velocity, u/u_c, of CH₄ flames ranged between about 10 and 30 for lifted flames, and between 0.3 and 3.6, at blow-off, for rim-attached flames. The latter comprise an important category, often intermediate between lifted flames and downwash-attached flames.     

Simultaneous stereo‐PIV and OH×CH2O PLIF measurements in turbulent ultra lean CH4/H2 swirling wall‐impinging flames

We present experimental results from turbulent low-swirl lean H₂/CH₄ flames impinging on an inclined, cooled iso-thermal wall, based on simultaneous stereo-PIV and OH×CH₂O PLIF measurements. By increasing the H₂ fraction in the fuel while keeping Karlovitz number (Ka) fixed in a first series of flames, a fuel dependent near-wall flame structure is identified. Although Ka is constant, flames with high H₂ fraction exhibit significantly more broken reaction zones. In addition, these high H₂ fraction flames interact significantly more with the wall, stabilizing through the inner shear layer and well inside the near-wall swirling flow due to a higher resistance to mean strain rate. This flame-wall interaction is argued to increase the effective local Ka due to heat loss to the wall, as similar flames with a (near adiabatic) ceramic wall instead of a cooled wall exhibit significantly less flame brokenness. A second series of leaner flames were investigated near blow-off limit and showed complete quenching in the inner shear layer, where the mean strain rate matches the extinction strain rate extracted from 1D flames. For pure CH₄ flames (Ka ≈ 30), the reaction zone remains thin up to the quenching point, while conversely for the 70% H₂ flames (Ka ≈ 1100), the reaction zone is highly fragmented. Remarkably, in all near blow-off cases with CH₄ in the fuel, a large cloud of CH₂O persists downstream the quenching point, suggesting incomplete combustion. Finally, ultra lean pure hydrogen flames were also studied for equivalence ratios as low as 0.22, and through OH imaging, exhibit a clear transition from a cellular flame structure to a highly fragmented flame structure near blow-off.     

Control of oscillating combustion and noise based on local flame structure

To control combustion oscillations, the characteristics of an oscillating swirl injection premixed flame have been investigated, and control of oscillating combustion and noise based on local flame structure has been conducted. The r.m.s. value of pressure fluctuations and noise level show significantly large values between = 0.8 and 1.1. The beating of pressure fluctuations is observed for the large oscillating flame conditions in this combustor. Relationship between beating of pressure fluctuations and local flame structure was observed by the simultaneous measurement of CH/OH planar laser induced fluorescence and pressure fluctuations. The local flame structure and beating of pressure fluctuations are related and the most complicated flame is formed in the middle pressure fluctuating region of beating. The beating of pressure fluctuations, which plays important roles in noise generation and nitric oxide emission in this combustor, could be controlled by injecting secondary fuel into the recirculating region of oscillating flames. Injecting secondary fuel prevented lean blowout, and low NO_x combustion was also achieved even for the case of pure methane injection as a secondary fuel. By injecting secondary fuel into the recirculating region near the swirl injector, the flame lifted from the swirl injector and its reaction region became uniform and widespread, hence resulting in low nitric oxide emission. Secondary mixture injection, fuel diluted with air, is not effective for control of combustion oscillations suppression and lean blowout prevention.     

Stabilization, propagation and instability of tribrachial triple flames

A tribrachial (or triple) flame is one kind of edge flame that can be encountered in nonpremixed mixing layers, consisting of a lean and a rich premixed flame wing together with a trailing diffusion flame all extending from a single point. The flame could play an important role on the characteristics of various flame behaviors including lifted flames in jets, flame propagation in two-dimensional mixing layers, and autoignition fronts. The structure of tribrachial flame suggests that the edge is located along the stoichiometric contour in a mixing layer due to the coexistence of all three different types of flames. Since the edge has a premixed nature, it has unique propagation characteristics. In this review, the propagation speed of tribrachial flames will be discussed for flames propagating in mixing layers, including the effects of concentration gradient, velocity gradient, and burnt gas expansion. Based on the tribrachial edge structure observed experimentally in laminar lifted flames in jets, the flame stabilization characteristics including liftoff height, reattachment, and blowout behaviors and their buoyancy-induced instability will be explained. Various effects on liftoff heights in both free and coflow jets including jet velocity, the Schmidt number of fuel, nozzle diameter, partial premixing of air to fuel, and inert dilution to fuel are discussed. Implications of edge flames in the modeling of turbulent nonpremixed flames and the stabilization of turbulent lifted flames in jets are covered.     

Decreasing liftoff height behavior in diluted laminar lifted methane jet flames

Stabilization of laminar lifted coflow jet flames of nitrogen-diluted methane was investigated experimentally and numerically. As the fuel jet velocity was increased, two distinct behaviors in liftoff height were observed depending on the initial fuel mole fraction; a monotonically increasing trend and a decreasing and then increasing trend (U-shaped behavior). The former was observed in the jet-developing region and the latter in the jet-developed region. Because the decreasing behavior of liftoff height with jet velocity has not been observed at ambient temperature, the present study focuses on decreasing liftoff height behavior. To elucidate the physical mechanism underlying the U-shaped behavior, numerical simulations of reacting jets were conducted by adopting a skeletal mechanism. The U-shaped behavior was related to the buoyancy. At small jet velocities, the relative importance of the buoyancy over convection was strong and the flow field was accelerated in the downstream region to stabilize the lifted flame. As the jet velocity increased, the relative importance of buoyancy decreased and the liftoff height decreased. As the jet velocity further increased, the flame stabilization was controlled by jet momentum and the liftoff height increased.     

Lagrangian investigation of local extinction,re-ignition and auto-ignition in turbulent flames

Lagrangian PDF investigations are performed of the Sandia piloted flame E and the Cabra H₂/N₂ lifted flame to help develop a deeper understanding of local extinction, re-ignition and auto-ignition in these flames, and of the PDF models' abilities to represent these phenomena. Lagrangian particle time series are extracted from the PDF model calculations and are analyzed. In the analysis of the results for flame E, the particle trajectories are divided into two groups: continuous burning and local extinction. For each group, the trajectories are further sub-divided based on the particles' origin: the fuel stream, the oxidizer stream, the pilot stream, and the intermediate region. The PDF calculations are performed using each of three commonly used models of molecular mixing, namely the EMST, IEM and modified Curl mixing models. The calculations with different mixing models reproduce the local extinction and re-ignition processes observed in flame E reasonably well. The particle behavior produced by the IEM and modified Curl models is different from that produced by the EMST model, i.e., the temperature drops prior to (and sometimes during) re-ignition. Two different re-ignition mechanisms are identified for flame E: auto-ignition and mixing-reaction. In the Cabra H₂/N₂ lifted flame, the particle trajectories are divided into different categories based on the particles' origin: the fuel stream, the oxidizer stream, and the intermediate region. The calculations reproduce the whole auto-ignition process reasonably well for the Cabra flame. Four stages of combustion in the Cabra flame are identified in the calculations by the different mixing models, i.e., pure mixing, auto-ignition, mixing-ignition, and fully burnt, although the individual particle behavior by the IEM and modified Curl models is different from that by the EMST model. The relative importance of mixing and reaction during re-ignition and auto-ignition are quantified for the IEM model.     

The differential diffusion effect of the intermediate species on the stability of premixed flames propagating in microchannels

The propagation of premixed flames in adiabatic and non-catalytic planar microchannels subject to an assisted or opposed Poiseuille flow is considered. The diffusive–thermal model and the well-known two-step chain-branching kinetics are used in order to investigate the role of the differential diffusion of the intermediate species on the spatial and temporal flame stability. This numerical study successfully compares steady-state and time-dependent computations to the linear stability analysis of the problem. Results show that for fuel Lewis numbers less than unity, Le_F < 1, and at sufficiently large values of the opposed Poiseuille flow rate, symmetry-breaking bifurcation arises. It is seen that small values of the radical Lewis number, Le_Z, stabilise the flame to symmetric shape solutions, but result in earlier flashback. For very lean flames, the effect of the radical on the flame stabilisation becomes less important due to the small radical concentration typically found in the reaction zone. Cellular flame structures were also identified in this regime. For Le_F > 1, flames propagating in adiabatic channels suffer from oscillatory instabilities. The Poiseuille flow stabilises the flame and the effect of Le_Z is opposite to that found for Le_F < 1. Small values of Le_Z further destabilise the flame to oscillating or pulsating instabilities.     

Stability characteristics and flowfields of turbulent non-premixed swirling flames

A simple, yet representative, burner geometry is used for the investigation of highly swirling turbulent unconfined, non-premixed, flames of natural gas. The burner configuration comprises a ceramic faced bluff-body with a central fuel jet. The bluff-body is surrounded by an annulus that delivers a swirling primary flow of air. The entire burner assembly is housed in a wind tunnel providing a secondary co-flowing stream of air. This hybrid bluff-body/swirl burner configuration stabilizes complex turbulent flames not unlike those found in practical combustors, yet is amenable to modelling because of its well-defined boundary conditions. Full stability characteristics including blow-off limits and comprehensive maps of flame shapes are presented for swirling flames of three different fuel mixtures: compressed natural gas (CNG), CNG–air (1:2 by volume) and CNG–H₂ (1:1 by volume).

It is found that with increased fuel flow, flame blow-off mode may change with swirl number, S_g. At low swirl, the flame remains stable at the base but blows off in the neck region further downstream. At higher swirl numbers, the flames peel off completely from the burner's base. Swirling CNG–air flames are distinct in that they only undergo base blow-off. In the low range of swirl number, increasing S_g causes limited improvement in the blow-off limits of the flames investigated and (for a few cases) can even lead to some deterioration over a small intermediate range of S_g. It is only above a certain threshold of swirl that significant improvements in blow-off limits appear. Six flames are selected for further detailed flowfield and composition measurements and these differ in the combination of swirl number, primary axial velocity through the annulus, U_s, and bulk fuel jet velocity, U_j. Only velocity field measurements are presented in this paper. A number of flow features are resolved in these flames, which resemble those already associated with non-reacting swirling flows of equivalent swirl obtained with the present burner configuration. Additionally, asymmetric flowfields inherent to some flames are revealed where the fluidic centreline of the flow (defined in the two-dimensional (U–W velocity pair) velocity field by the ?ω? = 0 tangential velocity contour), meanders strongly on either side of the geometric centreline downstream by about one bluff-body diameter. Flow structures revealed by the velocity data are correlated to flame shapes to yield a better understanding of how the velocity field influences the flames physical characteristics.     

LES of oxy-fuel jet flames using the Eulerian Stochastic Fields method with differential diffusion

The Eulerian Stochastic Fields (ESF) Monte Carlo method to solve the transported PDF (TPDF) equation is extended to account for differential diffusion effects by incorporating species individual molecular diffusivities. The method has been applied in Large Eddy simulation (LES) to non-piloted oxy-fuel jet flames at different Reynolds numbers experimentally investigated by Sevault et al. [1]. Due to the high H₂ content in the fuel stream and CO₂ in the oxidizer these flames pose new challenges to combustion modeling as the flame structures are different compared to CH₄/air flames. The simulations show very good agreement with the experiments in terms of mixture fraction conditional mean values for temperature and mayor species on the fuel lean side and the reaction zone, deviations on the fuel rich side are discussed. The trend and location of localized extinction is reproduced well in the simulations, as well as differential diffusion effects in the near field. Additionally, it is shown that a neglect of differential diffusion in the combustion model leads to a lifted flame.     

The diffusion flame structure of an ammonium perchlorate based composite propellant at elevated pressures

Examination of the surface behavior and flame structure of a bimodal ammonium perchlorate (AP) composite propellant at elevated pressure was performed using high speed (5 kHz) planar laser-induced fluorescence (PLIF) from 1 to 12 atm and visible surface imaging spanning 1–20 atm. The dynamics of the combustion of single, coarse AP crystals were resolved using these techniques. It was found that the ignition delay time for individual AP crystals contributed significant to the particle lifetime only at pressures below about 6 atm. In situ AP crystal burning rates were found to be higher than rates reported for pure AP deflagration studies. The flame structure was studied by exciting OH molecules in the gas phase. Two types of diffusion flames were observed above the composite propellant: jet-like flames and v-shaped, inverted, overventilated, flames (IOF) lifted off the surface. While jet-like diffusion flames have been imaged at low pressures and simulated by models, the lifted IOFs have not been previously reported or predicted. The causes for the observed flame structures are explained by drawing on an understanding of the surface topography and disparities in the burning rates of the fuel and oxidizer.     

Dynamic behavior of diffusion flame interacting with a large-scale vortex by laser imaging techniques

Instantaneous and simultaneous measurements of two-dimensional temperature and OH-LIF profiles by combining Rayleigh scattering with laser induced fluorescence (LIF) were demonstrated in a nitrogen-diluted hydrogen (H₂ 30% + N₂ 70%) laminar normal diffusion flame interacting with a large-scale vortex by oscillating central fuel flow or in an inverse diffusion flame by oscillating central airflow. The dynamic behavior of the diffusion flame extinction and reignition during the flame–vortex interaction processes was investigated. The results obtained are described as follows. (1) The width of the reaction zone decreases remarkably, and a decrease in flame temperature and OH-LIF is seen with increasing central airflow in an inverse diffusion flame. OH-LIF increases, and temperature does not change with increasing central fuel flow in a normal diffusion flame. The computations predict the experimental results well, and it is revealed that flame temperature characteristics result from the preferential diffusion of heat and species, which induces excess enthalpy or on enthalpy deficit, and an increase or decrease in H₂ mole fraction in the flame. (2) When a large velocity fluctuation is given to the central flow, the temperature and the OH-LIF at the reaction zone become thin at the convex and circumferential part of the vortex where a high temperature layer exists, and the temperature at the reaction zone is lowered in the inverse flame and the normal flame. (3) The width and temperature of the reaction zone interacting with the vortex recover quickly to that of the laminar steady flame after the vortex passing in the normal flame, but the recovery to that of the steady flame after the vortex passing is delayed in the inverse flame. (4) When a remarkably large velocity fluctuation is given to the central airflow in the inverse flame, thinning of temperature and reaction zone starts at the convex and circumferential part of the vortex, resulting in a and flame extinction completely occurs at the tail part of the vortex and makes the pair of edge flames. The outside edge flame reignites and connects with the upstream reaction zone. The inside edge flame finally extinguishes as the supply of fuel is interrupted by the outside edge flame.     

A review of chemical diffusion: Criticism and limits of simplified methods for diffusion coefficient calculation

A review of the physics and modelling of mass diffusion involving different gaseous chemical species is firstly proposed. Both accurate and simplified models for mass diffusion involve the calculation of individual species diffusion coefficients. Since these are computationally expensive, in CFD they are commonly estimated by assuming constant Lewis or Schmidt numbers for each chemical species. The constant Lewis number assumption is particularly used. As a matter of fact, these assumptions have never been theoretically justified nor verified in practical flames. The only published information are the first observations by Smooke and Giovangigli about the Lewis number against temperature distributions in methane–air premixed and counterflow diffusion one-dimensional flames. The aim of this work is to verify these assumptions. Functional dependences of molecular properties appearing in these numbers are made explicit to show that while Sc _i depends only on composition, Le _i depends also on temperature and therefore it certainly cannot be assumed constant in a flame. Then, accurately calculating molecular properties, distributions of these characteristic numbers against temperature are obtained a posteriori from numerical simulations of different flames, premixed and non-premixed, and burning different fuels. For non-premixed flames, individual species Lewis number distributions are broad for most of the species considered in this article, whilst they are tight for premixed flames. Some attention is focused on the particular shape of Lewis distributions in non-premixed flames: they are characterized by four or five (when extinction is experienced) branches associated to precise regions in the flame (basically, lean, rich and stoichiometric combusting zones). Instead, the Schmidt distributions are always tighter, also when extinctions take place: for many species they can be approximatively assumed constant. Finally, a simplified procedure to estimate individual species diffusion coefficients is suggested, assuming the median of non-premixed flame Schmidt distributions has a constant value for each chemical species.     

Cup-burner flame structure and extinguishment by CF3Br and C2HF5 in microgravity

The effects of fire-extinguishing agents CF₃Br and C₂HF₅ on the structure and extinguishing processes of microgravity cup-burner flames have been studied numerically. Propane and a propane–ethanol–water fuel mixture, prescribed for a Federal Aviation Administration (FAA) aerosol can explosion simulator test, were used as the fuel. The time-dependent, two-dimensional numerical code, which includes a detailed kinetic model (177 species and 2986 reactions), diffusive transport, and a gray-gas radiation model, revealed unique flame structure and predicted the minimum extinguishing concentration of agent when added to the air stream. The peak reactivity spot (i.e., reaction kernel) at the flame base stabilized a trailing flame. The calculated flame temperature along the trailing flame decreased downstream due to radiative cooling, causing local extinction at <1250 K and flame tip opening. As the mole fraction of agent in the coflow (X_a) was increased gradually: (1) the premixed-like reaction kernel weakened (i.e., lower heat release rate) (but nonetheless formed at higher temperature); (2) the flame base stabilized increasingly higher above the burner rim, parallel to the axis, until finally blowoff-type extinguishment occurred; (3) the calculated maximum flame temperature remained at nearly constant (≈1700 K) or mildly increased; and (4) the total heat release of the entire flame decreased (inhibited) for CF₃Br but increased (enhanced) for C₂HF₅. In the lifted flame base with added C₂HF₅, H₂O (formed from hydrocarbon-O₂ combustion) was converted further to HF and CF₂O through exothermic reactions, thus enhancing the heat-release rate peak. In the trailing flame, “two-zone” flame structure developed: CO₂ and CF₂O were formed primarily in the inner and outer zones, respectively, while HF was formed in both zones. As a result, the unusual (non-chain branching) reactions and the combustion enhancement (increased total heat release) due to the C₂HF₅ addition occurred primarily in the trailing diffusion flame.     

Lean or ultra-lean stretched planar methane/air flames

Lean premixed combustion has potential advantages of reducing pollutants and improving fuel economy. In some lean engine concepts, the fuel is directly injected into the combustion chamber resulting in a distribution of lean fuel/air mixtures. In this case, very lean mixtures can burn when supported by hot products from more strongly burning flames. This study examines the downstream interaction of opposed jets of a lean-limit CH₄/air mixture vs. a lean H₂/air flame. The CH₄ mixtures are near or below the lean flammability limit. The flame composition is measured by laser-induced Raman scattering and is compared to numerical simulations with detailed chemistry and molecular transport including the Soret effect. Several sub-limit lean CH₄/air flames supported by the products from the lean H₂/air flame are studied, and a small amount of CO₂ product (around 1% mole fraction) is formed in a “negative flame speed” flame where the weak CH₄/air mixture diffuses across the stagnation plane into the hot products from the H₂/air flame. Raman scattering measurements of temperature and species concentration are compared to detailed simulations using GRI-3.0, C₁, and C₂ chemical kinetic mechanisms, with good agreement obtained in the lean-limit or sub-limit flames. Stronger self-propagating CH₄/air mixtures result in a much higher concentration of product (around 6% CO₂ mole fraction), and the simulation results are sensitive to the specific chemical mechanism. These model-data comparisons for stronger CH₄/air flames improve when using either the C₂ or the Williams mechanisms.     

Fuel dependence of benzene pathways

The relative importance of formation pathways for benzene, an important precursor to soot formation, was determined from the simulation of 22 premixed flames for a wide range of equivalence ratios (1.0-3.06), fuels (C₁-C₁₂), and pressures (20-760 torr). The maximum benzene concentrations in 15 out of these flames were well reproduced within 30% of the experimental data. Fuel structural properties were found to be critical for benzene production. Cyclohexanes and C₃ and C₄ fuels were found to be among the most productive in benzene formation; and long-chain normal paraffins produce the least amount of benzene. Other properties, such as equivalence ratio and combustion temperatures, were also found to be important in determining the amount of benzene produced in flames. Reaction pathways for benzene formation were examined critically in four premixed flames of structurally different fuels of acetylene, n-decane, butadiene, and cyclohexane. Reactions involving precursors, such as C₃ and C₄ species, were examined. Combination reactions of C₃ species were identified to be the major benzene formation routes with the exception of the cyclohexane flame, in which benzene is formed exclusively from cascading fuel dehydrogenation via cyclohexene and cyclohexadiene intermediates. Acetylene addition makes a minor contribution to benzene formation, except in the butadiene flame where C₄H₅ radicals are produced directly from the fuel, and in the n-decane flame where C₄H₅ radicals are produced from large alkyl radical decomposition and H atom abstraction from the resulting large olefins.