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.     

The effect of preferential diffusion on soot formation in a laminar ethylene/air diffusion flame

The influence of preferential diffusion on soot formation in a laminar ethylene/air diffusion flame was investigated by numerical simulation using three different transport property calculation methods. One simulation included preferential diffusion and the other two neglected preferential diffusion. The results show that the neglect of preferential diffusion or the use of unity Lewis number for all species results in a significant underprediction of soot volume fraction. The peak soot volume fraction is reduced from 8.0 to 2.0 ppm for the studied flame when preferential diffusion is neglected in the simulation. Detailed examination of numerical results reveals that the underprediction of soot volume fraction in the simulation neglecting preferential diffusion is due to the slower diffusion of some species from main reaction zone to PAH and soot formation layer. The slower diffusion of these species causes lower PAH formation rate and thus results in lower soot inception rate and smaller particle surface area. The smaller surface area further leads to smaller surface growth rate. In addition, the neglect of preferential diffusion also leads to higher OH concentration in the flame, which causes the higher specific soot oxidation rate. The lower inception rate, smaller surface growth rate and higher specific oxidation rate results in the lower soot volume fraction when preferential diffusion is neglected. The finding of the paper implies the importance of preferential diffusion for the modeling of not only laminar but maybe also some turbulent flames.     

Modelling soot formation in a premixed flame using an aromatic-site soot model and an improved oxidation rate

An updated rate of O₂ oxidation of one to four ring polyaromatic hydrocarbons in premixed flames is presented based on density function theory simulations of oxygen attack at different radical sites on various PAHs. The rate is in agreement with other rates found in the literature; however, it is several orders of magnitude lower than the currently accepted oxidation rate of multi-ring aromatic species, including soot. Simulations are presented of a premixed flame using this improved rate and a new advanced soot particle model, which is developed in this paper. This model includes unprecedented detail of the particles in the ensemble, including the aromatic content, C/H composition and primary-particle aggregate structure. The O₂ oxidation rate calculated in this paper is shown to give a better prediction of particle number density and soot volume fraction for a premixed flame. The predicted particle size distributions are shown also to describe better the experimental data. Predicted C/H ratio and PAH size distributions are shown for the flame. Computed TEM-style images are compared to experimental TEM images, which show that the aggregate structure of the particles is well predicted.     

Addition of NO2 to a laminar premixed ethylene-air flame: Effect on soot formation

Experiments were conducted on a laminar premixed ethylene-air flame at equivalence ratios of 2.34 and 2.64. Comparisons were made between flames with 5% NO₂ added by volume. Soot volume fraction was measured using light extinction and light scattering and fluorescence measurements were also obtained to provide added insight into the soot formation process. The flame temperature profiles in these flames were measured using a spectral line reversal technique in the non-sooting region, while two-color pyrometry was used in the sooting region. Chemical kinetics modeling using the PREMIX 1-D laminar flame code was used to understand the chemical role of the NO₂ in the soot formation process. The modeling used kinetic mechanisms available in the literature. Experimental results indicated a reduction in the soot volume fraction in the flame with NO₂ added and a delay in the onset of soot as a function of height above the burner. In addition, fluorescence signals—often argued to be an indicator of PAH—were observed to be lower near the burner surface for the flames with NO₂ added as compared to the baseline flames. These trends were captured using a chemical kinetics model that was used to simulate the flame prior to soot inception. The reduction in soot is attributed to a decrease in the H-atom concentration induced by the reaction with NO₂ and a subsequent reduction in acetylene in the pre-soot inception region.     

Effects of cryogenic temperature on premixed hydrogen/air flame propagation in a closed channel

As a carbon-free fuel, hydrogen has received significant attention recently since it can help enable low-carbon-economy. Hydrogen has very broad flammability range and very low minimum ignition energy, and thereby there are severe safety concerns for hydrogen transportation and utilization. Cryo-compressed hydrogen is popularly used in practice. Therefore, it is necessary to investigate the combustion properties of hydrogen at extremely low or cryogenic temperatures. This study aims to assess and interpret the effects of cryogenic temperature on premixed hydrogen/air flame propagation and acceleration in a thin closed channel. Different initial temperatures ranging from normal temperature (T₀ = 300 K) to cryogenic temperature (T₀ = 100 K) are considered. Both one- and two-dimensional hydrogen/air flames are investigated through transient simulations considering detailed chemistry and transport. It is found that when the initial temperature decreases from T₀ = 300 K to T₀ = 100 K, the expansion ratio and equilibrium pressure both increase substantially while the laminar flame speeds relative to unburned and burned gasses decrease moderately. The one-dimensional flame propagation is determined by laminar flame speed and thereby the combustion duration increases as the initial temperature decreases. However, the opposite trend is found to happen to two-dimensional flame propagation, which is mainly controlled by the flame surface area increase due to the no-slip side wall constraint and flame instability. Based on the change in flame surface area, three stages including the initial acceleration, steady burning and rapid acceleration are identified and investigated. It is demonstrated that the large expansion ratio and high pressure rise at cryogenic temperatures can significantly increase the flame surface area in early stage and promote both Darrieus-Landau instability (hydrodynamic instability) and Rayleigh-Taylor instability in later stage. These two instabilities can substantially increase the flame surface area and thereby accelerate flame propagation in hydrogen/air mixtures at cryogenic temperatures. The present study provides useful insights into the fundamental physics of hydrogen flames at extremely low temperatures, and is closely related to hydrogen safety.     

Determination of spatial temperature distributions in a laminar premixed ethylene/air flame by laser-induced atomic fluorescence

A quick and simple detection system for spatially resolved temperature measurements in flames based on laser-induced thermally assisted atomic line fluorescence of seeded rubidium atoms is described. The fluorescence light from two atomic states is dispersed and simultaneously recorded by a CCD camera. The fluorescence ratio distributions lead directly to absolute temperature distributions. The practical use, the spatial and temperature resolution and error limits of the method are discussed and compared with other procedures for temperature measurements.     

Modeling of soot aggregate formation and size distribution in a laminar ethylene/air coflow diffusion flame with detailed PAH chemistry and an advanced sectional aerosol dynamics model

Soot aggregate formation and size distribution in a laminar ethylene/air coflow diffusion flame is modeled with a PAH-based soot model and an advanced sectional aerosol dynamics model. The mass range of solid soot phase is divided into 35 discrete sections and two variables are solved for in each section. The coagulation kernel of soot aggregates is calculated for the entire Knudsen number regime. Radiation from gaseous species and soot are calculated by a discrete-ordinate method with a statistical narrow-band correlated-k based band model. The discretized sectional soot equations are solved simultaneously to ensure convergence. Parallel computation with the domain decomposition method is used to save computational time. The flame temperature, soot volume fraction, primary particle size and number density are well reproduced. The number of primary particles per aggregate is overpredicted. This discrepancy is presumably associated with the unitary coagulation efficiency assumption in the current sectional model. Along the maximum soot volume fraction pathline, the number-based and mass-based aggregate size distribution functions are found to evolve from unimodal to bimodal and finally to unimodal again. The different shapes of these two aggregate size distribution functions indicate that the total number and mass of aggregates are dominated by aggregates of different sizes. The PAH-soot condensation efficiency γ is found to have a small effect on soot formation when γ is larger than 0.5. However, the soot level and primary particle number density are significantly overpredicted if the PAH-soot condensation process is neglected. Generally, larger γ predicts lower soot level and primary particle number density. Further study on soot aggregate coagulation efficiency should be pursued and more experimental data on soot aggregate structure and size distribution are needed for improving the current sectional soot model and for better understanding the complex soot aggregation phenomenon.     

Effect of fuel/air ratio and aromaticity on the molecular weight distribution of soot in premixed n-heptane flames

Soot growth from inception to mass-loading is studied in a wide range of molecular weights (MW) from 10⁵ to 10¹⁰u by means of size exclusion chromatography (SEC) coupled with on-line UV-visible spectroscopy. The evolution of MW distributions of soot is also numerically predicted by using a detailed kinetic model coupled with a discrete-sectional approach for the modeling of the gas-to-particle process. Two premixed flames burning n-heptane in slightly sooting and heavily sooting conditions are studied. The effect of aromatic addition to the fuel is studied by adding n-propylbenzene (10% by volume) to n-heptane in the heavily sooting condition. A progressive reduction of the MW distribution from multimodal to unimodal is observed along the flames testifying the occurrence of particle growth and agglomeration. These processes occur earlier in the aromatic-doped n-heptane flame due to the overriding role of benzene on soot formation which results in bigger young soot particles. Modeled MW distributions are in reasonable agreement with experimental data although the model predicts a slower coagulation process particularly in the slightly sooting n-heptane flame. Given the good agreement between model predictions and experiments, the model is used to explore the role of fuel chemistry on MW distributions. Two flames of n-heptane and n-heptane/n-propylbenzene in heavily sooting conditions with the same temperature profile and inert dilution are modeled. The formation of larger soot particles is still evident in the n-heptane/n-propylbenzene flame with respect to the n-heptane flame in the same operating conditions of temperature and dilution. In addition the model predicts a larger formation of molecular particles in the flame containing n-propylbenzene and shows that soot inception occurs in correspondence of their maximum formation thus indicating the importance of molecular growth in soot inception.     

Implementation of an advanced fixed sectional aerosol dynamics model with soot aggregate formation in a laminar methane/air coflow diffusion flame

An advanced fixed sectional aerosol dynamics model describing the evolution of soot particles under simultaneous nucleation, coagulation, surface growth and oxidation processes is successfully implemented to model soot formation in a two-dimensional laminar axisymmetric coflow methane/air diffusion flame. This fixed sectional model takes into account soot aggregate formation and is able to provide soot aggregate and primary particle size distributions. Soot nucleation, surface growth and oxidation steps are based on the model of Fairweather et al. Soot equations are solved simultaneously to ensure convergence. The numerically calculated flame temperature, species concentrations and soot volume fraction are in good agreement with the experimental data in the literature. The structures of soot aggregates are determined by the nucleation, coagulation, surface growth and oxidation processes. The result of the soot aggregate size distribution function shows that the aggregate number density is dominated by small aggregates while the aggregate mass density is generally dominated by aggregates of intermediate size. Parallel computation with the domain decomposition method is employed to speed up the calculation. Three different domain decomposition schemes are discussed and compared. Using 12 processors, a speed-up of almost 10 is achieved which makes it feasible to model soot formation in laminar coflow diffusion flames with detailed chemistry and detailed aerosol dynamics.     

Effect of ammonia addition on suppressing soot formation in methane co-flow diffusion flames

Due to issues surrounding carbon dioxide emissions from carbon-containing fuels, there is growing interest in ammonia (NH₃) as an alternative combustion fuel. One attractive method of burning NH₃ is to co-fire it with hydrocarbons, such as natural gas, and in this case soot formation is possible. To begin understanding the influence of NH₃ on soot formation when co-fired with hydrocarbons, soot volume fractions and mole fractions of gas-phase species were computationally and experimentally interrogated for CH₄ flames with up to 40% NH₃ by volumetric fuel fraction. Mole fractions of gas-phase species, including C₂H₂ and C₆H₆, were measured with on-line electron impact mass spectrometry, and soot volume fractions were obtained via color-ratio pyrometry. The simulations employed a detailed chemical mechanism developed for capturing nitrogen interactions with hydrocarbons during combustion. The results are compared to findings in N₂CH₄ flames, in order to separate thermal and dilution effects from the chemical influence of NH₃ on soot formation. Experimentally, C₂H₂ concentrations were found to decrease slightly for the NH₃CH₄ flames relative to N₂CH₄ flames, and a stronger suppression of C₆H₆ was found for NH₃ relative to N₂ additions. The measured results show a strong suppression of soot with the addition of NH₃, with soot concentrations reduced by over a factor of 10 with addition of up to 20% or more NH₃ by mole fraction. The model satisfactorily captured relative differences in maximum centerline C₂H₂, C₆H₆, and soot concentrations with addition of N₂, but was unable to match measured differences in NH₃CH₄ flames. These results highlight the need for an improved understanding of fuel-nitrogen interactions with higher hydrocarbons to enable accurate models for predicting particulate emissions from NH₃/hydrocarbon combustion.     

Effects of gas and soot radiation on soot formation in counterflow ethylene diffusion flames

Numerical study of soot formation in counterflow ethylene diffusion flames at atmospheric pressure was conducted using detailed chemistry and complex thermal and transport properties. Soot kinetics was modelled using a semi-empirical two-equation model. Radiation heat transfer was calculated using the discrete-ordinates method coupled with an accurate band model. The calculated soot volume fractions are in reasonably good agreement with the experimental results in the literature. The individual effects of gas and soot radiation on soot formation were also investigated.     

Understanding soot particle size evolution in laminar ethylene/air diffusion flames using novel soot coalescence models

Two coalescence models based on different merging mechanisms are introduced. The effects of the soot coalescence process on soot particle diameter predictions are studied using a detailed sectional aerosol dynamic model. The models are applied to a laminar ethylene/air diffusion flame, and comparisons are made with experimental data to validate the models. The implementation of coalescence models significantly improves the agreement of prediction of particle diameters with the experimental data. Sensitivity of the soot prediction to the coalescence parameters is analysed. Finally, an update to the coalescence model based on experimental observations of soot particles in the flame oxidation regions has been introduced to improve its predicting capabilities.