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.     

Analysis of uncertainties in instantaneous soot volume fraction measurements using two-dimensional, auto-compensating, laser-induced incandescence (2D-AC-LII)

Laser-induced incandescence (LII) is an optical measurement technique capable of measuring soot volume fraction over a wide range of conditions. However, development of two-dimensional auto-compensating LII (2D-AC-LII) in the literature has been limited and until now, instantaneous measurements have not been demonstrated. In this paper, we successfully demonstrate instantaneous 2D-AC-LII soot volume fraction (SVF) measurements in an ethylene-air co-annular diffusion flame. Results were then used to support a detailed uncertainty analysis based on a Monte-Carlo simulation. Agreement between both the instantaneous and average SVF measurements with published data from attenuation measurements under identical conditions was found to be good. Uncertainties are discussed both in terms of an overall accuracy of the SVF measurement, which is strongly dominated by uncertainty in the optical properties of soot, and the comparative uncertainties with optical properties fixed. The uncertainty in an instantaneous 2D determination of SVF for a comparative measurement is dominated by photon shot noise, and in regions of high soot volume fraction it is below 25% (95% confidence interval). Shot noise uncertainty could be further reduced with additional pixel averaging at the expense of spatial resolution. This diagnostic shows significant promise for quantitative planar soot concentration measurements within turbulent flames.     

Three-angle scattering/extinction versus TEM measurements on soot in premixed ethylene/air flame

In this work, a three-angle scattering and extinction technique has been applied in order to study soot formation and growth in a rich ethylene/air premixed flame (Φ=2.34). The Rayleigh–Debye–Gans theory together with the fractal-like approach has been applied to derive soot parameters, in terms of volume concentration and morphology. A mathematical procedure is presented to obtain the radius of gyration by considering scattering signals collected at two supplementary angles. TEM measurements, carried out at different locations on the flame axis, are used to derive some parameters, such as fractal prefactor, fractal dimension and size distribution, to be entered in the mathematical treatment of optical data. The radius of gyration and the primary particle size as obtained by TEM and by optical measurements are compared. Good agreement has been found in the upper part of the flame investigated. Discrepancies observed low in the flame are discussed.     

Numerical accuracy of “equivalent” spherical approximations for computing ensemble-averaged scattering properties of fractal soot aggregates

Numerical accuracy is quantitatively assessed in conjunction with the application of four “equivalent” spherical approximations in the computation of the optical properties of small aggregate soot particles. The approximations are based on equal volume, equal surface area, the radius of gyration, and a collection of independent spheres with the same volume and the same volume-to-projected area ratio as the original nonspherical particle. A diffusion-limited cluster-cluster aggregation algorithm is used to specify the geometries of soot particles. Furthermore, the Generalized Multi-particle Mie (GMM) method is utilized to compute the single-scattering properties of individual soot aggregate particles assumed to be randomly oriented in space. The ensemble-averaged single-scattering properties of the particles are obtained by accounting for the probability distribution functions (PDF) of the number of monomers per aggregate at two wavelengths, 0.628 and 1.1 μm. It is shown that all of the aforementioned equivalent-spherical approximations lead to large errors in the computation of the phase function.     

Laser induced incandescence measurements of soot volume fraction and effective particle size in a laminar co-annular non-premixed methane/air flame at pressures between 0.5–4.0 MPa

An auto-compensating laser-induced incandescence (AC-LII) technique was applied for the first time to measure soot volume fraction (SVF) and effective primary particle diameter (dp_eff) in a high pressure methane/air non-premixed flame. The measured dp_eff profiles had annular structures and radial symmetry, and the particle size increased with increasing pressure. LII-determined SVFs were lower than those measured by a line of sight attenuation (LOSA) technique. The LOSA measured soot volume fractions were corrected for light scattering using the Rayleigh–Debye–Gans polydisperse fractal aggregate (RDG-PFA) theory, the dp_eff data, and assumptions regarding the soot aggregate size distribution. The correction dramatically improved agreement between data obtained using these two measurement techniques. Qualitatively, soot volume distributions obtained using LII had more annular shapes than those obtained using LOSA. Nonetheless, it has been demonstrated that the AC-LII technique is very well suited for application in media where attenuation of the excitation laser pulse energy can exceed 45%. This paper also underlines the importance of correcting LOSA SVF measurements for light scattering in high pressure flames. PACS 07-60.-j; 47.70.Pq; 65.80.+n; 78.67.-n     

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.     

2D aggregate sizing by combining laser-induced incandescence (LII) and elastic light scattering (ELS)

The combination of laser-induced incandescence and elastic light scattering has been further developed to allow for a quantitative two-dimensional determination of characteristic properties of soot aggregates, namely radius of gyration R _g and number N _p of primary particles per aggregate. In demonstrating the principle of the method, we have in a first approach approximated the particle ensemble as monodisperse and used a structure factor with an exponential cut-off function. Nonetheless, experiments performed on a laminar premixed ethene flame demonstrate basically good agreement with observations from literature and data from electron microscopy on thermophoretically obtained samples.     

Relationship between soot volume fraction and LII signal in AC-LII: effect of primary soot particle diameter polydispersity

Theoretical analysis and numerical calculations were conducted to investigate the relationship between soot volume fraction and laser-induced incandescence (LII) signal within the context of the auto-compensating LII technique. The emphasis of this study lies in the effect of primary soot particle diameter polydispersity. The LII model was solved for a wide range of primary soot particle diameters from 2 to 80 nm. For a log-normally distributed soot particle ensemble encountered in a typical laminar diffusion flame at atmospheric pressure, the LII signals at 400 and 780 nm were calculated. To quantify the effects of sublimation and differential conduction cooling on the determined soot volume fraction in auto-compensating LII, two new quantities were introduced and demonstrated to be useful in LII study: an emission intensity distribution function and a scaled soot volume fraction. When the laser fluence is sufficiently low to avoid soot mass loss due to sublimation, accurate soot volume fraction can be obtained as long as the LII signals are detected within the first 200 ns after the onset of the laser pulse. When the laser fluence is in the high fluence regime to induce significant sublimation, however, the LII signals should be detected as early as possible even before the laser pulse reaches its peak when the laser fluence is sufficiently high. The analysis method is shown to be useful to provide guidance for soot volume fraction measurements using the auto-compensating LII technique.     

Laser-induced incandescence for soot particle size and volume fraction measurements using on-line extinction calibration

A novel technique for two-dimensional measurements of soot volume fraction and particle size has been developed. It is based on a combined measurement of extinction and laser-induced incandescence using Nd:YAG laser wavelengths of 532 nm and 1064 nm. A low-energy laser pulse at 532 nm was used for extinction measurements and was followed by a more intense pulse at 1064 nm, delayed by 15 ns, for LII measurements. The 532-nm beam was split into a signal beam passing the flame and a reference beam, both of which were directed to a dye cell. The resulting fluorescence signals, from which the extinction was deduced, together with the LII signal, were registered on a single CCD detector. Thus the two-dimensional LII image could be converted to a soot volume fraction map through a calibration procedure during the same laser shot. The soot particle sizes were evaluated from the ratio of the temporal LII signals at two gate time positions. The uncertainty in the particle sizing arose mainly from the low signal for small particles at long gate times and the uncertainty in the flame temperature. The technique was applied to a well-characterized premixed flat flame, the soot properties of which had been previously thoroughly investigated. Received: 21 June 2000 / Revised version: 11 September 2000 / Published online: 7 February 2001     

Characterization of Nano‐Particles Using Laser‐Induced Incandescence

Laser‐induced incandescence (LII) is introduced as a valuable tool for the characterization of nanoparticles. This optical measurement technique is based on the heating of the particles by a short laser pulse and the subsequent detection of the thermal radiation. It has been applied successfully for the investigation of soot in different fields of application, which is described here in the form of an overview with a focus on work done at the LTT‐Erlangen during the last 10 years. In laboratory flames the soot primary particle size, volume concentration, and relative aggregate size have been determined in combination with the number density of primary particles. Furthermore, the primary particle sizes of carbon blacks have been measured in situ and online under laboratory conditions and also in production reactors. Measurements with different types of commercially available carbon black powders, which were dispersed in a measurement chamber yielded a good correlation between LII results and the specified product properties. Particle diameters determined by LII in a furnace black reactor correlate very well with the CTAB‐absorption number, which is a measure for the specific surface area. It turned out that the LII method is not affected by variations of the aggregate structure of the investigated carbon blacks. The LII signal also contains information on the primary particle size distribution, which can be reconstructed by the evaluation of the signal decay time at, at least, two different time intervals. Additionally, soot mass concentrations have been determined inside diesel engines and online measurements were performed in the exhaust gas of such engines for various engine conditions simultaneously providing information about primary particle size, soot volume, and number concentration. The LII results exhibit good correlation with traditional measurement techniques, e.g., filter smoke number measurements. In addition to the soot measurements, primarily tests with other nanoparticles like TiO₂ or metal particles are encouraging regarding the applicability of the technique for the characterization of such different types of nanoparticles.     

Influence of polydisperse distributions of both primary particle and aggregate size on soot temperature in low-fluence LII

An improved aggregate-based low-fluence laser-induced incandescence (LII) model has been developed. The shielding effect in heat conduction between aggregated soot particles and the surrounding gas was modeled using the concept of the equivalent heat transfer sphere. The diameter of such an equivalent sphere was determined from direct simulation Monte Carlo calculations in the free molecular regime as functions of the aggregate size and the thermal accommodation coefficient of soot. Both the primary soot particle diameter and the aggregate size distributions are assumed to be lognormal. The effective temperature of a soot particle ensemble containing different primary particle diameters and aggregate sizes in the laser probe volume was calculated based on the ratio of the total thermal radiation intensities of soot particles at 400 and 780 nm to simulate the experimentally measured soot particle temperature using two-color optical pyrometry. The effect of primary particle diameter polydispersity is in general important and should be considered. The effect of aggregate size polydispersity is relatively unimportant when the heat conduction between the primary particles and the surrounding gas takes place in the free-molecular regime; however, it starts to become important when the heat conduction process occurs in the near transition regime. The model developed in this study was also applied to the re-determination of the thermal accommodation coefficient of soot in an atmospheric pressure laminar ethylene diffusion flame. PACS 44.05.+e; 61.46.Df; 65.80.+n     

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.     

Numerical investigation of the effect of signal trapping on soot measurements using LII in laminar coflow diffusion flames

Laser-induced incandescence has been rapidly developed into a powerful diagnostic technique for measurements of soot in many applications. The incandescence intensity generated by laser-heated soot particles at the measurement location suffers the signal trapping effect caused by absorption and scattering by soot particles present between the measurement location and the detector. The signal trapping effect was numerically investigated in soot measurements using both a 2D LII setup and the corresponding point LII setup at detection wavelengths of 400 and 780 nm in a laminar coflow ethylene/air flame. The radiative properties of aggregated soot particles were calculated using the Rayleigh–Debye–Gans polydisperse fractal aggregate theory. The radiative transfer equation in emitting, absorbing, and scattering media was solved using the discrete-ordinates method. The radiation intensity along an arbitrary direction was obtained using the infinitely small weight technique. The contribution of scattering to signal trapping was found to be negligible in atmospheric laminar diffusion flames. When uncorrected LII intensities are used to determine soot particle temperature and the soot volume fraction, the errors are smaller in 2D LII setup where soot particles are excited by a laser sheet. The simple Beer–Lambert exponential attenuation relationship holds in LII applications to axisymmetric flames as long as the effective extinction coefficient is adequately defined.     

Effect of aggregation on the absorption cross-section of fractal soot aggregates and its impact on LII modelling

This study concerns the effect of particle aggregation on laser heating rate of soot aggregates in laser-induced incandescence. Three aggregate absorption models were investigated: the Rayleigh-Debye-Gans approximation, the electrostatics approximation, and the numerically exact generalized multi-sphere Mie-solution method. Fractal aggregates containing 5-893 primary particles of 30 nm in diameter were generated numerically using a combined particle-cluster and cluster-cluster aggregation algorithm with specific fractal parameters typical of soot. The primary particle size parameters considered are 0.089, 0.177, and 0.354. The Rayleigh-Debye-Gans approximation neglects the effect of particle aggregation on absorption; so it underestimates the aggregate absorption cross-section area by approximately 10%, depending on the aggregate size and primary particle size parameter. The electrostatics approximation is somewhat better than the Rayleigh-Debye-Gans approximation, but cannot account for the effect of primary particle size parameter. The aggregate absorption submodel affects the calculated soot temperature in laser-induced incandescence mainly in the low laser fluence regime. At high laser fluences, the effect diminishes due to the enhanced importance of soot sublimation cooling and neglect of aggregation effect in the sublimation in the present numerical model of laser-induced incandescence.     

Measurements of Silica Aggregate Particle Growth Using Light Scattering and Thermophoretic Sampling in a Coflow Diffusion Flame

The evolution of silica aggregate particles in a coflow diffusion flame has been studied experimentally using light scattering and thermophoretic sampling techniques. An attempt has been made to calculate the aggregate number density and volume fraction using the measurements of scattering cross section from 90° light scattering with combination of measuring the particle size and morphology from the localized sampling and a TEM image analysis. Aggregate or particle number densities and volume fractions were calculated using Rayleigh–Debye–Gans and Mie theory for fractal aggregates and spherical particles, respectively. Using this technique, the effects of H₂ flow rates on the evolution of silica aggregate particles have been studied in a coflow diffusion flame burner. As the flow rate of H₂ increases, the primary particle diameters of silica aggregates have been first decreased, but, further increase of H₂ flow rate causes the diameter of primary particles to increase and for sufficiently larger flow rates, the fractal aggregates finally become spherical particles. For the cases of high flame temperatures, the particle sizes become larger and the number densities decrease by coagulation as the particles move up within the flame. For cases of low flame temperatures, the primary particle diameters of aggregates vary a little following the centerline of burner and for the case of the lowest flame temperature in the present experiments, the sizes of primary particles even decrease as particles move upward.     

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

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

Effects of soot absorption and scattering on LII intensities in laminar coflow diffusion flames

Absorption and scattering of laser-induced incandescence (LII) intensities by soot particles present between the measurement volume and the detector were numerically investigated at detection wavelengths of 400 and 780 nm in a laminar coflow ethylene/air flame. The radiative properties of aggregated soot particles were calculated using the Rayleigh-Debye-Gans polydisperse fractal aggregate theory. The radiative transfer equation in emitting, absorbing, and scattering media was solved using the discrete-ordinates method. The radiation intensity along an arbitrary direction was obtained using the infinitely small weight technique. The effects of absorption and scattering on LII intensities are found to be significant under the conditions of this study, especially at the shorter detection wavelength and when the soot volume fraction is higher. Such a wavelength-dependent signal-trapping effect leads to a lower soot particle temperature estimated from the ratio of uncorrected LII intensities at the two detection wavelengths. The corresponding soot volume fraction derived from the absolute LII intensity technique is overestimated. The Beer-Lambert relationship can be used to describe radiation attenuation in absorbing and scattering media with good accuracy provided the effective extinction coefficient is adequately.     

Single-shot measurement of soot aggregate sizes by wide-angle light scattering (WALS)

The wide-angle light scattering (WALS) approach has been utilized for the measurement of soot aggregate sizes (radii of gyration) in flames on a single-shot basis. Key elements are a pulsed laser and an ellipsoidal mirror, which images the light scattered within a plane onto an intensified CCD camera, thus allowing for an instantaneous acquisition of a full scattering diagram with high resolution. Results for a laminar premixed flame exhibit good agreement with averaged data and demonstrate the feasibility of the method. The applicability of the technique to unsteady combustion processes is demonstrated by measuring aggregate sizes in a weakly turbulent jet-diffusion flame. In both cases light scattering results are verified by data obtained from electron microscopy analysis of sampled soot.     

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.