Experimental and numerical investigation of low-pressure laminar premixed synthetic natural gas/O2/N2 and natural gas/H2/O2/N2 flames

The oxidation of laminar premixed natural gas flames has been studied experimentally and computationally with variable mole fractions of hydrogen (0, 20, and 60%) present in the fuel mixture. All flames were operated at low pressure (0.079 atm) and at variable overall equivalence ratios (0.74<?<1.0) with constant cold gas velocity. At the same global equivalence ratio, there is no significant effect of the replacement of natural gas by 20% of H₂. The small differences recorded for the intermediate species and combustion products are directly due to the decrease of the amount of initial carbon. However, in 60% H₂ flame, the reduction of hydrocarbon species is due both to kinetic effects and to the decrease of initial carbon mole fraction. The investigation of natural gas and natural gas/hydrogen flames at similar C/O enabled identification of the real effects of hydrogen. It was shown that the presence of hydrogen under lean conditions activated the H-abstraction reactions with H atoms rather than OH and O, as is customary in rich flames of neat hydrocarbons. It was also demonstrated that the presence of H₂ favors CO formation.     

Effect of differential diffusion on turbulent lean premixed hydrogen enriched flames through structure analysis

This study presents the flame structure influenced by the differential diffusion effects and evaluates the structural modifications induced by the turbulence, thus to understand the coupling effects of the diffusively unstable flame fronts and the turbulence distortion. Lean premixed CH₄/H₂/air flames were conducted using a piloted Bunsen burner. Three hydrogen fractions of 0, 30% and 60% were adopted and the laminar flame speed was kept constant. The turbulence was generated with a single-layer perforated plate, which was combined with different bulk velocities to obtain varied turbulence intensities. Quasi-laminar flames without the plate were also performed. Explicit flame morphology was obtained using the OH-PLIF. The curvature, flame surface density and turbulent burning velocity were measured. Results show that the preferential transport of hydrogen produces negatively curved cusps flanked with positively curved bulges, which are featured by skewed curvature pdfs and consistent with the typical structure caused by the Darrieus-Landau instability. Prevalent bulge-cusp like wrinkles remain with relatively weak turbulence. However, stronger turbulence can break the bulges to be finer, and induce random positively curved cusps, therefore to destroy the bulge-cusp structures. Evident positive curvatures are generated in this process modifying the skewed curvature pdfs to be more symmetric, while the negative curvatures are not affected seriously. From low to high turbulence intensities, the hydrogen addition always strengthens the flame wrinkling. The augmentation of flame surface density and turbulent burning velocity with hydrogen is even more obvious at higher turbulence intensity. It is suggested that the differential diffusion can persist and even be strengthened with strong turbulence.     

NOx emission characteristics of partially premixed laminar flames of H2/CO/CO2 mixtures

NO_x emission indices were experimentally measured for partially premixed laminar flames of five different H₂/CO/CO₂ fuels over a wide range of equivalence ratios. Through those fuels, the effects of H₂/CO ratio and CO₂ concentration on NO_x emissions, flame appearance, visible flame height and flame temperature are presented. EINO_x values increase when 1.0 ≤ Φ ≤ 1.6, then remain near the highest value, before decreasing slowly when 3.85 ≤ Φ ≤ ∞. The increase of the CO₂ concentration reduces the EINO_x for the whole range of equivalence ratios, while the increase in the H₂/CO ratio reduces the EINO_x when Φ ≤ 2.0 and is inconsequential for richer mixtures. The variation in flame temperatures approximates EINO_x trends. The variation of flame color from blue to orange when the H₂/CO ratio is increased might be explained by higher CO levels in by-product combustion.     

Quantitative HCN measurements in CH4/N2O/O2/N2 flames using mid-infrared polarization spectroscopy

Quantitative measurements of hydrogen cyanide (HCN) were nonintrusively performed using mid-infrared polarization spectroscopy (IRPS) in atmospheric pressure flames. The lifted flat, laminar, premixed CH₄/N₂O/O₂/N₂ flames were stabilized on a 7 cm diameter home-built McKenna-type burner with variable proportion of N₂O and O₂. The characteristic spectral structure of HCN molecules was identified in the rotational line-resolved IRPS spectra collected in flames at around 3248 cm⁻¹. The P20 line belonging to the CH stretching band was chosen for quantitative measurements and the line-integrated IRPS signal was recorded in a series of fuel-rich CH₄/O₂/N₂O/N₂ flames with equivalence ratios of 1.2, 1.4 and 1.6. Absolute mole fractions of HCN molecules in these flames were obtained through in situ calibration of the optical system with nonreactive gas flow of N₂ seeded with known amount of HCN on the same burner. Moreover, the experimental results were compared with calculations performed using the Konnov detailed C/H/N/O mechanism, which implements NCN prompt-NO reactions. Generally good agreement was found with some discrepancies indicating the need for further model improvement.     

Onset of cellular instability in adiabatic H2/O2/N2 premixed flames anchored to a flat-flame heat-flux burner

The onset of cellular instability in adiabatic H₂/O₂/N₂ premixed flames anchored to a heat-flux burner is investigated numerically. Both hydrodynamic instability and diffusional-thermal instability are shown to play an important role in the onset of cellular flames. The burner can effectively suppress cellular instability when the flames are close to the burner, otherwise the burner can suppress the instabilities only at large wavenumbers. Because of differential diffusion, local extinction can occur in lean H₂/O₂/N₂ flames. When the flames develop to take on cellular shapes, the surface length, the overall heat release rate and the mean burning velocity are all increased. For near stoichiometric fuel-rich flames the mean burning velocity can increase by as much as 20%–30%. For lean flames with an equivalence ratio of 0.56, the mean burning velocity can be 2–3 times of the burning velocity of the corresponding planar flame.     

Laminar flame speed studies of lean premixed H2/CO/air flames

Laminar flame speeds of lean premixed H₂/CO/air mixtures were measured in the counterflow configuration over a wide range of H₂ content at lean conditions. The values were determined by extrapolating the referenced flame speed to zero stretch rate using the non-linear extrapolation method to reduce the systematic error. Detailed calculation of laminar flame speed was also conducted using PREMIX code coupled with three different kinetic models. In general, simulation results agreed well with the experimental data. Both the experimental and calculation results revealed that the laminar flame speeds of lean premixed H₂/CO/air mixtures increased with H₂ content significantly when H₂ content was small (?15%) and gradually when H₂ content was large (>15%).     

Extinction of premixed H2/air flames: Chemical kinetics and molecular diffusion effects

Laminar flame speed has traditionally been used for the partial validation of flame kinetics. In most cases, however, its accurate determination requires extensive data processing and/or extrapolations, thus rendering the measurement of this fundamental flame property indirect. Additionally, the presence of flame front instabilities does not conform to the definition of laminar flame speed. This is the case for Le<1 flames, with the most notable example being ultralean H₂/air flames, which develop cellular structures at low strain rates so that determination of laminar flame speeds for such mixtures is not possible. Thus, this low-temperature regime of H₂ oxidation has not been validated systematically in flames. In the present investigation, an alternative/supplemental approach is proposed that includes the experimental determination of extinction strain rates for these flames, and these rates are compared with the predictions of direct numerical simulations. This approach is meaningful for two reasons: (1) Extinction strain rates can be measured directly, as opposed to laminar flame speeds, and (2) while the unstretched lean H₂/air flames are cellular, the stretched ones are not, thus making comparisons between experiment and simulations meaningful. Such comparisons revealed serious discrepancies between experiments and simulations for ultralean H₂/air flames by using four kinetic mechanisms. Additional studies were conducted for lean and near-stoichiometric H₂/air flames diluted with various amounts of N₂. Similarly to the ultralean flames, significant discrepancies between experimental and predicted extinction strain rates were also found. To identify the possible sources of such discrepancies, the effect of uncertainties on the diffusion coefficients was assessed and an improved treatment of diffusion coefficients was advanced and implemented. Under the conditions considered in this study, the sensitivity of diffusion coefficients to the extinction response was found to be significant and, for certain species, greater than that of the kinetic rate constants.     

Experimental and computational study on partially premixed flames in a centerbody burner

The centerbody burner was designed with the objective of understanding the coupled processes of soot formation, growth, and burnout. Fuel that issues from the center of the burner establishes two flame zones – one associated with the recirculation zone (RZ) and the other, with the trailing jet. The sooting characteristics in these two flame zones can be quite different because of variations in residence time and transport of reactants and products. Calculations performed for this burner operating under a partially premixed fuel jet suggested that soot in the RZ decreases and that soot in the trailing jet flame increases with the amount of premixing. An experimental and numerical study is performed to aid the understanding of these differences. A time-dependent, axisymmetric, detailed-chemistry computational-fluid-dynamics (CFD) model known as Unsteady Ignition and Combustion using ReactioNs (UNICORN) is used for simulating flames under different equivalence-ratio conditions. Combustion and PAH formation are modeled using the Wang–Frenklach (99 species and 1066 reactions) mechanism, and soot is simulated using a two-equation model of Lindstedt. A Lagrangian-based particle-tracking model is used for understanding the evolution of soot-like particles. Flame and recirculation-zone structures and soot in the experiments are identified using direct photographs taken with and without Mie scattering from soot particles as well as laser-induced-incandescence (LII) measurements. Calculations predict the structures of the partially premixed centerbody flames for various equivalence ratios reasonably well. Experiments confirm the predicted soot suppression in the RZs and enhancement of soot in the trailing jet flame when air is added to the fuel jet. It is found that flame movement in the RZ increases soot-particle burnout and, thereby, reduces the amount of soot within the RZ. As the flame moves closer to the fuel jet, more soot becomes entrained into the inner vortex. Motion of soot-like particles explained the spiral rings observed in the experiment. Increased particle burnout with partial premixing leads to shrinkage of soot spirals.     

Effect of high hydrogen enrichment on the outer-shear-layer flame of confined lean premixed CH4/H2/air swirl flames

In this study, we investigated the H₂-induced transition of confined swirl flames from the “V” to “M” shape. H₂-enriched lean premixed CH₄/H₂/air flames with H₂ fractions up to 80% were conducted. The flame structure was obtained with Planar Laser-Induced Fluorescence (PLIF) of the OH radical. Flow fields were measured with Particle Image Velocimetry (PIV). It was observed that the flame tip in the outer shear layer gradually propagated upstream and finally anchored to the injector with the hydrogen fractions increase, yielding the transition from the “V” to “M” flame. We examined the flame structures and the flame flow dynamics during the transition. The shape transition was directly related to the evolution of the corner flame along the outer shear layer. With H₂ addition, the outer recirculation zone first appeared downstream where the corner flame started to propagate upstream; then, the recirculation zone expanded upward to form a stable “M” flame gradually. The flow straining was observed to influence the stabilization of the outer shear layer flame significantly. This study can be useful for the understanding of recirculation-stabilized swirling flames with strong confinement. The flame structure and the flow characteristics of flames with a high H₂ content are also valuable for model validation.     

Effect of hydrogen enrichment on flame broadening of turbulent premixed flames in thin reaction regime

An experimental study to identify the effect of hydrogen enrichment and differential diffusion on the flame broadening is conducted. Turbulent lean premixed flames in the Broadened Preheat–Thin Reaction (BP-TR) regime are obtained. The flames are stabilized on a Bunsen burner and CH₄/H₂/air mixtures are adopted with three hydrogen fractions of 0, 30% and 60%. The preheat zone and heat release zone are captured with the multi-species Planar Laser-Induced Fluorescence (PLIF) of OH and CH₂O radicals. Flame thicknesses of the preheat and heat release layers are measured. Results show broadened preheat zone and thin heat release layers for the flames, as predicted by the BP-TR regime. The preheat zone thickness can be increased to about 3–6 times compared to the laminar preheat thickness. An apparently decreased preheat zone thickness with hydrogen addition is observed. The differential diffusion is anticipated to locally thicken the heat release zone along the flame front. The mean heat release thickness is nearly not affected by the turbulence or hydrogen addition.     

An experimental and numerical study on characteristics of laminar premixed H2/CO/CH4/air flames

The effects of variations in the fuel composition on the characteristics of H₂/CO/CH₄/air flames of gasified biomass are investigated experimentally and numerically. Experimental measurements and numerical simulations of the flame front position and temperature are performed in the premixed stoichiometric H₂/CO/CH₄/air opposed-jet flames with various H₂ and CO contents in the fuel. The adiabatic flame temperatures and laminar burning velocities are calculated using the EQUIL and PREMIX codes of Chemkin collection 3.5, respectively. Whereas the flame structures of the laminar premixed stoichiometric H₂/CO/CH₄/air opposed-jet flames are simulated using the OPPDIF package with the GRI-Mech 3.0 chemical kinetic mechanisms and detailed transport properties. The measured flame front position and temperature of the stoichiometric H₂/CO/CH₄/air opposed-jet flames are closely predicted by the numerical calculations. Detailed analysis of the calculated chemical kinetic structures reveals that the reaction rate of reactions (R38), (R46), and (R84) increase with increasing H₂ content in the fuel mixture. It is also found that the increase in the laminar flame speed with H₂ addition is most likely due to an increase in active radicals during combustion (chemical effect), rather than from changes in the adiabatic flame temperature (thermal effect). Chemical kinetic structure and sensitivity analyses indicate that for the stoichiometric H₂/CO/CH₄/air flames with fixed H₂ concentration in the fuel mixture, the reactions (R99) and (R46) play a dominant role in affecting the laminar burning velocity as the CO content in the fuel is increased.     

Laminar propagation of lean premixed flames ignited in stratified mixture

Combustion in stratified mixtures is envisaged in practical energy systems such as direct-injection spark-ignited (DISI) car engines, gas turbines, for reducing CO₂ and pollutant emissions while protecting their efficiency. The mixture gradients change the fundamental properties of the flame, especially by a difference in temperature and composition between the burnt gases and those of a flame consuming a homogeneous mixture. This paper presents an investigation of the properties of the flame propagating in a lean homogeneous mixture after ignition in a richer mixture according to the magnitude of the stratification. Three magnitudes of stratification are investigated. The local flame burning velocity is determined by an original PIV algorithm developed previously. The local equivalence ratio in the fresh gases is measured from anisole PLIF. From the simultaneous PIV–PLIF measurements, the flame burning velocities conditioned on the local stretch rate and equivalence ratio in fresh gases are measured. The flame propagating through the homogeneous lean mixture has properties depending on the ignition conditions in the stratified layer. The flame propagating in the lean mixture is back-supported longer for ignition under the richer condition. The change of stretch sensitivity and burning velocity of the flame in the lean mixture is measured over time for the three magnitudes of mixture stratification investigated. The ignition in richer mixtures compensates for the nonequidiffusion effect of lean propane flame and sustains its robustness to stretch. The flame propagation in the lean homogeneous mixture is enhanced by ignition in a richer stratified layer, as much by their robustness to stretch as by an increase in the flame speed or the burning velocity. The decay time of this influence of the stratification, called memory effect, is determined.     

Influence of H2 on the response of lean premixed CH4 flames to high strained flows

A combined experimental and numerical investigation on the effects of H₂ addition to lean-premixed CH₄ flames in highly strained counterflow fields (with strain rates up to 8000 s⁻¹) using preheated flows indicate significant enhancement of lean flammability limits and extinction strain rates for relatively small amounts of H₂ addition. Numerical modeling of the counterflow opposed jet configuration used in this study indicated extinction strain rates which were within 5% of experimentally measured values for equivalence ratios ranging from 0.75 to less than 0.4. Both experimental and numerical results indicate that increasing H₂ in the fuel significantly increases flame speeds and thus extinction strain rates. Furthermore, increasing H₂ decreases the dependency of extinction equivalence ratio on the strain rate of the flow. For all of the mixtures investigated, extinction temperatures depend primarily on equivalence ratio and not fuel composition for the range of H₂ content studied, which suggests that extinction can be correlated to flame temperature and O₂ concentration. Nonetheless, H₂ addition greatly increases the maximum allowable strain rate before extinction temperatures are reached. Inspection of the model-predicted species profiles suggest that the enhancement of CH₄ burning rates with H₂ addition is driven by early H₂ breakdown increasing radical production rates early in the flame zone to enhance CH₄ ignition under conditions where otherwise CH₄ combustion might be prone to undergo extinction.     

Experimental and modeling studies of C2H4/O2/Ar, C2H4/methylal/O2/Ar and C2H4/ethylal/O2/Ar rich flames and the effect of oxygenated additives

The addition of dimethoxymethane (DMM or methylal) and diethoxymethane (DEM or ethylal) to a rich ethylene/oxygen/argon flame has been investigated by measuring the depletion of soot precursors. Three rich premixed ethylene/oxygen/argon (with and without added methylal or ethylal) flat flames have been stabilized at low-pressure (50 mbar) on a Spalding–Botha type burner with the same equivalence ratio of 2.50. Identification and monitoring of signal intensity profiles of species within the flames have been carried out by using molecular beam mass spectrometry (M.B.M.S.). The replacement of some C₂H₄ by C₃H₈O₂ or C₅H₁₂O₂ is responsible for a decrease of the maximum mole fractions of the detected intermediate species. This phenomenon is noticeable for C₂–C₄ intermediates and becomes more effective for C₅–C₁₀ species, mainly when C₃H₈O₂ added.A new kinetic model has been elaborated and contains 546 reactions and 107 chemical species in order to simulate the three investigated flames: C₂H₄/O₂/Ar, C₂H₄/DMM/O₂/Ar and C₂H₄/DEM/O₂/Ar. The reaction mechanism well reproduces experimental mole fraction profiles of major and intermediate species, and underlines the effect of methylal and ethylal addition on species concentration profiles for these flames.     

Flame front structure and burning velocity of turbulent premixed CH4/H2/air flames

Flame front structure of turbulent premixed CH₄/H₂/air flames at various hydrogen fractions was investigated with OH-PLIF technique. A nozzle-type burner was used to achieve the stabilized turbulent premixed flames. Hot-wire anemometer measurement and OH-PLIF observation were performed to measure the turbulent flow and detect the instantaneous flame front structure, respectively. The hydrogen fractions of 0%, 5%, 10% and 20% were studied. Results show that the flame front structures of the turbulent premixed flames are wrinkled flame front with small scale convex and concave structures compared to that of the laminar-flame front. The wrinkle intensity of flame front is promoted with the increase of turbulence intensity as well as hydrogen fraction. Hydrogen addition promotes the flame intrinsic instability which leads to the active response of laminar flame to turbulence and results in the much more wrinkled flame front structure. The value of S_T/S_L increases monotonically with the increase of u′/S_L and hydrogen fraction. The increase of S_T/S_L with the increase of hydrogen fraction is mainly attributed to the diffusive-thermal instability effects represented by the effective Lewis number, Le_eff. A general correlation between S_T/S_L and u′/S_L is provided from the experimental data fitting in the form of S_T/S_L ∝ a(u′/S_L)ⁿ, and the exponent, n, gives the constant value of 0.35 for all conditions and at various hydrogen fractions.     

Laser-induced fluorescence measurements of NCN in low-pressure CH4/O2/N2 flames and its role in prompt NO formation

NCN profiles were measured for five rich and lean premixed, low-pressure methane flames using laser-induced fluorescence (LIF). A semiquantitative determination of the NCN mole fractions as a function of spatial height above the burner is made by calibrating the NCN LIF signals using highly accurate OH LIF measurements in an adjacent spectral region. The resulting calibration yields an uncertainty estimate of a factor of 3 for the absolute values, but only ±25% for the relative NCN profiles. For all flame conditions, the NCN profiles occur immediately downstream of previously measured CH profiles. In addition, high correlations are found between the peak CH and peak NCN concentrations and the peak NCN and postflame NO concentrations over all equivalence ratios. These observations are consistent with NCN being the primary product channel from the CH + N₂ reaction and the initial intermediate in the prompt NO formation. This is the first mechanistic study in hydrocarbon flames that provides such experimental evidence. The experimental profiles are compared to numerical calculations using modified versions of two well-established hydrocarbon kinetic mechanisms. Reasonable agreement between the calculations and experiment is found for NCN profile shape, location of peak NCN concentrations, and absolute mole fractions. However, the dependence on stoichiometry of the peak NCN concentration is overestimated. Further work is required on NCN kinetics for modeling prompt NO in laminar premixed flames.     

Extinction studies of near-limit lean premixed syngas/air flames

Extinction studies of weakly-stretched near-limit lean premixed syngas/air flames were conducted in a twin-flame counterflow configuration. Experiments showed that buoyancy-induced natural convection at normal gravity strongly disturbed these flames. In order to validate the simulation, accurate extinction data was obtained at micro-gravity. Experimental data obtained from the 3.6 s micro-gravity drop tower showed that the extinction equivalence ratio increased with the increasing global stretch rate and decreased with the increasing H₂ mole fraction in the fuel. Numerical simulation was conducted with CHEMKIN software using GRI 3.0 and USC-Mech II mechanisms. The predicted extinction limit trend was in agreement with the micro-gravity experimental data. Sensitivity analyses showed that the competition between the main branching reaction H + O₂ ⇔ O + OH and the main termination reaction H + O₂ + M ⇔ HO₂ + M in the H₂/O₂ chemistry determined the extinction limits of the flames. The dominant species for syngas/air flame extinction was the H radical. The key exothermal reaction changed from OH + CO ⇔ H + CO₂ to OH + H₂ ⇔ H + H₂O with the increasing H₂ mole fraction in the fuel. Also, the mass diffusion played a more important role than chemical kinetics in the flame extinction. When the H₂ mass diffusion was suppressed, the reaction zone was pushed to the stagnation plane and the flame became weaker; while H mass diffusion is suppressed, the reaction zone slightly shifted towards the upstream and the flame was slightly strengthened.     

Kinetics and modeling study of a Mg(BH4)2/Ca(BH4)2 destabilized system

Borohydrides such as Mg(BH₄)₂ and Ca(BH₄)₂ continue to attract attention as potential hydrogen storage materials because of their high hydrogen content. In this study the desorption kinetics of Mg(BH₄)₂, Ca(BH₄)₂ and a 5Mg(BH₄)₂/Ca(BH₄)₂ mixture of the two were compared in the two-phase region at the same temperature and at a constant pressure thermodynamic driving force. The rate of hydrogen desorption from the two-phase region in the 5Mg(BH₄)₂/Ca(BH₄)₂ mixture was faster than that from either of the constituents. This indicated that Ca(BH₄)₂ was able to serve as a destabilizing agent for Mg(BH₄)₂. The results also showed that the desorption rates from the two-phase region were much faster than those from the single phase region. Modeling studies showed that the rate of hydrogen release from Mg(BH₄)₂, during the first 80% of the reaction, is diffusion controlled while in Ca(BH₄)₂ the reaction rate is phase boundary controlled. In the mixture the rate appears to be under the mixed control of both processes.     

Evolution of soot size distribution in premixed ethylene/air and ethylene/benzene/air flames: Experimental and modeling study

The effect of benzene concentration in the initial fuel on the evolution of soot size distribution in ethylene/air and ethylene/benzene/air flat flames was characterized by experimental measurements and model predictions of size and number concentration within the flames. Experimentally, a scanning mobility particle sizer was used to allow spatially resolved and online measurements of particle concentration and sizes in the nanometer-size range. The model couples a detailed kinetic scheme with a discrete-sectional approach to follow the transition from gas-phase to nascent particles and their coagulation to larger soot particles. The evolution of soot size distribution (experimental and modeled) in pure ethylene and ethylene flames doped with benzene showed a typical nucleation-sized (since particles do not actually nucleate in the classical sense particle inception is often used in place of nucleation) mode close to the burner surface, and a bimodal behavior at greater height above burner (HAB). However, major features were distinguished between the data sets. The growth of nucleation and agglomeration-sized particles was faster for ethylene/benzene/air flames, evidenced by the earlier presence of bimodality in these flames. The most significant changes in size distribution were attributed to an increase in benzene concentration in the initial fuel. However, these changes were more evident for high temperature flames. In agreement with the experimental data, the model also predicted the decrease of nucleation-sized particles in the postflame region for ethylene flames doped with benzene. This behavior was associated with the decrease of soot precursors after the main oxidation zone of the flames.