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.     

An aerosol model to predict size and structure of soot particles

An aerosol model to simulate soot formation and growth was developed using moving- and fixed-sectional methods. The new model is composed of a set of subroutines that can be easily combined with the Chemkin package. Using the model, we have simulated soot formation and growth in plug flow reactors.

Our model was compared with a previously published method of moments model for a simulation of the plasma pyrolysis of methane in a plug flow reactor. Inclusion of the transition correction factor for the condensation coefficient led to the prediction of a smaller condensation rate compared with the method of moments model. The average coagulation rate calculated by the sectional model was much higher than that by the method of moments model for a broad particle size distribution. The two models predicted significantly different soot precursor concentration and rates of aerosol processes, but substantially similar particle mass and number for the pyrolysis process.

We have also simulated soot formation and growth in a jet-stirred/plug flow reactor (JSR/PFR) system for which soot size distribution measurements are available in the literature. It is shown that the adjusted-point fixed-sectional method can provide comparable accuracy to the moving-sectional model in a simulation of soot formation and growth. It is also shown that the measured surface growth rate could be much higher than the value used in this study. Soot mass concentrations and size distributions for particles larger than 10 nm were well predicted with a surface reaction enhancement. The primary particle size was underpredicted by only about 30% compared with the measurements, without any model adjustments. As the new model can predict both the particle size distribution and structure, and is suitable for application in complex flows, its application to diverse soot formation conditions will enhance our knowledge on the evolution of soot structures.     

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

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.     

In-situ monitoring of UF5 nanoparticle growth in a plug-flow reactor by a low-pressure differential mobility analyzer

The dependence of UF₅ nanoparticle size distribution on growth time was investigated for various initial concentrations of UF₅ monomers. UF₅ nanoparticles prepared by photodissociation of UF₆ in a mixed-flow reactor were allowed to pass through a plug-flow reactor (PFR), and their size distribution was measured by an in-situ size-monitoring system, namely a low-pressure differential mobility analyzer. By changing the length of the PFR, the growth time in the PFR was controlled from 0 to 30 s. An analysis using Lee’s coagulation model revealed that the growth mechanism of the UF₅ nanoparticles involved a Brownian coagulation process and that the sticking probability of the nanoparticles was approximately unity. Received: 27 August 2001 / Accepted: 6 September 2001 / Published online: 20 December 2001     

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.     

Two-color laser-induced incandescence and cavity ring-down spectroscopy for sensitive and quantitative imaging of soot and PAHs in flames

Laser-induced incandescence is a technique which enables the measurement of soot volume fractions. However, the laser-induced soot emission might be affected by a fluorescence background generally ascribed to the polycyclic aromatic hydrocarbon compounds (PAHs) present at the soot location. In this paper, spatially resolved distributions of PAH absorbance and soot are obtained in sooting diffusion flames. The original method developed here consists in comparing the emission distributions induced by two different laser wavelengths: (1) at 1064 nm emission signals are exempt from PAH fluorescence and (2) at 532 nm both soot incandescence and PAH emission contribute to the total signal. In addition, the absolute absorption coefficient of the PAH mixture is determined by comparing absorption measurements obtained by cavity ring-down spectroscopy (CRDS) at 1064 nm and 532 nm. The proposed method can provide highly sensitive 2D imaging of PAHs and soot using the fundamental and the second-harmonic frequencies of a single YAG laser. Finally, 2D distributions of PAH absorbance and soot volume fraction calibrated by CRDS are obtained in two diffusion flames, particularly in a very low-sooting flame exhibiting a maximum PAH absorbance of 6×10^-4 cm^-1 and a maximum soot volume fraction of 3 ppb only. The respective spatial distributions of PAHs and soot are shown to vary with the initial C/O ratio. PACS 33.20.Lg; 42.62.Fi; 44.40.+a     

燃油分级多点喷射低污染燃烧室的化学反应网络模型分析

本文采用基于详细化学反应机理的化学反应网络模型分析了航空发动机燃油径向分级多点喷射低污染燃烧室的NO_x排放特性。该分级燃烧室不同于传统燃烧室,头部由值班区和主燃区两个不同的燃烧区域,根据CFD得到的流场特性和当量比的分布特性对燃烧室进行分区构建化学反应器网络模型,研究了值班级当量比以及值班级和主燃级两级供油比例对排放的影响。同时,还分析了空气进口温度对NO_x排放的影响。得到了较为合理的变化趋势,为低污染燃烧室的初步设计提供了有益的指导。     

Electric-arc high-capacity reactor for the synthesis of carbon soot with a high content of endohedral metallofullerenes

An electric-arc reactor is designed for synthesizing carbon soot containing endohedral metallofullerenes during the sequential evaporation of five composite graphite electrodes. The preparation conditions of the composite graphite electrodes and the electric-arc evaporation parameters are optimized, which increases the reactor capacity and the content of endohedral metallofullerenes in carbon soot.     

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     

Fragmentation of pulverized coal in a laminar drop tube reactor: Experiments and model

Fragmentation during pulverized coal particles conversion shifts the particle size distribution of the fuel towards smaller particle sizes, affecting both conversion rates and heat release. After pyrolysis of a high volatiles Colombian coal in CO₂ atmosphere in a drop tube reactor at 1573?K, solid carbonaceous particles of different size, from 100?µm of the particle feed down to the nanometric size, have been observed. A fragmentation model has been used to predict the fate of Colombian coal particles under the experimental conditions of the drop tube experiment and predict the particle size distribution (PSD). Model and experimental results are in very good agreement and indicate that in the DTR experiment the coal underwent almost complete pyrolysis and that fragmentation generated a 36?wt% population of particles with size close to 30?µm. The close match between the PSDs obtained from experiments and from the fragmentation model is an important novelty. It demonstrates that fragmentation occurs not only under fluidized bed conditions but also under the conditions of pulverized coal combustion. Experimentalists are warned against the fact that the fine particulate sampled at the outlet of laminar flow reactors and boilers is not always composed of soot only. Char fragments can be misidentified as soot. The implementation of fragmentation submodels in pulverized fuel combustion and gasification codes is highly recommended.     

Characterization of high heating rate chars of biomass fuels

Data on biomass chars obtained under conditions similar to those of practical applications (high heating rate and low residence time) are required for co-combustion and gasification plants. A methodological procedure is developed and applied to two biomass fuels (cacao shells and olive cake) for producing high heating rate chars and characterizing their reactivity and morphology after the first steps of devolatilization. Different chars are produced in a drop tube reactor (rapid pyrolysis) by varying the nominal temperature and the residence time. Oxidation in air is performed to compare typical temperatures and kinetic parameters and evaluate the effect of the operating conditions on char reactivity. A detailed SEM analysis allows to assess the structural variations during the pyrolysis and detect the main phenomena (softening, swelling, melting, formation of bubbles). A quantitative morphological study is also performed to provide size and shape (important for biomasses) distributions of the parent fuel and the chars. These data are more significant than average values in advanced model to correctly simulate the fluid dynamic behaviour of each dimensional class of particles in large scale furnaces and gasifiers and predict a more reliable residence time of the particles.     

The impact of residence time on ignitability and time to ignition in a toroidal jet-stirred reactor

Understanding of ignition processes is central to design for reliable and safe aerospace combustor systems. Ignition is influenced by many factors including combustor geometry, flow conditions, fuel composition, turbulence intensity, ignition source, and energy deposition method. A toroidal jet-stirred reactor (TJSR) utilizes bulk fluid motion, presence of recirculation zones, a bulk residence time, and turbulence intensities which emulate characteristics relevant to cavity stabilized and swirl stabilized combustors. In this work, a TJSR was used to quantify ignitability and time-to-ignition of premixed ethylene and air. The effects of inlet temperature, residence time, and reactivity were studied on forced ignition processes. Experimental conditions ranged from residence times of 15–35?ms, mixture temperatures of 340–450?K, and equivalence ratios of 0.5–1 using capacitive spark-discharge ignition. The minimum equivalence ratio for ignition (MER), or the equivalence ratio at 50% probability, shows an inverse relationship with mixture temperature and residence time. Prior theory of real engine combustor performance for lean light off, proposed by Ballal and Lefebvre, was compared to the MER and displayed similar trends to the model. Spatially integrated OH* chemiluminescence was used to measure time to ignition within the reactor. Reduction in ignitibility was experienced as the time-to-ignition approached the residence time stressing the importance of device flow time scales in relation to kernel growth dynamics and ignition probability.     

Time-resolved two-color LII: size distributions of nano-particles from gas-to-particle synthesis

A two-color LII technique for in situ measurements of particle size distributions is described. The technique is based on the simultaneous detection of time-resolved LII signals at two different wavelengths with one-dimensional spatial resolution using a newly developed experimental setup. The ratio of both LII signals yields particle temperatures as a function of time and location. Measured particle temperature decays are numerically simulated based on a detailed cooling model for particle ensembles. Particle size distributions are obtained by fitting simulated particle temperature decays to measured ones using multi-dimensional non-linear regression. The two-color LII technique for particle sizing can be applied to a wide range of materials because it is independent of the optical properties of the particle material. Exemplarily, the measuring technique is applied to investigate the synthesis of nanoscaled metal oxide particle in a laser vaporization reactor.     

Modelling biodiesel–diesel spray combustion using multicomponent vaporization coupled with detailed fuel chemistry and soot models

A multicomponent vaporization model is integrated with detailed fuel chemistry and soot models for simulating biodiesel–diesel spray combustion. Biodiesel, a fuel mixture comprised of fatty-acid methyl esters, is an attractive alternative to diesel fuel for use in compression-ignition engines. Accurately modelling of the spray, vaporization, and combustion of the fuel mixture is critical to predicting engine performance using biodiesel. In this study, a discrete-component vaporization model was developed to simulate the vaporization of biodiesel drops. The model can predict differences in the vaporization rates of different fuel components. The model was validated by use of experimental data of the measured biodiesel drop size history and spray penetration data obtained from a constant-volume chamber. Gas phase chemical reactions were simulated using a detailed reaction mechanism that also includes PAH reactions leading to the production of soot precursors. A phenomenological multi-step soot model was utilized to predict soot emissions from biodiesel–diesel combustion. The soot model considered various steps of soot formation and destruction, such as soot inception, surface growth, coagulation, and PAH condensation, as well as oxidation by oxygen and hydroxyl-containing molecules. The overall numerical model was validated with experimental data on flame structure and soot distributions obtained from a constant-volume chamber. The model was also applied to predict combustion, soot and NOx emissions from a diesel engine using different biodiesel–diesel blends. The engine simulation results were further analysed to determine the soot emissions characteristics by use of biodiesel–diesel fuels.     

Experimental and computational study of diethyl sulfide pyrolysis and mechanism

The pyrolysis of diethyl sulfide (C₂H₅SC₂H₅), a simulant for mustard chemical warfare agents, was studied in a turbulent flow reactor with extractive gas composition analysis by GC/MS and FTIR. Experiments were performed at approximately atmospheric pressure for four different temperatures between 630 and 740 °C with maximum residence times between 0.06 and 0.08 s. Temperature and species profiles were obtained on the centerline of the reactor. The mixing characteristics in the reactor were determined by using carbon monoxide as a tracer. 80% destruction of diethyl sulfide was observed for the experiment at the temperature of 740 °C and the residence time of 0.06 s. The following species were quantified: diethyl sulfide, ethylene, methane, ethane, acetylene, carbon disulfide, and thiophene. In addition, ethanethiol, methyl thiirane CH₃-(Cy-CH-CH₂-S), ethyl methyl disulfide, and diethyl disulfide were identified but not quantified. A light yellow solid containing sulfur condensed in sampling probes. Thermochemical properties for all species and a detailed mechanism were developed for modeling the reaction system. Thermodynamic and kinetic parameters were based on density functional theory and ab initio calculations using isodesmic work reactions for enthalpies. Kinetic parameters for chemical activation and unimolecular dissociation reactions were determined with multi frequency quantum RRK analysis for k(E) and master equation for fall-off. Important reactions were identified by sensitivity analysis and reaction pathway analysis of the mechanism. Model predictions show overall good agreement with experiment.     

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.     

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.     

Unraveling the carbene chemistry of oxymethylene ethers: Experimental investigation and kinetic modeling of the high-temperature pyrolysis of OME-2

Oxymethylene ethers (OMEs) form an interesting family of synthetic compounds to replace fossil fuels. This alternative liquid energy carrier can contribute to a circular carbon economy when produced via carbon capture and utilization technology using renewable electricity. Despite the potential to reduce greenhouse gas and particulate matter emissions and their ideal ignition characteristics, little is known about the thermal decomposition behavior of OMEs. In this work, new insights are obtained in the pyrolysis chemistry of oxymethylene ether-2 (OME-2) and the role of carbenes by performing experiments at high temperatures (> 850 K) in a tubular quartz reactor. The used continuous bench-scale pyrolysis unit has a dedicated on-line analysis section including comprehensive two-dimensional gas chromatography (GC $\times$ GC) coupled with flame ionization detection (FID) and mass spectroscopy (MS) to identify and quantify the full product spectrum over the complete temperature range. The reactor temperature was varied between 850 and 1150 K at a fixed pressure of 0.15 MPa and residence times of 400 to 850 ms. The major products are dimethoxymethane, formaldehyde, methyl formate, methane, CO₂, CO and H₂. Minor intermediate compounds comprise dimethyl ether, formic anhydride, formic acid, methoxymethyl formate and methoxymethanol. The yields of compounds with carbon-carbon bonds are low since no such bonds originally occur in OME-2. Precursors of aromatic compounds and soot particles are absent in the reactor effluent. The experimental results are simulated with a new first principles-based kinetic model for pyrolysis and combustion of OME-2. This model can predict the experimental trends of major products on average within the experimental uncertainty margin of ± 10% relative for major product species. A reaction pathway and sensitivity analysis are presented to highlight the importance of the carbenes for initiation of the radical chemistry under pyrolysis conditions.