1. |
- Donev, Aleksandar, et al.
(författare)
-
A First-Passage Kinetic Monte Carlo algorithm for complex diffusion-reaction systems
- 2010
-
Ingår i: Journal of Computational Physics. - : Elsevier BV. - 0021-9991 .- 1090-2716. ; 229:9, s. 3214-3236
-
Tidskriftsartikel (refereegranskat)abstract
- We develop an asynchronous event-driven First-Passage Kinetic Monte Carlo (FPKMC) algorithm for continuous time and space systems involving multiple diffusing and reacting species of spherical particles in two and three dimensions. The FPKMC algorithm presented here is based on the method introduced in Oppelstrup et al. [10] and is implemented in a robust and flexible framework. Unlike standard KMC algorithms such as the n-fold algorithm, FPKMC is most efficient at low densities where it replaces the many small hops needed for reactants to find each other with large first-passage hops sampled from exact time-dependent Green's functions, without sacrificing accuracy. We describe in detail the key components of the algorithm, including the event-loop and the sampling of first-passage probability distributions, and demonstrate the accuracy of the new method. We apply the FPKMC algorithm to the challenging problem of simulation of long-term irradiation of metals, relevant to the performance and aging of nuclear materials in current and future nuclear power plants. The problem of radiation damage spans many decades of time-scales, from picosecond spikes caused by primary cascades, to years of slow damage annealing and microstructure evolution. Our implementation of the FPKMC algorithm has been able to simulate the irradiation of a metal sample for durations that are orders of magnitude longer than any previous simulations using the standard Object KMC or more recent asynchronous algorithms. (C) 2010 Elsevier Inc. All rights reserved.
|
|
2. |
- Oppelstrup, Tomas, 1977-, et al.
(författare)
-
First-Passage Kinetic Monte Carlo method
- 2009
-
Ingår i: Physical Review E. Statistical, Nonlinear, and Soft Matter Physics. - 1063-651X .- 1095-3787. ; 80:6, s. 066701-
-
Tidskriftsartikel (refereegranskat)abstract
- We present a new effcient method for Monte Carlo simulations of diusion-reaction processes. First introduced by us in [Phys. Rev. Lett., 97:230602, 2006], the new algorithm skips the traditionalsmall diusion hops and propagates the diusing particles over long distances through asequence of super-hops, one particle at a time. By partitioning the simulation space into nonoverlappingprotecting domains each containing only one or two particles, the algorithm factorizesthe N-body problem of collisions among multiple Brownian particles into a set of much simplersingle-body and two-body problems. Ecient propagation of particles inside their protective domainsis enabled through the use of time-dependent Green's functions (propagators) obtained assolutions for the rst-passage statistics of random walks. The resulting Monte Carlo algorithmis event-driven and asynchronous; each Brownian particle propagates inside its own protectivedomain and on its own time clock. The algorithm reproduces the statistics of the underlyingMonte-Carlo model exactly. Extensive numerical examples demonstrate that for an importantclass of diusion-reaction models the new algorithm is effcient at low particle densities, whereother existing algorithms slow down severely.
|
|
3. |
- Oppelstrup, Tomas, 1977-, et al.
(författare)
-
First-Passage Monte Carlo Algorithm: Diffusion without All the Hops
- 2006
-
Ingår i: Physical Review Letters. - 0031-9007 .- 1079-7114. ; 97, s. 230602-
-
Tidskriftsartikel (refereegranskat)abstract
- We present a novel Monte Carlo algorithm for N diffusing finite particles that react on collisions. Using the theory of first-passage processes and time dependent Green’s functions, we break the difficult N-body problem into independent single- and two-body propagations circumventing numerous diffusion hopsused in standard Monte Carlo simulations. The new algorithm is exact, extremely efficient, and applicablet o many important physical situations in arbitrary integer dimensions.
|
|
4. |
|
|
5. |
- Litvinenko, A., et al.
(författare)
-
Computation of electromagnetic fields scattered from objects with uncertain shapes using multilevel monte carlo method
- 2019
-
Ingår i: IEEE Journal on Multiscale and Multiphysics Computational Techniques. - : Institute of Electrical and Electronics Engineers Inc.. - 2379-8793. ; 4, s. 51-64
-
Tidskriftsartikel (refereegranskat)abstract
- Computational tools for characterizing electromagnetic scattering from objects with uncertain shapes are needed in various applications ranging from remote sensing at microwave frequencies to Raman spectroscopy at optical frequencies. Often, such computational tools use the Monte Carlo (MC) method to sample a parametric space describing geometric uncertainties. For each sample, which corresponds to a realization of the geometry, a deterministic electromagnetic solver computes the scattered fields. However, for an accurate statistical characterization, the number of MC samples has to be large. In this paper, to address this challenge, the continuation multilevel Monte Carlo (CMLMC) method is used together with a surface integral equation solver. The CMLMC method optimally balances statistical errors due to sampling of the parametric space, and numerical errors due to the discretization of the geometry using a hierarchy of discretizations, from coarse to fine. The number of realizations of finer discretizations can be kept low, with most samples computed on coarser discretizations to minimize computational cost. Consequently, the total execution time is significantly reduced, in comparison to the standard MC scheme.
|
|
6. |
- Zepeda-Ruiz, L. A., et al.
(författare)
-
Atomistic simulations of grain boundary pinning in CuFe alloys
- 2005
-
Ingår i: Applied Physics Letters. - : AIP Publishing. - 0003-6951 .- 1077-3118. ; 87:23, s. 231904-
-
Tidskriftsartikel (refereegranskat)abstract
- We apply a hybrid Monte Carlo-molecular dynamics code to the study of grain boundary motion upon annealing of pure Cu and Cu with low concentrations of Fe. The hybrid simulations account for segregation and precipitation of the low solubility Fe, together with curvature-driven grain boundary motion. Grain boundaries in two different systems, a Sigma 7+U-shaped half-loop grain and a nanocrystalline sample, were found to be pinned in the presence of Fe concentrations exceeding 3%.
|
|