skip to main content
DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Monte Carlo Perturbation Theory Estimates of Sensitivities to System Dimensions

Abstract

Here, Monte Carlo methods are developed using adjoint-based perturbation theory and the differential operator method to compute the sensitivities of the k-eigenvalue, linear functions of the flux (reaction rates), and bilinear functions of the forward and adjoint flux (kinetics parameters) to system dimensions for uniform expansions or contractions. The calculation of sensitivities to system dimensions requires computing scattering and fission sources at material interfaces using collisions occurring at the interface—which is a set of events with infinitesimal probability. Kernel density estimators are used to estimate the source at interfaces using collisions occurring near the interface. The methods for computing sensitivities of linear and bilinear ratios are derived using the differential operator method and adjoint-based perturbation theory and are shown to be equivalent to methods previously developed using a collision history–based approach. The methods for determining sensitivities to system dimensions are tested on a series of fast, intermediate, and thermal critical benchmarks as well as a pressurized water reactor benchmark problem with iterated fission probability used for adjoint-weighting. The estimators are shown to agree within 5% and 3σ of reference solutions obtained using direct perturbations with central differences for the majority of test problems.

Authors:
ORCiD logo [1]; ORCiD logo [1]
  1. Univ. of Michigan, Ann Arbor, MI (United States)
Publication Date:
Research Org.:
Univ. of Michigan, Ann Arbor, MI (United States)
Sponsoring Org.:
USDOE National Nuclear Security Administration (NNSA), Office of Defense Nuclear Nonproliferation (NA-20); US Nuclear Regulatory Commission (NRC), Rockville, MD (United States)
OSTI Identifier:
1438024
Grant/Contract Number:  
NA0002576; NRC-HQ-13-G-38- 0007
Resource Type:
Accepted Manuscript
Journal Name:
Nuclear Science and Engineering
Additional Journal Information:
Journal Volume: 189; Journal Issue: 3; Journal ID: ISSN 0029-5639
Publisher:
American Nuclear Society - Taylor & Francis
Country of Publication:
United States
Language:
English
Subject:
97 MATHEMATICS AND COMPUTING; 70 PLASMA PHYSICS AND FUSION TECHNOLOGY; Monte Carlo; system dimensions; sensitivities

Citation Formats

Burke, Timothy P., and Kiedrowski, Brian C. Monte Carlo Perturbation Theory Estimates of Sensitivities to System Dimensions. United States: N. p., 2017. Web. doi:10.1080/00295639.2017.1388093.
Burke, Timothy P., & Kiedrowski, Brian C. Monte Carlo Perturbation Theory Estimates of Sensitivities to System Dimensions. United States. doi:10.1080/00295639.2017.1388093.
Burke, Timothy P., and Kiedrowski, Brian C. Mon . "Monte Carlo Perturbation Theory Estimates of Sensitivities to System Dimensions". United States. doi:10.1080/00295639.2017.1388093. https://www.osti.gov/servlets/purl/1438024.
@article{osti_1438024,
title = {Monte Carlo Perturbation Theory Estimates of Sensitivities to System Dimensions},
author = {Burke, Timothy P. and Kiedrowski, Brian C.},
abstractNote = {Here, Monte Carlo methods are developed using adjoint-based perturbation theory and the differential operator method to compute the sensitivities of the k-eigenvalue, linear functions of the flux (reaction rates), and bilinear functions of the forward and adjoint flux (kinetics parameters) to system dimensions for uniform expansions or contractions. The calculation of sensitivities to system dimensions requires computing scattering and fission sources at material interfaces using collisions occurring at the interface—which is a set of events with infinitesimal probability. Kernel density estimators are used to estimate the source at interfaces using collisions occurring near the interface. The methods for computing sensitivities of linear and bilinear ratios are derived using the differential operator method and adjoint-based perturbation theory and are shown to be equivalent to methods previously developed using a collision history–based approach. The methods for determining sensitivities to system dimensions are tested on a series of fast, intermediate, and thermal critical benchmarks as well as a pressurized water reactor benchmark problem with iterated fission probability used for adjoint-weighting. The estimators are shown to agree within 5% and 3σ of reference solutions obtained using direct perturbations with central differences for the majority of test problems.},
doi = {10.1080/00295639.2017.1388093},
journal = {Nuclear Science and Engineering},
number = 3,
volume = 189,
place = {United States},
year = {2017},
month = {12}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

Citation Metrics:
Cited by: 3 works
Citation information provided by
Web of Science

Save / Share:

Works referenced in this record:

SCALE Continuous-Energy Eigenvalue Sensitivity Coefficient Calculations
journal, March 2016

  • Perfetti, Christopher M.; Rearden, Bradley T.; Martin, William R.
  • Nuclear Science and Engineering, Vol. 182, Issue 3
  • DOI: 10.13182/NSE15-12

Development of a Generalized Perturbation Theory Method for Sensitivity Analysis Using Continuous-Energy Monte Carlo Methods
journal, March 2016

  • Perfetti, Christopher M.; Rearden, Bradley T.
  • Nuclear Science and Engineering, Vol. 182, Issue 3
  • DOI: 10.13182/NSE15-13

Non-Parametric Estimation of a Multivariate Probability Density
journal, January 1969

  • Epanechnikov, V. A.
  • Theory of Probability & Its Applications, Vol. 14, Issue 1
  • DOI: 10.1137/1114019

Impact of Perturbed Fission Source on the Effective Multiplication Factor in Monte Carlo Perturbation Calculations
journal, May 2005


Development of Calculation Technique for Iterated Fission Probability and Reactor Kinetic Parameters Using Continuous-Energy Monte Carlo Method
journal, November 2010


Surface Perturbation Theory
journal, June 1960


The OpenMC Monte Carlo particle transport code
journal, January 2013


Monte Carlo Method for Geometrical Perturbation and its Application to the Pulsed Fast Reactor
journal, August 1970


XGPT: Extending Monte Carlo Generalized Perturbation Theory capabilities to continuous-energy sensitivity functions
journal, October 2016


Kernel Density Estimation of Reaction Rates in Neutron Transport Simulations of Nuclear Reactors
journal, August 2017

  • Burke, Timothy P.; Kiedrowski, Brian C.; Martin, William R.
  • Nuclear Science and Engineering, Vol. 188, Issue 2
  • DOI: 10.1080/00295639.2017.1350000

Adjoint-Weighted Tallies for k -Eigenvalue Calculations with Continuous-Energy Monte Carlo
journal, July 2011

  • Kiedrowski, Brian C.; Brown, Forrest B.; Wilson, Paul P. H.
  • Nuclear Science and Engineering, Vol. 168, Issue 3
  • DOI: 10.13182/NSE10-22

Calculation of effective point kinetics parameters in the Serpent 2 Monte Carlo code
journal, March 2014


Kernel Density Estimation Method for Monte Carlo Global Flux Tallies
journal, March 2012

  • Banerjee, Kaushik; Martin, William R.
  • Nuclear Science and Engineering, Vol. 170, Issue 3
  • DOI: 10.13182/NSE10-77

Adjoint-Based k -Eigenvalue Sensitivity Coefficients to Nuclear Data Using Continuous-Energy Monte Carlo
journal, July 2013

  • Kiedrowski, Brian C.; Brown, Forrest B.
  • Nuclear Science and Engineering, Vol. 174, Issue 3
  • DOI: 10.13182/NSE12-46

ENDF/B-VII.1 Nuclear Data for Science and Technology: Cross Sections, Covariances, Fission Product Yields and Decay Data
journal, December 2011


Internal Interface Perturbations in Neutron Transport Theory
journal, January 1984


Application of the virtual density theory in fast reactor analysis based on the neutron transport calculation
journal, August 2017


Eigenvalue sensitivity to system dimensions
journal, April 2010