Towards a highly efficient and unbiased population-control algorithm for kinetic Monte Carlo simulations

Population-control methods are key to non-stationary Monte Carlo simulations of multiplying systems: they prevent either the unbounded growth or the disappearance of neutrons, occurring respectively in supercritical and subcritical conditions; furthermore, they contribute to an efficient allocation of computational resources by addressing the unbalance between the neutron and the precursor populations. In this paper, we present two alternative population-control algorithms: the legacy implementation in TRIPOLI-4®, the Monte Carlo code developed at CEA, and an improved version that is currently under investigation, based on the use of a simplified point-kinetics solver. We assess the performance of these methods through the simulation of a $2.2 step reactivity insertion in a fast system (Flattop-Pu), leading to an increase of the neutron population by a factor 200, which is benchmarked against point kinetics. We show that the new implementation not only suppresses the slight bias that was present in the legacy method due to a stochastic normalization factor, but also outperforms the previous algorithm in terms of variance reduction and improvement of the figure of merit.