Stochastic solution of population balance equations for reactor networks


This work presents a sequential modular approach to solve a generic network of reactors with a population balance model using a stochastic numerical method. Full-coupling to the gas-phase is achieved through operator-splitting. The convergence of the stochastic particle algorithm in test networks is evaluated as a function of network size, recycle fraction and numerical parameters. These test cases are used to identify methods through which systematic and statistical error may be reduced, including by use of stochastic weighted algorithms. The optimal algorithm was subsequently used to solve a one-dimensional example of silicon nanoparticle synthesis using a multivariate particle model. This example demonstrated the power of stochastic methods in resolving particle structure by investigating the transient and spatial evolution of primary polydispersity, degree of sintering and TEM-style images.

Access options

Keywords: coagulation, population balance, weights, networks, CSTR, compartmental model

Associated Projects: Numerics and Nanoparticles

*Corresponding author:
Telephone:Department +44 (0)1223 762784 (Dept) 769010 (CHU)
Mobile +49 173 3045528 and +44 7944 237879
Address:Department of Chemical Engineering
University of Cambridge
West Cambridge Site
Philippa Fawcett Drive
United Kingdom
Website:Personal Homepage