Stochastic modeling of soot particle size and age distributions in laminar premixed flames

Abstract In the present study, the evolution of size distributions of soot particles in premixed laminar flames is investigated. A computationally efficient stochastic approach is used to model the dynamics of the soot particle ensemble. The stochastic method is enhanced to include a formalism for modeling coagulation at high pressures. A detailed kinetic model is used to describe gas phase reactions and the formation, growth, and oxidation of soot particles. Four fuel rich laminar premixed flames at different pressure are investigated. It is found that calculated size distributions vary markedly in shape between the different flames. The distributions calculated in the two flames at elevated pressure and in the sub-atmospheric flame exhibit log-normal shape in the post-flame zone. In contrast, a bimodal distribution is found for the atmospheric flame, which persists throughout the entire flame. The bimodal behavior could be attributed to continuous nucleation. Furthermore, the effect of surface ageing, i.e., the deactivation of sites on the soot particles’ surface available for reaction with gas phase species, is investigated. For this purpose, a definition of a particle’s age is introduced, and age distributions are calculated. Subsequently, it was investigated if the surface reactivity could be correlated with the particle age. Two different correlations were investigated: (a) a step-function attributing a high surface activity to young particles and low activity to old particles and (b) an exponential function giving a smooth transition of surface activity with particle age. Good agreement with measured soot volume fractions could be obtained with both approaches. The decay constant in the exponential correlation was found to be a linear function of maximum flame temperature for three of the four calculated flames. For these three flames, the experimentally established trend that surface deactivation proceeds faster at higher temperatures could thus be reproduced.

[1]  Michael Frenklach,et al.  Monte carlo simulation of soot aggregation with simultaneous surface growth-why primary particles appear spherical , 1998 .

[2]  N. Brown,et al.  Hydrogen migration in polyaromatic growth , 1998 .

[3]  H. Bockhorn,et al.  A detailed numerical study of the evolution of soot particle size distributions in laminar premixed flames. , 2001, Chemosphere.

[4]  H. Bockhorn,et al.  Some new observations concerning the mass growth of soot in premixed hydrocarbon-oxygen flames , 1989 .

[5]  B. Haynes,et al.  Soot surface growth at active sites , 1991 .

[6]  H. Bockhorn,et al.  Kinetic modeling of soot formation with detailed chemistry and physics: laminar premixed flames of C2 hydrocarbons , 2000 .

[7]  M. Frenklach,et al.  Detailed modeling of soot formation in laminar premixed ethylene flames at a pressure of 10 bar , 1995 .

[8]  A. M. Weiner,et al.  CHEMICAL KINETICS OF SOOT PARTICLE GROWTH , 1985 .

[9]  Daniel T. Gillespie,et al.  An Exact Method for Numerically Simulating the Stochastic Coalescence Process in a Cloud , 1975 .

[10]  H. Bockhorn,et al.  Soot Formation in Premixed Hydrocarbon Flames: Prediction of Temperature and Pressure Dependence , 1995 .

[11]  Adel F. Sarofim,et al.  A fully integrated Kinetic Monte Carlo/Molecular Dynamics approach for the simulation of soot precursor growth , 2002 .

[12]  S. Harris,et al.  Determination of the Rate Constant for Soot Surface Growth , 1983 .

[13]  B. Haynes,et al.  The Surface Growth Phenomenon in Soot Formation , 1982 .

[14]  Markus Kraft,et al.  An efficient stochastic algorithm for simulating Nano-particle dynamics , 2002 .

[15]  Markus Kraft,et al.  Measurement and numerical simulation of soot particle size distribution functions in a laminar premixed ethylene-oxygen-argon flame , 2003 .

[16]  D. E. Rosner,et al.  Multivariate Population Balances via Moment and Monte Carlo Simulation Methods: An Important Sol Reaction Engineering Bivariate Example and “Mixed” Moments for the Estimation of Deposition, Scavenging, and Optical Properties for Populations of Nonspherical Suspended Particles , 2003 .

[17]  Robert J. Kee,et al.  PREMIX :A F ORTRAN Program for Modeling Steady Laminar One-Dimensional Premixed Flames , 1998 .

[18]  Henning Bockhorn,et al.  Soot Formation in Combustion: Mechanisms and Models , 1994 .

[19]  M. Frenklach.  On surface growth mechanism of soot particles , 1996 .

[20]  Michael Frenklach,et al.  Aerosol dynamics modeling using the method of moments , 1987 .

[21]  M. Frenklach.  Method of moments with interpolative closure , 2002 .

[22]  M. Frenklach,et al.  On the relative contribution of acetylene and aromatics to soot particle surface growth , 1998 .

[23]  M. Frenklach,et al.  A detailed kinetic modeling study of aromatics formation in laminar premixed acetylene and ethylene flames , 1997 .

[24]  Mitchell D. Smooke,et al.  Computational and experimental study of soot formation in a coflow, laminar diffusion flame , 1999 .

[25]  Fabian Mauss,et al.  Inception and growth of soot particles in dependence on the surrounding gas phase , 1994 .

[26]  Markus Kraft,et al.  Modelling of a RDC using a combined CFD-population balance approach , 2004 .

[27]  M. Kraft,et al.  A stochastic approach to calculate the particle size distribution function of soot particles in laminar premixed flames , 2003 .

[28]  M. Frenklach,et al.  Detailed modeling of soot particle nucleation and growth , 1991 .

[29]  Michael Frenklach,et al.  Reaction mechanism of soot formation in flames , 2002 .