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.     

An experimental and kinetic modeling study of premixed nitroethane flames at low pressure

An experimental and kinetic modeling study is reported on three premixed nitroethane/oxygen/argon flames at low pressure (4.655 kPa) with the equivalence ratios (Φ) of 1.0, 1.5 and 2.0. Over 30 flame species were identified with tunable synchrotron vacuum ultraviolet photoionization mass spectrometry, with their mole fractions quantified as the function of the height above burner. The flame temperature profiles were measured with a Pt–6%Rh/Pt–30%Rh thermocouple. A detailed kinetic mechanism with 115 species and 730 reactions was proposed and validated against experimental results. The computed predictions have shown satisfactory agreement with the experimental results. Basing on the rate-of-production analysis, the reaction pathways that feature the combustion of nitroethane were revealed, including the primary decomposition of C–N bond fission, the oxidation of C2 and C1 hydrocarbons and the formation of nitrogenous species. The presence of NO₂ and NO has been proved to be important for these processes.     

Experimental and semi-detailed kinetic modeling study of decalin oxidation and pyrolysis over a wide range of conditions

Decalin is the simplest polycyclic alkane (polynaphtenic hydrocarbon) found in liquid fuels (jet fuels, Diesel). In order to better understand the combustion characteristics of decalin, this study provides new experimental data for its oxidation in a jet-stirred reactor. For the first time, stable species concentration profiles were measured in a jet-stirred reactor at a constant mean residence time of 0.1 s and 0.5 s at respectively 1 and 10 atm, over a range of equivalence ratios (? = 0.5–1.5) and temperatures (750–1350 K). The oxidation of decalin under these experimental conditions was modeled using a semi-detailed chemical kinetic reaction mechanism (11,000 reactions involving 360 species) derived from a previously proposed scheme for the ignition of the same fuel in a shock-tube. The proposed mechanism that includes both low- and high-temperature chemistry shows reasonably good agreement with the present experimental data set. It can also represent well decalin pyrolysis and oxidation data available in the literature. Reaction path analyses and sensitivity analyses were conducted to interpret the results.     

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.     

Modeling and measurements of size distributions in premixed ethylene and benzene flames

This study demonstrates the major differences in the evolution of the particle size distributions (PSDs), both measured and modeled, of soot in premixed benzene and ethylene flat flames. In the experiments, soot concentration and PSDs were measured by using a scanning mobility particle sizer (SMPS, over the size range of 3-80 nm). The model employed calculations of gas phase species coupled with a discrete sectional approach for the gas-to-particle conversion. The model includes reaction pathways leading to the formation of nano-sized particles and their coagulation to larger soot particles. The particle size distribution, both experimental and modeled, evolved from a single particle mode (the nucleation mode) to a bimodal size distribution. An important distinction between the results for the ethylene and benzene flames is the behavior of the nucleation mode which persists at all heights above the burner (HAB) for ethylene whereas it was greatly suppressed at greater HAB for the benzene flames. The explanation for the decreased nucleation mode at higher elevations in the benzene flame is that the aromatics are consumed in the oxidation zone of the flame. Fair predictions of particle-phase concentrations and particle sizes in the two flames were obtained with no adjustments to the kinetic scheme. In agreement with experimental data, the model predicts a higher formation of particulate in the benzene flame as compared with the ethylene flame.     

An experimental and kinetic modeling study of a premixed nitromethane flame at low pressure

A premixed nitromethane/oxygen/argon flame at low pressure (4.67 kPa) has been investigated using tunable vacuum ultraviolet (VUV) photoionization and molecular-beam mass spectrometry. About 30 flame species including hydrocarbons, oxygenated and nitrogenous intermediates have been identified by measurements of photoionization efficiency spectra. Mole fraction profiles of the flame species have been determined by scanning burner position at some selected photon energies. The results indicate that N₂ and NO are the major nitrogenous products in the nitromethane flame. Compared with previous studies on nitromethane combustion, a number of unreported intermediates, including C₃H₄, C₄H₆, C₄H₈, C₂H₂O, C₂H₄O, CH₃CN, H₂CNHO, C₃H₃N and C₃H₇N, are observed in this work. Based on our experimental results and previous modeling studies, a detailed oxidation mechanism including 69 species and 314 reactions has been developed to simulate the flame structure. Despite some small discrepancies, the predictions by the modeling study are in reasonable agreement with the experimental results.     

Experimental and kinetic modeling study of combustion of JP-8, its surrogates and components in laminar premixed flows

Experimental and kinetic modeling studies are carried out to characterize premixed combustion of jet fuels, their surrogates, and reference components in laminar nonuniform flows. In previous studies, it was established that the Aachen surrogate made up of 80 % n-decane and 20 % trimethylbenzene by weight, and surrogate C made up of 57 % n-dodecane, 21 % methylcyclohexane and 22 % o-xylene by weight, reproduce key aspects of combustion of jet fuels in laminar nonpremixed flows. Here, these surrogates and a jet fuel are tested in premixed, nonuniform flows. The counterflow configuration is employed, and critical conditions of extinction are measured. In addition, the reference components tested are n-heptane, n-decane, n-dodecane, methylcyclohexane, trimethylbenzene, and o-xylene. Measured critical conditions of extinction of the Aachen surrogate and surrogate C are compared with those for the jet fuel. In general the alkanes n-heptane, n-decane, and n-dodecane, and methylcyclohexane are found to be more reactive than the aromatics o-xylene and trimethylbenzene. Flame structure and critical conditions of extinction are predicted for the reference components and the surrogates using a semi-detailed kinetic model. The predicted values are compared with experimental data. Sensitivity analysis shows that the lower reactivity of the aromatic species arises from the formation of resonantly stabilized radicals. These radicals are found to have a scavenging effect. The present study on premixed flows together with previous studies on nonpremixed flows show that the Aachen surrogate and surrogate C reproduce many aspects of premixed and nonpremixed combustion of jet fuels.     

Prediction of electron and ion concentrations in low-pressure premixed acetylene and ethylene flames

Flame stabilisation and extinction in a number of different flows can be affected by application of electric fields. Electrons and ions are present in flames, and because of charge separation, weak electric fields can also be generated even when there is no externally applied electric field. In this work, a numerical model incorporating ambipolar diffusion and plasma kinetics has been developed to predict gas temperature, species, and ion and electron concentrations in laminar premixed flames without applied electric fields. This goal has been achieved by combining the existing CHEMKIN-based PREMIX code with a recently developed methodology for the solution of electron temperature and transport properties that uses a plasma kinetics model and a Boltzmann equation solver. A chemical reaction set has been compiled from seven sources and includes chemiionisation, ion-molecule, and dissociative–recombination reactions. The numerical results from the modified PREMIX code (such as peak number densities of positive ions) display good agreement with previously published experimental data for fuel-rich, non-sooting, low-pressure acetylene and ethylene flames without applied electric fields.     

Subgrid combustion modeling of 3-D premixed flames in the thin-reaction-zone regime

Large eddy simulation (LES) is used to investigate three-dimensional (3D) lean premixed turbulent methane–air flames in the thin-reaction-zone regime. In this regime, the Kolmogorov scale is smaller than the preheat zone thickness, but larger than the reaction zone thickness. Past numerical studies of similar flames were primarily direct numerical simulation either in two-dimensions or using the artificially thickened flame approach in 3D. For an LES the effect of small (unresolved) scales on the scalar field must be, modeled accurately to capture the correct flame structure. A subgrid combustion model based on the linear-eddy-mixing (LEM) model is used within an LES framework (called LEM–LES hereafter) to capture the 3D flame-structure of the highly stretched premixed flames. A finite-rate, one-step methane–air chemistry with a non-unity Lewis number formulation is used in this study. The simulated flame structure resembles flames experimentally studied in the thin-reaction-zone regime. Even though the preheat zone is broadened by the penetration of small eddies, the chemical reaction zone remains thin and localized. This feature is captured properly in the current LEM–LES approach. The flame structure and other statistics such as the flame area evolution, curvature, and strain-rate statistics computed using the LEM–LES are also in good agreement with the past DNS studies.     

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

Experimental and kinetic modeling study of oxidation of acetonitrile

Oxidation of acetonitrile has been studied in a flow reactor in the absence and presence of nitric oxide. The experiments were conducted at atmospheric pressure in the temperature range 1150–1450 K, varying the excess air ratio from slightly fuel-lean to very lean. Oxidation of CH₃CN was slow below 1300 K. Nitric oxide, hydrogen cyanide and nitrous oxide were detected as important products. A detailed chemical kinetic model for oxidation of acetonitrile was developed, based on a critical evaluation of data from literature. The rate coefficients for the reactions of CH₃CN and CH₂CN with O₂ were calculated from ab initio theory. Modeling predictions were in satisfactory agreement with experiments. Calculations were sensitive to thermal dissociation of CH₃CN and to the branching fraction for CH₃CN + OH to CH₂CN + H₂O and HOCN + CH₃, respectively. More work is desirable for these steps, as well as for reactions of CH₂CN and HCCN.     

Experimental and kinetic study of 2,4,4-trimethyl-1-pentene and iso-octane in laminar flames

Laminar flame speeds of 2,4,4-trimethyl-1-pentene are investigated at equivalence ratios of 0.7–1.6, initial temperatures of 298–453 K and initial pressures of 0.1–0.5 MPa. The comparison between 2,4,4-trimethyl-1-pentene and iso-octane is also performed. Results show that 2,4,4-trimethyl-1-pentene has faster laminar flame speed than iso-octane. Chemical kinetic models (Metcalfe model, Modified model I) were tested against the present experimental data. The laminar flame speeds are apparently over-estimated by the Metcalfe model and under-predicted by the Modified model I. Therefore, high-level quantum mechanical calculations were used to revise the Modified model I to obtain Modified model II and it can give fairly good prediction at various conditions on laminar flame speeds. In addition, the chemical kinetic analysis was conducted. The analysis indicates both thermal and kinetic effects result in the discrepancy of laminar flame speeds between 2,4,4-trimethyl-1-pentene and iso-octane. Furthermore, IC₄H₈ plays a dominant role in laminar flame speeds of 2,4,4-trimethyl-1-pentene and iso-octane.     

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.