An explicit adaptive grid approach for the numerical solution of the population balance equation

Abstract

Many engineering applications, such as the formation of soot in hydrocarbon combustion or the precipitation of nanoparticles from aqueous solutions, encompass a polydispersed particulate phase that is immersed in a reacting carrier flow. From a Eulerian perspective, the evolution of the particulate phase both in physical and in particle property space can be described by the population balance equation (PBE). In this article, we present an explicit solution-adaptive numerical scheme for discretizing the spatially inhomogeneous and unsteady PBE along a one-dimensional particle property space. This scheme is based on a space and time dependent coordinate transformation which redistributes resolution in particle property space according to the shapes of recent solutions for the particle property distribution. In particular, the coordinate transformation is marched in time explicitly. In comparison to many existing moving or adaptive grid approaches, this has the advantage that the semi-discrete PBE does not need to be solved in conjunction with an additional system of equations governing the movement of nodes in particle property space.By design, our adaptive grid technique is able to accurately capture sharp features such as peaks or near-discontinuities, while maintaining the semi-discrete system size and adhering to a uniform fixed grid discretization in transformed particle property space. This is particularly advantageous if the PBE is combined with a spatially and temporally fully resolved flow model and a standard Eulerian solution scheme is applied in physical space. In order to accommodate localized source terms and to control the grid stretching, we develop a robust scheme for modifying the coordinate transformation such that constraints on the resolution in physical particle property space are obeyed.As an example, we consider the precipitation of BaSO4 particles from an aqueous solution in a plug flow reactor. Our findings demonstrate that for a given accuracy of the numerical solution the explicit adaptive grid technique requires over an order of magnitude fewer grid points than a comparable fixed grid discretization scheme.