# Population balance equation

Population balance equations (PBEs) have been introduced in several branches of modern science, mainly in Chemical Engineering, to describe the evolution of a population of particles. This includes topics like crystallization, liquid-liquid extraction, gas-liquid dispersions, liquid-liquid reactions, comminution, aerosol engineering, biology (where the separate entities are cells), polymerization, etc. Population balance equations can be said to be derived as an extension of the Smoluchowski coagulation equation which describes only the coalescence of particles. PBEs, more generally, define how populations of separate entities develop in specific properties over time. They are a set of Integro-partial differential equations which gives the mean-field behavior of a population of particles from the analysis of behavior of single particle in local conditions.[1] Particulate systems are characterized by the birth and death of particles. For example, consider precipitation process (formation of solid from liquid solution) which has the subprocesses nucleation, agglomeration, breakage, etc., that result in the increase or decrease of the number of particles of a particular radius (assuming formation of spherical particles). Population balance is nothing but a balance on the number of particles of a particular state (in this example, size).

## Formulation of PBE

Consider the average number of particles with particle properties denoted by a particle state vector (x,r) (where x corresponds to particle properties like size, density, etc. also known as internal coordinates and, r corresponds to spatial position or external coordinates) dispersed in a continuous phase defined by a phase vector Y(r,t) (which again is a function of all such vectors which denote the phase properties at various locations) is denoted by f(x,r,t). Hence it gives the particle characteristics in property and space domains. Let h(x,r,Y,t) denote the birth rate of particles per unit volume of particle state space, so the number conservation can be written as[1]

${\displaystyle {\frac {d}{dt}}\int _{\Omega _{x}(t)}dV_{x}\int _{\Omega _{r}(t)}dV_{r}\,f({\mathbf {x}},{\mathbf {r}},t)=\int _{\Omega _{x}(t)}dV_{x}\int _{\Omega _{r}(t)}dV_{r}\,h({\mathbf {x}},{\mathbf {r}},{\mathbf {Y}},t)}$

This is a generalized form of PBE.[1]

## Solution to PBE

Monte Carlo methods [2] ,[3] discretization methods and moment methods [2][3][4][5][6] are mainly used to solve these equations. The choice depends on the application and computing infrastructure.

## References

1. ^ a b c Ramkrishna, D.: Population Balances: Theory and Applications to Particulate Systems in Engineering, Academic Press, 2000
2. ^ a b Hashemian, N.; Armaou, A. (2016). "Simulation, model-reduction and state estimation of a two-component coagulation process". AIChE Journal. 62: 1557–1567. doi:10.1002/aic.15146.
3. ^ a b Hashemian, N.; Ghanaatpishe, M.; Armaou, A. (2016). "Development of a reduced order model for bi-component granulation processes via laguerre polynomials". Proceedings of the American Control Conference: 3668–3673. doi:10.1109/ACC.2016.7525483.
4. ^ Description of Aerosol Dynamics by the Quadrature Method of Moments, Robert McGrawa, Aerosol Science and Technology, Volume 27, Issue 2, 1997, pages 255-265
5. ^ Yu, M., Lin, J., and Chan, T. (2008). A New Moment Method for Solving the Coagulation Equation for Particles in Brownian Motion. Aerosol Sci. Technol., 42(9):705–713.
6. ^ Marchisio, D. L., and Fox, R. O. (2005). Solution of Population Balance Equa- tions Using the Direct Quadrature Method of Moments. J. Aerosol Sci., 36(1):43–73.