skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: 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.

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., E-mail: Gregor.Kotalczyk@uni-due.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., E-mail: Gregor.Kotalczyk@uni-due.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., E-mail: Gregor.Kotalczyk@uni-due.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., E-mail: Gregor.Kotalczyk@uni-due.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 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.},
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 balance-Monte 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 coagulation-rule matrix of differentially-weighted 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 » a route to attain the compromise between accuracy and cost of PBMC methods. Further, in order to avoid double looping over all simulation particles when considering the two-particle events (typically, particle coagulation), the weighted majorant kernel is introduced to estimate the maximum coagulation rates being used for acceptance–rejection processes by single-looping over all particles, and meanwhile the mean time-step of coagulation event is estimated by summing the coagulation kernels of rejected and accepted particle pairs. The computational load of these fast differentially-weighted PBMC simulations (based on the Markov jump model) is reduced greatly to be proportional to the number of simulation particles in a zero-dimensional system (single cell). Finally, for a spatially inhomogeneous multi-dimensional (multi-cell) simulation, the proposed fast PBMC is performed in each cell, and multiple cells are parallel processed by multi-cores on a GPU that can implement the massively threaded data-parallel tasks to obtain remarkable speedup ratio (comparing with CPU computation, the speedup ratio of GPU parallel computing is as high as 200 in a case of 100 cells with 10 000 simulation particles per cell). These accelerating approaches of PBMC are demonstrated in a physically realistic Brownian coagulation case. The computational accuracy is validated with benchmark solution of discrete-sectional method. The simulation results show that the comprehensive approach can attain very favorable improvement in cost without sacrificing computational accuracy.« less
  • The initial phase of thin-film 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 Au-Pd and of Ag-Pd depositions were fitted to experimental data from earliermore » work, which allowed us to determine the differences of the atomic energy parameters of the components, e.g., of the energies of adsorption {ital E}{sub {ital a}} and of diffusion {ital E}{sub {ital d}} with high accuracy. These differences are ({ital E}{sub {ital a}}{minus}{ital E}{sub {ital d}}){sub Pd}{minus}({ital E}{sub {ital a}}{minus}{ital E}{sub {ital d}}){sub Au}=0.12{plus minus}0.03 eV, and ({ital E}{sub {ital a}}{minus}{ital E}{sub {ital d}}){sub Pd}{minus}({ital E}{sub {ital a}}{minus}{ital E}{sub {ital d}}){sub Ag}=0.25{plus minus}0.05 eV.« less
  • 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 two-point nonlinear boundary value problem, and, the secondmore » explores the application of the Volterra series in characterizing the controls. The relative merits of these schemes, vis-à-vis the method based on ideas from the FORM, are discussed. Illustrative examples, involving archetypal single degree of freedom (dof) nonlinear oscillators, and a multi-degree of freedom nonlinear dynamical system, are presented. The credentials of the proposed procedures are established by comparing the solutions with pertinent results from direct Monte Carlo simulations. - Highlights: • The distance minimizing control forces minimize a bound on the sampling variance. • Establishing Girsanov controls via solution of a two-point boundary value problem. • Girsanov controls via Volterra's series representation for the transfer functions.« less
  • Recently a considerable interest has been triggered in the investigation of the composition of individual particles by X-ray fluorescence microanalysis. The sources of these micro-samples 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 X-rays.more » Hence, a way out is to approximate the real sample shape with the other one or to use Monte Carlo (MC) simulation method. In current work authors utilized the iterative MC simulation to assess an elemental percentage of individual particles. The set of glass micro-spheres, made of NIST K3089 material of known chemical composition, with diameters in the range between 25 and 45 {mu}m was investigated. The microspheres were scanned with X-ray tube primary radiation. Results of MC simulation were compared with these of some analytical approaches based on particle shape approximation. An investigation showed that the low-Z elements (Si, Ca, Ti) were the most sensitive on changes of particle shape and sizes. For high-Z elements (Fe--Pb) concentrations were nearly equal regardless of method used. However, for the all elements considered, results of MC simulation were more accurate then these of analytical relationships taken into comparison.« less