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.     

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.     

Size distribution and morphology of nascent soot in premixed ethylene flames with and without benzene doping

Particle size distribution functions of nascent soot formed in four burner-stabilized, premixed ethylene-oxygen-argon flames were studied in a spatially resolved manner by online sampling/scanning mobility particle sizer. Particle morphology was analyzed by atomic force microscopy (AFM) of substrate-deposited samples. Two of the four flames were doped with benzene. An aerosol electrometer is introduced to extend the lower detection limit to around 1.5 nm in diameter. The results show that the bimodal behavior of particle size is applicable to all premixed ethylene flames studied. The variation of the size distribution from flame to flame is conclusively attributed to flame temperature variation. Under the condition of an equal carbon concentration, benzene doping leads to negligible changes in the characteristics of the size distribution. For all flames studied, AFM observations show that nascent soot is liquid-like and spreads extensively upon impact on a substrate surface.     

In-flame soot particle structure on the up- and down-swirl side of a wall-interacting jet in a small-bore diesel engine

This study shows how the structure of soot particles within the flame changes due to the relative direction of the swirl flow in a small-bore diesel engine in which significant flame–wall interactions cause about half of the flame travelling against the swirl flow while the other half penetrating in the same direction. The thermophoresis-based particle sampling method was used to collect soot from three different in-flame locations including the flame–wall impingement point near the jet axis and the two 60° off-axis locations on the up-swirl and down-swirl side of the wall-interacting jet. The sampled soot particle images were obtained using transmission electron microscopes and the image post-processing was conducted for statistical analysis of size distribution of soot primary particles and aggregates, fractal dimension, and sub-nanoscale parameters such as the carbon layer fringe length, tortuosity, and spacing. The results show that the jet-wall impingement region is dominated by many small immature particles with amorphous internal structure, which is very different to large, fractal-like soot aggregates sampled from 60° downstream location on the down-swirl side. This structure variation suggests that the small immature particles underwent surface growth, coagulation and aggregation as they travelled along the piston-bowl wall. During this soot growth, the particle internal structure exhibits the transformation from amorphous carbon segments to a typical core–shell structure. Compared to those on the down-swirl side, the soot particles sampled on the up-swirl side show much lower number counts and more compact aggregates composed of highly concentrated primary particles. This soot aggregate structure, together with much narrower carbon layer gap, indicates higher level of soot oxidation on the up-swirl side of the jet.     

The Formation of Carbon in Combustion and how to Quantify the Impact on Human Health

We summarize our current research on combustion aerosols. First, sampling devices for the analyses of flame gases are described. The flame gas samples are investigated by mass spectroscopy and by standard aerosol techniques. Time-of-flight mass spectroscopy is well suited to study formation and growth of soot precursor molecules. Fullerenes can also be seen in some mass spectra of flame gases. Presumably, the fullerenes are evaporated from small soot particles in the mass spectrometer by the ionizing laser. Size spectra of soot particles from the flame are presented. The flame is optionally seeded with palladium aerosol to demonstrate that the particle size distribution is not altered during the sampling procedure. It is found that soot particles are already present low in the flame where large molecules are absent.Photoemission is applied to study surface properties of soot particles from the flame. It is shown that the surface of the particles is covered with polycyclic aromatic hydrocarbons (PAH). The PAH can be removed by heating and the properties of the carbon core are revealed. One can thereby distinguish a soot growth from a soot burnout region in the flame. Time-resolved desorption experiments of perylene (a PAH) from model aerosol particles are presented. It is shown that they follow a first order rate law. The photoelectric PAH sensor is introduced as a personal air quality monitor. The danger from inhaling combustion aerosol can be expressed in units of standard cigarettes.     

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.     

Experimental and numerical study of the evolution of soot primary particles in a diffusion flame

The evolution of primary soot particles is studied experimentally and numerically along the centreline of a co-flow laminar diffusion flame. Soot samples from a flame fueled with C₂H₄ are taken thermophoretically at different heights above the burner (HAB), their size and nano-structure are analysed through TEM. The experimental results suggest that after inception, the nascent soot particles coagulate and coalesce to form larger primary particles (?～?5 to 15 nm). As these primary particles travel along the centreline, they grow mainly due coagulation and condensation and a layer of amorphous hydrocarbons (revealed by HRTEM) forms on their surface. This amorphous layer appears to promote the aggregation of primary particles to form fractal structures. Fast carbonisation of the amorphous layer leads to a graphitic-like shell around the particles. Further graphitization compacts the primary particles, resulting in a decrease of their size. Towards the flame tip the primary particles decrease in size due to rapid oxidation. A detailed population balance model is used to investigate the mechanisms that are important for prediction of primary particle size distributions. Suggestions are made regarding future model development efforts. Simulation results indicate that the primary particle size distributions are very sensitive to the parameterization of the coalescence and particle rounding processes. In contrast, the average primary particle size is less sensitive to these parameters. This demonstrates that achieving good predictions for the average primary particle size does not necessarily mean that the distribution has been accurately predicted.     

Instantaneous in-flame measurement of soot volume fraction, primary particle diameter, and aggregate radius of gyration via auto-compensating laser-induced incandescence and two-angle elastic light scattering

A new combination of soot diagnostics employing two-angle elastic light scattering and laser-induced incandescence is described that is capable of producing non-intrusive, instantaneous, and simultaneous, in situ measurements of soot volume fraction, primary particle size, and aggregate radius of gyration within flames. Controlled tests of the new apparatus on a well-characterized laminar flame show good agreement with existing measurements in the literature. From a detailed and comprehensive Monte Carlo uncertainty analysis of the results, it was found that the uncertainty in all three measured parameters is dominated by knowledge of soot properties and aggregation behavior. The soot volume fraction uncertainty is dominated by uncertainty in the soot refractive index light absorption function; the primary particle diameter uncertainty is dominated by uncertainty in the fractal prefactor; while the uncertainty in the aggregate radius of gyration is dominated by the uncertainty in the width of the distribution of aggregate sizes.     

Soot particle size distributions in a well-stirred reactor/plug flow reactor

A well-stirred reactor (WSR) followed by a plug flow reactor (PFR) is being used to study polycyclic aromatic hydrocarbon (PAH) growth and soot inception. Soot size distributions were measured using a dilution probe followed by a nano-differential mobility analyzer (Nano-DMA). A rapid insertion probe was fabricated to thermophoretically collect particles from the reactor for transmission electron microscopy (TEM) imaging. Results are presented on the effect of equivalence ratio on the soot size distributions obtained for fixed dilution ratio, the effect of dilution ratio on the soot size distributions obtained for fixed equivalence ratio, and the effect of temperature on the soot size distributions obtained for fixed equivalence ratio. In addition to particle sizing measurements, gas samples were analyzed by a gas chromatograph to determine the concentration of gaseous species in the PFR thought to be important in soot formation. Our soot size distribution measurements demonstrate that the mixing conditions in the flame zone affect whether or not a nucleation mode was detected in the size distribution.     

Modeling nongray gas-phase and soot radiation in luminous turbulent nonpremixed jet flames

Much progress has been made in radiative heat transfer modeling with respect to treatment of nongray radiation from both gas-phase species and soot particles, while radiation modeling in turbulent flame simulations is still in its infancy. Aiming at reducing this gap, this paper introduces state-of-the-art models of gas-phase and soot radiation to turbulent flame simulations. The full-spectrum k-distribution method (Modest, M.F., 2003, Journal of Quantitative Spectroscopy & Radiative Transfer, 76, 69–83) is implemented into a three-dimensional unstructured CFD code for nongray radiation modeling. The mixture full-spectrum k-distributions including nongray absorbing soot particles are constructed from a narrow-band k-distribution database created for individual gas-phase species, and an efficient scheme is employed for their construction in CFD simulations. A detailed reaction mechanism including NO _x and soot kinetics is used to predict flame structure, and a detailed soot model using a method of moments is employed to determine soot particle size distributions. A spherical-harmonic P₁ approximation is invoked to solve the radiative transfer equation. An oxygen-enriched, turbulent, nonpremixed jet flame is simulated, which features large concentrations of gas-phase radiating species and soot particles. Nongray soot modeling is shown to be of greater importance than nongray gas modeling in sooty flame simulations, with gray soot models producing large errors. The nongray treatment of soot strongly influences flame temperatures in the upstream and the flame-tip region and is essential for accurate predictions of NO. The nongray treatment of gases, however, weakly influences upstream flame temperatures and, therefore, has only a small effect on NO predictions. The effect of nongray soot radiation on flame temperature is also substantial in downstream regions where the soot concentration is small. Limitations of the P₁ approximation are discussed for the jet flame configuration; the P₁ approximation yields large errors in the spatial distribution of the computed radiative heat flux for highly anisotropic radiation fields such as those in flames with localized, near-opaque soot regions.     

Modelling nongrey gas-phase and soot radiation in luminous turbulent nonpremixed jet flames

Much progress has been made in radiative heat transfer modelling with respect to the treatment of nongrey radiation from both gas-phase species and soot particles, while radiation modelling in turbulent flame simulations is still in its infancy. Aiming at reducing this gap, this paper introduces state-of-the-art models of gas-phase and soot radiation to turbulent flame simulations. The full-spectrum k-distribution method (M.F. Modest, 2003, Journal of Quantitative Spectroscopy & Radiative Transfer, 76, 69–83) is implemented into a three-dimensional unstructured computational fluid dynamics (CFD) code for nongrey radiation modelling. The mixture full-spectrum k-distributions including nongrey absorbing soot particles are constructed from a narrow-band k-distribution database created for individual gas-phase species, and an efficient scheme is employed for their construction in CFD simulations. A detailed reaction mechanism including NO _x and soot kinetics is used to predict flame structure, and a detailed soot model using a method of moments is employed to determine soot particle size distributions. A spherical harmonic P₁ approximation is invoked to solve the radiative transfer equation. An oxygen-enriched, turbulent, nonpremixed jet flame is simulated, which features large concentrations of gas-phase radiating species and soot particles. Nongrey soot modelling is shown to be of greater importance than nongrey gas modelling in sooty flame simulations, with grey soot models producing large errors. The nongrey treatment of soot strongly influences flame temperatures in the upstream and the flame-tip region and is essential for accurate predictions of NO. The nongrey treatment of gases, however, weakly influences upstream flame temperatures and, therefore, has only a small effect on NO predictions. The effect of nongrey soot radiation on flame temperature is also substantial in downstream regions where the soot concentration is small. Limitations of the P₁ approximation are discussed for the jet flame configuration; the P₁ approximation yields large errors in the spatial distribution of the computed radiative heat flux for highly anisotropic radiation fields such as those in flames with localized, near-opaque soot regions.     

Investigation of the transition from lightly sooting towards heavily sooting co-flow ethylene diffusion flames

Laminar, sooting, ethylene-fuelled, co-flow diffusion flames at atmospheric pressure have been studied experimentally and theoretically as a function of fuel dilution by inert nitrogen. The flames have been investigated experimentally using a combination of laser diagnostics and thermocouple-gas sampling probe measurements. Numerical simulations have been based on a fully coupled solution of the flow conservation equations, gas-phase species conservation equations with complex chemistry and the dynamical equations for soot spheroid growth. Predicted flame heights, temperatures and the important soot growth species, acetylene, are in good agreement with experiment. Benzene simulations are less satisfactory and are significantly under-predicted at low dilution levels of ethylene. As ethylene dilution is decreased and soot levels increase, the experimental maximum in soot moves from the flame centreline toward the wings of the flame. Simulations of the soot field show similar trends with decreasing dilution of the fuel and predicted peak soot levels are in reasonable agreement with the data. Computations are also presented for modifications to the model that include: (i) use of a more comprehensive chemical kinetics model; (ii) a revised inception model; (iii) a maximum size limit to the primary particle size; and (iv) estimates of radiative optical thickness corrections to computed flame temperatures.     

热电偶沉积法测量乙烯非预混火焰的试验研究

烟黑容积份额的测量是研究烟黑生成的反应机理的额的薪方法.本文详细描述了采用热电偶沉积法测量烟黑容积份额的理论基础和数据处理过程,并将此方法应用于层流乙烯非预混火焰的测量中.测量结果表明,该火焰中烟黑容积份额的分布同火焰结构和火焰温度都有关.     

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.     

Surface deposition and coagulation efficiency of combustion generated nanoparticles in the size range from 1 to 10 nm

The size distribution of the nanoparticles formed in premixed ethylene–air flames and collected thermophoretically on mica cleaved substrates is obtained by atomic force microscopy (AFM). The distribution function extends from 1 to about 5 nm in non-sooting flames and in the soot pre-inception region of the richer flames, while it becomes bimodal and larger particles are formed in the soot inception region of the slightly sooting flames. The distribution is compared with the size distribution of nano-sized organic carbon (NOC) and soot particles, obtained by “in situ” multi-wavelength extinction and light scattering methods. The deposition efficiency is estimated from the differences between these two size distribution functions as a function of the equivalent diameter of the nanoparticles. Furthermore, the coagulation coefficient of particles in flame is obtained from the temporal evolution of the number concentration of the nanoparticles inside the flames. NOC particles, which are rapidly produced in locally rich combustion regions, have peculiar properties since their sticking coefficient both for coagulation and adhesion result to be orders of magnitudes lower than that expected by larger aerosols, like soot particles. The experimental results are interpreted by modelling the van der Waals interactions of the nanoparticles in terms of Lennard-Jones potentials and in the framework of the gas kinetic theory. The estimated adhesion and coagulation efficiencies are in good agreement with those calculated from AFM and optical data. The very low efficiency values observed for the smaller particles could be ascribed to the high energy of these particles due to their Brownian motion, which causes thermal rebound effects prevailing over adhesion mechanisms due to van der Waals forces.     

Influence of soot particle aggregation on time-resolved laser-induced incandescence signals

Laser-induced incandescence (LII) is a versatile technique for quantitative soot measurements in flames and exhausts. When used for particle sizing, the time-resolved signals are analysed as these will show a decay rate dependent on the soot particle size. Such an analysis has traditionally been based on the assumption of isolated primary particles. However, soot particles in flames and exhausts are usually aggregated, which implies loss of surface area, less heat conduction and hence errors in estimated particle sizes. In this work we present an experimental investigation aiming to quantify this effect. A soot generator, based on a propane diffusion flame, was used to produce a stable soot stream and the soot was characterised by transmission electron microscopy (TEM), a scanning mobility particle sizer (SMPS) and an aerosol particle mass analyzer coupled in series after a differential mobility analyzer (DMA-APM). Despite nearly identical primary particle size distributions for three selected operating conditions, LII measurements resulted in signal decays with significant differences in decay rate. However, the three cases were found to have quite different levels of aggregation as shown both in TEM images and mobility size distributions, and the results agree qualitatively with the expected effect of diminished heat conduction from aggregated particles resulting in longer LII signal decays. In an attempt to explain the differences quantitatively, the LII signal dependence on aggregation was modelled using a heat and mass transfer model for LII given the primary particle and aggregate size distribution data as input. Quantitative agreement was not reached and reasons for this discrepancy are discussed.     

基于光谱诊断的烃类火焰碳烟形成机理研究综述

碳烟主要是烃类燃料不完全燃烧生成的产物,其对人类健康、空气质量以及燃烧装置的使用寿命都会产生有害影响。碳烟生成是一个复杂的物理化学过程,控制碳烟排放,需要克服碳烟生成和燃烧过程中物理和化学演化的巨大差异,这些差异表现为对碳烟纳观结构和表面官能团随碳烟氧化活性反应变化的深入探索研究。近些年,研究人员对碳烟的生成机理开展了系列研究,对碳烟生成各个物理化学反应阶段有了一定认识。结合光谱诊断技术可深入了解燃烧系统碳烟形成过程,确定碳烟颗粒分子组成、精细结构、浓度分布等特征,也可从碳烟结构变化、黑体辐射强度等方面详细了解碳烟形成过程。该文旨在阐述光谱诊断技术对烃类火焰碳烟表征的研究进展和发展趋势,探讨LIBS, LII和LIF等作为诊断工具在包含背景辐射的火焰中检测碳烟生成过程产生辐射强度准确性等问题。主要介绍了烃类火焰碳烟的形成机理(从前驱体产生、生长到颗粒生成、凝聚,最后进行颗粒氧化)。总结了探测碳烟性质光谱诊断方法的应用以及光谱诊断技术对燃烧过程中碳烟表征的研究现状,包括对碳烟体积分数、温度和基于图像处理的碳烟结构表征,反应碳烟前驱体(多环芳烃)、反应气氛、温度等对碳烟颗粒物生成的影响。最...     

Comparison of LII derived soot temperature measurements with LII model predictions for soot in a laminar diffusion flame

Laser-induced incandescence (LII) was used to derive temperatures of pulsed laser heated soot particles from their thermal emission intensities detected at two wavelengths in a laminar ethylene/air co-annular diffusion flame. The results are compared to those of a numerical nanoscale heat and mass transfer model. Both aggregate and primary particle soot size distributions were measured using transmission electron microscopy (TEM). The model predictions were numerically averaged over these experimentally derived size distributions. The excitation laser wavelength was 532 nm, and the LII signal was detected at 445 nm and 780 nm. A wide range of laser fluence from very low to moderate (0.13 to 1.56 mJ/mm²) was used in the experiments. A large part of the temporal decay curve, beginning 12–15 nsec after the peak of the laser excitation pulse, is successfully described by the model, resulting in the determination of accommodation coefficients, which varies somewhat with soot temperature and is in the range of 0.36 to 0.46. However, in the soot evaporative regime, the model greatly overpredicts the cooling rate shortly after the laser pulse. At lower fluences, where evaporation is negligible, the initial experimental cooling rates, immediately following the laser pulse, are anomalously high. Potential physical processes that could account for these effects are discussed. From the present data the soot absorption function, E(m), of 0.4 at 532 nm is obtained. A procedure for correcting the measured signals for the flame radiation is presented. It is further shown that accounting for the local gas temperature increase due to heat transfer from soot particles to the gas significantly improves the agreement in the temperature dependence of soot cooling rates between model and experiments over a large range of laser fluences.     

A comparative study of nanoparticles in premixed flames by scanning mobility particle sizer, small angle neutron scattering, and transmission electron microscopy

Scanning mobility particle sizer (SMPS) and transmission electron microscopy (TEM) studies were conducted for TiO₂ and soot particles. The TiO₂ particles were produced from a premixed stagnation ethylene-oxygen-argon flame (? = 0.36) doped with titanium tetraisopropoxide. Soot was generated from a burner-stabilized premixed ethylene-oxygen-argon flame (? = 2.5). The close agreement among SMPS, TEM, and X-ray diffraction results for TiO₂ nanoparticles demonstrates that the probe sampling/mobility measurement technique is accurate for on-line analysis of the size distribution of particles as small as 3 nm in diameter. In the case of soot, notable disagreement between the SMPS and TEM sizes was found and attributable to the fact that the soot taken from the flame studied herein is liquid-like and that upon deposition on the TEM grid, the primary particles do not retain their sphericity. This interpretation is supported by measurements with photo ionization aerosol mass spectrometry, small angle neutron scattering, and thermocouple particle densitometry.