Numerical approximations of a population balance model for coupled batch preferential crystallizers

  • Authors:
  • S. Qamar;I. Angelov;M. P. Elsner;A. Ashfaq;A. Seidel-Morgenstern;G. Warnecke

  • Affiliations:
  • Institute for Analysis and Numerics, Otto-von-Guericke University, PSF 4120, 39106 Magdeburg, Germany;Max Planck Institute for Dynamics of Complex Technical Systems, Sandtorstrasse 1, 39106 Magdeburg, Germany;Max Planck Institute for Dynamics of Complex Technical Systems, Sandtorstrasse 1, 39106 Magdeburg, Germany;Institute for Analysis and Numerics, Otto-von-Guericke University, PSF 4120, 39106 Magdeburg, Germany;Max Planck Institute for Dynamics of Complex Technical Systems, Sandtorstrasse 1, 39106 Magdeburg, Germany;Institute for Analysis and Numerics, Otto-von-Guericke University, PSF 4120, 39106 Magdeburg, Germany

  • Venue:
  • Applied Numerical Mathematics
  • Year:
  • 2009

Quantified Score

Hi-index 0.00

Visualization

Abstract

This article is concerned with the numerical approximations of population balance equations for modeling coupled batch preferential crystallization processes. The current setup consists of two batch crystallizers interconnected with two fines dissolution pipes. The crystallization of both enantiomers is assumed to take place in separate crystallizers after seeding with their corresponding crystals. The withdrawn fines are assumed to be dissolved in the dissolution unit after heating. crystallizer, the crystallizer temperature before entering to the opposite crystallizer. Two types of numerical methods are proposed for the simulation of this process. The first method uses high resolution finite volume schemes, while the second method is the so-called method of characteristics. On the one hand, the finite volume schemes which were derived for general system in divergence form, are computationally efficient, give desired accuracy on coarse grids, and are robust. On the other hand, the method of characteristics is in general a powerful tool for solving growth processes, has capability to overcome numerical diffusion and dispersion, gives highly resolved solutions, as well as being computationally efficient. A numerical test problem with both isothermal and non-isothermal conditions is considered here. The numerical results show clear advantages of the proposed schemes.