Predicting Correlation Coefficients for Monte Carlo Eigenvalue Simulations

Monte Carlo methods are most often considered as a reference for neutron transport simulations since very limited approximations are made abount nuclear data and system geometry. To report uncertainty of any tally evaluated as generation averages, the sample variance is divided by the number of acti...

Full description

Bibliographic Details
Main Authors: Miao, Jilang (Contributor), Forget, Benoit Robert Yves (Contributor), Smith, Kord S. (Contributor)
Other Authors: Massachusetts Institute of Technology. Department of Nuclear Science and Engineering (Contributor)
Format: Article
Language:English
Published: American Nuclear Society, 2017-04-10T18:08:46Z.
Subjects:
Online Access:Get fulltext
Description
Summary:Monte Carlo methods are most often considered as a reference for neutron transport simulations since very limited approximations are made abount nuclear data and system geometry. To report uncertainty of any tally evaluated as generation averages, the sample variance is divided by the number of active generations, which is based on the assumption that the neutron generations are independent. Correlation effects between neutrons in multiplying systems, particularly when performing power iteration to evaluate eigenvalues have been observed in previous work. Neglecting the correlation effect results in an underestimate of uncertainty reported by Monte Carlo calculations. Previous work has also proposed methods to predict the underestimation ratio. Yamamoto et al expanded the fission source distribution with diffusion equation modes, performed numerical simulation of the AR(autoregressive) process of the expansion coefficients and used the correlation of the AR process to predict that of the Monte Carlo eigenvalue simulation. Sutton applied the discretized phase space (DPS) approach to predict the underestimation ratio but the method cannot predict the ratio when one neutron generates offspring in different phase space regions or generates a random number of offspring. This paper presents a method to predict the correlation effect with the model of multitype branching processes (MBP). The method requires simulations for one generation of neutrons without knowing the source distribution and can predict the underestimation ratio for the cases where the traditional DPS approach does not work. The generation-to-generation correlation determines the convergence rate of active generations, the bias of variance estimator for each generation and the underestimation ratio of variance estimator for tallies averaged over active generations. The generation-to-generation correlation is characterized by the Auto-Correlation Coefficients (ACC) between tallies from different generations.
United States. Dept. of Energy (Consortium for Advanced Simulation of Light Water Reactors. Contract DE-AC05-00OR22725)