High resolution algorithms for multidimensional population balance equations

Population balance equations have been used to model a wide range of processes including polymerization, crystallization, cloud formation, and cell dynamics. Rather than developing new algorithms specific to population balance equations, it is proposed to adapt the high-resolution finite volume methods developed for compressible gas dynamics, which have been applied to aerodynamics, astrophysics, detonation waves, and related fields where shock waves occur. High-resolution algorithms are presented for simulating multidimensional population balance equations with nucleation and size-dependent growth rates. For sharp distributions, these high-resolution algorithms can achieve improved numerical accuracy with orders-of-magnitude lower computational cost than other finite difference and finite volume algorithms. The algorithms are implemented in the ParticleSolver software package, which is applied to batch and continuous processes with one and multiple internal coordinates.