Numerical analysis of crystallization of high aspect ratio crystals with breakage

Abstract Two-dimensional population balance model for continuous cooling crystallization of high aspect ratio crystals, including nucleation, size-dependent growth of the two characteristic crystal facets and binary breakage along the crystal length is developed and analysed. The population balance equation is solved numerically using the 2D quadrature method of moments. The behaviour of the crystallizer and the effects of kinetic and process parameters on the characteristics of the crystal size distribution are analysed by detailed numerical experimentation. The binary breakage is described by four breakage rate functions comparing their effects on the mean length and mean aspect ratio of crystals. The simulation results revealed several interesting phenomena caused by interactions of the nonlinear sub-processes playing a relevant role in formulation of high aspect ratio crystals. Interactions of nucleation and size-dependent breakage may induce decrease in production of crystals. Size-dependent crystal breakage, decreasing the mean length of crystals induces some increase in their mean width.

[1]  Influence of breakage on crystal size distribution in a continuous cooling crystallizer , 2012 .

[3]  Chao Y. Ma,et al.  Model identification of crystal facet growth kinetics in morphological population balance modeling of l-glutamic acid crystallization and experimental validation , 2012 .

[4]  Alexander Reinhold,et al.  High dimensional population balances for the growth of faceted crystals: Combining Monte Carlo integral estimates and the method of characteristics , 2015 .

[5]  Jayanta Chakraborty,et al.  A new framework for solution of multidimensional population balance equations , 2007 .

[6]  Daniele Marchisio,et al.  Multivariate Quadrature-Based Moments Methods for turbulent polydisperse gas–liquid systems , 2013 .

[7]  K. Piotrowski,et al.  Monte Carlo modelling of continuous adiabatic MSMPR crystallizer , 2005 .

[8]  Ioannis G. Kevrekidis,et al.  Coarse-grained computation for particle coagulation and sintering processes by linking Quadrature Method of Moments with Monte-Carlo , 2010, J. Comput. Phys..

[9]  Marco Mazzotti,et al.  Measurement of size and shape distributions of particles through image analysis , 2008 .

[10]  Roberto Irizarry,et al.  Fast Monte Carlo methodology for multivariate particulate systems—I: Point ensemble Monte Carlo , 2008 .

[11]  Kai Sundmacher,et al.  Dimension reduction of bivariate population balances using the quadrature method of moments , 2011, Comput. Chem. Eng..

[12]  K. Sundmacher,et al.  On the Prediction of Crystal Shape Distributions in a Steady State Continuous Crystallizer , 2009 .

[13]  J. Favero,et al.  Solution of the population balance equation by the direct dual quadrature method of generalized moments , 2013 .

[14]  Peter J. Jansens,et al.  Two-dimensional population balance model with breakage of high aspect ratio crystals for batch crystallization , 2008 .

[15]  Carl M Schoellhammer,et al.  Computational and experimental investigation of needle-shaped crystal breakage. , 2011, International journal of pharmaceutics.

[16]  Robert McGraw,et al.  Representation of generally mixed multivariate aerosols by the quadrature method of moments: I. Statistical foundation , 2004 .

[17]  Robert McGraw,et al.  Description of Aerosol Dynamics by the Quadrature Method of Moments , 1997 .

[18]  Wolfgang Wagner,et al.  Numerical study of a stochastic particle algorithm solving a multidimensional population balance model for high shear granulation , 2010, J. Comput. Phys..

[19]  Stefan Heinrich,et al.  An efficient numerical technique for solving population balance equation involving aggregation, breakage, growth and nucleation , 2008 .

[20]  Robert McGraw,et al.  Representation of generally mixed multivariate aerosols by the quadrature method of moments: II. Aerosol dynamics , 2004 .

[21]  R. Fox,et al.  Application of the direct quadrature method of moments to polydisperse gas–solid fluidized beds , 2004 .

[22]  Kevin J. Roberts,et al.  Multi-dimensional population balance modeling of the growth of rod-like L-glutamic acid crystals using growth rates estimated from in-process imaging , 2007 .

[23]  Roberto Irizarry,et al.  Stochastic simulation of population balance models with disparate time scales: Hybrid strategies , 2011 .

[24]  D. E. Rosner,et al.  Bivariate Extension of the Quadrature Method of Moments for Modeling Simultaneous Coagulation and Sintering of Particle Populations. , 2001, Journal of colloid and interface science.

[25]  François Puel,et al.  Simulation and analysis of industrial crystallization processes through multidimensional population balance equations. Part 1: a resolution algorithm based on the method of classes , 2003 .

[26]  Béla G. Lakatos,et al.  Stability and dynamics of isothermal CMSMPR crystallizers , 2007 .

[27]  Ernst Dieter Gilles,et al.  A population model for crystallization processes using two independent particle properties , 2001 .

[28]  D. L. Ma,et al.  High-Resolution Simulation of Multidimensional Crystal Growth , 2002 .

[29]  B. Lakatos,et al.  Stability and Bifurcation of Continuous Cooling Crystallizers , 2005 .

[30]  Akos Borsos,et al.  Investigation and simulation of crystallization of high aspect ratio crystals with fragmentation , 2014 .

[31]  François Puel,et al.  Simulation and analysis of industrial crystallization processes through multidimensional population balance equations. Part 2: a study of semi-batch crystallization , 2003 .

[32]  Rodney O. Fox,et al.  Bivariate direct quadrature method of moments for coagulation and sintering of particle populations , 2006 .

[33]  Q. Yin,et al.  Determination of growth and breakage kinetics of l-threonine crystals , 2006 .

[34]  Jitendra Kumar,et al.  Accurate and efficient solution of bivariate population balance equations using unstructured grids , 2013 .

[35]  Margaritis Kostoglou,et al.  Optimal low order methods of moments for solving the fragmentation equation , 2004 .

[36]  Heiko Briesen,et al.  Simulation of crystal size and shape by means of a reduced two-dimensional population balance model , 2006 .

[37]  Zoltan K. Nagy,et al.  Simultaneous Quadrature Method of Moments for the Solution of Population Balance Equations, Using a Differential Algebraic Equation Framework , 2009 .

[38]  B. Biscans Impact attrition in crystallization processes. Analysis of repeated impacts events of individual crystals , 2004 .

[39]  Béla G. Lakatos,et al.  Modelling and simulation of suspension polymerization of vinyl chloride via population balance model , 2013, Comput. Chem. Eng..