DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Singular value decomposition of adjoint flux distributions for Monte Carlo variance reduction

Abstract

Monte Carlo (MC) shielding calculations often use weight windows (WWs) and biased sources formed from a deterministic estimate of the adjoint flux to improve the convergence rate of tallies. This requires a significant amount of computer memory, which can limit the memory available for high-resolution tally output. A new method is proposed for reducing these memory requirements by using singular value decomposition (SVD) in linear or logarithmic space to approximate the adjoint flux. This method’s performance is evaluated using the Shift and Denovo codes for streaming and diffusion base case problems, followed by problems using the Westinghouse AP1000 and the Joint European Torus. The log SVD reduced WW memory requirements by an order of magnitude in all cases without a significant performance penalty. Additionally, the linear SVD reduced biased source memory requirements by an order of magnitude, but further investigation is needed to account for observed limitations.

Authors:
ORCiD logo [1]; ORCiD logo [1]; ORCiD logo [1]; ORCiD logo [1]
  1. Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States)
Publication Date:
Research Org.:
Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States)
Sponsoring Org.:
USDOE
OSTI Identifier:
1606986
Grant/Contract Number:  
AC05-00OR22725
Resource Type:
Accepted Manuscript
Journal Name:
Annals of Nuclear Energy (Oxford)
Additional Journal Information:
Journal Name: Annals of Nuclear Energy (Oxford); Journal Volume: 141; Journal Issue: C; Journal ID: ISSN 0306-4549
Publisher:
Elsevier
Country of Publication:
United States
Language:
English
Subject:
73 NUCLEAR PHYSICS AND RADIATION PHYSICS; Monte Carlo radiation transport; Singular value decomposition; Variance reduction

Citation Formats

Biondo, Elliott, Evans, Thomas, Davidson, Gregory, and Hamilton, Steven P. Singular value decomposition of adjoint flux distributions for Monte Carlo variance reduction. United States: N. p., 2020. Web. doi:10.1016/j.anucene.2020.107327.
Biondo, Elliott, Evans, Thomas, Davidson, Gregory, & Hamilton, Steven P. Singular value decomposition of adjoint flux distributions for Monte Carlo variance reduction. United States. https://doi.org/10.1016/j.anucene.2020.107327
Biondo, Elliott, Evans, Thomas, Davidson, Gregory, and Hamilton, Steven P. Tue . "Singular value decomposition of adjoint flux distributions for Monte Carlo variance reduction". United States. https://doi.org/10.1016/j.anucene.2020.107327. https://www.osti.gov/servlets/purl/1606986.
@article{osti_1606986,
title = {Singular value decomposition of adjoint flux distributions for Monte Carlo variance reduction},
author = {Biondo, Elliott and Evans, Thomas and Davidson, Gregory and Hamilton, Steven P.},
abstractNote = {Monte Carlo (MC) shielding calculations often use weight windows (WWs) and biased sources formed from a deterministic estimate of the adjoint flux to improve the convergence rate of tallies. This requires a significant amount of computer memory, which can limit the memory available for high-resolution tally output. A new method is proposed for reducing these memory requirements by using singular value decomposition (SVD) in linear or logarithmic space to approximate the adjoint flux. This method’s performance is evaluated using the Shift and Denovo codes for streaming and diffusion base case problems, followed by problems using the Westinghouse AP1000 and the Joint European Torus. The log SVD reduced WW memory requirements by an order of magnitude in all cases without a significant performance penalty. Additionally, the linear SVD reduced biased source memory requirements by an order of magnitude, but further investigation is needed to account for observed limitations.},
doi = {10.1016/j.anucene.2020.107327},
journal = {Annals of Nuclear Energy (Oxford)},
number = C,
volume = 141,
place = {United States},
year = {Tue Jan 21 00:00:00 EST 2020},
month = {Tue Jan 21 00:00:00 EST 2020}
}

Works referenced in this record:

Monte Carlo variance reduction with deterministic importance functions
journal, January 2003


Implementation, capabilities, and benchmarking of Shift, a massively parallel Monte Carlo radiation transport code
journal, March 2016

  • Pandya, Tara M.; Johnson, Seth R.; Evans, Thomas M.
  • Journal of Computational Physics, Vol. 308
  • DOI: 10.1016/j.jcp.2015.12.037

Efficient solution of the simplified P N equations
journal, March 2015


Reactor noise analysis based on the singular value decomposition (SVD)
journal, August 1998


Efficient Subspace Methods-Based Algorithms for Performing Sensitivity, Uncertainty, and Adaptive Simulation of Large-Scale Computational Models
journal, July 2008

  • Abdel-Khalik, Hany S.; Turinsky, Paul J.; Jessee, Matthew A.
  • Nuclear Science and Engineering, Vol. 159, Issue 3
  • DOI: 10.13182/NSE159-256

Singular value decomposition for genome-wide expression data processing and modeling
journal, August 2000

  • Alter, O.; Brown, P. O.; Botstein, D.
  • Proceedings of the National Academy of Sciences, Vol. 97, Issue 18, p. 10101-10106
  • DOI: 10.1073/pnas.97.18.10101

Comparison of global variance reduction techniques for Monte Carlo radiation transport simulations of ITER
journal, October 2011


Review of Hybrid Methods for Deep-Penetration Neutron Transport
journal, April 2019


Nuclide depletion capabilities in the Shift Monte Carlo code
journal, April 2018


Denovo: A New Three-Dimensional Parallel Discrete Ordinates Code in SCALE
journal, August 2010

  • Evans, Thomas M.; Stafford, Alissa S.; Slaybaugh, Rachel N.
  • Nuclear Technology, Vol. 171, Issue 2
  • DOI: 10.13182/NT171-171

Optimization of processor allocation for domain decomposed Monte Carlo calculations
journal, September 2019


Shutdown dose rate analysis with CAD geometry, Cartesian/tetrahedral mesh, and advanced variance reduction
journal, May 2016


Importance Estimation in Forward Monte Carlo Calculations
journal, January 1984

  • Booth, Thomas E.; Hendricks, John S.
  • Nuclear Technology - Fusion, Vol. 5, Issue 1
  • DOI: 10.13182/FST84-A23082

MC21 v.6.0 – A continuous-energy Monte Carlo particle transport code with integrated reactor feedback capabilities
journal, August 2015


ITER oriented neutronics benchmark experiments on neutron streaming and shutdown dose rate at JET
journal, November 2017


Singular Value Decomposition (SVD) Image Coding
journal, April 1976


Automated Weight Windows for Global Monte Carlo Particle Transport Calculations
journal, January 2001

  • Cooper, Marc A.; Larsen, Edward W.
  • Nuclear Science and Engineering, Vol. 137, Issue 1
  • DOI: 10.13182/NSE00-34

3D radiation transport benchmark problems and results for simple geometries with void region
journal, January 2001


Estimating the parameters of exponentially damped sinusoids and pole-zero modeling in noise
journal, December 1982

  • Kumaresan, R.; Tufts, D.
  • IEEE Transactions on Acoustics, Speech, and Signal Processing, Vol. 30, Issue 6
  • DOI: 10.1109/TASSP.1982.1163974

FW-CADIS Method for Global and Regional Variance Reduction of Monte Carlo Radiation Transport Calculations
journal, January 2014

  • Wagner, John C.; Peplow, Douglas E.; Mosher, Scott W.
  • Nuclear Science and Engineering, Vol. 176, Issue 1
  • DOI: 10.13182/NSE12-33

Analysis for the Effect of Spatial Discretization Method on AP1000 Reactor Pressure Vessel Fluence Calculation
journal, January 2016

  • Zheng, Junxiao; Zhang, Bin; Shi, Shengchun
  • Science and Technology of Nuclear Installations, Vol. 2016
  • DOI: 10.1155/2016/3461290

Continuous-energy Monte Carlo neutron transport on GPUs in the Shift code
journal, June 2019


Integration of the Full Tokamak Reference Model with the Complex Model for ITER Neutronic Analysis
journal, August 2018