A differentially weighted Monte Carlo method for two-component coagulation

The direct simulation Monte Carlo (DSMC) method for population balance modeling is capable of retaining the history of each simulation particle and is thus able to deal with multivariate properties in a simple and straightforward manner. As opposed to conventional DSMC approaches that track equally weighted simulation particles, a differentially weighted Monte Carlo method is extended to simulate two-component coagulation processes and is thereby able to simulate the micromixing of the components. A new feature of the method for this bivariate population balance modeling is that it is possible to specify how the simulation particles are distributed over the compositional axis. This allows us to obtain information about particles in those regions of the size and composition distribution functions where the non-weighted MC methods place insufficient simulation particles to obtain an inaccurate solution. The new feature results in lower statistical noise for simulating two-component coagulation, which is validated by using two-component coagulation cases for which analytical solutions exist (a discrete process with sum kernel for initial monodisperse populations and a process with constant kernel for initial polydisperse populations).

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

[2]  T. Mexia,et al.  Author ' s personal copy , 2009 .

[3]  Roberto Irizarry,et al.  Fast Monte Carlo methodology for multivariate particulate systems—I: Point ensemble Monte Carlo , 2008 .

[4]  Scott L. Diamond,et al.  A General Algorithm for Exact Simulation of Multicomponent Aggregation Processes , 2002 .

[5]  Gerald Warnecke,et al.  Solving population balance equations for two-component aggregation by a finite volume scheme , 2007 .

[6]  Kangtaek Lee,et al.  Compositional distributions in multicomponent aggregation , 2008 .

[7]  S. Friedlander.  Smoke, Dust, and Haze: Fundamentals of Aerosol Dynamics , 2000 .

[8]  J. Seinfeld.  Atmospheric Chemistry and Physics of Air Pollution , 1986 .

[9]  John H. Seinfeld,et al.  Simulation of multicomponent aerosol dynamics , 1992 .

[10]  John H. Seinfeld,et al.  Simulation of multicomponent aerosol condensation by the moving sectional method , 1990 .

[12]  J. M. Fernández-Díaz,et al.  Exact solution of Smoluchowski's continuous multi-component equation with an additive kernel , 2007 .

[13]  John H. Seinfeld,et al.  Coagulation and growth of a multicomponent aerosol , 1978 .

[14]  Chuguang Zheng,et al.  Reducing Statistical Noise and Extending the Size Spectrum by Applying Weighted Simulation Particles in Monte Carlo Simulation of Coagulation , 2009 .

[15]  Ian J. Laurenzi.  Erratum to "A general algorithm for exact simulation of multicomponent aggregation processes" [J. Comput. Phys. 177(2002) 418-449] , 2006, J. Comput. Phys..

[16]  D. Gillespie.  A General Method for Numerically Simulating the Stochastic Time Evolution of Coupled Chemical Reactions , 1976 .

[17]  D. Ramkrishna,et al.  ON THE SOLUTION OF POPULATION BALANCE EQUATIONS BY DISCRETIZATION--I . A FIXED PIVOT TECHNIQUE , 2003 .

[18]  A. A Lushnikov,et al.  Evolution of coagulating systems: III. Coagulating mixtures , 1976 .

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

[20]  Jost O.L. Wendt,et al.  Toxic metal emissions from incineration: Mechanisms and control , 1993 .

[21]  F. Einar Kruis,et al.  Direct simulation Monte Carlo for simultaneous nucleation, coagulation, and surface growth in dispersed systems , 2004 .

[22]  Kangtaek Lee,et al.  Mixing of components in two‐component aggregation , 2006 .

[23]  F. Einar Kruis,et al.  Direct simulation Monte Carlo method for particle coagulation and aggregation , 2000 .

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

[25]  R. Axelbaum,et al.  Monte Carlo Simulation of Multicomponent Aerosols Undergoing Simultaneous Coagulation and Condensation , 2004 .

[26]  Costas Kiparissides,et al.  Dynamic evolution of the particle size distribution in suspension polymerization reactors: A comparative study on Monte Carlo and sectional grid methods , 2008 .

[27]  Robert McGraw,et al.  Chemically resolved aerosol dynamics for internal mixtures by the quadrature method of moments , 2003 .

[28]  Francis Filbet,et al.  Numerical Simulation of the Smoluchowski Coagulation Equation , 2004, SIAM J. Sci. Comput..

[29]  C. Kiparissides.  Polymerization reactor modeling: A review of recent developments and future directions , 1996 .

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

[31]  Mixing selectivity in bicomponent, bipolar aggregation , 2002 .

[32]  Timothy F. L. McKenna,et al.  Solution of the Population Balance Equation for Two-Component Aggregation by an Extended Fixed Pivot Technique , 2005 .

[33]  Chuguang Zheng,et al.  Correcting the multi-Monte Carlo method for particle coagulation , 2009 .

[34]  G. Raga,et al.  Monte Carlo simulations of two-component drop growth by stochastic coalescence , 2008 .

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

[36]  James D. Litster,et al.  Scaleup of wet granulation processes: science not art , 2003 .