Stationary combustion regimes and extinction limits of one-dimensional stretched premixed flames in a gap between two heat conducting plates

Stationary combustion regimes, their linear stability and extinction limits of stretched premixed flames in a narrow gap between two heat conducting plates are studied by means of numerical simulations in the framework of one-dimensional thermal-diffusion model with overall one-step reaction. Various stationary combustion modes including normal flame (NF), near-stagnation plane flame (NSF), weak flame (WF) and distant flame (DF) are detected and found to be analogous to the same-named regimes of conventional counterflow flames. For the flames stabilized in the vicinity of stagnation plane at moderate and large stretch rates (which are NF, NSF and WF) the effect of channel walls is basically reduced to additional heat loss. For distant flame characterized by large flame separation distance and small stretch rates intensive interphase heat transfer and heat recirculation are typical. It is shown that in mixture content / stretch rate plane the extinction limit curve has ε-shape, while for conventional counterflow flames it is known to be C-shaped. This result is quite in line with recent experimental findings and is explained by extension of extinction limits at small stretch rates at the expense of heat recirculation. Analysis of the numerical results makes possible to reveal prime mechanisms of flame quenching on different branches of ε-shaped extinction limit curve. Namely, two upper limits are caused by stretch and heat loss. These limits are direct analogs of the upper and lower limits on conventional C-shaped curve. Two other limits are related with weakening of heat recirculation and heat dissipation to the burner. Thus, the present study provides a satisfactory explanation for the recent experimental observations of stretched flames in narrow channel.     

Kinetics and extinction of non-premixed cool and warm flames of dimethyl ether at elevated pressure

The growing demand of clean and efficient propulsion and energy systems has sparked an interest in understanding low-temperature combustion at high pressure. Cool flame transition and extinction limits as well as oxygen concentration dependence at elevated pressures provide insights of the low-temperature and high-pressure fuel reactivity. A new experimental high-pressure counterflow burner platform was designed and developed to achieve the studies of high-pressure cool flames. Dimethyl ether (DME) was chosen to study its non-premixed cool flame in high-pressure counterflow burner at pressure up to 5 atm, perhaps for the first time. This paper investigates the effects of pressure on cool flame structure, extinction and transition limits, and oxygen concentration dependence as well as ozone assisted warm flames of DME in experiments and numerical simulations. The results show that the reignition transition from cool flame to hot flame occurs either with the decrease of the strain rate at a given fuel concentration and pressure or with the increase of fuel mole fraction or pressure at a given strain rate. Furthermore, it is shown that the higher pressure shifts the cool flame to higher strain rates and results in higher cool flame extinction strain rates. However, the existing kinetic model of DME fails in predicting the cool flame extinction limit at elevated pressures. Besides, the cool flame extinction limits are proportional to nth power of the oxygen concentration, [O₂]ⁿ, and the increase of pressure leads to stronger extinction limit dependence (larger n) on oxygen concentration. The present experiment and detailed kinetic analysis show clearly that increasing pressure promotes the low-temperature chemistry including the oxygen addition reactions. In addition, stable warm flame was first experimentally observed by using DME at elevated pressure with ozone sensitization.     

Structure and dynamic response of radiative diffusion flames

Nitrogen-diluted hydrogen burning in air is modeled numerically using a constant density and one-step reaction model in a plane two-dimensional counterflow configuration. An optically thin assumption is used to investigate the effects of radiation on the dynamics, structure, and extinction of diffusion flames. While there exist dual steady-state extinction limits for the 1D radiative flame response, it is found that as the 1D radiative extinction point is approached the 1D low-stretch diffusion flame exhibits oscillatory response, even with sub-unity Lewis number fuel. These radiation-induced limit cycle oscillations are found to have increasing amplitude and decreasing frequency as the stretch rate is reduced. Flame oscillation eventually leads to permanent extinction at the stretch rate which is larger than the steady-state radiative extinction value. Along the 1D radiative response curve, the transition from 1D flame to 2D structure and the differences in the resulting 2D flame patterns are also examined using a variety of initial profiles, with special emphasis on the comparison of using the initial profiles with and without a flame edge. Similar to the previous studies on the high-stretch adiabatic edge flames using the same configuration, the high-stretch radiative flames are found to resist 1D blow-off quenching through various 2D structures, including propagating front and steady cellular flames for initial profiles with and without flame edges. For all initial profiles studied, the low-stretch radiative flames are also found to exhibit different 2D flame phenomena near the 1D radiative extinction limit, such as transient cellular structures, steady cellular structures, and pulsating ignition fronts. Although the results demonstrate the presence of low-stretch and high-stretch 2D bifurcation branches close to the corresponding 1D extinction limits irrespective of the initial profile used, particular 2D flame structures in certain stretch rate range are initial profile dependent. The existence of two-dimensional flame structures beyond the 1D steady-state radiative extinction limit suggests that the flammable range is expanded as compared to that predicted by the 1D model. Hence, multi-dimensional flame patterns need to be accounted for when determining the flammability limits for a given system.     

Kinetic effects of NO addition on n-dodecane cool and warm diffusion flames

The kinetic effects of NO addition on the flame dynamics and burning limits of n-dodecane cool and warm diffusion flames are investigated experimentally and computationally using a counterflow system. The results show that NO plays different roles in cool and warm flames due to their different reaction pathway sensitivities to the flame temperature and interactions with NO. We observe that NO addition decreases the cool flame extinction limit, delays the extinction transition from warm flame to cool flame, and promotes the ignition transition from warm flame to hot flame. In addition, jet-stirred reactor (JSR) experiments of n-dodecane oxidation with and without NO addition are also performed to develop and validate a n-dodecane/NO_x kinetic model. Reaction pathway and sensitivity analyses reveal that, for cool flames, NO addition inhibits the low-temperature oxidation of n-dodecane and reduces the flame temperature due to the consumption of RO₂ via NO+RO₂?NO₂+RO, which competes with the isomerization reaction that continues the peroxy radical branching sequence. The model prediction captures well the experimental trend of the inhibiting effect of NO on the cool flame extinction limit. For warm flames, two different kinds of warm flame transitions, the warm flame extinction transition to cool flame and the warm flame reignition transition to hot flame, were observed. The results suggest that warm extinction transition to cool flame is suppressed by NO addition while the warm flame reignition transition to hot flame is promoted. The kinetic model developed captures well the experimentally observed warm flame transitions to cool flame but fails to predict the warm flame reignition to hot flame at similar experimental conditions.     

Transient interactions between a premixed double flame and a vortex

The interaction between a laminar flame and a vortex is an important study for understanding the fundamentals of turbulent combustion. In the past, however, flame-vortex interactions have been investigated only for high-temperature flames. In this study, the impact of a vortex on a premixed double flame, which consists of a coupled cool flame and a hot flame, is examined experimentally and computationally using dimethyl ether/oxygen/ozone mixtures. The double flame is first shown to occur near the extinction limit of the hot flame. The differences between steady-state cool flames, double flames, and hot flames are explored in a one-dimensional counterflow configuration. The transient interactions between double flames and impinging vortices are then investigated experimentally using a micro-jet and numerically in two-dimensional transient modeling. It is seen that the vortex can extinguish the near-limit hot flame locally, resulting in a lone cool flame. At higher vortex intensities, the cool flame may also be extinguished after the extinction of the hot flame. It is found that there can be three different transient flame structures coexisting at the same time: an extinguished flame hole, a cool flame, and a double flame. Moreover, flame curvature is shown to play an important role in determining whether the vortex weakens or strengthens the cool flame and double flame.     

The two-dimensional structure of low strain rate counterflow nonpremixed-methane flames in normal and microgravity

The structure and extinction of low strain rate nonpremixed methane–air flames was studied numerically and experimentally. A time-dependent axisymmetric two-dimensional (2D) model considering buoyancy effects and radiative heat transfer was developed to capture the structure and extinction limits of normal gravity (1-g) and zero gravity (0-g) flames. For comparison with the 2D modelling results, a one-dimensional (1D) flamelet computation using a previously developed numerical code was exercised to provide information on the 0-g flames. A 3-step global reaction mechanism was used in both the 1D and 2D computations to predict the measured extinction limit and flame temperature. Photographic images of flames undergoing the process of extinction were compared with model calculations. The axisymmetric numerical model was validated by comparing flame shapes, temperature profiles, and extinction limits with experiments and with the 1D computational results. The 2D computations yielded insight into the extinction mode and flame structure. A specific maximum heat release rate was introduced to quantify the local flame strength and to elucidate the extinction mechanism. The contribution by each term in the energy equation to the heat release rate was evaluated to investigate the multi-dimensional structure and radiative extinction of the 1-g flames. Two combustion regimes depending on the extinction mode were identified. Lateral heat loss effects and multi-dimensional flame and flow structure were also found. At low strain rates in 1-g flames (‘regime A’), the flame is extinguished from the weak outer edge of the flame, which is attributed to a multi-dimensional flame structure and flow field. At high strain rates, (‘regime B’), the flame extinction initiates near the flame centreline owing to an increased diluent concentration in the reaction zone, similar to the extinction mode of 1D flames. These two extinction modes can be clearly explained by consideration of the specific maximum heat release rate.     

Studies of the dynamics of autoignition assisted outwardly propagating spherical cool and double flames under shock-tube conditions

The initiation, propagation, and transition of the autoignition assisted spherical cool flame and double flame are studied numerically and experimentally using n-heptane/air/He mixtures under shock-tube experimental conditions over a wide range of temperatures. The primary goal of the current study is to understand the effects of the ignition Damkohler number, ignition energy, flame curvature, and autoignition-induced flow compression on the propagation of spherical flames to ensure the proper interpretation of shock-tube flame speed measurements at engine-relevant conditions. The results show that at high ignition Damkohler number, there are three different flame regimes, cool flame, double flame, and hot flame. The cool flame speed accelerates dramatically with the increase of ignition Damkohler number. In addition, it is found that the change of flame regime, low-temperature autoignition, flame stretch, and autoignition-induced flow compression result in a complicated non-linear dependence of flame speed on stretch. The results also reveal that the spherical cool flame has much lower Markstein length compared to the hot flame at T > 600 K. Moreover, it is found that both the autoignition assisted cool flame and the trailing hot flame front in the double flame can propagate much faster that the hot flame alone at the same mixture conditions, leading to a nonlinear dependence of flame speed on the mixture initial temperature. The simulated flame trajectories and the flame speed dependence on temperature agree qualitatively well with the shock-tube experiments. A quantitative criterion to ensure the accurate speed measurement of the cool and hot flame is proposed. The present study provides important physical insight and guidance for the flame speed measurement using a shock-tube at engine relevant conditions.     

Extinction and flame bifurcations of stretched dimethyl ether premixed flames

Extinction limits and flame bifurcation of lean premixed dimethyl ether–air flames are numerically investigated using the counterflow flame with a reduced chemistry. Emphasis is paid to the combined effect of radiation and flame stretch on the extinction and flammability limits. A method based on the reaction front is presented to predict the Markstein length. The predicted positive Markstein length agrees well with the experimental data. The results show that flow stretch significantly reduces the flame speed and narrows the flammability limit of the stretched dimethyl ether–air flame. It is found that the combined effect of radiation and flow stretch results in a new flame bifurcation and multiple flame regimes. At an equivalence ratio slightly higher than the flammability limit of the planar flame, the distant flame regime appears at low stretch rates. With an increase in the equivalence ratio, in addition to the distant flame, a weak flame isola emerges at moderate stretch rates. With a further increase in the equivalence ratio, the distant flame and the weak flame branches merge together, resulting in the splitting of the weak flame branch into two weak flame branches, one at low stretch and the other at high stretch. Flame stability analysis demonstrates that the high stretch weak flame is also stable. Furthermore, a K-shaped flammability limit diagram showing various flame regimes and their extinction limits is obtained.     

Extinction characteristics of isolated n-alkane fuel droplets during low temperature cool flame burning in air

Large carbon number n-alkanes are a notable component in all real transportation fuels, and their chemical structure fosters substantial low temperature kinetic reactivity. Normal alkanes have been studied in various canonical configurations but rarely in systems with strong coupling between low temperature chemistry and transport for pure as well as for multi-component n-alkane mixtures. The Flame Extinguishment (FLEX) experiments onboard the International Space Station provided a unique platform for investigating low temperature multi-phase n-alkane and iso-alkane combustion. Among the many interesting phenomena experimentally observed, cool flame extinction can occur, accompanied by the concurrent formation of a surrounding cloud of condensed vapor. In this work we conduct numerical simulations of high and low temperature combustion of large, initially single-component n-heptane, n-decane and n-dodecane droplets. The role of initial droplet diameter, operating pressure, and n-alkyl carbon number on the extinction of hot and low temperature flames is investigated and compared against the available experimental data. While all three fuels exhibit similar hot flame behavior, cool flame activity increases with the carbon number, resulting in an increased cool flame temperature and decreased extinction diameter. Multi-cyclic “hot/cool flame transitions” are found in air as pressure is slightly increased above one atmosphere. The cyclic behaviors correspond to continuously varying hot and cool flame transitions across the high, low, and negative temperature coefficient (NTC) kinetic regimes. Further increase in pressure results in a second stage steady “Warm flame” transition. The extinction of hot and cool flame has a strong non-linear dependence on ambient pressure but as the hot flame extinction diameter increases with pressure the extinction diameter of the cool flame decreases. The computational results are compared with a recent asymptotic analysis of FLEX n-alkane cool flames.     

Minimum ignition energy and propagation dynamics of laminar premixed cool flames

Laminar premixed cool flames, induced by the coupling of low-temperature chemistry and convective-diffusive transport process, have recently attracted extensive interest in combustion and engine research. In this work, numerical simulations have been conducted using a recently developed open-source reacting flow platform reactingFOAM-SCT, to investigate the minimum ignition energy (MIE) and propagation dynamics of premixed cool flames in a 1D spherical coordinate. Results have shown that when ignition energy is below the MIE of regular hot flames, a class of cool flames could be initiated, which allow much wider flammability limits, both lean and rich, compared to hot flames. Furthermore, the overall cool flame propagation dynamics exhibit intrinsic similarity to those of hot flames, in that, they begin with an ignition kernel propagation regime, followed by two transition regimes, and eventually reach a normal flame propagation regime. However, a spherical expanding cool flame responds completely differently to stretch. Specifically, a regular outwardly propagating hot spherical flame accelerates with increasing stretch rate when the mixture Le < 1 and decelerates when Le > 1. However, it is found that a cool flame always tends to decelerate with increasing stretch rate regardless of mixture composition, exhibiting unique flame aerodynamic characteristic. This research discovers novel features of premixed cool flame initiation and propagation dynamics and sheds light on flame transition, spark-ignition system design, and advanced engine combustion control.     

Effectiveness of flame suppressants on cool flames and hot flames

While the effectiveness of various flame suppressants such as bromotrifluoromethane and trimethylphosphate on hot flames has been relatively well studied over the years, such suppressants have not been examined in the context of low-temperature cool flames. This investigation solves this issue by exploring the extinction limits of six suppressants on both hot flames and cool flames in the counterflow geometry using n-dodecane as the fuel. In contrast to hot flames, it is found both experimentally and numerically that cool flames are relatively impervious to chemically based suppressants such as bromotrifluoromethane; these suppressants are essentially diluents at low temperatures. Detailed examination of the computed flame structure reveals that the reactions composing the catalytic cycles that interfere with hydrogen radical and hydroxyl radical production in hot flames are orders of magnitude lower in cool flames. Furthermore, mildly flammable suppressants such as trimethylphosphate and 2‑bromo-3,3,3-trifluoropropene are observed to ignite under the conditions necessary to initiate cool flames, which limits measurements of the cool flame extinction limits. This premature oxidation is not predicted by kinetic models describing the suppressant chemistry.     

An apparatus-independent extinction strain rate in counterflow flames

Resistance to extinction by stretch is a key property of any flame, and recent work has shown that this property controls the overall structure of several important types of turbulent flames. Multiple definitions of the critical strain rate at extinction (ESR) have been presented in the literature. However, even if the same definition is used, different experiments report different extinction strain rates for flames burning the same fuel-air mixture at very similar temperatures using similarly constructed opposed-flow instruments. Here we show that at extinction, all these flames are essentially identical, so one would expect that each would be assigned the same value of a parameter representing its intrinsic resistance-to-stretch-induced-extinction, regardless of the specifics of the experimental apparatus. A similar situation arises in laminar flame speed measurements since different apparatuses could result in different strain rate distributions. In that instance, the community has agreed to report the unstretched laminar flame speed, and methods have been developed to translate the experimental (stretched) flame speed into a universal unstretched laminar flame speed. We propose an analogous method for translating experimental measurements for stretch-induced extinction into an unambiguous and apparatus-independent quantity (ESR_∞) by extrapolating to infinite opposing burner separation distance. The uniqueness of the flame at extinction is shown numerically and supported experimentally for twin premixed, single premixed, and diffusion flames at Lewis numbers greater than and less than one. A method for deriving ESR_∞ from finite-boundary experimental studies is proposed and demonstrated for methane and propane experimental diffusion and premixed single flame data. The two values agree within the range of ESR differences typically observed between experimental measurements and simulation results for the traditional ESR definition.     

Three stage cool flame droplet burning behavior of n-alkane droplets at elevated pressure conditions: Hot,warm and cool flame

Transient, isolated n-alkane droplet combustion is simulated at elevated pressure for helium-diluent substituted-air mixtures. We report the presence of unique quasi-steady, three-stage burning behavior of large sphero-symmetric n-alkane droplets at these elevated pressures and helium substituted ambient fractions. Upon initiation of reaction, hot-flame diffusive burning of large droplets is initiated that radiatively extinguishes to establish cool flame burning conditions in nitrogen/oxygen “air” at atmospheric and elevated pressures. However, at elevated pressure and moderate helium substitution for nitrogen (X_He?>?20%), the initiated cool flame burning proceeds through two distinct, quasi-steady-state, cool flame burning conditions. The classical “Hot flame” (～1500?K) radiatively extinguishes into a “Warm flame” burning mode at a moderate maximum reaction zone temperature (～ 970?K), followed by a transition to a lower temperature (～765?K), quasi-steady “Cool flame” burning condition. The reaction zone (“flame”) temperatures are associated with distinctly different yields in intermediate reaction products within the reaction zones and surrounding near-field, and the flame-standoff ratios characterizing each burning mode progressively decrease. The presence of all three stages first appears with helium substitution near 20%, and the duration of each stage is observed to be strongly dependent on helium substitutions level between 20–60%. For helium substitution greater than 60%, the hot flame extinction is followed by only the lower temperature cool flame burning mode. In addition to the strong coupling between the diffusive loss of both energy and species and the slowly evolving degenerate branching in the low and negative temperature coefficient (NTC) kinetic regimes, the competition between the low-temperature chain branching and intermediate-temperature chain termination reactions control the “Warm” and “Cool” flame quasi-steady conditions and transitioning dynamics. Experiments onboard the International Space Station with n-dodecane droplets confirm the existence of these combustion characteristics and predictions agree favorably with these observations.     

Studies of autoignition-assisted nonpremixed cool flames

Autoignition-assisted nonpremixed cool flames of diethyl ether (DEE) are investigated in both laminar counterflow and turbulent jet flame configurations. First, the ignition and extinction limits of laminar nonpremixed cool flames of diluted DEE are measured and simulated using detailed kinetic models. The laminar flame measurements are used to validate the kinetic models and guide the turbulent flame measurements. The results show that, below a critical mixture condition, for elevated temperature and dilute mixtures, the cool flame extinction limit and the low-temperature ignition limit merge, leading to autoignition-assisted cool flame stabilization without hysteresis. Based on the findings from the laminar flame experiments, autoignition-assisted turbulent lifted cool flames are established using a Co-flow Axisymmetric Reactor-Assisted Turbulent (CARAT) burner. The lift-off heights of the turbulent cool flames are quantified using formaldehyde planar laser-induced fluorescence. Based on an analogy with autoignition-assisted lifted hot flames, a correlation is proposed such that the autoignition-assisted cool flame lift-off height scales with the product of the flow velocity and the square of the first-stage ignition delay time. Using this scaling, we demonstrate that the kinetic mechanism that most accurately predicts the laminar flame ignition and extinction limits also best predicts the turbulent cool flame lift-off height.     

Bifurcation and extinction limit of stretched premixed flames with chain-branching intermediate kinetics and radiative loss

Premixed counterflow flames with thermally sensitive intermediate kinetics and radiation heat loss are analysed within the framework of large activation energy. Unlike previous studies considering one-step global reaction, two-step chemistry consisting of a chain branching reaction and a recombination reaction is considered here. The correlation between the flame front location and stretch rate is derived. Based on this correlation, the extinction limit and bifurcation characteristics of the strained premixed flame are studied, and the effects of fuel and radical Lewis numbers as well as radiation heat loss are examined. Different flame regimes and their extinction characteristics can be predicted by the present theory. It is found that fuel Lewis number affects the flame bifurcation qualitatively and quantitatively, whereas radical Lewis number only has a quantitative influence. Stretch rates at the stretch and radiation extinction limits respectively decrease and increase with fuel Lewis number before the flammability limit is reached, while the radical Lewis number shows the opposite tendency. In addition, the relation between the standard flammability limit and the limit derived from the strained near stagnation flame is affected by the fuel Lewis number, but not by the radical Lewis number. Meanwhile, the flammability limit increases with decreased fuel Lewis number, but with increased radical Lewis number. Radical behaviours at flame front corresponding to flame bifurcation and extinction are also analysed in this work. It is shown that radical concentration at the flame front, under extinction stretch rate condition, increases with radical Lewis number but decreases with fuel Lewis number. It decreases with increased radiation loss.     

Influences of heat flux on extinction characteristics of steady/unsteady premixed stagnation flames

This paper investigates the extinction characteristics of premixed stagnation flames (PSFs) with controlled heat losses and flow disturbances. The low-frequency air flow pulsations that imitate the operational transients in practical combustors were specially introduced. The tunable diode-laser absorption spectroscopy (TDLAS) measurement was applied to obtain the temperature profile and wall heat flux. It is found that, for steady flame with a fixed equivalence ratio, the extinction stretch rate dramatically increases as the wall heat flux decreases. The extinction criterion is summarized as a global Karlovitz number of 0.57 by establishing a relationship between the global and local stretch rates. Numerical simulations reveal that the local extinction Karlovitz number of steady PSFs is approximately 1.0 regardless of the conditions such as heat flux and equivalence ratio. Further experiments present that the air pulsations with a repetition of ～5 Hz significantly deteriorate the flame stability. Particularly, for unsteady perturbed flames, the extinction stretch rate exhibits a nonlinear trend, yielding two regimes with discrepant sensitivities to wall heat flux. The unsteady simulation then highlights a local stretch rate overshoot in the presence of pulsation. It is caused by the time delay between the inlet velocity and flame front movement that eventually leads to poor flame stability. Moreover, in the high heat-flux regime, a smaller local stretch rate overshoot results in the weak dependence of extinction limits on heat fluxes.     

Critical temperature and reactant mass flux for radiative extinction of ethylene microgravity spherical diffusion flames at 1 bar

In microgravity combustion, where buoyancy is not present to accelerate the flow field and strain the flame, radiative extinction is of fundamental importance, and has implications for spacecraft fire safety. In this work, the critical point for radiative extinction is identified for normal and inverse ethylene spherical diffusion flames via atmospheric pressure experiments conducted aboard the International Space Station, as well as with a transient numerical model. The fuel is ethylene with nitrogen diluent, and the oxidizer is an oxygen/nitrogen mixture. The burner is a porous stainless-steel sphere. All experiments are conducted at constant reactant flow rate. For normal flames, the ambient oxygen mole fraction was varied from 0.2 to 0.38, burner supply fuel mole fraction from 0.13 to 1, total mass flow rate, ṁ_total, from 0.6 to 12.2 mg/s, and adiabatic flame temperature, T_ad, from 2000 to 2800 K. For inverse flames, the ambient fuel mole fraction was varied from 0.08 to 0.12, burner supply oxygen mole fraction from 0.4 to 0.85, ṁ_total from 2.3 to 11.3 mg/s, and T_ad from 2080 to 2590 K. Despite this broad range of conditions, all flames extinguish at a critical extinction temperature of 1130 K, and a fuel-based mass flux of 0.2 g/m²-s for normal flames, and an oxygen-based mass flux of 0.68 g/m²-s for inverse flames. With this information, a simple equation is developed to estimate the flame size (i.e., location of peak temperature) at extinction for any atmospheric-pressure ethylene spherical diffusion flame given only the reactant mass flow rate. Flame growth, which ultimately leads to radiative extinction if the critical extinction point is reached, is attributed to the natural development of the diffusion-limited system as it approaches steady state and the reduction in the transport properties as the flame temperature drops due to increasing flame radiation with time (radiation-induced growth.)     

Effect of addition of a non-equidiffusional reactant to an equidiffusional diffusion flame

A Burke–Schumann (flame-sheet) formulation is developed for diffusion flames between a fuel and oxidiser with Lewis numbers of unity, subject to addition to the fuel and/or oxidiser stream of a different reactant for which the Lewis number differs from unity. This formulation is applied to laminar counterflow diffusion-flame experiments, reported here, in which hydrogen was added to either methane–nitrogen mixtures or oxygen–nitrogen mixtures at normal atmospheric pressure, with both feed streams at normal room temperature. Experimental conditions were adjusted to fix selected values of the stoichiometric mixture fraction and the adiabatic flame temperature, and the strain rate was increased gradually, maintaining the momentum balance of the two streams, until extinction occurred. At the selected sets of values, the strain rate at extinction was measured as a function of the hydrogen concentration in the fuel or oxidiser stream. The ratio of the fraction of the oxidiser flux that consumes hydrogen to the fraction that consumes fuel was calculated from the new Burke–Schumann formulation, and it was found that, within experimental uncertainty, the ratio of the extinction strain rate with hydrogen addition to that without was the same at any given value of this oxygen flux ratio, irrespective of whether the hydrogen was added on the fuel or oxidiser side. This experimental result was also in close agreement with computational predictions employing detailed chemistry. These results imply that differences in detailed hydrogen concentration profiles within the reaction zone have little or no influence on the chemical kinetics of extinction when the stoichiometric mixture fraction, the adiabatic flame temperature, and the proportion of oxygen that consumes the added fuel are fixed. This same correspondence may be expected to apply for other fuels and additives.     

Oscillatory cool flame combustion behavior of submillimeter sized n-alkane droplet under near limit conditions

This paper reports simulation results of oscillatory cool flame burning of an isolated, submillimeter sized n-heptane (n-C₇H₁₆) droplet in a selectively ozone (O₃) seeded nitrogen-oxygen (N₂-O₂) environments at atmospheric pressure. An evolutionary one-dimensional droplet combustion code encompassing relevant physics and detailed chemistry was employed to explore the roles of low-temperature chemistry, O₃ seeding, and dynamic flame structure on burning behaviors. For X_O2= 21% and a range of selective ozone seeding, near-quasi-steady cool flame burning is achieved directly (without requiring hot flame initiation and radiative extinction). Under low oxygen index conditions, but with significant O₃ seeding (X_O3 = 5%), a nearly quasi-steady cool flame is initially established that then transitions to a dynamically oscillating cool flame burning mode which continues until the droplet is completely consumed. It is found that the oscillation occurs as result of a initial depletion of fuel vapor-oxidizer layer evolving near the droplet surface and its dynamic re-establishment through liquid vaporization and vapor/oxidizer transport. A kinetic analysis indicates that the dynamic competition between the reaction classes- (a) degenerate chain branching and (b) chain termination/propagation - along with continuous fuel and oxygen leakage through the flame location contributes to an oscillatory burning phenomena of ever-increasing amplitude. Analysis based on single full-cycle of oscillatory burning shows that the reaction progression matrices (evolution of heat and species) for QOOH➔chain propagation/termination reactions (here, Q = C₇H₁₄-) directly scales with the gas phase temperature field. On the contrary, the QOOH➔degenerate branching reactions undergoes three distinct stages within the same oscillatory cycle. The coupled flame dynamics and kinetics suggest that in the oscillatory burning mode, kinetic processes dynamically cross through conditions characterizing the negative temperature coefficient (NTC) turnover temperature, separating low temperature and NTC kinetic regimes. In addition, a parametric study is conducted to determine the role of O₃ seeding level on the observed oscillation phenomena.