Lattice Boltzmann method for multi-dimensional population balance models in crystallization

Abstract In this work, lattice Boltzmann method (LBM) is developed for efficient and accurate solution of multi-dimensional population balance equations (PBEs) used to model crystallization processes with growth and nucleation. Detailed derivation of LBM for multi-dimensional advection equation is presented, where the velocity is a function of space coordinates. The developed scheme is subsequently applied to solve multi-dimensional PBEs with size-dependent growth rate by drawing an analogy between the advection equation and the PBE. The computational advantage of LBM is shown by solving several benchmark examples taken from the literature and comparing the results with those obtained using well-established high resolution (HR) method. It is found that LBM provides at least as accurate solution as HR method, while requiring lower computation time.

[1]  R. LeVeque.  Finite Volume Methods for Hyperbolic Problems: Characteristics and Riemann Problems for Linear Hyperbolic Equations , 2002 .

[2]  Monte Carlo simulation for the solution of the bi-variate dynamic population balance equation in batch particulate systems , 2007 .

[3]  Peter V. Coveney,et al.  HemeLB: A high performance parallel lattice-Boltzmann code for large scale fluid flow in complex geometries , 2008, Comput. Phys. Commun..

[4]  S. Ansumali,et al.  Hydrodynamics beyond Navier-Stokes: exact solution to the lattice Boltzmann hierarchy. , 2007, Physical review letters.

[5]  Peter M. A. Sloot,et al.  Lattice-Boltzmann hydrodynamics on parallel systems , 1998 .

[6]  Shamsul Qamar,et al.  Adaptive high-resolution schemes for multidimensional population balances in crystallization processes , 2007, Comput. Chem. Eng..

[7]  D. L. Ma,et al.  Optimal control and simulation of multidimensional crystallization processes , 2002 .

[8]  A. Seidel-Morgenstern,et al.  Numerical solutions of population balance models in preferential crystallization , 2008 .

[9]  Richard D. Braatz,et al.  High‐order simulation of polymorphic crystallization using weighted essentially nonoscillatory methods , 2009 .

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

[11]  Jitendra Kumar,et al.  The cell average technique for solving multi-dimensional aggregation population balance equations , 2008, Comput. Chem. Eng..

[12]  I. Karlin,et al.  Stabilization of the lattice boltzmann method by the H theorem: A numerical test , 2000, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[13]  X. Wang,et al.  Crystal Growth Rate Dispersion Modeling Using Morphological Population Balance , 2008 .

[14]  Doraiswami Ramkrishna,et al.  Efficient solution of population balance equations with discontinuities by finite elements , 2002 .

[15]  Iliya V. Karlin,et al.  Perfect entropy functions of the Lattice Boltzmann method , 1999 .

[16]  Arvind Rajendran,et al.  Entropic lattice Boltzmann method for crystallization processes , 2010 .

[17]  Rosner,et al.  Monte Carlo Simulation of Particle Aggregation and Simultaneous Restructuring. , 1999, Journal of colloid and interface science.

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

[19]  Heinz Pitsch,et al.  Accuracy of higher-order lattice Boltzmann methods for microscale flows with finite Knudsen numbers , 2007, J. Comput. Phys..

[20]  Michael J. Hounslow,et al.  Finite-element methods for steady-state population balance equations , 1998 .

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

[22]  Shiyi Chen,et al.  LATTICE BOLTZMANN METHOD FOR FLUID FLOWS , 2001 .

[23]  S. Succi,et al.  Mesoscopic Models of Liquid/Solid Phase Transitions , 1998 .

[24]  Taslima Khanam,et al.  Three dimensional digital holographic profiling of micro-fibers. , 2009, Optics express.

[25]  Iliya V. Karlin,et al.  Elements of the lattice Boltzmann method I: Linear advection equation , 2006 .

[26]  Sanjeev Kumar,et al.  A new discretization of space for the solution of multi-dimensional population balance equations , 2008 .

[27]  Zoltan K. Nagy,et al.  Automatic differentiation‐based quadrature method of moments for solving population balance equations , 2012 .

[28]  Richard D. Braatz,et al.  Parallel high‐resolution finite volume simulation of particulate processes , 2008 .

[29]  Y. Pomeau,et al.  Lattice-gas automata for the Navier-Stokes equation. , 1986, Physical review letters.

[30]  James B. Rawlings,et al.  The Potential of Current High-Resolution Imaging-Based Particle Size Distribution Measurements for Crystallization Monitoring , 2009 .

[31]  D. Ramkrishna,et al.  On the solution of population balance equations by discretization - III. Nucleation, growth and aggregation of particles , 1997 .

[32]  K. Binder,et al.  A Guide to Monte Carlo Simulations in Statistical Physics , 2000 .

[33]  S. Katz,et al.  Some problems in particle technology: A statistical mechanical formulation , 1964 .

[34]  Haibo Zhao,et al.  Analysis of four Monte Carlo methods for the solution of population balances in dispersed systems , 2007 .

[35]  Caliane Bastos Borba Costa,et al.  Considerations on the crystallization modeling: Population balance solution , 2007, Comput. Chem. Eng..

[36]  M. Mazzotti,et al.  Model-Based Optimization of Particle Size Distribution in Batch-Cooling Crystallization of Paracetamol , 2004 .

[37]  D. Patil,et al.  Higher-order Galilean-invariant lattice Boltzmann model for microflows: single-component gas. , 2010, Physical review. E, Statistical, nonlinear, and soft matter physics.

[38]  Jonas Latt,et al.  Hydrodynamic limit of lattice Boltzmann equations , 2007 .

[39]  Thomas J. Naughton,et al.  Microparticle characterization using digital holography , 2010 .

[40]  R. Benzi,et al.  The lattice Boltzmann equation: theory and applications , 1992 .

[41]  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 .

[42]  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 .

[43]  Richard D. Braatz,et al.  Modelling and control of combined cooling and antisolvent crystallization processes , 2008 .

[44]  Costas Kiparissides,et al.  Part V: Dynamic evolution of the multivariate particle size distribution undergoing combined particle growth and aggregation , 2009 .

[45]  Wei Shyy,et al.  Lattice Boltzmann Method for 3-D Flows with Curved Boundary , 2000 .

[46]  S. Engell,et al.  WENO scheme with static grid adaptation for tracking steep moving fronts , 2009 .

[47]  M. Larson,et al.  Crystal size distributions in continuous crystallizers when growth rate is size dependent , 1968 .

[48]  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.

[49]  Zoltan K. Nagy,et al.  Combined Quadrature Method of Moments and Method of Characteristics Approach for Efficient Solution of Population Balance Models for Dynamic Modeling and Crystal Size Distribution Control of Crystallization Processes , 2009 .

[50]  Francis J. Doyle,et al.  A two-level discretisation algorithm for the efficient solution of higher-dimensional population balance models , 2008 .

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

[52]  M. Mazzotti,et al.  Modeling and Experimental Analysis of PSD Measurements through FBRM , 2000 .

[53]  Lynn F. Gladden,et al.  Single- and two-phase flow in fixed-bed reactors: MRI flow visualisation and lattice-Boltzmann simulations , 2001 .

[54]  Zhaolin Gu,et al.  Solution of population balance equation using quadrature method of moments with an adjustable factor , 2007 .

[55]  Kevin J. Roberts,et al.  Crystal growth measurement using 2D and 3D imaging and the perspectives for shape control , 2008 .

[56]  Sauro Succi,et al.  Massively Parallel Lattice-Boltzmann Simulation of Turbulent Channel Flow , 1997 .

[57]  Themis Matsoukas,et al.  Constant-number Monte Carlo simulation of population balances , 1998 .

[58]  S Succi,et al.  Lattice Boltzmann model for anisotropic liquid-solid phase transition. , 2001, Physical review letters.

[59]  S. Orszag,et al.  Extended Boltzmann Kinetic Equation for Turbulent Flows , 2003, Science.

[60]  V. Alopaeus,et al.  Solution of population balances with growth and nucleation by high order moment-conserving method of classes , 2007 .

[61]  Arvind Rajendran,et al.  Fast High-Resolution Method for Solving Multidimensional Population Balances in Crystallization , 2010 .

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

[63]  S. Succi,et al.  Phase-field lattice kinetic scheme for the numerical simulation of dendritic growth. , 2005, Physical review. E, Statistical, nonlinear, and soft matter physics.

[64]  H. C. Ottinger,et al.  Minimal entropic kinetic models for hydrodynamics , 2002, cond-mat/0205510.

[65]  Sauro Succi,et al.  Galilean-invariant lattice-Boltzmann models with H theorem. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[66]  J. Boon.  The Lattice Boltzmann Equation for Fluid Dynamics and Beyond , 2003 .

[67]  Doraiswami Ramkrishna,et al.  Analysis of population balance—IV: The precise connection between Monte Carlo simulation and population balances , 1981 .

[68]  R. Braatz,et al.  High resolution algorithms for multidimensional population balance equations , 2004 .

[69]  En Sup Yoon,et al.  On the solution of population balance equations (PBE) with accurate front tracking methods in practical crystallization processes , 2002 .

[70]  Qinjun Kang,et al.  Lattice Boltzmann model for crystal growth from supersaturated solution , 2004 .

[71]  D. Ramkrishna,et al.  Solution of population balance equations by MWR , 1977 .