Abstract
Iteratively solved Monte Carlo (MC) codes are frequently used for plasma edge simulations. However, their accuracy and convergence assessment are still unresolved issues. In analogy with the error classification recently developed for coupled finite-volume/Monte Carlo (FV-MC) codes, we define different error contributions and analyse them separately in a simplified non-linear MC code. Three iterative procedures are examined: Random Noise (RN), where different seeds are used in each iteration; Correlated Sampling, where particle trajectories remain correlated between iterations; and Robbins Monro, where averaging is used during the simulation. We show that, as in FV-MC codes, RN is the most efficient iterative procedure provided averaging is used to decrease the statistical error. In addition, we conclude that the accuracy can be assessed using the same techniques as in FV-MC codes.
| Original language | English |
|---|---|
| Pages (from-to) | 652-658 |
| Number of pages | 7 |
| Journal | Contributions to Plasma Physics |
| Volume | 58 |
| Issue number | 6-8 |
| DOIs | |
| State | Published - Jul 1 2018 |
| Externally published | Yes |
Funding
This work was sponsored by Flanders Innovation and Entrepreneurship (IWT.141064). 1Department of Mechanical Engineering, KU Leuven, Leuven, Belgium 2Department of Computer Science, KU Leuven, Leuven, Belgium *Correspondence Kristel Ghoos, Department of Mechanical Engineering, KU Leuven, Celestijnenlaan 300A, 3001 Leuven, Belgium. Email: [email protected] Funding Information This research was supported by the Flanders Innovation and Entrepreneurship, IWT.141064.
Keywords
- Monte Carlo
- accuracy
- convergence
- numerical errors