Optimisation of criticality and burnup calculations in monk®

The primary goal of this paper is to increase the efficiency of criticality and burnup calculations in the ANSWERS MONK® Monte Carlo code [1]. Two ways of achieving this goal are investigated as part of the H2020 McSAFE Project: creating a unified energy grid for all materials in the model, and reducing the spread in variances of fluxes for depletable materials using a generated optimised importance map. The average tracking speedup factor across all cycles of all burnup calculations ran using the unified energy grid, at base temperature, was found to be 1.96. For criticality calculations at 400K with runtime Doppler broadening, the unified grid approach gave a total speedup factor of 7.32. This demonstrates the potential importance of this method to reduce the calculation time with models with runtime Doppler broadening. The use of the generated optimised importance map has been demonstrated to significantly reduce the variance in the standard deviations on the fluxes in the fuel pins across two different test cases. If a solution is required in which the standard deviation in none of the fuel pins exceeds 5% it was found that the number of scoring stages required was more than halved, highlighting the potential for the outlined methodology to speedup burnup credit calculations.