Population balance breakage probability

The main contribution from the work of Luo [95, 96] was a closure model for binary breakage of fluid particles in fully developed turbulence flows based on isotropic turbulence - and probability theories. The author(s) also claimed that this model contains no adjustable parameters, a better phrase may be no additional adjustable parameters as both the isotropic turbulence - and the probability theories involved contain adjustable parameters and distribution functions. Hagesaether et al [49, 50, 51, 52] continued the population balance model development of Luo within the framework of an idealized plug flow model, whereas Bertola et al [13] combined the extended population balance module with a 2D algebraic slip mixture model for the flow pattern. Bertola et al [13] studied the effect of the bubble size distribution on the flow fields in bubble columns. An extended k-e model was used describing turbulence of the mixture flow. Two sets of simulations were performed, i.e., both with and without the population balance involved. Four different superficial gas velocities, i.e., 2,4,6 and 8 (cm/s) were used, and the superficial liquid velocity was set to 1 (cm/s) in all the cases. The population balance contained six prescribed bubble classes with diameters set to = 0.0038 (m), d = 0.0048 (m), di = 0.0060 (m), di = 0.0076 (m), di = 0.0095 (m) and di = 0.0120 (m). 

Lee et al [66] and Prince and Blanch [92] adopted the basic ideas of Coulaloglou and Tavlarides [16] formulating the population balance source terms directly on the averaging scales performing analysis of bubble breakage and coalescence in turbulent gas-liquid dispersions. The source term closures were completely integrated parts of the discrete numerical scheme adopted. The number densities of the bubbles were thus defined as the number of bubbles per unit mixture volume and not as a probability density in accordance with the kinetic theory of gases. 

The first term of the population balance, PBjSip, relates to the formation of droplets with a volume in the interval [F, F+dV] by breakage of droplets with a larger volume (maximally F ax)- In this term, the mesoscale parameter g U) is the breakage coefficient for a droplet with a volume U, u(U) is the number of droplets formed upon breakage of a droplet with a volume U (typically two), and P U,V) reflects the probability that a droplet with a volume U breaks into a droplet with a volume F. The second term represents the formation of droplets in the volume interval [F,... 

---