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 solutionframework 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 nonweighted Monte Carlo particles. The tuning of the accuracy (named ‘stochastic resolution’ in this paper) of those transformations allows the construction of a constantnumber 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 constantnumber scheme: the lowweight 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 GPUbased simulation method which is validated by comparison with the discretesectional simulation technique. Two test models describing a constantrate nucleation coupled to a simultaneous coagulation in 1) the freemolecular regime or 2) the continuum regime are simulated for this purpose.
 Authors:
 Publication Date:
 OSTI Identifier:
 22622302
 Resource Type:
 Journal Article
 Resource Relation:
 Journal Name: Journal of Computational Physics; Journal Volume: 340; Other Information: Copyright (c) 2017 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved.; Country of input: International Atomic Energy Agency (IAEA)
 Country of Publication:
 United States
 Language:
 English
 Subject:
 71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; ACCURACY; ALGORITHMS; BALANCES; COMPARATIVE EVALUATIONS; COMPUTERIZED SIMULATION; INCLUSIONS; MATHEMATICAL SOLUTIONS; MONTE CARLO METHOD; NOISE; NUCLEATION; PARTICLES; RANDOMNESS; RESOLUTION; STOCHASTIC PROCESSES; TRANSFORMATIONS; WEIGHT
Citation Formats
Kotalczyk, G., Email: Gregor.Kotalczyk@unidue.de, and Kruis, F.E.. A Monte Carlo method for the simulation of coagulation and nucleation based on weighted particles and the concepts of stochastic resolution and merging. United States: N. p., 2017.
Web. doi:10.1016/J.JCP.2017.03.041.
Kotalczyk, G., Email: Gregor.Kotalczyk@unidue.de, & Kruis, F.E.. A Monte Carlo method for the simulation of coagulation and nucleation based on weighted particles and the concepts of stochastic resolution and merging. United States. doi:10.1016/J.JCP.2017.03.041.
Kotalczyk, G., Email: Gregor.Kotalczyk@unidue.de, and Kruis, F.E.. 2017.
"A Monte Carlo method for the simulation of coagulation and nucleation based on weighted particles and the concepts of stochastic resolution and merging". United States.
doi:10.1016/J.JCP.2017.03.041.
@article{osti_22622302,
title = {A Monte Carlo method for the simulation of coagulation and nucleation based on weighted particles and the concepts of stochastic resolution and merging},
author = {Kotalczyk, G., Email: Gregor.Kotalczyk@unidue.de and Kruis, F.E.},
abstractNote = {Monte Carlo simulations based on weighted simulation particles can solve a variety of population balance problems and allow thus to formulate a solutionframework 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 nonweighted Monte Carlo particles. The tuning of the accuracy (named ‘stochastic resolution’ in this paper) of those transformations allows the construction of a constantnumber 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 constantnumber scheme: the lowweight 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 GPUbased simulation method which is validated by comparison with the discretesectional simulation technique. Two test models describing a constantrate nucleation coupled to a simultaneous coagulation in 1) the freemolecular regime or 2) the continuum regime are simulated for this purpose.},
doi = {10.1016/J.JCP.2017.03.041},
journal = {Journal of Computational Physics},
number = ,
volume = 340,
place = {United States},
year = 2017,
month = 7
}

This paper proposes a comprehensive framework for accelerating population balanceMonte Carlo (PBMC) simulation of particle coagulation dynamics. By combining Markov jump model, weighted majorant kernel and GPU (graphics processing unit) parallel computing, a significant gain in computational efficiency is achieved. The Markov jump model constructs a coagulationrule matrix of differentiallyweighted simulation particles, so as to capture the time evolution of particle size distribution with low statistical noise over the full size range and as far as possible to reduce the number of time loopings. Here three coagulation rules are highlighted and it is found that constructing appropriate coagulation rule providesmore »

Monte Carlo simulation of the nucleation and growth of binaryalloy particles of Au, Ag, and Pd on NaCl(100) substrates
The initial phase of thinfilm formation from two components that are simultaneously deposited on a substrate was simulated by Monte Carlo calculations of adatom processes, e.g., surface diffusion, reevaporation, nucleation, and cluster growth by collisions and capture. Experiments had revealed that the condensation of Pd is nearly complete on NaCl at 300 {degree}C, whereas Au and Ag exhibit very low initial condensation coefficients. Therefore, the composition of the nuclei differs strongly from that of the vapor beam. Monte Carlo calculations of the composition at the early stage of AuPd and of AgPd depositions were fitted to experimental data from earliermore » 
Girsanov's transformation based variance reduced Monte Carlo simulation schemes for reliability estimation in nonlinear stochastic dynamics
The study considers the problem of simulation based time variant reliability analysis of nonlinear randomly excited dynamical systems. Attention is focused on importance sampling strategies based on the application of Girsanov's transformation method. Controls which minimize the distance function, as in the first order reliability method (FORM), are shown to minimize a bound on the sampling variance of the estimator for the probability of failure. Two schemes based on the application of calculus of variations for selecting control signals are proposed: the first obtains the control force as the solution of a twopoint nonlinear boundary value problem, and, the secondmore » 
Quantitative Elemental Microanalysis Of Individual Particles With Use Of XRay Fluorescence Method And Monte Carlo Simulation
Recently a considerable interest has been triggered in the investigation of the composition of individual particles by Xray fluorescence microanalysis. The sources of these microsamples are mostly diversified. These samples come from space dust, air and ash, soil as well as environment and take the shape of a sphere or an oval. In analysis this kind of samples the geometrical effects caused by different sizes and shapes influence on accuracy of results. This fact arises from the matrix effect. For these samples it is not possible to find analytically a solution of equation taking into account an absorption of Xrays.more »