Technical Report 26, c4e-Preprint Series, Cambridge

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

ref: Technical Report 26, c4e-Preprint Series, Cambridge

Associated Themes: Numerics and Nanoparticles


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 (DSA) and various splitting techniques. LPDA is designed for all kinds of population balance problems where non-linear 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.

Material from this preprint has been published in: SIAM Journal on Scientific Computing 28 (1), 303-320, (2006)


PDF (398.5 KB)