Size Distribution

In order to account for a size distribution of the dispersed phase, the mass and momentum transport equations have to be combined with a population balance equation (PBE). The central object of a PBE is the particle number density n(dp), which gives the number of particles (droplets or bubbles) whose diameters range from dp to dp+d(dp).

In its general form, the PBE reads:

Figure 1. EQUATION_DISPLAY
n(dp)t+∇⋅(vn(dp))=BD
(2169)

where B and D are birth and death terms due to breakup, coalescence, nucleation, dissolution, and so on [531].

The two most popular numerical methods for solving the PBE are the method of moments and the method of classes. Instead of solving for the number density, the method of moments solves for a (usually small) number of its moments:
Figure 2. EQUATION_DISPLAY
Sγ=dpγn(dp)(dp)
(2170)

In Simcenter STAR-CCM+, the method of moments based S-Gamma model is implemented. This model implementation relies on transporting one or two moments of the size distribution. The S-Gamma method is a computationally efficient tool suitable for engineering estimates.

The method of classes implies that the polydisperse phase is split into a number of size classes. Each class has its own mass and momentum balance. The standard method of classes, which is also known as multiple size-group (MUSIG) method [504], prescribes the sizes of the size-groups a priory. The sizes are constant throughout the computational domain, which makes the method time-consuming. In Simcenter STAR-CCM+, the partitioning of the particles ensemble into size classes is adaptive in space and time [565], therefore the name adaptive multiple size-group (AMUSIG) method. The AMUSIG method lets you work with a smaller number of size-groups than other class methods.