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.     

Understanding in-cylinder soot reduction in the use of high pressure fuel injection in a small-bore diesel engine

This study shows how soot particles inside the cylinder of the engine are reduced due to high pressure fuel injection used in a light-duty single-cylinder optical diesel engine fuelled with methyl decanoate, a selected surrogate fuel for the diagnostics. For various injection pressures, planar laser induced incandescence (PLII) imaging and planar laser-induced fluorescence of hydroxyl (OH-PLIF) imaging were performed to understand the temporal and spatial development of soot and high-temperature flames. In addition, a thermophoresis-based particle sampling technique was used to obtain transmission electron microscope (TEM) images of soot aggregates and primary particles for detailed morphology analysis. The OH-PLIF images suggest that an increase in the injection pressure leads to wider distribution of high-temperature flames likely due to better mixing. The enhanced high-temperature reaction can promote soot formation evidenced by both a faster increase of LII signals and larger soot aggregates on the TEM images. However, the increased OH radicals at higher injection pressure accelerates the soot oxidation as shown in a higher decreasing rate of LII signals as well as dramatic reduction of the sampled soot aggregates at later crank angles. The analysis of nanoscale carbon layer fringe structures also shows a consistent trend that, at higher injection pressure, the soot particles are more oxidized to form more graphitic carbon layer structures. Therefore, it is concluded that the in-cylinder soot reduction at higher injection pressure conditions is due to enhanced soot oxidation despite increased soot formation.     

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.     

Soot particles in piston-top pool fires and exhaust at 5 and 15 MPa injection pressure in a gasoline direct-injection engine

This study shows the structure of soot particles sampled directly from wall wetting-induced pool fires formed on the piston top in a spark-ignition direct-injection (SIDI) engine. Of particular interest is its variation with injection pressure considering the current trend of high-pressure DI system development to reduce engine-out particulate emissions. Thermophoretic particle sampling was performed for transmission electron microscope (TEM) imaging, which was post-processed for statistical analysis of key morphology and internal structure parameters. These include the size distribution of soot aggregates and primary particles as well as carbon-layer fringe-to-fringe gap and concentricity. With the fixed engine speed and load conditions, in-flame soot particles are compared to the exhaust particles sampled simultaneously at selected 5 and 15 MPa injection pressures corresponding to low-speed/low-load and high-speed/high-load in practical engine operation. From the TEM images and statistical analysis, it was found that more concentrated and taller pool fire developed for 5-MPa injection leads to smaller soot aggregates composed of smaller soot primary particles due to suppressed soot growth in fuel-rich flames. However, the soot particles formed in 15-MPa injection-induced pool fires are at a more reactive status evidenced by less defined core-shell boundaries and higher fringe separation. The higher soot reactivity results in enhanced soot oxidation, which explains smaller soot aggregates and primary particles found for the 15-MPa injection in the exhaust sample.     

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.     

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.     

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.     

Prediction of soot and thermal radiation in a model gas turbine combustor burning kerosene fuel spray at different swirl levels

Combustion of kerosene fuel spray has been numerically simulated in a laboratory scale combustor geometry to predict soot and the effects of thermal radiation at different swirl levels of primary air flow. The two-phase motion in the combustor is simulated using an Eulerian–Lagragian formulation considering the stochastic separated flow model. The Favre-averaged governing equations are solved for the gas phase with the turbulent quantities simulated by realisable k–? model. The injection of the fuel is considered through a pressure swirl atomiser and the combustion is simulated by a laminar flamelet model with detailed kinetics of kerosene combustion. Soot formation in the flame is predicted using an empirical model with the model parameters adjusted for kerosene fuel. Contributions of gas phase and soot towards thermal radiation have been considered to predict the incident heat flux on the combustor wall and fuel injector. Swirl in the primary flow significantly influences the flow and flame structures in the combustor. The stronger recirculation at high swirl draws more air into the flame region, reduces the flame length and peak flame temperature and also brings the soot laden zone closer to the inlet plane. As a result, the radiative heat flux on the peripheral wall decreases at high swirl and also shifts closer to the inlet plane. However, increased swirl increases the combustor wall temperature due to radial spreading of the flame. The high incident radiative heat flux and the high surface temperature make the fuel injector a critical item in the combustor. The injector peak temperature increases with the increase in swirl flow mainly because the flame is located closer to the inlet plane. On the other hand, a more uniform temperature distribution in the exhaust gas can be attained at the combustor exit at high swirl condition.     

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.     

Soot-velocity measurements by particle vaporization velocimetry

We demonstrate a new imaging technique for velocity measurements in particle-laden flows. The technique, particle vaporization velocimetry, is a form of flow tagging based on laser vaporization of absorbing particles at defined locations in the flow. The locations of these tagged regions are then interrogated after a known delay to determine the convective velocity. Results are presented for vaporization of carbonaceous (soot) particles in a nonreacting gas jet and a hydrocarbon flame, with interrogation provided by either elastic scattering or laser-induced incandescence from the soot. The long lifetime of the tagged soot regions (>2 ms) allows measurements to be made over a wide range of velocities.     

Detailed investigation of soot formation from jet fuel in a diffusion flame with comprehensive and hybrid chemical mechanisms

A numerical and experimental study is performed to investigate soot formation from jet fuel in a laminar coflow diffusion flame. The combustion chemistry of the fuel is simulated using (1) the MURI jet fuel surrogate (Dooley et al. 2012) with a modestly reduced Ranzi mechanism (Ranzi et al. 2012), and (2) the recently proposed HyChem model (Xu et al. 2018) combined with the KAUST PAH mechanism 2 (Wang et al. 2013). The two reaction mechanisms are coupled with a sectional soot model to simulate a coflow diffusion flame of methane doped with the MURI jet fuel surrogate. The combined laser extinction and two-angle elastic light scattering method is used to perform non-intrusive in situ measurements of soot volume fraction, primary particle diameter and number density. The good agreement including soot particle size and number density between the experimental data and the simulation results computed with the reduced Ranzi mechanism demonstrate the robustness of the soot model to changes in fuel composition, as the model parameters are unchanged with a previous numerical study of soot formation of n-propylbenzene/n-dodecane mixtures (Zhang and Thomson, 2018). The computation with the combined HyChem/KAUST mechanism predicts similar results as the computation with the detailed chemistry of the reduced Ranzi mechanism for fuel breakdown, thus the basic premise of the HyChem model that the fuel decomposition process can be greatly simplified with the lumped reaction steps is supported. The results also show that by adding a PAH growth scheme to the HyChem model, the approach can be used to predict soot formation from jet fuel combustion in a laminar coflow diffusion flame. Finally, the dependency of the soot prediction on PAH chemistry is discussed and it is suggested that more experimental data is needed to validate the PAH mechanism and improve the predictive accuracy of the model.     

LES study of diesel flame/wall interaction and mixing mechanisms at different wall distances

In this paper, the flame-wall interaction of reacting diesel spray under engine like conditions is investigated using large eddy simulations. The aim of this study is to understand the influence of the distance between the wall and the spray nozzle on the air entrainment rate, which is a key variable in formation/oxidation process of soot. Three experimental cases are investigated, a free jet case and two wall impingement cases with a distance from nozzle to wall of 30 mm and 50 mm, which are considered as characteristic wall impingement distances for light- and heavy-duty bores in diesel engines, respectively. The optical soot measurements imply a positive influence of wall on the rate of soot oxidation. Numerical simulations are employed to elucidate importance of different mechanisms for the air entrainment, i.e., air entrainment prior to flame lift-off position, enhanced mixing due to the wall impingement and enhanced mixing by the entrainment wave. The results show that oxidation process after the end of injection is driven by a different mixing mechanism depending on the distance to the wall. The 30 mm case resulted in a “mixing boost”, where the dominant mixing mechanism is the wall impingement vortex mixing, which gives rise to the fastest soot decay among the cases. The mixing in the 50 mm case is governed by a late wall impingement vortex mixing, giving rise to a low, but a constant air entrainment rate, i.e., a “mixing plateau”. The free jet case resulted in mixing governed by the entrainment wave mechanism. Both wall impingement cases have faster soot oxidation rate compared with the free jet case, but due to a different underlying mixing process. LES is shown to be able to replicate the line-of-sight measurements of natural OH* chemiluminescence and distribution of soot region from the optical soot diagnostics.     

Numerical simulation aided relative optical density analysis of TEM images for soot morphology determination

The relative optical density (ROD) method provides a means to measure three-dimensional information about soot aggregates from two-dimensional transmission electron microscopy (TEM) micrographs of soot. The method is dependent on accurate calibration of the relationship between the measured soot ROD in TEM images and the actual soot thickness perpendicular to the imaging plane. A novel calibration method based on the comparison between probability distributions of measured soot ROD in TEM images and that of virtual soot thickness of numerically simulated soot is introduced. Soot aggregates of various prescribed fractal structure parameters were numerically generated using a tunable cluster-cluster aggregation model. The probability histograms of the local soot thickness for the simulated soot aggregates and ROD of the TEM images of flame generated soot aggregates were found to be quite similar and were used as a basis to establish a quantitative relationship between ROD and the local soot thickness. The calibration constant obtained from the analysis of the simulated soot was found to be insensitive to the fractal structure parameters over a wide range. The calibrated ROD method is successfully applied to the morphology analysis of soot aggregates generated in an atmospheric laminar co-flow ethylene-air diffusion flame based on thermophoretic sampling (TS) and TEM analysis techniques. With the ROD method, an overlap coefficient is introduced to identify and eliminate non-soot-aggregate structures and the selection of a cut-off overlap coefficient was found to have little influence on the final results over a relatively wide range. ROD is independent of empirical constants and human judgments and has been found to be an accurate and reliable TEM image analysis method for studying the morphology of soot aggregates.     

Soot formation and temperature field structure in co-flow laminar methane-air diffusion flames at pressures from 10 to 60 atm

The effects of pressure on soot formation and the structure of the temperature field were studied in co-flow methane-air laminar diffusion flames over a wide pressure range, from 10 to 60 atm in a high-pressure combustion chamber. The selected fuel mass flow rate provided diffusion flames in which the soot was completely oxidized within the visible flame envelope and the flame was stable at all pressures considered. The spatially resolved soot volume fraction and soot temperature were measured by spectral soot emission as a function of pressure. The visible (luminous) flame height remained almost unchanged from 10 to 100 atm. Peak soot concentrations showed a strong dependence on pressure at relatively lower pressures; but this dependence got weaker as the pressure is increased. The maximum conversion of the fuel’s carbon to soot, 12.6%, was observed at 60 atm at approximately the mid-height of the flame. Radial temperature gradients within the flame increased with pressure and decreased with flame height above the burner rim. Higher radial temperature gradients near the burner exit at higher pressures mean that the thermal diffusion from the hot regions of the flame towards the flame centerline is enhanced. This leads to higher fuel pyrolysis rates causing accelerated soot nucleation and growth as the pressure increases.     

Soot studies of laminar diffusion flames with recirculation zones

This paper describes the unusual sooting structure of three flames established by the laminar recirculation zones of a centerbody burner. The vertically mounted burner consists of an annular air jet and a central fuel jet separated by a bluff-body. The three ethylene fueled flames are identified as: fully sooting, donut-shape, and ring-shape sooting flames. Different shapes of the soot structures are obtained by varying the N₂ dilution in the fuel and air jets while maintaining a constant air and fuel velocity of 1.2 m/s. All three flames have the unusual characteristic that the soot, entrained into the recirculation zone, follows discrete spiral trajectories that terminate at the center of the vortex. The questions are what cause: (1) the unusual sooting structures and (2) the spiral trajectories of the soot? Flame photographs, laser sheet visualizations, and calculations with a 2D CFD-based code (UNICORN) are used to answer these questions. The different sooting structures are related to the spiral transport of the soot, the spatial location of the stoichiometric flame surface with respect to the vortex center, and the burnout of the soot particles. Computations indicate that the spiral trajectories of the soot particles are due to thermophoresis.     

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

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

Single-walled carbon nanotube formation on iron oxide catalysts in diffusion flames

Single-walled carbon nanotubes (SWCNTs) are shown to grow rapidly on iron oxide catalysts on the fuel side of an inverse ethylene diffusion flame. The pathway of carbon in the flame is controlled by the flame structure, leading to formation of SWCNTs free of polycyclic aromatic hydrocarbons (PAH) or soot. By using a combination of oxygen-enrichment and fuel dilution, fuel oxidation is favored over pyrolysis, PAH growth, and subsequent soot formation. The inverse configuration of the flame prevents burnout of the SWCNTs while providing a long carbon-rich region for nanotube formation. Furthermore, flame structure is used to control oxidation of the catalyst particles. Iron sub-oxide catalysts are highly active toward SWCNT formation while Fe and Fe₂O₃ catalysts are less active. This can be understood by considering the effects of particle oxidation on the dissociative adsorption of gas-phase hydrocarbons. The optimum catalyst particle composition and flame conditions were determined in near real-time using a scanning mobility particle sizer (SMPS) to measure the catalyst and SWCNT size distributions. In addition, SMPS results were combined with flame velocity measurement to measure SWCNT growth rates. SWCNTs were found to grow at rates of over 100 μm/s.     

Prediction of soot formation characteristics in a pulverized-coal combustion field by large eddy simulations with the TDP model

In this study, the soot formation characteristics in a pulverized-coal combustion field formed by a 4 kW Central Research Institute of Electric Power Industry (CRIEPI) jet burner were predicted by large eddy simulation (LES) employing a tabulated-devolatilization-process model (TDP model) [N. Hashimoto et al., Combust. Flame 159 (2012) 353–366]. This model enables to take into account the effect of coal particle heating rate on coal pyrolysis. The coal-derived soot formation model proposed by Brown and Fletcher [A. L. Brown and T. H. Fletcher, Energy Fuels 12 (1998) 745–757] was employed in the LES. A comparison between the data predicted by LES and the soot volume fraction distribution data measured by laser induced incandescence confirmed that the soot formation characteristics in the coal combustion field of the CRIEPI burner can be accurately predicted by LES. A detailed analysis of the data predicted by LES showed that the soot particle distribution in this burner is narrow because the net soot formation rate is negative on both sides of the base of the soot volume fraction. At these positions, soot particles diffused from the peak position of soot volume fraction are oxidized due to a relatively high oxygen concentration. Finally, the effect of soot radiation on the predicted gas temperature distribution was examined by comparing the simulation results obtained with and without soot radiation. This comparison showed that the maximum gas temperature predicted by the simulation performed with soot radiation was over 100 K lower than that predicted by the simulation performed without soot radiation. From result strongly suggests the importance of considering a soot formation model for performing numerical simulations of a pulverized-coal combustion filed.