A control oriented study on the numerical solution of the population balance equation for crystallization processes

The population balance equation provides a well established mathematical framework for dynamic modeling of numerous particulate processes. Numerical solution of the population balance equation is often complicated due to the occurrence of steep moving fronts and/or sharp discontinuities. This study aims to give a comprehensive analysis of the most widely used population balance solution methods, namely the method of characteristics, the finite volume methods and the finite element methods, in terms of the performance requirements essential for on-line control applications. The numerical techniques are used to solve the dynamic population balance equation of various test problems as well as industrial crystallization processes undergoing simultaneous nucleation and growth. The time-varying supersaturation profiles in the latter real-life case studies provide more realistic scenarios to identify the advantages and pitfalls of a particular numerical technique.The simulation results demonstrate that the method of characteristics gives the most accurate numerical predictions, whereas high computational burden limits its use for complex real crystallization processes. It is shown that the high order finite volume methods in combination with flux limiting functions are well capable of capturing sharp discontinuities and steep moving fronts at a reasonable computational cost, which facilitates their use for on-line control applications. The finite element methods, namely the orthogonal collocation and the Galerkin's techniques, on the other hand may severely suffer from numerical problems. This shortcoming, in addition to their complex implementation and low computational efficiency, makes the finite element methods less appealing for the intended application.