Chebyshev Spectral Collocation Method for Population Balance Equation in Crystallization

The population balance equation (PBE) is the main governing equation for modeling dynamic crystallization behavior. In the view of mathematics, PBE is a convection–reaction equation whose strong hyperbolic property may challenge numerical methods. In order to weaken the hyperbolic property of PBE, a diffusive term was added in this work. Here, the Chebyshev spectral collocation method was introduced to solve the PBE and to achieve accurate crystal size distribution (CSD). Three numerical examples are presented, namely size-independent growth, size-dependent growth in a batch process, and with nucleation, and size-dependent growth in a continuous process. Through comparing the results with the numerical results obtained via the second-order upwind method and the HR-van method, the high accuracy of Chebyshev spectral collocation method was proven. Moreover, the diffusive term is also discussed in three numerical examples. The results show that, in the case of size-independent growth (PBE is a convection equation), the diffusive term should be added, and the coefficient of the diffusive term is recommended as 2 G × 10 −3 to G × 10 −2 , where G is the crystal growth rate.