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.
10.1051_epjconf_202430209006.pdf
Main Document
http://purl.org/coar/version/c_970fb48d4fbd8a85
openaccess
CC BY
1.3 MB
Adobe PDF
62a6ce7e7fc748f4f3994931b354f3ad