A time-driven constant-number Monte Carlo method for the GPU-simulation of particle breakage based on weighted simulation particles

Abstract Monte Carlo (MC) simulations based on weighted particles offer novel and more precise techniques for the solution of the population balance equation for particulate systems. A recent constant-number approach named stochastic weighted algorithm (SWA) (Lee et al. (2015), J. Comput. Phys. (303) 1–18) has been developed, which renders the breakage of a simulation particle by an alteration of its properties, without the creation of novel simulation particles. The theoretic justification of the general formulation for all possible SWAs is limited to binary breakage kernels. We present a novel approach for the derivation of the properties of the MC particles representing fragments, which is applicable for all sorts of breakage kernels. This general scheme encompasses the already introduced SWA schemes, especially a number-based (SWA1, named NBS in this paper) and volume-based (SWA2, named VBS in this paper) breakage scheme, and it makes novel formulations possible: the low volume scheme (LVS), which renders preferably lower fragment sizes, and the combination of LVS with the NBS (LVS-NBS) or VBS (LVS-VBS). The implementation of these breakage schemes in the context of a GPU-based time-driven method is presented and the gained results are validated by comparison with results of the analytic solutions of a homogeneous binary breakage kernel. It is found, that the SWA methods (NBS and VBS) are only able to render large particle sizes, and that LVS, NBS-LVS and VBS-LVS are able to render the whole spectrum of particle sizes. Smaller noise levels are found for VBS and specific VBS-LVS schemes, making both more suitable for prolonged simulations than the other presented methods.

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

[2]  P. Das.  Monte Carlo simulation of drop breakage on the basis of drop volume , 1996 .

[3]  Transient and asymptotic behaviour of the binary breakage problem , 2005 .

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

[5]  Bijay K. Mishra.  Monte Carlo simulation of particle breakage process during grinding , 2000 .

[6]  Markus Kraft,et al.  Simulation of coalescence and breakage: an assessment of two stochastic methods suitable for simulating liquid–liquid extraction , 2004 .

[7]  Luc Devroye,et al.  Chapter 4 Nonuniform Random Variate Generation , 2006, Simulation.

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

[9]  Robert M. Ziff,et al.  The kinetics of cluster fragmentation and depolymerisation , 1985 .

[10]  Kangtaek Lee,et al.  Simultaneous coagulation and break-up using constant-N Monte Carlo , 2000 .

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

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

[13]  Doraiswami Ramkrishna,et al.  Simulation of particulate systems using the concept of the interval of quiescence , 1977 .

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

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

[16]  Benjamin Jourdain,et al.  A stochastic approach for the numerical simulation of the general dynamics equation for aerosols , 2003 .

[17]  Gregor Kotalczyk,et al.  A Monte Carlo method for the simulation of coagulation and nucleation based on weighted particles and the concepts of stochastic resolution and merging , 2017, J. Comput. Phys..

[18]  Cheng,et al.  Scaling theory of fragmentation. , 1988, Physical review letters.

[19]  Jitendra Kumar,et al.  Development and Convergence Analysis of a Finite Volume Scheme for Solving Breakage Equation , 2015, SIAM J. Numer. Anal..

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

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

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

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

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

[25]  Jitendra Kumar,et al.  Numerical simulation and convergence analysis of a finite volume scheme for solving general breakage population balance equations , 2013, Appl. Math. Comput..

[26]  Roberto Irizarry,et al.  Fast Monte Carlo methodology for multivariate particulate systems-II: τ-PEMC , 2008 .

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

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

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

[30]  Ka Ming Ng,et al.  New discretization procedure for the breakage equation , 1995 .

[31]  D. Gillespie.  Approximate accelerated stochastic simulation of chemically reacting systems , 2001 .

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