Inhibition of premixed and nonpremixed flames with phosphorus-containing compounds

The inhibition/extinction of various flames—premixed stoichiometric C₃H₈/air, nonpremixed counterflow CH₄/O₂/N₂, and nonpremixed coflow n-heptane/air cup-burner flames doped with a number of phosphorus-containing compounds (PCCs)—has been investigated experimentally. More than 20 PCCs (organic phosphates, phosphonates, phosphates) and their fluorinated derivatives were studied. All PCCs exhibited similar dependencies in burning velocities, extinction strain rates, and extinction volume fractions of CO₂ upon PCC loading in the range of mole fractions of 0–7000 ppm within an experimental deviation of ± 5%. This confirms that the inhibition effectiveness of the PCCs is influenced by the phosphorus content in the PCC molecule rather than by the structure of the molecule. The burning velocity of a stoichiometric C₃H₈/air mixture and the extinction strain rate of a nonpremixed counterflow CH₄/O₂/N₂ flame doped with trimethylphosphate were calculated. Satisfactory agreement between experimental and modeling results confirms the conclusion that the reactions of phosphorus oxyacids with radicals are responsible for flame inhibition.     

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

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.     

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.     

Quantitative NH measurements by using laser-based diagnostics in low-pressure flames

NH is a key short-lived radical involved in the prompt-NO formation. Quantification of NH is thus particularly important for testing the NO kinetic mechanisms. However, quantitative measurements of native NH in hydrocarbon/oxygen/nitrogen flames remain very scarce. Therefore, in this work, the mole fractions of native NH were obtained using a combination of laser-based diagnostics; Laser Induced Fluorescence (LIF) and Cavity Ring-Down Spectroscopy (CRDS). The NH species was probed after exciting the transition R₁(6) in the A³Π-X³Σ^? (0-0) system at 333.9?nm. The mole fraction profiles of NH were successfully obtained in premixed low-pressure flames of CH₄/O₂/N₂ and C₂H₂/O₂/N₂ at two equivalence ratios of 1.00 and 1.25. The estimated detection limit for the NH radical was around 4.5?×?10⁸ molecule cm^?3 (i.e. 2 ppb in mole fraction at 1600?K), which is nearly 2 orders of magnitude lower than previous values reported in the literature. These new experimental results were compared with predictions by a recently developed NO model (namely NOMecha2.0). In the case of the CH₄ flames, a satisfying agreement between the experiment and model was observed. However, in the case of the C₂H₂ flames, some discrepancies were observed. Model analysis has highlighted the importance of the HCCO radicals in the NH formation through the HCNO→HNCO→NH₂ reactions pathway. Modification of the rate constant values of the reactions C₂H₂+?O and HCCO?+?O₂, which are key reactions for both the acetylene laminar flame speed and the HCCO predictions, has enabled the model to satisfactorily predict the experimental NH and NO profiles also in the C₂H₂ flames.     

Effects of Soret diffusion on the laminar flame speed and Markstein length of syngas/air mixtures

The effects of Soret diffusion on premixed syngas/air flames at normal and elevated temperatures and pressures are investigated numerically including detailed chemistry and transport. The emphasis is placed on assessing and interpreting the influence of Soret diffusion on the unstretched and stretched laminar flame speed and Markstein length of syngas/air mixtures. The laminar flame speed and Markstein length are obtained by simulating the unstretched planar flame and positively-stretched spherical flame, respectively. The results indicate that at atmospheric pressure the laminar flame speed of syngas/air is mainly reduced by Soret diffusion of H radical while the influence of H₂ Soret diffusion is negligible. This is due to the facts that the main reaction zone and the Soret diffusion for H radical (H₂) are strongly (weakly) coupled, and that Soret diffusion reduces the H concentration in the reaction zone. Because of the enhancement in the Soret diffusion flux of H radical, the influence of Soret diffusion on the laminar burning flux increases with the initial temperature and pressure. Unlike the results at atmospheric pressure, at elevated pressures the laminar flame speed is shown to be affected by the Soret diffusion of H₂ as well as H radical. For stretched spherical flame, it is shown that the Soret diffusion of both H and H₂ should be included so that the stretched flame speed can be accurately predicted. Similar to the laminar flame speed, the Markstein length is also reduced by Soret diffusion. However, the reduction is found to be mainly caused by Soret diffusion of H₂ rather than that of H radical. Moreover, the influence of Soret diffusion on the Markstein length is demonstrated to decrease with the initial temperature and pressure.     

Termolecular chemistry facilitated by radical-radical recombinations and its impact on flame speed predictions

Recent theoretical studies have shown that termolecular chemistry can be facilitated through reactions of flame radicals (H, O, and OH) or O₂ with highly-energized collision complexes (either radical or stable species) formed in exothermic reactions. In this work, radical-radical recombination reaction induced termolecular chemistry and its impact on combustion modeling was studied. Two recombination reactions, H + CH₃ + M → CH₄ + M and H + OH + M → H₂O + M, were analyzed using ab-initio master equation analyses guided by quasiclassical trajectory results. The dynamics results and the master equation calculations indicate that CH₄^? and H₂O^? (formed in the two radical-radical reactions outlined above) react rapidly with flame radicals and O₂ at rates that are competitive with collisional cooling. The addition of these processes into conventional combustion modeling requires two modifications: the inclusion of the new nonthermal termolecular reaction rates and the simultaneous reduction of the competing recombination reaction rates. The former is described with newly derived Arrhenius expressions based on quasiclassical trajectories, and the latter is achieved by perturbing the recombination reaction rate during the simulation. Kinetic modeling was used to gauge the impact of including this nonthermal chemistry for H₂/CH₄-air laminar flames speeds. Inclusion of this nonthermal chemistry has a noticeable impact on simulated flame speeds. The procedure developed here can be utilized to properly quantify the effects of such nonthermal reactions in macroscopic kinetic models.     

Inhibition of hydrogen–oxygen flames by iron pentacarbonyl at atmospheric pressure

The chemistry of inhibition of laminar premixed hydrogen–oxygen flames by iron pentacarbonyl at atmospheric pressure was studied experimentally and by numerical simulation. Flame speed and chemical structure were analyzed. Flame burning velocities and inhibition effectiveness were measured and simulated for various equivalence ratios. The concentration profiles of a number of Fe-containing products of Fe(CO)₅ combustion, including Fe, FeO₂, FeOH, and Fe(OH)₂, were first measured using probing molecular beam mass spectrometry in an atmospheric-pressure H₂/O₂/N₂ flame. A comparison of the experimental and modeling results shows that they are in satisfactory agreement with each other, indicating that the reaction mechanism proposed previously for flame inhibition by iron pentacarbonyl is adequate for predicting the chemical structure of flames. The key recombination stages of active species catalyzed by Fe-containing species for flames of various stoichiometries can be determined by calculations of the production rates of H and O atoms and OH radicals as well as by analysis of the kinetic model.     

Strain rate and fuel composition dependence of chemiluminescent species profiles in non-premixed counterflow flames: comparison with model results

A detailed comparison has been conducted between chemiluminescence (CL) species profiles of OH^?, CH^?, and C₂ ^?, obtained experimentally and from detailed flame kinetics modeling, respectively, of atmospheric pressure non-premixed flames formed in the forward stagnation region of a fuel flow ejected from a porous cylinder and an air counterflow. Both pure methane and mixtures of methane with hydrogen (between 10 and 30 % by volume) were used as fuels. By varying the air-flow velocities methane flames were operated at strain rates between 100 and 350 s^?1, while for methane/hydrogen flames the strain rate was fixed at 200 s^?1. Spatial profiles perpendicular to the flame front were extracted from spectrograms recorded with a spectrometer/CCD camera system and evaluating each spectral band individually. Flame kinetics modeling was accomplished with an in-house chemical mechanism including C₁–C₄ chemistry, as well as elementary steps for the formation, removal, and electronic quenching of all measured active species. In the CH₄/air flames, experiments and model results agree with respect to trends in profile peak intensity and position. For the CH₄/H₂/air flames, with increasing H₂ content in the fuel the experimental CL peak intensities decrease slightly and their peak positions shift towards the fuel side, while for the model the drop in mole fraction is much stronger and the peak positions move closer to the fuel side. For both fuel compositions the modeled profiles peak closer to the fuel side than in the experiments. The discrepancies can only partly be attributed to the limited attainable spatial resolution but may also necessitate revised reaction mechanisms for predicting CL species in this type of flame.     

An experimental study on the formation of polycyclic aromatic hydrocarbons in laminar coflow non-premixed methane/air flames doped with four isomeric butanols

Experimental measurements were conducted for temperatures and mole fractions of C1–C16 combustion intermediates in laminar coflow non-premixed methane/air flames doped with 3.9% (in volume) 1-butanol, 2-butanol, iso-butanol and tert-butanol, respectively. Synchrotron vacuum ultraviolet photoionization mass spectrometry (SVUV-PIMS) technique was utilized in the measurements of species mole fractions. The results show that the variant molecular structures of butyl alcohols have led to different efficiencies in the formation of polycyclic aromatic hydrocarbons (PAHs) that may cause the variations in sooting tendency. Detailed species information suggests that the presence of allene and propyne promotes benzene formation through the C₃H₃ + C₃H₄ reactions and consequently PAH formation through the additions of C2 and C3 species to benzyl or phenyl radicals. As a matter of fact, PAHs formed from the 1-butanol doped flame are the lowest among the four investigated flames, because 1-butanol mainly decomposes to ethylene and oxygenates rather than C3 hydrocarbon species. Meanwhile, the tert-butanol doped flame generates the largest quantities of allene and propyne among the four flames and therefore is the sootiest one.     

Laminar burning velocity and ignition delay time for premixed isooctane–air flames with syngas addition

The effects of blending syngas in different proportions to isooctane on the laminar burning velocity and ignition delay time of the fuel–air mixture have been studied in SI engine relevant conditions. The syngas is assumed to be composed of 50% H₂ and 50% CO. Simulations have been carried out using a skeletal mechanism containing 143 species and 643 reaction steps. It has been found that the blending of syngas augments the laminar burning velocity of isooctane due to increase of the thermal diffusivity of the reactant mixture and alteration in the chemistry of the flame reactions. For the mixture of 30% isooctane/70% syngas, the laminar burning velocity and the ignition delay time values are very close to those corresponding to pure isooctane. Additionally, the effects of exhaust gas recirculation have been explored for the 30% isooctane/70% syngas–air flame. It is seen that the reduction in laminar burning velocity due to the dilution by the recirculated exhaust gas can be compensated by an increase in the unburnt gas temperature. The effect of the exhaust gas dilution on the ignition delay time of 30% isooctane/70% syngas–air mixture has been found to be negligible.     

Numerical investigation of flame–vortex interactions in laminar cross-flow non-premixed flames in the presence of bluff bodies

Flame stabilisation in a combustor having vortices generated by flame holding devices constitutes an interesting fundamental problem. The presence of vortices in many practical combustors ranging from industrial burners to high speed propulsion systems induces vortex–flame interactions and complex stabilisation conditions. The scenario becomes more complex if the flame sustains after separating itself from the flame holder. In a recent study [P.K. Shijin, S.S. Sundaram, V. Raghavan, and V. Babu, Numerical investigation of laminar cross-flow non-premixed flames in the presence of a bluff-body, Combust. Theory Model. 18, 2014, pp. 692–710], the authors reported details of the regimes of flame stabilisation of non-premixed laminar flames established in a cross-flow combustor in the presence of a square cylinder. In that, the separated flame has been shown to be three dimensional and highly unsteady. Such separated flames are investigated further in the present study. Flame–vortex interactions in separated methane–air cross flow flames established behind three bluff bodies, namely a square cylinder, an isosceles triangular cylinder and a half V-gutter, have been analysed in detail. The mixing process in the reactive flow has been explained using streamlines of species velocities of CH₄ and O₂. The time histories of z-vorticity, net heat release rate and temperature are analysed to reveal the close relationship between z-vorticity and net heat release rate spectra. Two distinct fluctuating layers are visible in the proper orthogonal decomposition and discrete Fourier transform of OH mass fraction data. The upper fluctuating layer observed in the OH field correlates well with that of temperature. A detailed investigation of the characteristics of OH transport has also been carried out to show the interactions between factors affecting fluid dynamics and chemical kinetics that cause multiple fluctuating layers in the OH.     

Prediction of NO in turbulent diffusion flames using Eulerian particle flamelet model

The combustion characteristics for the turbulent diffusion flames using the unsteady flamelet concept have been numerically investigated. The Favre-averaged Navier–Stokes equations are solved by a finite volume method of SIMPLE type that incorporates the laminar flamelet concept with a modified k ? ε turbulence model. The NO formation is estimated by solving the Eulerian particle transport equations in a postprocessing mode. Two test problems are considered: CH₄/H₂/N₂ jet flame and CH₄/H₂ stabilised bluff body flame. The temperature and species profiles are well captured by the flamelet model. Two different chemical mechanisms (GRI 2.11 and 3.0) give nearly identical results for temperature and species except NO. The GRI 3.0 gives significantly higher NO levels compared to the GRI 2.11. This is mainly attributed to the difference in NO formation by the prompt mechanism. The NO formation is sensitive to the number of flamelet particles. The NO levels for two test flames do not change when the flamelet particle number exceeds six.     

Further studies of the reaction kernel structure and stabilization of jet diffusion flames

The structure of axisymmetric laminar jet diffusion flames of ethane, ethylene, acetylene, and propane in quasi-quiescent air has been studied numerically in normal earth gravity (1g) and zero gravity (0g). The time-dependent full Navier–Stokes equations with buoyancy were solved using an implicit, third-order accurate numerical scheme, including a C₃-chemistry model and an optically thin-media radiation model for heat losses. Observations of the flames were also made at the NASA Glenn 2.2-Second Drop Tower. For all cases of the fuels and gravity levels investigated, a peak reactivity spot, i.e., reaction kernel, was formed in the flame base, thereby holding a trailing diffusion flame. The location of the reaction kernel with respect to the burner rim depended inversely on the reaction-kernel reactivity or velocity. In the C₂ and C₃ hydrocarbon flames, the H₂–O₂ chain reactions were important at the reaction kernel, yet the CH₃ + O → CH₂O + H reaction, a dominant contributor to the heat-release rate in methane flames studied previously, did not outweigh other exothermic reactions. Instead of the C₁-route oxidation pathway in methane flames, the C₂ and C₃ hydrocarbon fuels dehydrogenated on the fuel side and acetylene was a major hydrocarbon fragment burning at the reaction kernel. The reaction-kernel correlations between the reactivity (the heat-release or oxygen-consumption rate) and the velocity, obtained previously for methane, were developed further for various fuels in more universal forms using variables related to local Damköhler numbers and Peclet numbers.     

Multi-scalar measurements in a premixed swirl burner using 1D Raman/Rayleigh scattering

Instantaneous measurements of temperature, equivalence ratio, and major species were performed along a one-dimensional probe volume using simultaneous Raman/Rayleigh scattering in an unconfined turbulent lean-premixed swirling methane/air flame. Temperature was determined from Rayleigh scattering and the major species, CO₂, O₂, N₂, CH₄, H₂O, and H₂ from Raman scattering. Effective Rayleigh cross-sections were corrected using the local chemical composition obtained from Raman scattering. These experiments were conducted to investigate the compositional structure of a lean-premixed swirling flame in detail and to complement previous measurements of the underlying flow field. The flame was classified within a revised regime diagram at the cross-over between corrugated flames and thin reaction zones. Instantaneous temperature profiles varied significantly showing shapes ranging from laminar-like flamelets to mixing between reacted fluid elements and secondary air. Different thermo-kinetic states could be assigned to the inner and outer recirculation zones and to the inner and outer mixing layers. Linked to published velocity data of this flame, the present multi-scalar data are useful for validation of numerical simulations.     

Soot emission reduction in oxygenated co-flow jet flames

A computational study was performed for ethylene/air non-premixed laminar co-flow jet flames using an axisymmetric CFD code to explore the effect of oxygenation on PAH and soot emissions. Oxygenated flames were established using N₂ diluted fuel stream along with O₂ enriched air stream such that the stoichiometric mixture fraction (Ζ_st) is varied but the adiabatic flame temperature is not materially changed. Simulations were carried out using a spatially and temporally accurate algorithm with detailed chemistry and transport. A detailed kinetic model involving 111 species and 784 reactions and a fairly detailed soot model were incorporated into the code. Two different approaches, one with constant flame height and other with constant inlet velocity are comprehensively examined to bring out the effects of changes in flame structure and residence time on soot emissions with respect to Z_st. With increase in Ζ_st, a drastic reduction in the formation of soot precursors (acetylene and benzene) and thus in soot emissions are observed. In the present study, oxygenated flames with Ζ_st ≥ 0.424 are considered as blue flames or completely soot free. For various oxygenated flames a C/O ratio between 0.45 and 0.6 is found to be most favorable for soot formation.     

Nitrous oxide conversion in laminar premixed flames of CH4 + O2 + Ar

Experimental measurements of the adiabatic burning velocity in neat and NO formation in CH₄ + O₂ + Ar flames doped with small amounts of N₂O are presented. The oxygen content in the oxidizer was varied from 15 to 17%. Non-stretched flames were stabilized on a perforated plate burner at 1 atm. The Heat Flux method was used to determine burning velocities under conditions when the net heat loss of the flame is zero. Adiabatic burning velocities of methane + oxygen + argon mixtures were found in satisfactory agreement with the modeling. The NO concentrations in the flames doped with N₂O (100 ppm in the argon stream before mixing) were measured in the burnt gases at a fixed distance from the burner using probe sampling. Axial profiles of [NO] were found insensitive to the downstream heat losses. Experimental dependencies of [NO] versus equivalence ratio had a maximum between φ = 1.1 and 1.2. Calculated concentrations of NO were in good agreement with the measurements. In lean flames calculated concentrations of NO strongly depends on the rate constant of reaction N₂O + O=NO + NO if too high values proposed in the literature are employed. These new experimental data thus allowed for validation of the key reactions of the nitrous oxide mechanism of NO formation in flames.     

Experimental and kinetic modeling study of the positive ions in premixed ethylene flames over a range of equivalence ratios

Understanding the ion chemistry in flames is crucial for developing ion sensitive technologies for controlling combustion processes. In this work, we measured the spatial distributions of positive ions in atmospheric-pressure burner-stabilized premixed flames of ethylene/oxygen/argon mixtures in a wide range of equivalence ratios ϕ = 0.4÷1.5. A flame sampling molecular beam system coupled with a quadrupole mass spectrometer was used to obtain the spatial distributions of cations in the flames, and a high mass resolution time-of-flight mass spectrometer was utilized for the identification of the cations having similar m/z ratios. The measured profiles of the flame ions were corrected for the contribution of hydrates formed during sampling in the flames slightly upstream the flame reaction zone. We also proposed an updated ion chemistry model and verified it against the experimental profiles of the most abundant cations in the flames. Our model is based on the kinetic mechanism available in the literature extended with the reactions for C₃H₅⁺ cation. Highly accurate W2-F12 quantum chemical calculations were used to obtain a reliable formation enthalpy of C₃H₅⁺. The model was found to reproduce properly the measured relative abundance of the key oxygenated cations (viz., CH₅O⁺, C₂H₃O⁺) in the whole range of equivalence ratios employed, and the C₃H₅⁺ cation abundance in the richest flame with ϕ=1.5, but significantly underpredicts the relative mole fraction of C₃H₃⁺, which becomes a key species under fuel-rich conditions. Apart from this, several aromatic and cyclic C_xH_y cations dominating under fuel-rich conditions were identified. We also considered the most important directions for the further refinement of the mechanism.     

Elemental surface enrichment in flame synthesis: A mechanism associated with particle melting-solidification

Elemental surface enrichment is important for functionalities of flame-synthesized particle materials, but its mechanism is poorly understood. In this paper, a mechanism associated with particle melting-solidification is proposed based on an experimental study. Y₂O₃:Eu particles were generated by flame assisted spray pyrolysis (FASP), using H₂/O₂ flames or H₂/air flames. The particles were analyzed by transmission electron microscopy (TEM), X-ray diffraction (XRD), inductively coupled plasma mass spectrometry (ICP-MS), and X-ray photoelectron spectroscopy (XPS). H₂/O₂ flames resulted in particles with Eu surface enrichment, i.e. the surface Eu concentration was several times higher than the overall Eu concentration; there was no elemental surface enrichment in particles from H₂/air flames. The Eu surface enrichment in H₂/O₂ flames was attributed to elemental partitioning during solidification of molten Y₂O₃:Eu particles; in H₂/air flames the particles did not melt and hence there was no elemental surface enrichment. The findings of this study suggest that elemental surface enrichment may be a common phenomenon for binary metal oxide particles that experience melting-solidification. Such particles should be examined for elemental surface enrichment, both for understanding their functionality and for their potential biological effects.