Skip to main content
U.S. Department of Energy
Office of Scientific and Technical Information

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

Journal Article · · Annals of Nuclear Energy (Oxford)
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.
Research Organization:
Oak Ridge National Laboratory (ORNL), Oak Ridge, TN (United States)
Sponsoring Organization:
USDOE
Grant/Contract Number:
AC05-00OR22725
OSTI ID:
1606986
Journal Information:
Annals of Nuclear Energy (Oxford), Journal Name: Annals of Nuclear Energy (Oxford) Journal Issue: C Vol. 141; ISSN 0306-4549
Publisher:
ElsevierCopyright Statement
Country of Publication:
United States
Language:
English

References (24)

The approximation of one matrix by another of lower rank journal September 1936
3D radiation transport benchmark problems and results for simple geometries with void region journal January 2001
Monte Carlo variance reduction with deterministic importance functions journal January 2003
Reactor noise analysis based on the singular value decomposition (SVD) journal August 1998
MC21 v.6.0 – A continuous-energy Monte Carlo particle transport code with integrated reactor feedback capabilities journal August 2015
Nuclide depletion capabilities in the Shift Monte Carlo code journal April 2018
Continuous-energy Monte Carlo neutron transport on GPUs in the Shift code journal June 2019
Comparison of global variance reduction techniques for Monte Carlo radiation transport simulations of ITER journal October 2011
Shutdown dose rate analysis with CAD geometry, Cartesian/tetrahedral mesh, and advanced variance reduction journal May 2016
ITER oriented neutronics benchmark experiments on neutron streaming and shutdown dose rate at JET journal November 2017
Efficient solution of the simplified P N equations journal March 2015
Implementation, capabilities, and benchmarking of Shift, a massively parallel Monte Carlo radiation transport code journal March 2016
Optimization of processor allocation for domain decomposed Monte Carlo calculations journal September 2019
Singular value decomposition for genome-wide expression data processing and modeling journal August 2000
Review of Hybrid Methods for Deep-Penetration Neutron Transport journal April 2019
Integration of the Full Tokamak Reference Model with the Complex Model for ITER Neutronic Analysis journal August 2018
Estimating the parameters of exponentially damped sinusoids and pole-zero modeling in noise journal December 1982
Singular Value Decomposition (SVD) Image Coding journal April 1976
Analysis for the Effect of Spatial Discretization Method on AP1000 Reactor Pressure Vessel Fluence Calculation journal January 2016
Importance Estimation in Forward Monte Carlo Calculations journal January 1984
Automated Weight Windows for Global Monte Carlo Particle Transport Calculations journal January 2001
FW-CADIS Method for Global and Regional Variance Reduction of Monte Carlo Radiation Transport Calculations journal January 2014
Efficient Subspace Methods-Based Algorithms for Performing Sensitivity, Uncertainty, and Adaptive Simulation of Large-Scale Computational Models journal July 2008
Denovo: A New Three-Dimensional Parallel Discrete Ordinates Code in SCALE journal August 2010