A General Algorithm for Exact Simulation of Multicomponent Aggregation Processes

详细信息    查看全文

• 作者：Laurenzi ; Ian J. ; Bartels ; John D. ; Diamond ; Scott L.
• 关键词：coagulation ; clustering ; collection ; moments ; gel transition
• 刊名：Journal of Computational Physics
• 出版年：2002
• 出版时间：April 10, 2002
• 年：2002
• 卷：177
• 期：2
• 页码：418-449
• 全文大小：935 K

文摘

A Monte Carlo (MC) algorithm is presented for the simulation of the time evolution of aggregation processes featuring multiple components, properties, or conservation laws. Instead of using deterministic differential population balance equations, the MC algorithm utilizes a stochastic approach to aggregation kinetics. As a result, exact simulation of spatially independent aggregation processes is possible without the need for numerical approximations. Furthermore, simulations exactly predict all moments of the size and composition distributions of aggregating particles for both nongelling and gelling kernels and extend these results to the postgelation period. The algorithm is shown to require at most O((Ω₁Ω₂…Ω_κ)^1/(κ+1)) rate-limiting operations per time step for a κ-component aggregation process featuring Ω_i monomers of each component i—a substantial performance improvement over the potential of previous methods. Simulation results are presented for bivariate sum, product, and constant kernels, and for the perikinetic (Brownian) kernel.