Unbiasedness and optimization of regional weight cancellation.

ABSTRACT

The Monte Carlo method is often used to simulate systems which can be modeled by random walks. To calculate observables, in many implementations the "walkers" carry a statistical weight which is generally assumed to be positive. Some random walk simulations, however, may require walkers to have positive or negative weights it has been shown that the presence of a mixture of positive and negative weights can impede the statistical convergence, and special weight-cancellation techniques must be adopted to overcome these issues. In a recent work we demonstrated the usefulness of one such method, exact regional weight cancellation, to solve eigenvalue problems in nuclear reactor physics in three spatial dimensions. The method previously exhibited had several limitations (including multigroup transport and isotropic scattering) and needed homogeneous cuboid cancellation regions. In this paper we lift the previous limitations, in view of applying exact regional cancellation to more realistic continuous-energy neutron transport problems. This extended regional cancellation framework is used to optimize the efficiency of the weight cancellation. Our findings are illustrated on a benchmark configuration for reactor physics.