A method of characteristics for solving population balance equations (PBE) describing the adsorption of impurities during crystallization processes

• Published in: Chemical engineering science Vol. 65; no. 10; pp. 3191 - 3198
• Main Authors: FEVOTTE, François, FEVOTTE, Gilles
• Format: Journal Article
• Language: English
• Published: Kidlington Elsevier Ltd 15-05-2010; Elsevier
• Subjects: Adsorption; Applied sciences; Chemical and Process Engineering; Chemical engineering; Crystallization; Crystallization, leaching, miscellaneous separations; Dynamic simulation; Engineering Sciences; Exact sciences and technology; Impurities; Nucleation; Particulate processes; Population balance; Dynamic simulation; Particulate processes; Nucleation; Impurities; Population balance; Crystallization; Crystal growth; Method of characteristics; Impurity; Crystals; Unsteady state; Supersaturation; Algorithm; Modeling; Contamination; Kinetic model; Particle size distribution; Pollution; Adsorption; Dynamic model
• ISSN: 0009-2509; 1873-4405
• DOI: 10.1016/j.ces.2010.02.009

Possible hindering effects of impurities on the crystal growth were shown to take place because of the adsorption of impurity species on the crystal surface. Transient features of this adsorption were observed, such that the growth of a given crystal does not depend on supersaturation only, but also on the time a given particle spent in contact with impurities present in the mother liquor. Meanwhile, few kinetic models describe the effect of impurities on the growth of crystals in solution, and published models are usually derived from data obtained, thanks to specific experiments based on the evaluation of the growth rate of single crystals. Such models are obviously questionable because, in the industrial practice, distributed properties of crystals are actually involved. Considering the “time of contamination” of particles as a new internal variable is thus made necessary. This is the reason why a specific PBE resolution algorithm is presented in this paper. The numerical scheme for the resolution of PBEs is based on the method of characteristics and shown to allow fast and accurate simulation of transient features of the crystal size distribution in the particular case when the growth or nucleation rates are assumed to exhibit unsteady-state dynamics. The algorithm is finally used to simulate the isothermal desupersaturation crystallization of citric acid in water.