A framework for the simulation of mass crystallization considering the effect of fluid dynamics

The behavior of large-scale crystallizers is strongly affected by the fluid dynamic characteristics of the apparatus like e.g. local supersaturation and velocity profiles. The simulation of this effect is a complex multi-scale and multi-phenomena problem. This contribution presents an approach to solve the coupled problem of crystallization and fluid dynamics by means of software integration. Existing specialized software tools (Fluent and Parsival) are employed for the solution of specific subproblems, namely the solution of population balance models and fluid dynamics. To reflect the phenomena on their characteristic scales different grids are used during the simulation of the respective subproblems. The population balance and crystallization kinetics are formulated in the coarse scale compartments while the fluid dynamics are solved on the fine CFD grid. The problem decomposition needed for the formulation of the subproblems and the proper selection of the compartments are discussed. The approach is validated with a simple model of a crystallization in a tubular reactor which can be solved using reduced methods without introducing a systematic error.

[1]  Alex Pothen,et al.  Computing the block triangular form of a sparse matrix , 1990, TOMS.

[2]  H. Kramer,et al.  Towards on-scale crystalliser design using compartmental models , 1998 .

[3]  John K. Reid,et al.  An Implementation of Tarjan's Algorithm for the Block Triangularization of a Matrix , 1978, TOMS.

[4]  Wolfgang Marquardt,et al.  CHEOPS: A tool-integration platform for chemical process modelling and simulation , 2004, International Journal on Software Tools for Technology Transfer.

[5]  Benjamin Kuipers,et al.  Model Decomposition and Simulation , 2004 .

[6]  Costas Kiparissides,et al.  CFD analysis of turbulence non-homogeneity in mixing vessels: A two-compartment model , 2002 .

[7]  G. R. Jerauld,et al.  Simple conditions for the appearance of sustained oscillations in continuous crystallizers , 1983 .

[8]  M. Hounslow,et al.  A discretized population balance for nucleation, growth, and aggregation , 1988 .

[9]  Daniele Marchisio,et al.  Solution of population balance equations using the direct quadrature method of moments , 2005 .

[10]  Paul I. Barton,et al.  Index and Characteristic Analysis of Linear PDAE Systems , 2002, SIAM J. Sci. Comput..

[11]  Daniele Marchisio,et al.  Implementation of the Quadrature Method of Moments in CFD codes for aggregation-breakage problems , 2003 .

[12]  Joel H. Ferziger,et al.  Computational methods for fluid dynamics , 1996 .

[13]  M. Wulkow,et al.  The simulation of molecular weight distributions in polyreaction kinetics by discrete Galerkin methods , 1996 .

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

[15]  Jesse T. Pikturna,et al.  Quadrature method of moments for population‐balance equations , 2003 .

[16]  Rodney O. Fox,et al.  Comparison of micromixing models for CFD simulation of nanoparticle formation , 2004 .

[17]  Alfons Mersmann,et al.  Brittle fracture in crystallization processes Part A. Attrition and abrasion of brittle solids , 1999 .

[18]  I. Duff,et al.  Direct Methods for Sparse Matrices , 1987 .

[19]  Heiko Briesen,et al.  Scale Integration for the Coupled Simulation of Crystallization and Fluid Dynamics , 2005 .

[20]  Tuomas Koiranen,et al.  Secondary nucleation due to crystal-impeller and crystal- vessel collisions by population balances in CFD-modelling , 2002 .

[21]  Michael J. Hounslow,et al.  A micro-mechanical model for the rate of aggregation during precipitation from solution , 2001 .

[22]  D. Ramkrishna,et al.  On the solution of population balance equations by discretization—II. A moving pivot technique , 1996 .

[23]  William R. Paterson,et al.  Orthokinetic Aggregation During Precipitation: A Computational Model for Calcium Oxalate Monohydrate , 1997 .

[24]  Ulrich Nieken,et al.  Modeling and simulation of crystallization processes using parsival , 2001 .

[25]  W. Marquardt,et al.  Structural analysis of differential-algebraic equation systems—theory and applications , 1995 .

[26]  Sandro Macchietto,et al.  A general methodology for hybrid multizonal/CFD models: Part I. Theoretical framework , 2004, Comput. Chem. Eng..

[27]  Peter J.T. Verheijen,et al.  Modeling of industrial crystallizers for control and design purposes , 2000 .

[28]  Sandro Macchietto,et al.  A general methodology for hybrid multizonal/CFD models: Part II. Automatic zoning , 2004, Comput. Chem. Eng..

[29]  Wolfgang Marquardt,et al.  Modular dynamic simulation for integrated particulate processes by means of tool integration , 2005 .

[30]  T. Serra,et al.  Effect of the shear and volume fraction on the aggregation and breakup of particles , 1998 .

[31]  Z. Sha,et al.  Mixing and crystallization in suspensions , 2000 .

[32]  Panagiotis D. Christofides,et al.  Nonlinear control of particulate processes , 1999 .

[33]  Jerzy Bałdyga,et al.  Turbulent Mixing and Chemical Reactions , 1999 .

[34]  P. Ayazi Shamlou,et al.  Crystal break-up in dilute turbulently agitated suspensions , 1993 .