Abstract
Mathematical modeling is essential for the effective control of many chemical engineering processes, including crystallization. However, most existing crystallization models used in industry and academia assume ideal mixing. As a result, the unclear effects of imperfect mixing on crystallization, reported in experimental studies, remain largely unexplained. In this work we aim to address this gap in understanding by examining antisolvent crystallization processes on a general theoretical level, using a novel dimensionless model. To address the impact of mixing on crystallization, we employ the Engulfment model coupled with a population balance, and we nondimensionalize the model equations. Using this model, we explore the dependence of the mean particle size on the homogenization rate, represented by the Damköhler number for crystallization. Moreover, we study the impact of mixing at various values of the model's kinetic parameters to simulate difference in properties of individual products. We show that we are able to explain the complex interaction between crystallization and mixing, proving our model can serve as a tool for achieving a better understanding of the processes involved. Finally, due to its efficiency and reduced number of parameters, the model is suitable for direct fitting to experimental data.