Experimental low-stretch gaseous diffusion flames in buoyancy-induced flowfields

The present study experimentally investigates the structure and instabilities associated with extremely low-stretch (1 s⁻¹) gaseous diffusion flames. Ultra-low-stretch flames are established in normal gravity by bottom burning of a methane/nitrogen mixture discharged from a porous spherically symmetric burner of large radius of curvature. OH-PLIF and IR imaging techniques are used to characterize the reaction zone and the burner surface temperature, respectively. A flame stability diagram mapping the response of the ultra-low-stretch diffusion flame to varying fuel injection rate and nitrogen dilution is explored. In this diagram, two main boundaries are identified. These boundaries separate the stability diagram into three regions: sooting flame, non-sooting flame, and extinction. Two distinct extinction mechanisms are noted. For low fuel injection rates, flame extinction is caused by heat loss to the burner surface. For relatively high injection rates, at which the heat loss to burner surface is negligible, flame radiative heat loss is the dominant extinction mechanism. There also exists a critical inert dilution level beyond which the flame cannot be sustained. The existence of multi-dimensional flame phenomena near the extinction limits is also identified. Various multi-dimensional flame patterns are observed, and their evolutions are studied using direct chemiluminescence and OH-PLIF imaging. The results demonstrate the usefulness of the present burner configuration for the study of low-stretch gaseous diffusion flames.     

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.     

Dynamics and burning limits of near-limit hot,warm, and cool diffusion flames of dimethyl ether at elevated pressures

The near-limit diffusion flame regimes and extinction limits of dimethyl ether at elevated pressures and temperatures are examined numerically in the counterflow geometry with and without radiation at different oxygen concentrations. It is found that there are three different flame regimes—hot flame, warm flame, and cool flame—which exist, respectively, at high, intermediate, and low temperatures. Furthermore, they are governed by three distinct chain-branching reaction pathways. The results demonstrate that the warm flame has a double reaction zone structure and plays a critical role in the transition between cool and hot flames. It is also shown that the cool flame can be formed in several different ways: by either radiative extinction or stretch extinction of a hot flame or by stretch extinction of a warm flame. A warm flame can also be formed by radiative extinction of a hot flame or ignition of a cool flame. A general €-shaped flammability diagram showing the burning limits of all three flame regimes at different oxygen mole fractions is obtained. The results show that thermal radiation, reactant concentration, temperature, and pressure all have significant impacts on the flammable regions of the three flame regimes. Increases in oxidizer temperature, oxygen concentration, and pressure shift the cool flame regime to higher stretch rates and cause the warm flame to have two extinction limits. At elevated temperatures, it is found that there is a direct transition between the hot flame and warm flame at low stretch rates. The results also show that, unlike the hot flame, the cool flame structure cannot be scaled by using pressure-weighted stretch rates due to the its significant reactant leakage and strong dependence of reactivity on pressure. The present results advance the understanding of near-limit flame dynamics and provide guidance for experimental observation of different flame regimes.     

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.     

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.     

辐射传热对微重力环境下预混火焰传播特性的影响

1引言在火焰中，辐射过程是一种重要的传热方式。对该过程尽可能精确的计算，对于改进燃烧设备的设计、改善设备的运行性能十分有益。在正常重力环境下，与其它的释热现象相比，预混火焰中的辐射热损失十分微弱，因而，过去对预混火焰的分析中，往往忽略了辐射热损失的影响。近年来，对微重力（ug）环境下的预混火焰的研究结果表明，可燃极限与＃s最小点火能无关，自媳灭火焰（SEFs）发生时；其释放的能量比通常观察到的点火极限时的能量大几个数量级山，因此火焰伸张并不能解释“g环境下观察到的实验结果，辐射热损失可能是影响＃g火焰可…     

Modeled quenching limits of spherical hydrogen diffusion flames

Hydrogen–air diffusion flames were modeled with an emphasis on kinetic extinction. The flames were one-dimensional spherical laminar diffusion flames supported by adiabatic porous burners of various diameters. Behavior of normal (H₂ flowing into quiescent air) and inverse (air flowing into quiescent H₂) configurations were considered using detailed H₂/O₂ chemistry and transport properties with updated light component diffusivities. For the same heat release rate, inverse flames were found to be smaller and 290 K hotter than normal flames. The weakest normal flame that could be achieved before quenching has an overall heat release rate of 0.25 W, compared to 1.4 W for the weakest inverse flame. There is extensive leakage of the ambient reactant for both normal and inverse flames near extinction, which results in a premixed flame regime for diffusion flames except for the smallest burners with radii on the order of 1 μm. At high flow rates H + OH(+M) → H₂O(+M) contributes nearly 50% of the net heat release. However at flow rates approaching quenching limits, H + O₂(+M) → HO₂(+M) is the elementary reaction with the largest heat release rate.     

On the existence of steady-state gaseous microgravity spherical diffusion flames in the presence of radiation heat loss

Whether steady-state gaseous microgravity spherical diffusion exist in the presence of radiation heat loss is an important fundamental question and has important implications for spacecraft fire safety. In this work, experiments aboard the International Space Station and a transient numerical model are used to investigate the existence of steady-state microgravity spherical diffusion flames. Gaseous spherical diffusion flames stabilized on a porous spherical burner are employed in normal (i.e., fuel flowing into an ambient oxidizer) and inverse (i.e., oxidizer flowing into an ambient fuel) flame configurations. The fuel is ethylene and the oxidizer oxygen, both diluted with nitrogen. The flow rate of the reactant gas from the burner is held constant. It is found that steady-state gaseous microgravity spherical diffusion flames can exist in the presence of radiation heat loss, provided that the steady-state flame size is less than the flame size for radiative extinction, and the flame develops fast enough that radiation heat loss does not drop the flame temperature below the critical temperature for radiative extinction (1130 K). A simple model is provided that allows for the identification of initial conditions that can lead to steady-state spherical diffusion flames. In the spherical, infinite domain configuration, the characteristic time for the diffusion-controlled system to effectively reach steady-state is found to be on the order of 100,000 s. Despite a narrow range of attainable conditions, flames that exhibit steady-state behavior are observed aboard the ISS for up to 870 s, even with the constraint of a finite boundary. Steady-state flames are simulated using the numerical model for over 100,000 s.     

2D computations of FREI with cool flames for n-heptane/air mixture

Two-dimensional axisymmetric numerical simulation reproduced flames with repetitive extinction and ignition (FREI) in a micro flow reactor with a controlled temperature profile with a stoichiometric n-heptane/air mixture, which have been observed in the experiment. The ignition of hot flame occurred from consumption reactions of CO that was remained in the previous cycle of FREI. Between extinction and ignition locations of hot flames, several other heat release rate peaks related to cool and blue flames were observed for the first time. After the extinction of the hot flame, cool flame by the low-temperature oxidation of n-heptane appeared first and was stabilized in a low wall temperature region. In the downstream of the stable cool flame, a blue flame by the consumption reactions of cool flame products of CH₂O and H₂O₂ appeared. After that, the hot flame ignition occurred from the remaining CO in the downstream of the blue flame. Then after the next hot flame ignition, the blue flame was swept away by the propagating hot flame. Soon before the hot flame merged with the stable cool flame, the hot flame propagation was intensified by the cool flame. After the hot flame merged with the stable cool flame, the hot flame reacted with the incoming fresh mixture of n-C₇H₁₆ and O₂.     

Effects of radiation heat loss on laminar premixed ammonia/air flames

Ammonia (NH₃) direct combustion is attracting attention for energy utilization without CO₂ emissions, but fundamental knowledge related to ammonia combustion is still insufficient. This study was designed to examine effects of radiation heat loss on laminar ammonia/air premixed flames because of their very low flame speeds. After numerical simulations for 1-D planar flames with and without radiation heat loss modeled by the optically thin model were conducted, effects of radiation heat loss on flame speeds, flame structure and emissions were investigated. Simulations were also conducted for methane/air mixtures as a reference. Effects of radiation heat loss on flame speeds were strong only near the flammability limits for methane, but were strong over widely diverse equivalence ratios for ammonia. The lower radiative flame temperature suppressed the thermal decomposition of unburned ammonia to hydrogen (H₂) at rich conditions. The equivalence ratio for a low emission window of ammonia and nitric oxide (NO) in the radiative condition shifted to a lower value than that in the adiabatic condition.     

Numerical study of buoyancy effects on the structure and propagation of triple flames

The structure and propagation properties of diffusion neutral triple flames subject to buoyancy effects are studied numerically using a high-accuracy scheme. A wide range of gravity conditions, heat release, and mixing widths for a scalar mixing layer are computed for downward-propagating (in the same direction as the gravity vector) and upward-propagating (in the opposite direction to the gravity vector) triple flames. These results are used to identify non-dimensional quantities, which parametrize the triple flame responses. Results show that buoyancy acts primarily to modify the overall span of the premixed branches in response to gas acceleration across the triple flame. The impact of buoyancy on the structure of triple flame is less pronounced than its impact on the topology of the branches. The trailing diffusion branch is affected by buoyancy primarily as a result of the changes in the overall flame size, which consequently modifies the rates of diffusion of excess fuel and oxidizer from the premixed branches to the diffusion branch. A simple analytical model for the triple flame speed, which accounts for both buoyancy and heat release is developed. Comparisons of the proposed model with the numerical results for a wide range of gravity, heat release and mixing width conditions, yield very good agreement. The analysis shows that under neutral diffusion, downward propagation reduces the triple flame speed, while upward propagation enhances it. For the former condition, a critical Froude number may be evaluated, which corresponds to a vanishing triple flame speed.     

微尺度扩散火焰特性的数值解析

本文以均匀空气流中圆管形成的甲烷射流扩散火焰为对象,用数值解析的方法研究了微尺度扩散火焰的火焰结构和燃烧特性。燃烧反应采用甲烷/空气一步总包反应,喷管壁面采用绝热条件。在Re一定情况下,改变喷口尺寸和喷口流速考察了微扩散火焰的结构和火焰熄灭的尺度效应。数值结果表明,随着喷口直径的增大,微火焰的上方出现回流; Re=12条件下,在喷口直径=0.07 mm时存在熄灭极限;稳定燃烧区的最小发热率约为0.5 W;微尺度条件下,Da数对火焰结构和火焰的熄灭有一定的影响。     

圆柱火焰的可燃极限及稳定性分析

本文从理论上分析了有辐射热损失和曲率的圆柱火焰,推导出了关于火焰位置、火焰温度同热损失和来流速度之间的关系式。并在此基础上对圆柱火焰的可燃极限进行了研究,结果表明热损失对可燃极限产生很大的影响。另外,作者运用线性稳定分析法对有辐射热损失的圆柱火焰作了稳定性分析,得出了判断圆柱火焰稳定与否的通用表达式。     

A comparison of computational and experimental lift-off heights of coflow laminar diffusion flames

As a sensitive marker of changes in flame structure, the number densities of excited-state CH (denoted CH^*), and excited-state OH (denoted OH^*) are imaged in coflow laminar diffusion flames. Measurements are made both in normal gravity and on the NASA KC-135 reduced-gravity aircraft. The spatial distribution of these radicals provides information about flame structure and lift-off heights that can be directly compared with computational predictions. Measurements and computations are compared over a range of buoyancy and fuel dilution levels. Results indicate that the lift-off heights and flame shapes predicted by the computations are in excellent agreement with measurement for both normal gravity (1g) and reduced gravity flames at low dilution levels. As the fuel mixture is increasingly diluted, however, the 1g lift-off heights become underpredicted. This trend continues until the computations predict stable flames at highly dilute fuel mixtures beyond the 1g experimental blow-off limit. To better understand this behavior, an analysis was performed, which indicates that the lift-off height is sensitive to the laminar flame speed of the corresponding premixed mixture at the flame edge. By varying the rates of two key “flame speed” controlling reactions, we were able to modify the predicted lift-off heights so as to be in closer agreement with the experiments. The results indicate that reaction sets that work well in low dilution systems may need to be modified to accommodate high dilution flames.     

Unsteady effects on flame extinction limits during gaseous and two-phase flame/vortex interactions

In highly fluctuating flows, it happens that high values of the strain-rate do not induce extinction of the flame front. Unsteady effects minimize the flame response to rapidly varying strain fields. In the present study, the effects of time-dependent flows on non-premixed flames are investigated during flame/vortex interactions. Gaseous flames and spray flames in the external sheath combustion regime are considered. To analyse the flame/vortex interaction process, the velocity field and the flame geometry are simultaneously determined using particle imaging velocimetry and laser-induced fluorescence of the CH radical. The influence of vortex flows on the extinction limits for different vortex parameters and for different gaseous and two-phase flames is examined. If the external perturbation is applied over an extended period of time, the extinction strain-rate is that corresponding to the steady-state flame, and this critical value mainly depends on the fuel and oxidizer compositions and the injection temperature. If the external perturbation is applied during a short period of time, extinction occurs at strain-rates above the steady-state extinction strain-rate. This deviation appears for flow fluctuation timescales below steady flame diffusion timescales. This behaviour is induced by diffusive processes, limiting the ability of the flame to respond to highly fluctuating flows. With respect to unsteady effects, the spray flames investigated in this article behave essentially like gaseous flames, because evaporation takes place in a thin layer before the flame front. Extinction limits are only slightly modified by the spray, controlling process being the competition between aerodynamic and diffusive timescales.     

Theoretical and experimental studies on mesoscale flame propagation and extinction

Mesoscale flame propagation and extinction of premixed flames in channels are investigated theoretically and experimentally. Emphasis is placed on the effect of wall heat loss and the wall–flame interaction via heat recirculation. At first, an analytical solution of flame speed in mesoscale channels is obtained. The results showed that channel width, flow velocity, and wall thermal properties have dramatic effects on the flame propagation and lead to multiple flame regimes and extinction limits. With the decrease in channel width, there exist two distinct flame regimes, a fast burning regime and a slow burning regime. The existence of the new flame regime and its extended flammability limit render the classical quenching diameter inapplicable. Furthermore, the results showed that at optimum conditions of flow velocity and wall thermal properties, mesoscale flames can propagate faster than the adiabatic flame. Second, numerical simulation with detailed chemistry demonstrated the existence of multiple flame regimes. The results also showed that there is a non-linear dependence of the flame speed on equivalence ratio. Moreover, it is shown that the Nusselt number has a significant impact on this non-linear dependence. Finally, the non-linear dependence of flame speed on equivalence ratio for both flame regimes is measured using a C₃H₈–air mixture. The results are in good agreement with the theory and numerical simulation.     

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

Flame spread over solids in buoyant and forced concurrent flows: Model computations and comparison with experiments

A detailed three-dimensional concurrent flame spread model is tested and compared with three sets of experiments. The parameters varied include: gravity, flow velocity, pressure, oxygen mole fraction, and sample width. In buoyant flows (normal and partial gravities), the computed steady spread rate and flame profiles agree favorably with experiment. The predicted extinction limits are lower but can be improved. Comparison in forced concurrent flow in microgravity shows correct trends. The predicted steady spread rates are lower than the experimental ones if the flames are short but higher than the experimental ones if the flames are long. It is believed that the experimental flames may not have fully reached steady state at the end of the 5-s microgravity drops. Longer duration microgravity experiments in future will be needed to substantiate this belief.     

Gravity effects on partially premixed flames: an experimental-numerical investigation

While premixed and nonpremixed microgravity flames have been extensively investigated, the corresponding literature regarding partially premixed flames (PPFs) is sparse. We report the first experimental investigation of burner-stabilized microgravity PPFs. Partially premixed flames with multiple reaction zones are established in microgravity on a Wolfhard–Parker slot burner in the 2.2 s drop tower at the NASA Glenn Research Center. Microgravity measurements include flame imaging, and thermocouple and radiometer data. Detailed simulations are also used to provide further insight into the steady and transient response of these flames to variations in g. The flame topology and interactions between the various reaction zones are strongly influenced by gravity. The flames widen substantially in microgravity. During the transition from normal to microgravity, the flame structure experiences a fast change and another relatively slower transient change. The fast response is due to the altered advection as the value of g is reduced, while the slow response is due to the changes in the diffusive fluxes. The radiative heat loss from the flames increases in microgravity. A scaling analysis based on a radiation Damköhler number is able to characterize the radiation heat loss.     

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.