Coarse mesh finite difference acceleration of the random ray method

Abstract

This paper presents the development and evaluation of the Coarse Mesh Finite Difference (CMFD) acceleration method to The Random Ray Method (TRRM). We demonstrate its effectiveness in accelerating the convergence of the fission sources and reducing the real statistical error in neutron transport calculations. TRRM treats the angular variable of the neutron flux stochastically and performs batch sampling of characteristic rays to transform the Method of Characteristics (MOC) into a stochastic process, that is capable of making significant strides in memory efficiency and computational performance for some problems. Despite its advantages, TRRM exhibits a potential challenge with a large number of inactive cycles and inherent inter-cycle correlation much like the Monte Carlo method. To address this, the CMFD acceleration method is explored and demonstrated to dramatically reduce the number of required inactive cycles and diminish inter-cycle correlation. Results from the numerical analysis of a 2D C5G7 core problem indicate that the application of CMFD leads to enhanced convergence, with the integration of a CMFD acceleration step every cycle offering the most substantial reduction in statistical noise and error. The study reveals that applying CMFD with every cycle effectively resolves the issue of needing inactive cycles and significantly lowers the inter-cycle correlation, thereby providing a more accurate estimation of standard deviation for pin power distributions. We conclude that using CMFD not only minimizes the number of inactive cycles of TRRM – much like normal Monte Carlo transport – but also lowers real statistical error effectively. For a targeted maximum standard deviation of 0.1% in the pin power, the addition of CMFD can decrease the number of necessary active cycles by 41% compared to standard TRRM, as demonstrated by the 2D C5G7 benchmark analysis.