Kinetic and radiative extinctions of spherical burner-stabilized diffusion flames

Extinction of steady, spherical diffusion flames stabilized by a spherical porous burner was investigated by activation energy asymptotics. An optically-thin radiation model was employed to study the effect of radiation on flame extinction. Four model flames with the same adiabatic flame temperature and fuel consumption rate but different stoichiometric mixture fraction and flow direction, namely the flames with fuel issuing into air, diluted fuel issuing into oxygen, air issuing into fuel, and oxygen issuing into diluted fuel, were adopted to understand the relative importance of residence time and radiation intensity. Results show that for a specified flow rate emerging from the burner, only the kinetic extinction limit at low Damköhler numbers (low residence times) exists. In the presence of radiative heat loss, extinction is promoted so that it occurs at a larger Damköhler number. By keeping the radiation intensity constant while varying the flow rate, both the kinetic and radiative extinction limits, representing the smallest and largest flow rates, between which steady burning is possible, are exhibited. For flames with low radiation intensity, extinction is primarily dominated by residence time such that the high-flow rate flames are easier to be extinguished. The opposite is found for flames suffering strong radiative heat loss. The kinetic extinction limit might occur at mass flow rates lower than what is needed to keep the flame outside of the burner and not observable. An extinction state on the radiative extinction branch can be either kinetic or radiative depending on the process.     

Study on stretch extinction limits of CH4/CO2 versus high temperature O2/CO2 counterflow non-premixed flames

Extinction limits of counterflow non-premixed flames with normal and high temperature oxidizers were studied experimentally and numerically for development of new-type oxygen-enriched mild combustion furnace. Extinction stretch rates of CH₄/CO₂ (at 300 K) versus O₂/CO₂ flames at oxygen mole fractions of 0.35 and 0.40 and oxidizer temperatures of 300 K, 500 K, 700 K and 1000 K were obtained. Investigation was also conducted for CH₄/N₂ (at 300 K) versus air (O₂/N₂) flames at the same oxidizer temperatures. An effect of radiative heat loss on stretch extinction limits of oxygen-enriched flames and air flames was investigated by computations with optical thin model (OTM) and adiabatic flame model (ADI). The results show influence of radiative heat loss on stretch extinction limits was not significant in relative high fuel mole fraction regions. The extinction curve of the oxygen-enriched flames with oxygen mole fraction of 0.35 was close to that of the air flames at the oxidizer temperature of 300 K. However, the extinction curve of air flames with high temperature oxidizer was comparable with that of oxygen-enriched flames with oxygen mole fraction of 0.40. Scaling analysis based on asymptotic solution of stretch extinction was applied and it was found that stretch extinction limits can be expressed by two terms. The first term is total enthalpy flux of fuel stream based on thermo-physical parameters. The second term is a kinetic term which reflects an effect of the chemical reaction rate on stretch extinction limits. OH radicals which play important roles in chain propagating and main endothermic reactions were used to represent the kinetic term of both oxygen-enriched and air flames. The global rates of OH formation in these two cases were compared to understand the contribution of kinetic term to stretch extinction limits. Variation of extinction curves of oxygen-enriched flames and air flames was well explained by the present scaling analysis. This offers an effective approach to estimate stretch extinction limits of oxygen-enriched flames based on those of air flames at the same oxidizer temperature.     

Effects of structure and hydrodynamics on the sooting behavior of spherical microgravity diffusion flames

This study is an examination of the sooting behavior of spherical microgravity diffusion flames burning ethylene at atmospheric pressure in a 2.2-s drop tower. In a novel application of microgravity, spherical flames were employed to allow convection across the flame to be either from fuel to oxidizer or from oxidizer to fuel. Thus, spherical microgravity flames are capable of allowing stoichiometric mixture fraction, Z_st, and direction of convection across the flame to be controlled independently. This allowed for a study of the phenomenon of permanently blue diffusion flames—flames that remain blue as strain rate approaches zero. Z_st was varied by changing inert concentrations such that adiabatic flame temperature did not change. At low Z_st, nitrogen was supplied with the oxidizer, and at high Z_st, it was provided with the fuel. Flame structure, quantified by Z_st, was found to have a profound effect on soot production. Soot-free conditions were observed at high Z_st and sooting conditions were observed at low Z_st regardless of convection direction. Convection direction was found to have a smaller impact on soot inception, suppressing formation when convection at the flame sheet was directed towards the oxidizer. A numerical analysis was developed to simulate steady state conditions and aided the interpretation of the results. The analysis revealed that steady state was not achieved for any of the flames, but particularly for those with pure ethylene or oxygen flowing from the porous burner. Furthermore, despite the fact that all flames had the same adiabatic flame temperature, the actual peak temperatures differed considerably. While transient burner heating and burner radiation reduced flame temperature, gas-phase radiative heat loss was the dominant mechanism accounting for these differences.     

Edge flame instability in low-strain-rate counterflow diffusion flames

Experiments in low-strain-rate methane-air counterflow diffusion flames diluted with nitrogen have been conducted to study flame extinction behavior and edge flame oscillation in which flame length is less than the burner diameter and thus lateral conductive heat loss, in addition to radiative loss, could be high at low global strain rates. The critical mole fraction at flame extinction is examined in terms of velocity ratio and global strain rate. Onset conditions of the edge flame oscillation and the relevant modes are also provided with global strain rate and nitrogen mole fraction in the fuel stream or in terms of fuel Lewis number. It is observed that flame length is intimately relevant to lateral heat loss, and this affects flame extinction and edge flame oscillation considerably. Lateral heat loss causes flame oscillation even at fuel Lewis number less than unity. Edge flame oscillations, which result from the advancing and retreating edge flame motion of the outer flame edge of low-strain-rate flames, are categorized into three modes: a growing, a decaying, and a harmonic-oscillation mode. A flame stability map based on the flame oscillation modes is also provided for low-strain-rate flames. The important contribution of lateral heat loss even to edge flame oscillation is clarified finally.     

A theoretical analysis of the extinction limits of a methane-air opposed-jet diffusion flame

A theoretical analysis is described for a methane-air diffusion flame stabilized in the forward stagnation region of a porous metal cylinder in a forced convective flow. The analysis includes effects of radiative heat loss from the porous metal surface and finite rate kinetics but neglects the effects of gravity. The theoretically predicted extinction limits compare well with experimentally observed extinction limits from the literature.After the predicted limits compared well with the experimental limits, a parametric study of the effect of fuel surface emissivity and Lewis number was conducted with the numerical model. It was found that the computed blowoff limit is independent of radiative heat loss for high fuel blowing velocities but is a strong function of Lewis number. At low fuel blowing velocities, the extinction limit varies with both radiative heat loss and Lewis number. It is discovered, however, that even if thermal losses from the fuel surface are absent, the flame can extinguish at the fuel surface independently of Lewis number due to excessive reaction zone thinning.     

Dynamics of an edge flame in a mixing layer

The dynamics of an edge flame in a mixing layer is considered. The flame, which stands at a well-defined distance from the plate separating the fuel and oxidizer streams, is stabilized by heat conducted back to the relatively cold plate. It has a tribrachial structure that consists of a lean premixed segment leaning toward the oxidizer side, a rich premixed segment leaning toward the fuel side, and a diffusion flame trailing behind. Within the context of a diffusive-thermal model, we describe numerically the steady and unsteady behavior of the flame. The primary objective is to systematically identify the conditions for the onset of oscillations examining, in particular, the influence of differential diffusion, mixture strength, flow rate, and radiative heat losses.     

A computational study of oscillatory extinction of spherical diffusion flames

The transient behavior of burner-supported spherical diffusion flames was studied in the transport-induced limit of low mass flow rate and the radiation-induced limit of high mass flow rate which characterize the isola response of flame extinction. Oscillatory instability was observed near both steady-state extinction limits. The oscillation typically grows in amplitude until it becomes large enough to extinguish the flame. The oscillatory behavior was numerically observed using detailed chemistry and transport for methane (50%CH₄/50%He into 21%O₂/79%He) and hydrogen (100% H₂ into 21%O₂/79%He) diffusion flames where the fuel was issued from a point source, and helium was selected as an inert to increase the Lewis number, facilitating the onset of oscillation. In both methane and hydrogen flames, the oscillation always leads to extinction, and no limit cycle behavior was found. The growth rate of the oscillation was found to be slow enough under certain conditions to allow the flame to oscillate for over 450 s, suggesting that such oscillations can possibly be observed experimentally. For the hydrogen flames, however, the frequency of oscillation near the transport-induced limit is much larger, approximately 60 Hz as compared to 0.35 Hz for the methane flame, and the maximum amplitude of temperature oscillations was about 5 K. The distinctively different structures of the hydrogen and methane flames suggest that while both instabilities are thermal-diffusive in origin, oscillations in the hydrogen flames resemble those of premixed flames, while oscillations in the methane flames are non-premixed in character.     

Suppression limits of low strain rate non-premixed methane flames

The suppression of low strain rate non-premixed flames was investigated experimentally in a counterflow configuration for laminar flames with minimal conductive heat losses. This was accomplished by varying the velocity ratio of fuel to oxidizer to adjust the flame position such that conductive losses to the burner were reduced and was confirmed by temperature measurements using thermocouples near the reactant ducts. Thin filament pyrometry was used to measure the flame temperature field for a curved diluted methane-air flame near extinction at a global strain rate of 20 s⁻¹. The maximum flame temperature did not change as a function of position along the curved flame surface, suggesting that the local agent concentration required for suppression will not differ significantly along the flame sheet. The concentration of N₂, CO_2, and CF₃Br added to the fuel and the oxidizer streams required to obtain extinction was measured as a function of the global strain rate. In agreement with previous measurements performed under microgravity conditions, limiting non-premixed flame extinction behavior in which the agent concentration obtained a value that insures suppression for all global strain rates was observed. A series of extinction measurements varying the air:fuel velocity ratio showed that the critical N₂ concentration was invariant with this ratio, unless conductive losses were present. In terms of fire safety, the measurements demonstrate the existence of a fundamental limit for suppressant requirements in normal gravity flames, analogous to agent flammability limits in premixed flames. The critical agent volume fraction in the methane fuel stream assuring suppression for all global strain rates was measured to be 0.841 ± 0.01 for N₂, 0.773 ± 0.009 for CO₂, and 0.437 ± 0.005 for CF₃Br. The critical agent volume fraction in the oxidizer stream assuring suppression for all global strain rates was measured as 0.299 ± 0.004 for N₂, 0.187 ± 0.002 for CO₂, and 0.043 ± 0.001 for CF₃Br.     

Effects of Damköhler number on flame extinction and reignition in turbulent non-premixed flames using DNS

Results from a parametric study of flame extinction and reignition with varying Damköhler number using direct numerical simulation are presented. Three planar, non-premixed ethylene jet flames were simulated at a constant Reynolds number of 5120. The fuel and oxidizer stream compositions were varied to adjust the steady laminar extinction scalar dissipation rate, while maintaining constant flow and geometric conditions. Peak flame extinction varies from approximately 40% to nearly global blowout as the Damköhler number decreases. The degree of extinction significantly affects the development of the jets and the degree of mixing of fuel, oxidizer, and combustion products prior to reignition. The global characteristics of the flames are presented along with an analysis of the modes of reignition. It is found that the initially non-premixed flame undergoing nearly global extinction reignites through premixed flame propagation in a highly stratified mixture. A progress variable is defined and a budget of key terms in its transport equation is presented.     

Numerical investigation of spherical diffusion flames at their sooting limits

Detailed numerical simulations are presented of laminar microgravity spherical diffusion flames at their experimentally observed sooting limits. Ten normal and inverse flames fueled by ethylene are considered. Observed in a drop tower, these flames were initially sooty but reached their sooting limits 2 s after ignition (or slightly before). The flames span broad ranges of stoichiometric mixture fraction (0.065–0.692), adiabatic flame temperature (2226–2670 K), and stoichiometric scalar dissipation rate (0.013–0.384 s^?1). They were modeled using a one-dimensional, transient diffusion flame code with detailed chemistry (up to toluene) and transport. Radiative losses from products were modeled using a detailed absorption/emission statistical narrow-band model coupled with a discrete-ordinates method. Flame structure at the sooting limits was examined, emphasizing the behavior of carbon to oxygen atom ratio, temperature, and scalar dissipation rate. For ethylene flames with sufficiently long flow times it was found that soot formation coincides with regions where the C/O atom ratio and temperature exceed critical values, specifically 0.53 and 1305 K, respectively. The scatter about these critical values is small, which is noteworthy considering the wide range of flame conditions. These observations are consistent with the expected effects of H radicals on the propargyl soot pathway.     

Investigation of soot transport and radiative heat transfer in an ethylene jet diffusion flame

Numerical and experimental investigations highlighting the heat and mass transfer phenomena in a laminar co-flowing jet diffusion flame have been carried out. The fuel under consideration is ethylene, with ambient air as the co-flowing oxidizer. The diffusion flame is modeled using a 17-step reduced reaction mechanism with finite rate chemistry and the effects of soot on the radiative heat transfer of the flame have been demonstrated. Soot growth and oxidation processes are studied using a two-equation transport model, while the radiative heat transfer is modeled using the P1 approximation. An in-house finite volume code has been developed to solve the axi-symmetric Navier–Stokes equations in cylindrical coordinates, along with the soot mass fraction, soot number density, energy and species conservation equations. Comparison of predictions with experimental results shows reasonable agreement with regard to the flame height and temperature distribution. A parametric study is also presented, which illustrates the effects of the fuel jet Reynolds number and the flow rate of co-flow air.     

A radical index for the determination of the chemical kinetic contribution to diffusion flame extinction of large hydrocarbon fuels

The extinction limits of diffusion flames have been measured experimentally and computed numerically for fuels of three different molecular structures pertinent to surrogate fuel formulation: n-alkanes, alkyl benzenes, and iso-octane. The focus of this study is to isolate the thermal and mass transport effects from chemical kinetic contributions to diffusion flame extinction, allowing for a universal correlation of extinction limit to molecular structure. A scaling analysis has been performed and reveals that the thermal and mass transport effects on the extinction limit can be normalized by consideration of the enthalpy flux to the flame via the diffusion process. The transport-weighted enthalpy is defined as the product of the enthalpy of combustion per unit mole of fuel and the inverse of the square root of fuel molecular weight. The chemical kinetic contribution provided by the specific fuel chemistry has thus been elucidated for tested individual component and multi-component surrogate fuels. A chemical kinetic flux analysis for n-decane flames shows that the production/consumption rates of the hydroxyl (OH) radical govern the heat release rate in these flames and therefore play significant roles in defining the extinction limit. The rate of OH formation has been defined by considering the OH concentration, flame thickness, and flow strain rate. A fuel-specific radical index has been introduced as a concept to represent and quantify the kinetic contribution to the extinction limit owing to the fuel-specific chemistry. A relative radical index scale, centered on the radical index of a series of n-alkanes which are observed and fundamentally explained to be common, is established. A universal correlation of the observed extinction limits of all tested fuels has been obtained through a combined metric of radical index and transport-weighted enthalpy. Finally, evidence as to the validity of the fundamental arguments presented is provided by the success of the universal correlation in predicting the extinction limits of the multi-component mixtures typical of surrogate fuels.     

Extinction of counterflow diffusion flames with radiative heat loss and nonunity Lewis numbers

The structure and extinction characteristics of counterflow diffusion flames with flame radiation and nonunity Lewis numbers of the fuel and oxidant are examined using multiscale asymptotic theory, and a model expressed in terms of the jump relations and reactant leakages with the proper consideration of the excess enthalpy overlooked in previous analyses is developed. The existence of the dual extinction limits in the presence of radiative heat loss, namely the kinetic limit at small Damköhler number (high stretch rate) and the radiative limit at large Damköhler number (low stretch rate), are identified. It is found that the former is minimally affected by radiative loss, while a substantial amount of heat loss is associated with the radiative limit. Reactant leakage, however, is the root cause for both limits. The influence of radiative loss on the extinction Damköhler numbers is found to be through its effects on the flame temperature, the excess enthalpy, and the reduced extinction Damköhler number. At both extinction limits, the contribution from the flame temperature is always important and dominant. The contributions from the other two, however, could be important in some special cases. At small Le_F, the contribution from the reduced extinction Damköhler number is large and even dominant under small radiative loss. The contribution from the excess enthalpy is important for small Le_O and it may be comparable to the contribution from the flame temperature when radiative loss is small. Thus, overlooking the excess enthalpy in previous analyses may have resulted in rather large error in the predicted extinction Damköhler numbers, especially the kinetic one.     

Flame structures and behaviors of opposed flow non-premixed flames in mesoscale channels

An opposed flow non-premixed flame (OFNPF) in a narrow channel was chosen as a model of a non-premixed flame in a mesoscale combustion space or micro-combustor. The stabilization limits and behaviors of methane-air flames and propane-air flames were compared for various experimental parameters such as flow velocity, nozzle distance, nozzle width, channel gap, and fuel dilution. Flames could be stabilized in a wide range of strain rates (0.9–150 s⁻¹) and dilution ratios (∼80% nitrogen at the fuel side). The flame extinction limits were classified into three types and their mechanisms were investigated: higher-strain-rate (HSR) extinction limit determined by the flame stretch, lower-strain-rate (LSR) extinction limit determined by the conductive or convective heat loss from the flame, and fuel-dilution-ratio (FDR) extinction limit determined by the decrease in the heat release rate from the flames. The HSR extinction limits in mesoscale channels could be explained with a modified strain rate, and the LSR extinction limits could be explained by employing a premixed quenching theory in which the heat loss through the dead space near the wall was considered as a major extinction mechanism. Finally, the variation of the extinction limits with the FDR in both the HSR and the LSR conditions could be explained with a modified global reaction rate in which the variations in flame temperature and species concentrations were reflected. This study provides an essential model for the stabilization and extinction of non-premixed flames in mesoscale combustion spaces.     

Experimental investigation of the morphology and stability of diffusion and edge flames in an opposed jet burner

The stability of methane/air and hydrogen/air flames in an axisymmetric counterflow burner was investigated experimentally for different burner geometries, degrees of fuel dilution, and combinations of flow velocities. Both planar diffusion flames and edge flames were observed, and the transitions between these flame types were studied. The experimental results confirmed previously published numerical predictions on diluted hydrogen/air flames: the existence of two distinct stable flame types; the possibility of switching between the two flame types by perturbing the flames, e.g., by suitably changing a flow velocity; and the strong hysteresis for the transition from one flame type to the other. Flame stability diagrams were compiled which delineate the range of fuel and air flow velocities for which the planar diffusion flame and the toroidal edge flame are stable. The lower boundary curve for the edge flame stability exhibits a characteristic minimum at a well-defined value of the fuel velocity. For fuel velocities lower than this value, the transition between the edge and the diffusion structure is reversible, and the flames exhibit bistable behavior. For higher fuel velocities, the decrease of air velocity leads to the extinction of the edge flame. An investigation of both the cold and the reactive flow field identified bistable behavior for the flow field as well. Except for very low flow rates, the stagnation plane stabilizes in two positions, close to either of the two nozzles. Detailed numerical simulations of hydrogen flames capture the essentials of this behavior. The observed flame extinction results from the interaction of the flame dynamics with the dynamics of the flow field.     

An analysis of sub-limit flame dynamics using opposite propagating flames in mesoscale channels

The excess enthalpy flames and their dynamics below the flammability limit are studied by considering two flames that propagate in opposite directions in parallel channels. The model enables the coupling between the external heat loss, convection preheating, diffusion transport and finite rate chemistry. Analytical expressions for the flame temperature, separation distance, and extinction limit are obtained. The results show that flame extinction can be caused by the external heat loss without heat conduction of inner wall in the streamwise direction. The heat recirculation across the separating wall dramatically increases the flame speed and extends the flammability limit. It is shown that the maximum and minimum flame speeds corresponding respectively to the fast and slow flame modes exist at all separation distances between the two flames. It is found that the flame can adjust its separation distance to adapt to the variation of heat loss, heat recirculation and fuel concentration. There exists a maximum flame separation distance beyond which sub-limit flame does not exist. The results also showed that heat recirculation significantly extends the flammability limit. Furthermore, at low fuel concentrations, the flame can be stabilized in a narrow range of separation distance. The present study not only generalized the previous analyses of the heat recirculation flames but also provided a model for the study and control of sub-limit flames in micro power devices and reactors.     

Soot and NO formation in methane–oxygen enriched diffusion flames

NO and soot formation were investigated both numerically and experimentally in oxygen-enriched counterflow diffusion flames. Two sets of experiments were conducted. In the first set, the soot volume fraction was measured as a function of oxygen content in the oxidizer jet at constant strain rate (20 s⁻¹). In the second set of experiments, the soot volume fraction was measured as a function of strain rate variation from 10 to 60 s⁻¹ and at constant oxygen content on the oxidizer side. A soot model was developed based on a detailed C₆ gas phase chemistry. The soot and molecular radiation were taken into account. Numerical results were verified against experimental data. The soot volume fraction was predicted with the maximum discrepancy less than 30% for all cases considered. It was found that oxygen variation significantly modified the diffusion flame structure and the flame temperature, resulting in a substantial increase of soot. The temperature increase promotes aromatics production in the fuel pyrolysis zone and changes the relative contributions of the thermal and Fenimore mechanisms into NO formation. As the strain rate increases, the residence time of incipient soot particles in the high temperature zone is reduced and the total amount of soot decreases. High concentration of soot in the flame leads to enhancement of radiant heat exchange: the reduction of temperature due to radiation was found to be between 10 and 50 K. This caused a reduction of peak NO concentrations by 20%–25%. The increase of oxygen content in the oxidizer stream resulted in a reduction of the distance between the plane of the maximum temperature and the stagnation plane.     

Effect of unequal fuel and oxidizer Lewis numbers on flame dynamics

The interaction of non-unity Lewis number (due to preferential diffusion and/or unequal rates of heat and mass transfer) with the coupled effect of radiation, chemistry and unsteadiness alters several characteristics of a flame. The present study numerically investigates this interaction with a particular emphasis on the effect of unequal and non-unity fuel and oxidizer Lewis numbers in a transient diffusion flame. The unsteadiness is simulated by considering the flame subjected to modulations in reactant concentration. Flames with different Lewis numbers (ranging from 0.5 to 2) and subjected to different modulating frequencies are considered. The results show that the coupled effect of Lewis number and unsteadiness strongly influences the flame dynamics. The impact is stronger at high modulating frequencies and strain rates, particularly for large values of Lewis numbers. Compared to the oxidizer side Lewis number, the fuel side Lewis number has greater influence on flame dynamics.     

喷管内传递特性对微尺度扩散火焰的影响

对均匀空气流中微尺度甲烷扩散燃烧进行了数值模拟,重点考察微喷管内的流动和传热传质对微尺度燃烧特性的影响.结果表明,在低流速下,内径为0.3 mm的微喷管内进气速度为1.0 m/s时燃料与空气的混合已经发生,混合气被管外的热量预热,同时火焰的热损失增加.在喷管直径一定时,减小燃料喷出速度,传热传质现象对微尺度甲烷扩散火焰特性的影响增强;当进气速度为0.5 m/s时,甲烷在微喷管内开始燃烧,放出热量.在进行微尺度解析计算时,必须包含一定的喷管区域.     

From diffusion to premixed flames in an H2/air opposed-jet burner: the role of edge flames

Edge flames obtained on a hydrogen/air non-premixed opposed-jet burner after the local extinction of the disk-shaped diffusion flame are investigated with 2-D direct numerical simulations using detailed chemical kinetics and transport. Over a large range of flowrates, edge flames were found to coexist with the well-known strongly burning diffusion flames corresponding to the upper branch of the S-shaped curve. The critical flowrates of the strong hysteresis associated with the transitions between the two solution branches were identified: re-establishment of the diffusion flame is controlled by the propagation of the edge flame and cannot be represented simply by the extinction scalar dissipation rate. It was also observed that in all the flow conditions simulated, the edge flame was able to consume all the supplied fuel by re-orienting itself, varying its flame surface area, or changing its structure. The latter was found to depend on the flow conditions (which strongly affects the degree of mixing ahead of the edge flame) and can take on different configurations ranging from a triple flame to an essentially premixed flame. Because of flame curvature and the preferential diffusion of hydrogen, the propagation speed of the edge flames was found to be higher than that of the corresponding planar premixed flames.