TOWARDS ZERO-VARIANCE SCHEMES FOR KINETIC MONTE-CARLO SIMULATIONS

DEN-Service d’études des réacteurs et de mathématiques appliquées (SERMA), CEA, Université Paris-Saclay, F-91191, Gif-sur-Yvette, France

davide.mancusi@cea.fr
andrea.zoia@cea.fr

Published online: 22 February 2021

Abstract

The solution of the time-dependent transport problem for neutrons and precursors in a nuclear reactor is hard to treat in a naive Monte-Carlo framework because of the largely different time scales associated to the prompt-fission chains and to the decay of precursors. The increasing computer power and the development of variance-reduction techniques specific for reactor kinetics have recently unlocked the possibility to calculate reference solutions to the time-dependent transport problem. However, the application of time-dependent Monte Carlo to large systems (i.e., a full reactor core) is still stifled by the enormous computational requirements. In this paper, we formulate the construction of an optimal Monte-Carlo strategy (in the sense that it results in a zero-variance estimator) for a specific observable in time-dependent transport, in analogy with the existing schemes for stationary problems. As far as we are aware, zero-variance Monte-Carlo schemes for neutron-precursor kinetics have never been proposed before. We verify our construction with numerical calculations for a benchmark transport problem.