The Linear Process Deferment Algorithm: A new technique for solving population balance equations

In this paper a new stochastic algorithm for the solution of population balance equations is presented. The population balance equations have the form of extended Smoluchowski equations which include linear and source terms. The new algorithm, called the linear process deferment algorithm (LPDA), is used for solving a detailed model describing the formation of soot in premixed laminar flames. A measure theoretic formulation of a stochastic jump process is developed and the corresponding generator presented. The numerical properties of the algorithm are studied in detail and compared to the direct simulation algorithm and various splitting techniques. LPDA is designed for all kinds of population balance problems where nonlinear processes cannot be neglected but are dominated in rate by linear ones. In those cases the LPDA is seen to reduce run times for a population balance simulation by a factor of up to 1000 with a negligible loss of accuracy.

[1]  Andreas Eibeck,et al.  Approximative solution of the coagulation–fragmentation equation by stochastic particle systems , 2000 .

[2]  Karl Sabelfeld,et al.  Stochastic particle methods for Smoluchowski coagulation equation: variance reduction and error estimations , 2003, Monte Carlo Methods Appl..

[3]  N. N. Yanenko.  Splitting Methods for Partial Differential Equations , 1971, IFIP Congress.

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

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

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

[7]  F. Einar Kruis,et al.  Direct simulation Monte Carlo for simultaneous nucleation, coagulation, and surface growth in dispersed systems , 2004 .

[8]  Themis Matsoukas,et al.  Constant-number Monte Carlo simulation of population balances , 1998 .

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

[10]  Karl K. Sabelfeld,et al.  Stochastic algorithms for solving Smolouchovsky coagulation equation and applications to aerosol growth simulation. , 1996, Monte Carlo Methods Appl..

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

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

[13]  Graeme A. Bird,et al.  Molecular Gas Dynamics , 1976 .

[14]  Benjamin Jourdain,et al.  A stochastic approach for the numerical simulation of the general dynamics equation for aerosols , 2003 .

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

[16]  Andreas Eibeck,et al.  Stochastic interacting particle systems and nonlinear kinetic equations , 2003 .