Population Balance Modelling
Guide to CFD for Polydisperse Flows
- Disperse Multiphase Flows
- Polydisperse Multiphase Flows
- Population Balance Modelling
- Population Balance Modelling in OpenFOAM
This article provides the basic concepts behind the population balance model based on the method of classes, used for the simulation of polydisperse multiphase flows. A detailed theoretical description can be found in the associated publication by Lehnigk et al., AIChE J, 2022, Vol. 68(3), e17539.
Number density function
The method solves a population balance equation, an “integro-partial differential equation” for the number density function, nv, for the distribution of particle volume, v. The number density function is the number of particles per unit volume with a property (here v, in the case of nv) within a range of values (here, ∆v), normalised by the range. For example, if there are 24 air bubbles per m3 with bubble volumes in the range ∆v = 3-5 mm3, the average number density corresponding to v = 4 mm3 (the midpoint) would be nv = 24/((5-3) m3·mm3) = 12 particles/(m3·mm3).
The same integration is carried out for the entire population balance equation, which leads to a set of transport equations for the individual number concentrations. The arising source terms are written in a way that the conservation of the total volume and number of particles is ensured, irrespective of how the particle size domain is discretized.
Size group fractions
Each number concentration is then converted into a volume concentration, αi , by multiplying it with the representative particle volume vi, by:
The sum of volume concentrations yields the phase fraction α of the dispersed phase:
The phase fraction α of the dispersed phase is thereby known from the number density function, so, in theory, there is no need to solve a transport equation for α . However, from a numerical standpoint it is critical that α is strictly bounded between 0 and 1 and that mass conservation is obeyed. Therefore in practice multiphaseEulerFoam solves a transport equation for α which ensures conservation and, by using MULES (multidimensional universal limited explicit solver), boundedness is enforced.
With the phase fractions already solved for, the sole task of the population balance model is to determine how the dispersed phase elements are distributed with respect to size. This is achieved by transforming the transport equations for the number concentrations into equations for size group fractions fi which are bounded between zero and one and sum up to unity:
The interfacial models in the Euler-Euler model of multiphase flow describing the exchange of mass, momentum and energy rely on a surface area averaged diameter, also referred to as the Sauter diameter ds. From the size group fractions, d is calculated from the equivalent diameters di of each size group, e.g. di = (6vi /π)1/3 for spherical particles.
The next article describes the basics of setting up a polydisperse multiphase flow simulation within multiphaseEulerFoam.