Approximation of Coalescence Integrals in Population Balance Models with Local Mass Conservation

Wolfgang Hackbusch

Contact the author: Please use for correspondence this email.
Submission date: 09. Oct. 2006 (revised version: April 2007)
Pages: 28
published in: Numerische Mathematik, 106 (2007) 4, p. 627-657 
DOI number (of the published article): 10.1007/s00211-007-0077-y
Bibtex
MSC-Numbers: 44A35, 42A55, 45E99, 45K05, 92D25
Keywords and phrases: convolution integral, non-uniform grids, conservation of mass
Download full preprint: PDF (326 kB), PS ziped (260 kB)

Abstract:
The solution of population balance equations is a function f(t,r,x) describing a population density of particles of the property x at time t and space r. For instance, the additional independent variable x may denote the mass of the particle. The describing equation contains additional sink and source terms involving integral operators. Since the coordinate x adds at least one further dimension to the spatial directions and time coordinate, an efficient numerical treatment of the integral terms is crucial. One of the more involved integral terms appearing in population balance models is the coalescence integral, which is of the form The discretisation may use a locally refined grid. In this paper we describe an algorithm which (i) is efficient (the cost is N: data size) and (ii) ensures local mass conservation.