Abstract

In the Monte Carlo simulation of particle transport, and especially for shielding applications, variance reduction techniques are widely used to help simulate realisations of rare events and reduce the relative errors on the estimated scores for a given computation time. Adaptive Multilevel Splitting is one of these variance reduction techniques that has recently appeared in the literature. In the present paper, we propose an alternative version of the AMS algortihm, adapted for the first time to the field of particle tranport. Within this context, it can be used to build an unbiased estimator of any quantity associated with particle tracks, such as flux, reaction rates or even non-Boltzmann tallies. Furthermore, the effciency of the AMS algorithm is shown not to be very sensitive to variations of its input parameters, which makes it capable of significant variance reduction without requiring extended user effort.

Highlights

  • The challenge in using Monte Carlo particle transport simulations for shielding applications is to minimize the computation time required to attain a reasonable variance on the quantity of interest, called score

  • We present a specific mathematical setting of the AMS algorithm, which will be afterwards fairly easy to adapt to the context of particle transport

  • These results are from photon-only simulations, as the coupled photons/electrons/positrons AMS is not yet fully operational. These preliminary results show that the use of the Adaptive Multilevel Splitting algorithm in TRIPOLI-4 R reduces significatively the variance

Read more

Summary

Introduction

The challenge in using Monte Carlo particle transport simulations for shielding applications is to minimize the computation time required to attain a reasonable variance on the quantity of interest, called score. The basic approach of variance reduction techniques is to modify the simulation behaviour so as to increase rare events occurrence while keeping an unbiased estimator of the score. In this view, multilevel splitting techniques were introduced to the field of particle transport by Kahn and Harris [1]. The core of this work is presented, in which we introduce for the first time a pratical implementation of AMS within a Monte Carlo particle transport simulation. In the last Section of this paper, we present some of the results obtained using AMS within TRIPOLI-4 R These examples illustrate the validity of the algorithm as described, as well as the efficiency of AMS as a variance reduction technique These examples illustrate the validity of the algorithm as described in Section 4, as well as the efficiency of AMS as a variance reduction technique

Mathematical setting
Objective and setup
Importance function
Initialization
Resampling process
Final step
Interpretation of the replicas weights
About the number of resampled replicas
Neutron transport theory
Monte Carlo simulation of neutral particle transport
Definitions and notations
Importance map
The AMS algorithm
Final step and scoring
Number of contributions to the score
Robustness with regard to the importance function
Usability
Parallelization
Validation and results
Figure of merit
Importance map calculation
Bypass simulation
Neutron flux attenuation in water
Gamma spectroscopy
Findings
Conclusion
Full Text
Paper version not known

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call