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.     

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

Soot particle size distribution measurements in a turbulent ethylene swirl flame

There is a need to better understand particle size distributions (PSDs) from turbulent flames from a theoretical, practical and even regulatory perspective. Experiments were conducted on a sooting turbulent non-premixed swirled ethylene flame with secondary (dilution) air injection to investigate exhaust and in-burner PSDs measured with a Scanning Mobility Particle Sizer (SMPS) and soot volume fractions (f_v) using extinction measurements. The focus was to understand the effect of systematically changing the amount and location of dilution air injection on the PSDs and f_v inside the burner and at the exhaust. The PSDs were also compared with planar Laser Induced Incandescence (LII) calibrated against the average f_v. LII provides some supplemental information on the relative soot amounts and spatial distribution among the various flow conditions that helps interpret the results. For the flame with no air dilution, f_v drops gradually along the centreline of the burner towards the exhaust and the PSD shows a shift from larger particles to smaller. However, with dilution air f_v reduces sharply where the dilution jets meet the burner axis. Downstream of the dilution jets f_v reduces gradually and the PSDs remain unchanged until the exhaust. At the exhaust, the flame with no air dilution shows significantly more particles with an f_v one to two orders of magnitude greater compared to the Cases with dilution. This dataset provides insights into soot spatial and particle size distributions within turbulent flames of relevance to gas turbine combustion with differing dilution parameters and the effect dilution has on the particle size. Additionally, this work measures f_v using both ex situ and in situ techniques, and highlights the difficulties associated with comparing results across the two. The results are useful for validating advanced models for turbulent combustion.     

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.     

Transported JPDF modelling and measurements of soot at elevated pressures

Accurate measurements and modelling of soot formation in turbulent flames at elevated pressures form a crucial step towards design methods that can support the development of practical combustion devices. A mass and number density preserving sectional model is here combined with a transported joint-scalar probability density function (JDPF) method that enables a fully coupled scalar space of soot, gas-phase species and enthalpy. The approach is extended to the KAUST turbulent non-premixed ethylene-nitrogen flames at pressures from 1 to 5 bar via an updated global bimolecular (second order) nucleation step from acetylene to pyrene. The latter accounts for pressure-induced density effects with the rate fitted using comparisons with full detailed chemistry up to 20 bar pressure and with experimental data from a WSR/PFR configuration and laminar premixed flames. Soot surface growth is treated via a PAH analogy and soot oxidation is considered via O, OH and O₂ using a Hertz-Knudsen approach. The impact of differential diffusion between soot and gas-phase particles is included by a gradual decline of diffusivity among soot sections. Comparisons with normalised experimental OH-PLIF and PAH-PLIF signals suggest good predictions of the evolution of the flame structure. Good agreement was also found for predicted soot volume statistics at all pressures. The importance of differential diffusion between soot and gas-phase species intensifies with pressure with the impact on PSDs more evident for larger particles which tend to be transported towards the fuel rich centreline leading to reduced soot oxidation.     

Modelling of aromatics and soot formation from large fuel molecules

There is a need for prediction models of soot particles and polycyclic aromatic hydrocarbons (PAHs) formation in parametric conditions prevailing in automotive engines: large fuel molecules and high pressure. A detailed kinetic mechanism able to predict the formation of benzene and PAHs up to four rings from C₂ fuels, recently complemented by consumption reactions of decane, was extended in this work to heptane and iso-octane oxidation. Species concentrations measured in rich, premixed flat flames and in a jet stirred reactor (JSR) were used to check the ability of the mechanism to accurately predict the formation of C₂ and C₃ intermediates and benzene at pressures ranging from 0.1 to 2.0 MPa. Pathways analyses show that propargyl recombination is the only significant route to benzene in rich heptane and iso-octane flames. When included as the first step of a soot particle formation model, the gas-phase kinetic mechanism predicts very accurately the final soot volume fraction measured in a rich decane flame at 0.1 MPa and in rich ethylene flames at 1.0 and 2.0 MPa.     

Numerical simulations of soot aggregation in premixed laminar flames

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

Nano organic carbon and soot in turbulent non-premixed ethylene flames

Spectral optical techniques are combined to characterise the distribution of large-molecule soot precursors, nanoparticles of organic carbon, and soot in two turbulent non-premixed ethylene flames with differing residence times. Laser-induced fluorescence, laser-induced incandescence and light scattering are used to define distributions across the particle size distribution. From the scattering and laser-induced emission measurements it appears that two classes of particles are formed. The first ones are preferentially formed in the fuel-rich region of the flame closer to the nozzle, have sizes of the order of few nanometers but are not fully solid particles, because the constituent molecules still maintain their individual identity exhibiting strong broadband fluorescence in the UV. The second class of particles constituted by solid particles, with sizes of the order of tens of nanometers are able to absorb a sufficient number of photons to be heated to incandescent temperatures. These larger particles are formed at larger residence times in the flame since they are the result of slow growth processes such as coagulation or carbonization. The flames are also modeled in order to produce mixture fraction maps. A new discovery is that nanoparticles of organic carbon concentration, unlike soot, does correlate well with mixture fraction, independent of position in the flame. This is likely to be a significant benefit to future modelling of soot inception processes in turbulent non-premixed flames.     

Experimental and kinetic modeling study of sooting atmospheric-pressure cyclohexane flame

Experimental data and modelling results of the main products and intermediates from a fuel-rich sooting premixed cyclohexane flame were presented in this work. Model predictions well agree with experimental data both in sooting and non-sooting flames. Major and minor species are properly predicted, together with the soot yield. The initial benzene peak was demonstrated to be due to the fast dehydrogenation reactions of the cycloalkane, which gives rise to cyclohexene and cyclohexadiene both via molecular and radical pathways. Once formed cyclohexadiene quickly forms benzene whereas in the postflame zone, benzene comes from the recombination and addition reactions of small radicals, with C₃H₃ + C₃H₃ playing the most important role in these conditions. An earlier soot inception was detected in the cyclohexane flame with respect to a n-hexane flame and this feature is not reproduced by the model that foresees soot formation significant only in the second part of the flame. The model insensitivity of soot to the reactant hydrocarbon was also observed comparing the predictions of three flames of cyclohexane, 1-hexene and n-hexane with the same temperature profile. A sensitivity analysis revealed that soot primarily comes from the HACA mechanism for the three flames, acetylene being the key species in the nucleation. Experimental data on soot inception seem to indicate the importance of the early formation of benzene, that depends on the fuel structure. It is thus important to further investigate the role of benzene and aromatics in order to explain this discrepancy.     

A computational framework for predicting laminar reactive flows with soot formation

Numerical modeling is an attractive option for cost-effective development of new high-efficiency, soot-free combustion devices. However, the inherent complexities of hydrocarbon combustion require that combustion models rely heavily on engineering approximations to remain computationally tractable. More efficient numerical algorithms for reacting flows are needed so that more realistic physics models can be used to provide quantitative soot predictions. A new, highly-scalable combustion modeling tool has been developed specifically for use on large multiprocessor computer architectures. The tool is capable of capturing complex processes such as detailed chemistry, molecular transport, radiation, and soot formation/destruction in laminar diffusion flames. The proposed algorithm represents the current state of the art in combustion modeling, making use of a second-order accurate finite-volume scheme and a parallel adaptive mesh refinement (AMR) algorithm on body-fitted, multiblock meshes. Radiation is modeled using the discrete ordinates method (DOM) to solve the radiative transfer equation and the statistical narrow-band correlated-k (SNBCK) method to quantify gas band absorption. At present, a semi-empirical model is used to predict the nucleation, growth, and oxidation of soot particles. The framework is applied to two laminar coflow diffusion flames which were previously studied numerically and experimentally. Both a weakly-sooting methane–air flame and a heavily-sooting ethylene–air flame are considered for validation purposes. Numerical predictions for these flames are verified with published experimental results and the parallel performance of the algorithm analyzed. The effects of grid resolution and gas-phase reaction mechanism on the overall flame solutions were also assessed. Reasonable agreement with experimental measurements was obtained for both flames for predictions of flame height, temperature and soot volume fraction. Overall, the algorithm displayed excellent strong scaling performance by achieving a parallel efficiency of 70% on 384 processors. The proposed algorithm proved to be a robust, highly-scalable solution method for sooting laminar flames.     

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.     

Particle formation in opposed-flow diffusion flames of ethylene: An experimental and numerical study

An experimental and numerical study on particles inception and growth is performed in opposed-flow diffusion flames of ethylene and air characterized by different sooting tendencies. Spectrally resolved UV-visible laser induced fluorescence, laser induced incandescence and laser light scattering measurements are used to characterize different classes of combustion-generated compounds. A detailed kinetic model accounting for both gas-phase and particle formation is used. Comparison between experimental results and numerical predictions gives a qualitative view of the mechanism of particle formation in opposed-flow flames.Particle inception is the result of both chemical growth and coagulation of aromatic compounds. In the region close to the flame front where the temperature is relatively high and radicals are abundant, the particle inception is due to a chemical growth mechanism by which aromatic molecules add aromatic radicals leading to the formation of biphenyl-like structures. The growth process continues as high-molecular mass aromatics are moved away from the flame zone towards the stagnation plane by the addition of acetylene and other aromatics forming particles of increasing sizes. Graphitization of these particles and thermal annealing lead to the formation of soot particles. At relatively lower temperatures, found across the stagnation plane, particles growth still occurs and it is mainly due to a process of physical coagulation of PAHs.The experimental and numerical results obtained in this work demonstrate and explain the sensibility of inception and growth of particles to radical concentration and temperature in opposed-flow flame configurations.     

Laser-induced incandescence technique to identify soot nucleation and very small particles in low-pressure methane flames

This paper presents the study we carried out on the formation of soot particles in low-pressure premixed CH₄/O₂/N₂ flames by using Laser-Induced Incandescence (LII). Flames were stabilised at 26.6 kPa (200 torr). Four different equivalence ratios were tested (Φ = 1.95, 205, 2.15 and 2.32), Φ = 1.95 corresponding to the equivalence ratio for which LII signals begin to be measurable along the flame. The evolution of the LII signals with laser fluence (fluence curve), time (temporal decay) and emission wavelength is reported at different heights above the burner. We specifically took advantage of the low-pressure conditions to probe with a good spatial resolution the soot inception zone of the flames. Significant different behaviours of the fluence curves are observed according to the probed region of the flames and Φ. In addition, while the surface growth process is accompanied by an increase in the LII decay-times (indicator of the primary particle diameter) at higher Φ, decay-times become increasingly short at lower Φ reaching a constant value along the flame at Φ = 1.95. These behaviours are consistent with the detection of the smallest incandescent particles in the investigated flames, these particles having experienced very weak surface growth. Flame modelling including soot formation has been implemented in flames Φ = 2.05 and 2.32. Experimental quantitative soot volume fraction profiles were satisfactorily reproduced by adjusting the fraction of reactive soot surface available for reactions. The qualitative variation of the computed soot particle diameter and the relative weight of surface growth versus nucleation were consistent with the experimental observations.     

Measurement of radiative gas and particle emissions in biomass flames

Radiation is the dominant mode of heat transfer near the burner of coal and biomass-fired boilers. Predicting and measuring heat transfer is critical to the design and operation of new boiler concepts. The individual contributions of gas and particle phases are dependent on gas and particle concentration, particle size, and gas and particle temperature which vary with location relative to the flame. A method for measuring the contributions of both gas and particle radiation capable of being applied in harsh high temperature and pressure environments has been demonstrated using emission from particles and water vapor using an optical fiber probe transmitting a signal to a Fourier Transform Infrared (FTIR) spectrometer. The method was demonstrated in four environments of varying gas and particle loading using natural gas and pulverized wood flames in a down-fired 130?kW_th cylindrical reactor. The method generates a gas and particle temperature, gas concentrations (H₂O and CO₂), total gas and particle intensities, and gas and particle total effective emissivity from line-of-sight emission measurements. For the conditions measured, downstream of the luminous flame zone, water vapor and CO₂ radiation were the dominant modes of heat transfer (effective emissivity 0.13–0.19) with particles making a minor contribution (effective emissivity 0.01–0.02). Within a lean natural gas flame, soot emission was low (effective emissivity 0.02) compared to gas (0.14) but within a luminous flame of burning wood particles (500?µm mean diameter) the particles (soot and burning wood) produced a higher effective emissivity (0.17) than the gas (0.12). The measurement technique was therefore found to be effective for several types of combustion environments.     

Characterization of nanoparticles of organic carbon (NOC) produced in rich premixed flames by differential mobility analysis

Differential mobility analysis (DMA) is used to measure on-line the size distributions of inception particles in atmospheric pressure premixed ethylene air flames ranging from C/O = 0.61 to 0.69, just at the onset of soot formation. DMA is also used, in combination with electrospray, to measure the size distributions of suspended flame products captured in water samples. The DMA systems used for this work employ detectors sensitive to the full range of molecular clusters/nanoparticles in gas-to-particle conversion processes (as small as about 1 nm) and they have much larger sheath gas flow rates than is typically used to reduce losses and peak broadening by diffusion. The measured size distributions show that the first particles observed in flames have a size of 2 nm, consistent with previous in situ measurements by light scattering and extinction (LSE) and the off-line measurements of material captured in water samples from the same flames. For richer flames, the quantity of the 2 nm particles measured increases, and the width of its size distribution shifts asymmetrically toward larger sizes. A numerical coagulation model assuming size-dependent coagulation efficiency predicts well the experimentally measured size distributions in the flames examined. Similarly, the slightly larger size distributions measured by atomic force microscopy of inception particles deposited on surfaces can also be attributed to the size-dependent coagulation/adhesion efficiency. The results imply that the smaller nanoparticles formed in combustion processes have a longer lifetime than those larger than 6-7 nm and may play an important role in the formation of fine organic carbon particulate in the atmosphere.     

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.     

Calculating the soot particle size distribution function in turbulent diffusion flames using a sectional method

Soot formation in a turbulent jet diffusion flame is modeled using an unsteady flamelet approach in post-process. In the present work, we apply a detailed kinetic soot model with a sectional method, and study the evolution of the particle size distribution. Detailed information on the evolution of the soot particle size distribution function is acquired. It is found that the particle size distribution function is bimodal throughout the flame. The transition from the small to large particle size distributions is strongly influenced by surface growth and oxidation reactions. We find that large particles are most likely to be emitted from the flame.     

Similarities and dissimilarities in n-hexane and benzene sooting premixed flames

The flame structure of atmospheric-pressure sooting premixed flames of aliphatic and aromatic hydrocarbons with the same carbon atom number (hexane and benzene) were studied at similar temperatures and C/O ratios by sampling and chemical and spectroscopic analysis. The differences in the oxidation mechanism of hexane and benzene in fuel-rich conditions were found to produce a different chemical environment in the yield of light hydrocarbons and their relative compositions where soot inception occurs. The predominance of acetylene and simple aromatic reactants in the oxidation region of the benzene flame favoured the early appearance and steep rise of soot particles. Large formations of saturated and unsaturated hydrocarbons were observed in the main oxidation region of the hexane flame whereas a delayed formation of aromatics (mainly PAH) was observed at soot inception only after complete oxygen consumption. There are differences in soot inception mechanisms reflected by the soot structure from UV-vis spectral shapes and mass specific absorption coefficients. In the benzene flame, they appeared to be more ordered and aromatic with a narrower size of aromatic systems and/or more curved aromatic structures. By contrast, less ordering with a more complex aliphatic/aromatic structure and a larger variety of aromatic systems were found to characterize soot formed in the hexane flame.