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.     

An experimental and numerical study on the adequacy of CH as a flame marker in premixed methane flames

The CH radical is frequently used as a flame marker because it is relatively short-lived and is present over a narrow region in flames. Discontinuities in the CH field are thus often interpreted as localized extinction of the flame. Recently, however, the adequacy of CH laser-induced fluorescence (LIF) as a flame marker was questioned by an experimental study of flame–vortex interactions in highly N₂-diluted premixed methane flames. We demonstrate both experimentally and numerically that anomalies in the transient response of CH in this earlier study were due to reactant composition variations in the vortex. In addition, we evaluate the adequacy of CH LIF as a flame marker over a much broader range of conditions. Previous numerical studies showed that heat release rate correlates reasonably well with peak [HCO] and the concentration product [OH][CH₂O], but poorly with [CH], in highly N₂-diluted premixed methane flames. Here, the correlation between heat release rate and CH is investigated both experimentally, by performing simultaneous measurements of CH, OH, and CH₂O LIF, and numerically. We consider undiluted and N₂-diluted premixed methane flames over a range of strain rates and stoichiometries. Results are reported for flames subjected to unsteady stretch and reactant composition variations. For all N₂-dilution levels considered, the peak CH LIF signal correlates poorly with heat release rate when the stoichiometry of the reactant mixture changes from rich to lean. However, when flames are subjected to stretch, the correlation between CH and heat release rate improves as the N₂-dilution level decreases. The correlation is reasonably good for undiluted flames with equivalence ratios of 0.8 < Φ < 1.2. This result is particularly encouraging, given the relevance of undiluted flames to practical applications, and it motivates further investigation of the parameter space for which difficulties may exist in using CH as a flame marker.     

Simultaneous measurement of localized heat release with OH/CH2O-LIF imaging and spatially integrated OH∗ chemiluminescence in turbulent swirl flames

The in-situ and localized observation of heat release in turbulent flames is important for the validation of computational modeling of turbulent flows with combustion. In the present work we obtain localized information on heat release rate (HRR) by the commonly accepted technique of the simultaneous and single-shot planar imaging of OH and CH₂O concentrations by laser-induced fluorescence (LIF). Additionally, we combine this with the simultaneous line-of-sight and temporally resolved chemiluminescence detection of OH^?, spatially integrated within the flame volume, interrogated by the laser sheets used for the HRR imaging technique. The combined diagnostic methods are demonstrated for a swirl-stabilized, premixed turbulent methane/air flame of 30-kW thermal power, and they show the existence of correlations between both HRR-sensitive diagnostic techniques.     

Simulations of laminar non-premixed flames of kerosene with hot combustion products as oxidiser

The effects of hot combustion product dilution in a pressurised kerosene-burning system at gas turbine conditions were investigated with laminar counterflow flame simulations. Hot combustion products from a lean (φ = 0.6) premixed flame were used as an oxidiser with kerosene surrogate as fuel in a non-premixed counterflow flame at 5, 7, 9 and 11 bar. Kerosene-hot product flames, referred to as ‘MILD’, exhibit a flame structure similar to that of kerosene–air flames, referred to as ‘conventional’, at low strain rates. The Heat Release Rate (HRR) of both conventional and MILD flames reflects the pyrolysis of the primary and intermediate fuels on the rich side of the reaction zone. Positive HRR and OH regions in mixture fraction space are of similar width to conventional kerosene flames, suggesting that MILD flames are thin fronts. MILD flames do not exhibit typical extinction behaviour, but gradually transition to a mixing solution at very high rates of strain (above A = 160, 000 s^?1 for all pressures). This is in agreement with literature that suggests heavily preheated and diluted flames have a monotonic S-shaped curve. Despite these differences in comparison with kerosene–air flames, MILD flames follow typical trends as a function of both strain and pressure. Further still, the peak locations of the overlap of OH and CH₂O mass fractions in comparison with the peak HRR indicate that the pixel-by-pixel product of OH- and CH₂O-PLIF signals is a valid experimental marker for non-premixed kerosene MILD and conventional flames.     

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

Extinction strain rates of premixed ammonia/hydrogen/nitrogen-air counterflow flames

Chemical energy vectors will play a crucial role in the transition of the global energy system, due to their essential advantages in storing energy in form of gaseous, liquid, or solid fuels. Ammonia (NH₃) has been identified as a highly promising candidate, as it is carbon-free, can be stored at moderate pressures, and already has a developed distribution infrastructure. As a fuel NH₃ has poor combustion properties that can be improved by the addition of hydrogen, which can be obtained energy-efficiently by partially cracking ammonia into hydrogen (H₂) and nitrogen (N₂) prior to the combustion process. The resulting NH₃/H₂/N₂ blend leads to significantly improved flame stability and resilience to strain-induced blow-out, despite similar laminar flame properties compared to equivalent methane/air flames. This study reports the first measurements of extinction strain rates, measured using the premixed twin-flame configuration in a laminar opposed jet burner, for two NH₃/H₂/N₂ blends over a range of equivalence ratios. Local strain rates are measured using particle tracking velocimetry (PTV) and are related to the inflow conditions, such that the local strain rate at the extinction point can be approximated. The results are compared with 1D-simulations using three recent kinetic mechanisms for ammonia oxidation. By relating the extinction strain rates to laminar flame properties of the unstretched flame, a comparison of the extinction behaviour of CH₄ and NH₃/H₂/N₂ blends can be made. For lean mixtures, NH₃/H₂/N₂-air flames show a significant higher extinction resistance in comparison to CH₄/air. In addition, a strong non-linear dependence between the resistance to extinction and equivalence ratio for NH₃/H₂/N₂ blends is observed.     

Effect of Soret diffusion on lean hydrogen/air flames at normal and elevated pressure and temperature

The influence of Soret diffusion on lean premixed flames propagating in hydrogen/air mixtures is numerically investigated with a detailed chemical and transport models at normal and elevated pressure and temperature. The Soret diffusion influence on the one-dimensional (1D) flame mass burning rate and two-dimensional (2D) flame propagating characteristics is analysed, revealing a strong dependency on flame stretch rate, pressure and temperature. For 1D flames, at normal pressure and temperature, with an increase of Karlovitz number from 0 to 0.4, the mass burning rate is first reduced and then enhanced by Soret diffusion of H₂ while it is reduced by Soret diffusion of H. The influence of Soret diffusion of H₂ is enhanced by pressure and reduced by temperature. On the contrary, the influence of Soret diffusion of H is reduced by pressure and enhanced by temperature. For 2D flames, at normal pressure and temperature, during the early phase of flame evolution, flames with Soret diffusion display more curved flame cells. Pressure enhances this effect, while temperature reduces it. The influence of Soret diffusion of H₂ on the global consumption speed is enhanced at elevated pressure. The influence of Soret diffusion of H on the global consumption speed is enhanced at elevated temperature. The flame evolution is more affected by Soret diffusion in the early phase of propagation than in the long run due to the local enrichment of H₂ caused by flame curvature effects. The present study provides new insights into the Soret diffusion effect on the characteristics of lean hydrogen/air flames at conditions that are relevant to practical applications, e.g. gas engines and turbines.     

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.     

Effect of H-atom concentration on soot formation in premixed ethylene/air flames

Soot formation is a major challenge in the development of clean and efficient combustion systems based on hydrocarbon fuels. Fundamental understanding of the reaction mechanism leading to soot formation can be obtained by investigating the role of key reactive species such as atomic hydrogen taking part in soot formation pathways. In this study, two-dimensional laser induced incandescence (LII) measurements using λ?=?1064?nm laser have been used to measure soot volume fraction (f_V) in a series of rich ethylene (C₂H₄)/air flames, stabilized over a McKenna burner fitted with a flame stabilizing metal disc. Moreover, a comparison of UV (λ?=?283?nm), visible (λ?=?532?nm) and IR (λ?=?1064?nm) laser excited LII measurements of soot is discussed. Recently developed, femtosecond two-photon laser-induced fluorescence (fs-TPLIF) technique has been applied for obtaining spatially resolved H-atom concentration ([H]) profiles under the same flame conditions. The structure of the flames has also been determined using hydroxyl radical (OH) planar laser induced fluorescence (PLIF) imaging. The results indicate an inverse dependence of f_V on [H] for a range of C₂H₄/air rich flames up to an equivalence ratio, Φ?=?3.0. Although an absolute relationship between [H] and f_V cannot be easily derived owing to the multiple steps involving H and other intermediate species in soot formation pathways, the present study demonstrates the feasibility to couple [H] and f_V obtained using advanced optical techniques for soot formation studies.     

Dilution effects analysis of opposed-jet H2/CO syngas diffusion flames

This paper reported the analysis of dilution effects on the opposed-jet H₂/CO syngas diffusion flames. A computational model, OPPDIF coupled with narrowband radiation calculation, was used to study one-dimensional counterflow syngas diffusion flames with fuel side dilution from CO₂, H₂O and N₂. To distinguish the contributing effects from inert, thermal/diffusion, chemical, and radiation effects, five artificial and chemically inert species XH₂, XCO, XCO₂, XH₂O and XN₂ with the same physical properties as their counterparts were assumed. By comparing the realistic and hypothetical flames, the individual dilution effects on the syngas flames were revealed. Results show, for equal-molar syngas (H₂/CO = 1) at strain rate of 10 s^?1, the maximum flame temperature decreases the most by CO₂ dilution, followed by H₂O and N₂. The inert effect, which reduces the chemical reaction rates by behaving as the inert part of mixtures, drops flame temperature the most. The thermal/diffusion effect of N₂ and the chemical effect of H₂O actually contribute the increase of flame temperature. However, the chemical effect of CO₂ and the radiation effect always decreases flame temperature. For flame extinction by adding diluents, CO₂ dilution favours flame extinction from all contributing effects, while thermal/diffusion effects of H₂O and N₂ extend the flammability. Therefore, extinction dilution percentage is the least for CO₂. The dilution effects on chemical kinetics are also examined. Due to the inert effect, the reaction rate of R84 (OH+H₂ = H+H₂O) is decreasing greatly with increasing dilution percentage while R99 (CO+OH→CO₂+H) is less affected. When the diluents participate chemically, reaction R99 is promoted and R84 is inhibited with H₂O addition, but the trend reverses with CO₂ dilution. Besides, the main chain-branching reaction of R38 (H+O₂→O+OH) is enhanced by the chemical effect of H₂O dilution, but suppressed by CO₂ dilution. Relatively, the influences of thermal/diffusion and radiation effects on the reaction kinetics are then small.     

An investigation of flame zones and burning velocities of laminar unconfined methane-oxygen premixed flames

Numerical and experimental investigations of unconfined methane-oxygen laminar premixed flames are presented. In a lab-scale burner, premixed flame experiments have been conducted using pure methane and pure oxygen mixtures having different equivalence ratios. Digital photographs of the flames have been captured and the radial temperature profiles at different axial locations have been measured using a thermocouple. Numerical simulations have been carried out with a C2 chemical mechanism having 25 species and 121 reactions and with an optically thin radiation sub-model. The numerical results are validated against the experimental and numerical results for methane-air premixed flames reported in literature. Further, the numerical results are validated against the results from the present methane-oxygen flame experiments. Visible regions in digital flame photographs have been compared with OH isopleths predicted by the numerical model. Parametric studies have been carried out for a range of equivalence ratios, varying from 0.24 to 1.55. The contours of OH, temperature and mass fractions of product species such as CO, CO₂ and H₂O, are presented and discussed for various cases. By using the net methane consumption rate, an estimate of the laminar flame speed has been obtained as a function of equivalence ratio.     

Chemical structure of atmospheric pressure premixed laminar formic acid/hydrogen flames

The work presents an experimental and kinetic modeling study of laminar premixed formic acid [HC(O)OH]/H₂/O₂/Ar flames at different equivalence ratios (φ=0.85, 1.1 and 1.3) stabilized on a flat burner at atmospheric pressure, as well as laminar flame speed of HC(O)OH/O₂/Ar flames (φ=0.5–1.5) at 1 atm. Flame structure as well as laminar flame speed were simulated using three different detailed chemical kinetic mechanisms proposed for formic acid oxidation. The components in the fuel blends show different consumption profiles, namely, hydrogen is consumed slower than formic acid. According to kinetic analysis, the reason of the observed phenomenon is that the studied flames have hydrogen as a fuel but also as an intermediate product formed from HC(O)OH decomposition. Comparison of the measured and simulated flame structure shows that all the mechanisms satisfactorily predict the mole fraction profiles of the reactants, main products, and intermediates. It is noteworthy that the mechanisms proposed by Glarborg et al., Konnov et al. and the updated AramcoMech2.0 adequately predict the spatial variations in the mole fractions of free radicals, such as H, OH O and HO₂. However, some drawbacks of the mechanisms used were identified; in particular, they predict different concentrations of CH₂O. As for laminar flame speed simulations, the Konnov et al. mechanism predicts around two times higher values than in experiment, while the Glarborg et al. and updated AramcoMech2.0 show good agreement with the experimental data.     

Numerical simulation of premixed H2–air cellular tubular flames

The detailed flame structure of laminar premixed cellular flames in the tubular domain is simulated in 2D using a fully-implicit primitive variable finite difference formulation that includes multicomponent transport and detailed chemical kinetics. Numerical results for H₂/air flames are presented and compared against spatially resolved experimental measurements of temperature and chemical species including atomic H and OH. The experimental results compare well for flame structure and cell number, despite the numerical model under-predicting the peak temperature by 200 K. Numerical experiments were performed to assess the ability for cellular tubular flames to impact experimental and numerical investigations of practical flames. The cellular flame structure is found to provide a highly sensitive geometry that is useful for validating diffusive transport modelling approximations. This capability is exemplified through the development of a simple and accurate approximation for thermal diffusion (i.e. the Soret effect) that is suitable for practical combustion codes.     

稀释剂对H2/CH4预混火焰熄灭特性影响

掺氢天然气在稀释气体作用下的熄灭特性研究对实际燃烧设备的设计和优化具有重要的指导意义。本文利用对冲火焰法测量了掺氢天然气层流火焰在N₂和CO₂作用下的熄灭拉伸率,并采用数值模拟耦合详细化学反应机理对N₂,CO₂和He的稀释剂效应展开研究。结果表明,Li、GRI Mech 3.0和FFCM-1机理均能定性反映燃料熄灭拉伸率随当量比的变化规律,且FFCM-1机理综合预测精度最高。实验和模拟发现,不同稀释剂气体对掺氢天然气熄灭拉伸率降低幅度满足:He22。进一步研究发现,CO₂由于热容大,在反应体系中会降低火焰温度,同时增强了链终止反应强度,通过热效应和化学效应两方面对火焰熄灭特性起作用。He则能显著改变燃料混合物的平均摩尔质量,从而改变体系中重要反应物和自由基的扩散特性,从扩散效应方面影响火焰的熄灭特性。     

Structure of partially premixed n-heptane–air counterflow flames

To avoid the complexities associated with the droplet/vapor transport and nonuniform evaporation processes, a fundamental investigation of liquid fuel combustion in idealized configurations is very useful. An experimental–computational investigation of prevaporized n-heptane nonpremixed and partially premixed flames established in a counterflow burner is described. There is a general agreement between various facets of our nonpremixed flame measurements and the literature data. The partially premixed flames are characterized by a double flame structure. This becomes more distinct as the strain rate decreases and partial premixing increases, which also increases the separation distance between the two reaction zones. The peak partially premixed flame temperature increases with increasing premixing of the fuel stream. The peak CO₂ and H₂O concentrations are relatively insensitive to partial premixing. The CO and H₂ peak concentrations on the premixed side increase as the fuel-side equivalence ratio decreases. These species are transported to the nonpremixed reaction zone where they oxidize. The C₂ species have peaks in the premixed reaction zone. The concentrations of olefins are ten times larger than those of the corresponding paraffins. The oxidizer is present in partially premixed flames throughout the combustion system and there are no regions characterized by simultaneous high temperature and high fuel concentration. As a result, pyrolysis reactions leading to soot formation are greatly diminished.     

Laser-induced fluorescence determination of flame temperatures in comparison with CARS measurements

Temperature profiles in several premixed low pressure H₂/O₂/N₂ flames and in an atmospheric pressure CH₄/air flame were determined by laser-induced fluorescence (LIF) and by CARS experiments. In the LIF study, temperatures were derived from OH excitation spectra, CARS temperatures were deduced from N₂ Q-branch spectra. The present study is the first quantitative comparison of these two methods for temperature determination in flames burning at pressures up to 1 bar. The resulting temperatures showed good agreement.     

Lean or ultra-lean stretched planar methane/air flames

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

Attenuation of combustion instability in a fuel-staged dual-nozzle gas turbine combustor with asymmetric hydrogen composition

The instability attenuation mechanism of fuel staging was investigated in a CH₄/H₂ fueled dual-nozzle gas turbine combustor. Fuel staging was implemented using an asymmetry in fuel composition between the two nozzles. The fuel composition of the upper nozzle was varied while keeping that of the lower nozzle constant. Under these conditions, the self-excited and forced responses of fuel-staged flames were analyzed using OH* chemiluminescence imaging, OH planar laser-induced fluorescence, and particle image velocimetry. In the self-excited measurements, although strong combustion instability was exhibited in the symmetric condition, it weakened gradually with increasing asymmetry in fuel composition. The symmetric flame exhibited significant fluctuations in the heat release rate around the flame tip, which acted as the primary cause of driving combustion instability. However, in asymmetric flames, the H₂ addition induced phase leads in heat release rate fluctuations at the upper region, which damped combustion instability. Thus, our observations revealed a high correlation between the phase leads and the attenuation of combustion instability. Analyses of the forced responses showed that the heat release rate fluctuations were induced by interactions between the flame and the shedding vortex released from the nozzle tip into the downstream. Although these characteristics of shedding vortices did not depend on the H₂ addition, the change in the axial position of the flame caused by the H₂ addition induced the relocation of the site, at which the flame interacted with the vortex. Subsequently, it induced phase leads in the heat release rate fluctuations. The phase difference of heat release rate fluctuations between the two flames due to this phase leads enlarged progressively with increasing asymmetry in fuel composition, leading to the attenuation of combustion instability in asymmetric conditions.