tailieunhanh - Adaptive multilevel splitting for Monte Carlo particle transport

In the present paper, we propose an alternative version of the AMS algorithm, adapted for the first time to the field of particle transport. | Adaptive multilevel splitting for Monte Carlo particle transport EPJ Nuclear Sci. Technol. 3 29 2017 Nuclear Sciences H. Louvin et al. published by EDP Sciences 2017 amp Technologies DOI epjn 2017022 Available online at http REGULAR ARTICLE Adaptive multilevel splitting for Monte Carlo particle transport Henri Louvin1 Eric Dumonteil2 Tony Lelièvre3 Mathias Rousset3 and Cheikh M. Diop1 1 CEA Saclay DEN DM2S SERMA LTSD 91191 Gif-sur-Yvette France 2 IRSN PSN-EXP SNC 92262 Fontenay-aux-Roses France 3 Université Paris-Est CERMICS ENPC INRIA 77455 Marne-la-Vallée France Received 6 January 2017 Received in final form 18 May 2017 Accepted 30 August 2017 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 AMS 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 algorithm adapted for the first time to the field of particle transport. 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 like pulse-height tallies and other spectra. Furthermore the efficiency 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. 1 Introduction general setting by Bréhier et al. 3 . This method also aims to duplicate the interesting particles of the simulation but The challenge in using Monte Carlo particle transport does not use an a priori definition of importance regions. simulations for shielding applications is to minimize the Instead the splitting levels are .