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.     

An experimental and kinetic modeling study on nitric oxide formation in premixed C3 alcohols flames

This study provides new quantitative NO concentrations measurements in n-propanol + air and i-propanol + air flames together with a new combustion kinetic model. The heat flux method was employed to stabilize propyl alcohols flames and the initial gas conditions were set to 323 K, 1 atm, and Φ=0.7–1.4. Saturated laser-induced fluorescence was employed to measure NO concentration in the post-combustion region. The presented and literature models, namely the POLIMI and Bohon et al. (2018) kinetic mechanisms, were assessed against new experimental data. Experimental results showed a higher NO formation in the thermal zone for n-propanol flames, whereas i-propanol flames indicate a higher amount of NO formed at fuel-rich conditions. Overall among the tested models, the present mechanism exhibited the best agreement in emulating NO experimental profiles; conversely, numerical simulations from the POLIMI model showed significant inconsistencies at fuel-rich conditions and the Bohon et al. (2018) model was unable to reproduce the measured data, notably underpredicting experimental values at all investigated conditions. However, the present model manifested some uncertainties in reproducing NO formation in the prompt region; therefore, in connection with this important aspect, the new experimental data obtained in this work will provide a valid support to further develop more reliable kinetic models.     

Measurements of burning velocities of dimethyl ether and air premixed flames at elevated pressures

Laminar burning velocities of dimethyl ether (DME) and air premixed flames at elevated pressures up to 10 atm were measured by using a newly developed pressure-release type spherical bomb. The measurement system was validated using laminar burning velocities of methane–air flames. A comparison with the previous experimental data shows an excellent agreement and demonstrates the accuracy and reliability of the present experimental system. The measured flame speeds of DME–air flames were compared with the previous experimental data and the predictions using the full and reduced mechanisms. At atmospheric pressure, the measured laminar burning velocities of DME–air flames are in reasonable agreement with the previous data from spherical bomb method, but are much lower than both predictions and the experimental data of the PIV based counterflow flame measurements. The laminar burning velocities of DME–air flames at 2, 6, and 10 atm were also measured. It was found that flame speed decreases considerably with the increase of pressure. Moreover, the measured flame speeds are also lower than the predictions at high pressures. In addition, experiments showed that at high pressures the rich DME–air flames are strongly affected by the hydrodynamic and thermal-diffusive instabilities. Markstein lengths and the overall reaction order at different equivalence ratios were extracted from the flame speed data at elevated pressures. Sensitivity analysis showed that reactions involving methyl and formyl radicals play an important role in DME–air flame propagation and suggested that systematic modification of the reactions rates associated with methyl and formyl formations are necessary to reduce the discrepancies between predictions and measurements.     

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.     

Experimental and numerical study of premixed, lean ethylene flames

Ethylene is a key intermediate in the combustion mechanisms of most practical fuels. It plays also an important role in the formation of aromatic hydrocarbons and soot particules. The latter has motivated many experimental and numerical studies carried out on rich ethylene-air mixtures. Less studies have been devoted to lean mixtures, and the development of strategies based on lean, premixed flames to reduce soot and NO_x production requires additional experimental data in lean conditions. In this work, the chemical structure of lean premixed ethylene-oxygen-nitrogen flames stabilized on a flat-flame burner at atmospheric pressure was determined experimentally. The species mole fraction profiles were also computed by the Premix code (Chemkin II version) and four detailed reaction mechanisms. A very good agreement was observed for the main flame properties: reactants consumption, final products (CO₂, H₂O) and the main intermediates: CO and H₂. Marked differences occurred in the prediction of active intermediate species present in small concentrations. Pathways analyses were performed to identify the origins of these discrepancies. It was shown that the same reactions were involved in the four mechanisms to describe the consumption of ethylene, but with marked differences in their relative importance. C₂H₃ and CH₂HCO are the main radicals formed in this first step and their consumption increases the differences between the mechanisms either by the use of different kinetic data for common reactions or by differences in the nature of the consumption reactions.     

Experimental and numerical study of product gas characteristics of ammonia/air premixed laminar flames stabilized in a stagnation flow

Ammonia has widely attracted interest as a potential candidate not only as a hydrogen energy carrier but also as a carbon free fuel for internal combustion engines, such as gas turbines. Because ammonia contains a nitrogen atom in its molecule, nitrogen oxides (NO_x) and other pollutants may be formed when it burns. Therefore, understanding the fundamental product gas characteristics of ammonia/air laminar flames is important for the design of ammonia-fueled combustors to meet stringent emission regulations. In this study, the product gas characteristics of ammonia/air premixed laminar flames for various equivalence ratios were experimentally and numerically investigated up to elevated pressure conditions. In the experiments, a stagnation flame configuration was employed because an ammonia flame can be stabilized by using such a configuration without a pilot flame. The experimental results showed that the maximum NO mole fraction was about 3,500 ppmv, at an equivalence ratio of 0.9 at 0.1 MPa. The NO mole fraction decreased as the equivalence ratio increased. In addition, the maximum value of the NO mole fraction decreased with an increase in mixture pressure. Furthermore, it was experimentally clarified that the simultaneous reduction of NO and unburnt ammonia can be achieved at an equivalence ratio of about 1.06, which is the target equivalence ratio for emission control in rich-lean two-stage ammonia combustors. Comparison of experimental and numerical results showed that even though the reaction mechanisms employed have been optimized for predicting the laminar burning velocity of ammonia/air flames, they failed to satisfactorily predict the measured species in this study. Sensitivity analysis was used to identify elementary reactions that control the species profiles but have negligible effects on the burning velocity. It is considered that these reaction models need to be updated for accurate prediction of product gas characteristics of ammonia/air flames.     

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 the effects of gravity on soot formation in laminar coflow methane/air diffusion flames under different air stream velocities

Numerical simulations of laminar coflow methane/air diffusion flames at atmospheric pressure and different gravity levels were conducted to gain a better understanding of the effects of gravity on soot formation by using relatively detailed gas-phase chemistry and complex thermal and transport properties coupled with a semi-empirical two-equation soot model. Thermal radiation was calculated using the discrete-ordinates method coupled with a non-grey model for the radiative properties of CO, CO₂, H₂O, and soot. Calculations were conducted for three coflow air velocities of 77.6, 30, and 5 cm/s to investigate how the coflowing air velocity affects the flame structure and soot formation at different levels of gravity. The coflow air velocity has a rather significant effect on the streamwise velocity and the fluid parcel residence time, especially at reduced gravity levels. The flame height and the visible flame height in general increase with decreasing the gravity level. The peak flame temperature decreases with decreasing either the coflow air stream velocity or the gravity level. The peak soot volume fraction of the flame at microgravity can either be greater or less than that of its normal gravity counterpart, depending on the coflow air velocity. At sufficiently high coflow air velocity, the peak soot volume fraction increases with decreasing the gravity level. When the coflow air velocity is low enough, soot formation is greatly suppressed at microgravity and extinguishment occurs in the upper portion of the flame with soot emission from the tip of the flame owing to incomplete oxidation. The numerical results provide further insights into the intimate coupling between flame size, residence time, thermal radiation, and soot formation at reduced gravity level. The importance of thermal radiation heat transfer and coflow air velocity to the flame structure and soot formation at microgravity is demonstrated for the first time.     

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.     

Experimental and numerical determination of heat release in counterflow premixed laminar flames

This paper presents an experimental and numerical study of heat release in atmospheric laminar counterflow premixed flames. The measurements are based on simultaneous planar laser-induced fluorescence (PLIF) of OH and HCHO. These measurements are compared to numerical results obtained using detailed chemistry and multicomponent transport properties. A low Mach number formulation along the stagnation streamline is employed to describe the reactive flow. The conservation equations are completed with CHEMKIN and EGLIB packages. They are solved using finite differences, Newton iterations, and an adaptive gridding technique. The comparison is done along the burner axis for both, maximum heat release location and heat release profile width. It is shown that the product of OH and HCHO concentrations yields a result closely related to the heat release. These comparisons lead to the conclusion that the experimental method used seems to be a good tool for the determination of heat release in flames.     

Temperature and CH* measurements and simulations of laminar premixed ethylene jet-wall stagnation flames

New experimental 2D measurements are reported to characterise the flame location, shape and temperature of laminar premixed ethylene jet-wall stagnation flames when the equivalence ratio, exit gas velocity and burner-plate separation distance are varied. Bandpass-filtered optical measurements of the CH* chemiluminescence were used to provide information about the shape and location of the flames. Thin filament pyrometry (TFP) using a 14 µm diameter SiC filament was used to make line measurements of the temperature to reconstruct the full 2D temperature field for the first time in premixed, jet-wall stagnation flames. The comparison of CH* measurements with (intrusive) and without (non-intrusive) the presence of the SiC filament showed that the filament resulted in minimal disturbance of the flame when the filament was placed downstream of the flame front. However, the flame was observed to attach to the filament, resulting in more significant disturbance, when it was placed upstream of the flame front. The flames were simulated using both 1D and 2D models. The 2D simulations were used to provide estimates of the velocity, kinematic viscosity and thermal conductivity required to calculate the gas temperature from the TFP data. The 1D simulations showed excellent agreement with the experimentally observed centreline quantities, but required the strain boundary condition to be fitted in order to match the experimentally observed flame location. The 2D simulations showed excellent agreement without the need for any fitting, and correctly predicted the flame shape, location and temperature as the experimental conditions were varied. A comparison of the set of simulated temperature-residence times along different streamlines showed relatively uniform distributions within each flame. However, the most uniform set of temperature-residence time distributions did not correlate with the flattest flame.     

The effects of strain and curvature on the mass burning rate of premixed laminar flames

The Markstein number characterizes the effect that flame stretch has on the burning velocity. Different expressions for this number are deduced from integral analysis. According to a phenomenological law, the Markstein number can be separated into a part for the curvature of the flame and a part for the straining of the flow. This separation is analysed here. It appears that the Markstein number for curvature and the combined one for both curvature and strain are unique. It is, however, not possible to introduce a separate and unique Markstein number for the flow straining that can be used to describe its influence in different combustion situations. The theoretical and numerical analysis is applied to flat steady counterflow flames as well as to steady, imploding and expanding spherical flames.     

Experimental and numerical study of product gas and N2O emission characteristics of ammonia/hydrogen/air premixed laminar flames stabilized in a stagnation flow

In order to achieve carbon neutrality, the use of ammonia as a fuel for power generation is highly anticipated. The utilization of a binary fuel consisting of ammonia and hydrogen can address the weak flame characteristics of ammonia. In this study, the product gas characteristics of ammonia/hydrogen/air premixed laminar flames stabilized in a stagnation flow were experimentally and numerically investigated for various equivalence ratios for the first time. A trade-off relationship between NO and unburnt ammonia was observed at slightly rich conditions. At lean conditions, NO reached a maximum value of 8,700 ppm, which was larger than that of pure ammonia/air flames. The mole fraction of nitrous oxide (N₂O) which has large global warming potential rapidly increased around the equivalence ratio of 0.6, which was attributed to the effect of a decrease in flame temperature downstream of the reaction zone owing to heat loss to the stagnation wall. To understand this effect further, numerical simulations of ammonia/hydrogen/air flames were conducted using the stagnation flame model for various equivalence ratios and stagnation wall temperatures. The results show that the important reactions for N₂O production and reductions are NH +NO = N₂O + H, N₂O + H = N₂ + OH, and N₂O (+M) = N₂ + O (+M). A decrease in flame temperature in the post flame region inhibited N₂O reduction through N₂O (+M) = N₂ + O (+M) because this reaction has a large temperature dependence, and thus N₂O was detected as a product gas. N₂O is reduced through N₂O (+M) = N₂ + O (+M) in the post flame region if the stagnation wall temperature is sufficiently high. On the other hand, it was clarified that an increase in equivalence ratio enhances H radical production and promotes N₂O reduction by H radical through the reaction of N₂O + H = N₂ + OH.     

Measurement of unstable burning velocities of iso-octane–air mixtures at high pressure and the derivation of laminar burning velocities

A new technique is reported for measuring burning velocities at high pressures in the final stages of two inwardly propagating flame kernels in an explosion bomb. The flames were initiated at diametrically opposite spark electrodes, close to the wall, in quiescent mixtures. Measurements of pressure and flame kernel propagation speeds by high-speed photography showed the burning velocities to be elevated above the corresponding laminar burning velocities as a result of the developing flame instabilities. The enhancement increased with increase in pressure and decreased with increase in Markstein number. When the Markstein number was negative, instabilities could be appreciable, as could the enhancement. For the iso-octane–air mixtures investigated, where the mixtures had well-characterised Markstein numbers or critical Peclet numbers at the relevant pressures and temperatures, it was possible to explain the enhancement quantitatively by the spherical explosion flame instability theory of Bechtold and Matalon, provided the critical Peclet number was that observed experimentally, and allowance was made for the changing pressure. With this theoretical procedure, it was possible to derive values of laminar burning velocity from the measured values of burning velocity over a wide range of equivalence ratios, pressures, and temperatures. The values became less reliable at the higher temperatures and pressures as the data on Markstein and critical Peclet numbers became less certain. It was found that with iso-octane as the fuel the laminar burning velocity decreased during isentropic compression.     

Numerical simulations of soot aggregation in premixed laminar flames

In this paper we make use of a detailed particle model and stochastic numerical methods to simulate the particle size distributions of soot particles formed in laminar premixed flames. The model is able to capture the evolution of mass and surface area along with the full structural detail of the particles. The model is validated against previous models for consistency and then used to simulate flames with bimodal and unimodal soot particle distributions. The change in morphology between the particles from these two types of flames provides further evidence of the interplay among nucleation, coagulation, and surface rates. The results confirm the previously proposed role of the strength of the particle nucleation source in defining the instant of transition from coalescent to fractal growth of soot particles.