Variance Estimation in Monte Carlo Eigenvalue Simulations Using Spectral Analysis Method
Monte Carlo (MC) simulation is widely used to solve the eigenvalue form of the Boltzmann transport equation that mathematically represents the neutron transport process through complex multiplying (fissionable) systems. Monte Carlo eigenvalue simulation starts with an assumed fission source distribution and uses the fission sites from the previous iteration (cycle) as the starting source in the current iteration. Important system parameters (MC tallies) such as fuel pinpower distribution are estimated over several cycles after the convergence of the fission source distribution to a stationary distribution. However, the MC fission source iteration algorithm that uses fission source sites from the previous cycle introduces a cycletocycle correlation. Monte Carlo simulations that do not account for the cycletocycle correlation systematically underestimate the variance of the estimated system parameters (sample mean). This paper presents the relationship between the spectral density in the frequency domain at frequency zero and the variance of the sample mean. This paper introduces a novel method in the frequency domain for the MC variance estimation. For the three test problems used in this paper, researchers have observed that the new method results in an improvement of more than one order of magnitude to the standard deviation of the sample mean.more »
 Authors:

^{[1]};
^{[2]}
 Texas A&M Univ.–Corpus Christi, Corpus Christi, TX (United States)
 Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States)
 Publication Date:
 Grant/Contract Number:
 AC0500OR22725
 Type:
 Accepted Manuscript
 Journal Name:
 Nuclear Science and Engineering
 Additional Journal Information:
 Journal Volume: 191; Journal Issue: 3; Journal ID: ISSN 00295639
 Publisher:
 American Nuclear Society  Taylor & Francis
 Research Org:
 Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States)
 Sponsoring Org:
 USDOE
 Country of Publication:
 United States
 Language:
 English
 Subject:
 73 NUCLEAR PHYSICS AND RADIATION PHYSICS; Monte Carlo; variance estimation; periodogram; spectral density
 OSTI Identifier:
 1459299
Jin, Lei, and Banerjee, Kaushik. Variance Estimation in Monte Carlo Eigenvalue Simulations Using Spectral Analysis Method. United States: N. p.,
Web. doi:10.1080/00295639.2018.1471269.
Jin, Lei, & Banerjee, Kaushik. Variance Estimation in Monte Carlo Eigenvalue Simulations Using Spectral Analysis Method. United States. doi:10.1080/00295639.2018.1471269.
Jin, Lei, and Banerjee, Kaushik. 2018.
"Variance Estimation in Monte Carlo Eigenvalue Simulations Using Spectral Analysis Method". United States.
doi:10.1080/00295639.2018.1471269.
@article{osti_1459299,
title = {Variance Estimation in Monte Carlo Eigenvalue Simulations Using Spectral Analysis Method},
author = {Jin, Lei and Banerjee, Kaushik},
abstractNote = {Monte Carlo (MC) simulation is widely used to solve the eigenvalue form of the Boltzmann transport equation that mathematically represents the neutron transport process through complex multiplying (fissionable) systems. Monte Carlo eigenvalue simulation starts with an assumed fission source distribution and uses the fission sites from the previous iteration (cycle) as the starting source in the current iteration. Important system parameters (MC tallies) such as fuel pinpower distribution are estimated over several cycles after the convergence of the fission source distribution to a stationary distribution. However, the MC fission source iteration algorithm that uses fission source sites from the previous cycle introduces a cycletocycle correlation. Monte Carlo simulations that do not account for the cycletocycle correlation systematically underestimate the variance of the estimated system parameters (sample mean). This paper presents the relationship between the spectral density in the frequency domain at frequency zero and the variance of the sample mean. This paper introduces a novel method in the frequency domain for the MC variance estimation. For the three test problems used in this paper, researchers have observed that the new method results in an improvement of more than one order of magnitude to the standard deviation of the sample mean. The new method also compares favorably with the previously introduced batch, bootstrap, and covarianceadjusted methods when applied to the three test problems investigated in this paper. This new method does not require modification of the MC eigenvalue algorithm (power iteration), is code agnostic, and is therefore easy to use when implementing in any existing MC code. In conclusion, the new estimate can be calculated without saving tally results of all active/stationary cycles.},
doi = {10.1080/00295639.2018.1471269},
journal = {Nuclear Science and Engineering},
number = 3,
volume = 191,
place = {United States},
year = {2018},
month = {6}
}