A Monte Carlo method for the simulation of coagulation and nucleation based on weighted particles and the concepts of stochastic resolution and merging

Abstract Monte Carlo simulations based on weighted simulation particles can solve a variety of population balance problems and allow thus to formulate a solution-framework for many chemical engineering processes. This study presents a novel concept for the calculation of coagulation rates of weighted Monte Carlo particles by introducing a family of transformations to non-weighted Monte Carlo particles. The tuning of the accuracy (named ‘stochastic resolution’ in this paper) of those transformations allows the construction of a constant-number coagulation scheme. Furthermore, a parallel algorithm for the inclusion of newly formed Monte Carlo particles due to nucleation is presented in the scope of a constant-number scheme: the low-weight merging. This technique is found to create significantly less statistical simulation noise than the conventional technique (named ‘random removal’ in this paper). Both concepts are combined into a single GPU-based simulation method which is validated by comparison with the discrete-sectional simulation technique. Two test models describing a constant-rate nucleation coupled to a simultaneous coagulation in 1) the free-molecular regime or 2) the continuum regime are simulated for this purpose.

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

[2]  Toivo T. Kodas,et al.  Aerosol Processing of Materials , 1998 .

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

[4]  Wolfgang Wagner,et al.  On a multivariate population balance model to describe the structure and composition of silica nanoparticles , 2012, Comput. Chem. Eng..

[5]  Vanni,et al.  Approximate Population Balance Equations for Aggregation-Breakage Processes. , 2000, Journal of colloid and interface science.

[6]  D. Oxtoby.  Homogeneous nucleation: theory and experiment , 1992 .

[7]  Chuguang Zheng,et al.  A new event-driven constant-volume method for solution of the time evolution of particle size distribution , 2009, J. Comput. Phys..

[8]  Stefan Heinrich,et al.  A generic population balance model for simultaneous agglomeration and drying in fluidized beds , 2007 .

[9]  Benjamin J. McCoy,et al.  Time evolution to similarity solutions for polymer degradation , 1998 .

[10]  Sotiris E. Pratsinis,et al.  Formation of agglomerate particles by coagulation and sintering—Part II. The evolution of the morphology of aerosol-made titania, silica and silica-doped titania powders , 1993 .

[11]  Jason Wittenberg,et al.  Clarify: Software for Interpreting and Presenting Statistical Results , 2003 .

[12]  Jianming Wei,et al.  A Fast Monte Carlo Method Based on an Acceptance-Rejection Scheme for Particle Coagulation , 2013 .

[13]  Andreas Eibeck,et al.  Stochastic Particle Approximations for Smoluchoski’s Coagualtion Equation , 2001 .

[14]  Antonios Armaou,et al.  Constant number Monte Carlo simulation of population balances with multiple growth mechanisms , 2010 .

[15]  T. Tsang,et al.  Comparison of Different Numerical Schemes for Condensational Growth of Aerosols , 1988 .

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

[17]  Sotiris E. Pratsinis,et al.  A discrete-sectional model for particulate production by gas-phase chemical reaction and aerosol coagulation in the free-molecular regime , 1990 .

[18]  Costas Kiparissides,et al.  Part III: Dynamic evolution of the particle size distribution in batch and continuous particulate processes: A Galerkin on finite elements approach , 2005 .

[19]  D. Ramkrishna,et al.  Population balance modeling: current status and future prospects. , 2014, Annual review of chemical and biomolecular engineering.

[20]  Daniel T. Gillespie,et al.  An Exact Method for Numerically Simulating the Stochastic Coalescence Process in a Cloud , 1975 .

[21]  Chuguang Zheng,et al.  Accelerating population balance-Monte Carlo simulation for coagulation dynamics from the Markov jump model, stochastic algorithm and GPU parallel computing , 2015, J. Comput. Phys..

[22]  Kurt Liffman,et al.  A direct simulation Monte-Carlo method for cluster coagulation , 1992 .

[23]  Markus Kraft,et al.  Stochastic solution of population balance equations for reactor networks , 2014, J. Comput. Phys..

[24]  David A. H. Jacobs,et al.  The State of the Art in Numerical Analysis. , 1978 .

[25]  Robert I. A. Patterson,et al.  Stochastic weighted particle methods for population balance equations with coagulation, fragmentation and spatial inhomogeneity , 2015, J. Comput. Phys..

[26]  Roberto Irizarry,et al.  Fast Compartmental Monte Carlo Simulation of Population Balance Models: Application to Nanoparticle Formation in Nonhomogeneous Conditions , 2012 .

[27]  Alejandro L. Garcia,et al.  A Monte Carlo simulation of coagulation , 1987 .

[28]  S. George,et al.  Nanocoating individual cohesive boron nitride particles in a fluidized bed by ALD , 2004 .

[29]  Doraiswami Ramkrishna,et al.  Population Balances: Theory and Applications to Particulate Systems in Engineering , 2000 .

[30]  Jianming Wei,et al.  A GPU-based parallelized Monte-Carlo method for particle coagulation using an acceptance–rejection strategy , 2013 .

[31]  Kangtaek Lee,et al.  Solution of the population balance equation using constant-number Monte Carlo , 2002 .

[32]  Markus Kraft,et al.  Coupling a stochastic soot population balance to gas-phase chemistry using operator splitting , 2007 .

[33]  Zhao Haibo,et al.  Multi-Monte Carlo approach for general dynamic equation considering simultaneous particle coagulation and breakage , 2005 .

[34]  Population Balance-Monte Carlo Simulation for Gas-to-Particle Synthesis of Nanoparticles , 2013 .

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

[36]  J. Seinfeld,et al.  Atmospheric Chemistry and Physics: From Air Pollution to Climate Change , 1997 .

[37]  Shih-Yuan Lu.  Collision integrals of discrete‐sectional model in simulating powder production , 1994 .

[38]  J. Wei,et al.  GPU-accelerated Monte Carlo simulation of particle coagulation based on the inverse method , 2013, J. Comput. Phys..

[39]  Ernst Dieter Gilles,et al.  Development, analysis and validation of population models for continuous and batch crystallizers , 2002 .

[40]  Sarma V. Pisupati,et al.  Models of agglomerate growth in fluidized bed reactors: Critical review, status and applications , 2014 .