TY - JOUR

T1 - High dimensional population balances for the growth of faceted crystals

T2 - Combining Monte Carlo integral estimates and the method of characteristics

AU - Reinhold, Alexander

AU - Briesen, Heiko

N1 - Publisher Copyright:
© 2015 Elsevier Ltd.

PY - 2015/5/4

Y1 - 2015/5/4

N2 - A numerical scheme is proposed for the solution of population balance equations that describe the growth of faceted crystals. Therefore, two issues are discussed. First, a numerical scheme for solution of high dimensional population balance equations is proposed. The scheme is based on Monte Carlo integration for the integral distribution properties (number, area, volume) combined with the method of characteristics. Samples' points are initially obtained from an arbitrary probability distribution via the Metropolis-Hastings algorithm. The evolution of these sample points is subsequently computed using the method of characteristics, while the sample points themselves are the basis of the Monte Carlo integration of the population to obtain e.g. total crystal volume to close the mass balance. The theoretical background of the Monte Carlo integration allows us to estimate the error introduced in the mass balance depending on a particular choice of sampling points. The effectiveness of the scheme is verified for an analytically solvable 3-dimensional population. That case study is also employed to demonstrate that an adequate choice of the probability distribution for the generation of the initial sample points can improve the accuracy without increasing the computational effort. The second issue that must simultaneously be considered comprises the disappearance of faces during growth. Such conditions create a constraint for each face distance, measured from a chosen origin, which, prior to the disappearance of that face, had freely evolved according to its growth rate. This constraint is defined for general faceted crystals, and a procedure is proposed to correct the growth rate so that the constraint can be maintained in dynamic simulations. The application of this procedure is demonstrated in a second case study that simulates a 7-dimensional population balance equation.

AB - A numerical scheme is proposed for the solution of population balance equations that describe the growth of faceted crystals. Therefore, two issues are discussed. First, a numerical scheme for solution of high dimensional population balance equations is proposed. The scheme is based on Monte Carlo integration for the integral distribution properties (number, area, volume) combined with the method of characteristics. Samples' points are initially obtained from an arbitrary probability distribution via the Metropolis-Hastings algorithm. The evolution of these sample points is subsequently computed using the method of characteristics, while the sample points themselves are the basis of the Monte Carlo integration of the population to obtain e.g. total crystal volume to close the mass balance. The theoretical background of the Monte Carlo integration allows us to estimate the error introduced in the mass balance depending on a particular choice of sampling points. The effectiveness of the scheme is verified for an analytically solvable 3-dimensional population. That case study is also employed to demonstrate that an adequate choice of the probability distribution for the generation of the initial sample points can improve the accuracy without increasing the computational effort. The second issue that must simultaneously be considered comprises the disappearance of faces during growth. Such conditions create a constraint for each face distance, measured from a chosen origin, which, prior to the disappearance of that face, had freely evolved according to its growth rate. This constraint is defined for general faceted crystals, and a procedure is proposed to correct the growth rate so that the constraint can be maintained in dynamic simulations. The application of this procedure is demonstrated in a second case study that simulates a 7-dimensional population balance equation.

KW - Crystal growth

KW - Disappearing faces

KW - Morphological changes

KW - Morphological population balance

KW - Numerical solution

UR - http://www.scopus.com/inward/record.url?scp=84939208220&partnerID=8YFLogxK

U2 - 10.1016/j.ces.2015.01.035

DO - 10.1016/j.ces.2015.01.035

M3 - Article

AN - SCOPUS:84939208220

SN - 0009-2509

VL - 127

SP - 220

EP - 229

JO - Chemical Engineering Science

JF - Chemical Engineering Science

ER -