skip to main content

DOE PAGESDOE PAGES

This content will become publicly available on December 11, 2018

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

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:
Grant/Contract Number:
NA0002576; NRC-HQ-13-G-38- 0007
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
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)
Country of Publication:
United States
Language:
English
Subject:
97 MATHEMATICS AND COMPUTING; 70 PLASMA PHYSICS AND FUSION TECHNOLOGY; Monte Carlo; system dimensions; sensitivities
OSTI Identifier:
1438024

Burke, Timothy P., and Kiedrowski, Brian C.. Monte Carlo Perturbation Theory Estimates of Sensitivities to System Dimensions. United States: N. p., 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.. 2017. "Monte Carlo Perturbation Theory Estimates of Sensitivities to System Dimensions". United States. doi:10.1080/00295639.2017.1388093.
@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}
}