Search

MiS Preprint Repository

We have decided to discontinue the publication of preprints on our preprint server as of 1 March 2024. The publication culture within mathematics has changed so much due to the rise of repositories such as ArXiV (www.arxiv.org) that we are encouraging all institute members to make their preprints available there. An institute's repository in its previous form is, therefore, unnecessary. The preprints published to date will remain available here, but we will not add any new preprints here.

MiS Preprint
111/2006

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

Wolfgang Hackbusch

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 $\int_{0}^{x}\kappa(x-y,y)f(y)f(x-y)\mathrm{d}y.$ The discretisation may use a locally refined grid. In this paper we describe an algorithm which (i) is efficient (the cost is $\mathcal{O}(N\log N),$ $N$: data size) and (ii) ensures local mass conservation.

Received:
Oct 9, 2006
Published:
Oct 9, 2006
MSC Codes:
44A35, 42A55, 45E99, 45K05, 92D25
Keywords:
convolution integral, non-uniform grids, conservation of mass

Related publications

inJournal
2007 Repository Open Access
Wolfgang Hackbusch

Approximation of coalescence integrals in population balance models with local mass conservation

In: Numerische Mathematik, 106 (2007) 4, pp. 627-657