Capturing and Utilizing the Random Feature in Monte Carlo Fission Source Distributions

Pengfei Shen,Xiaodong Huo,Shanfang Huang,Yuchuan Guo,Zeng Shao,Haifeng Yang,Kan Wang

DOI: https://doi.org/10.1016/j.anucene.2022.109468

IF: 1.9

2023-01-01

Annals of Nuclear Energy

Abstract:The Monte Carlo algorithm commonly adopts the power iteration (PI) process in criticality calculations. Previous studies have shown that random error term in the PI process significantly influences the source convergence diagnosis and variance underestimation phenomenon. This paper uses the Sliced Wasserstein (SW) distance of the fission source distributions (FSDs) to estimate the error terms in the PI process. Three mesh-free SW-based methods are proposed to capture the random error term's 1/root m feature for the OECD source convergence fissile slab model, sphere array model, and the BEAVRS model. Then, a mesh-free source convergence auto-diagnosis algorithm is introduced based on the methods. The algorithm sets the batch size at five stages. The 1/root m feature of the random error term is utilized to achieve auto-diagnosis at the target batch size. Numerical calculation results show that the auto-diagnosis algorithm is practical and efficient for source convergence diagnosis and acceleration, which can automatically save 58 % to 81 % of the computational time to reach a converged fission source in models with the dominance ratio of around 0.99.