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

Title: Calculating Alpha Eigenvalues and Eigenfunctions with a Markov Transition Rate Matrix Monte Carlo Method

Abstract

For a nuclear system in which the entire α-eigenvalue spectrum is known, eigenfunction expansion yields the time-dependent flux response to any arbitrary source. Applications in which this response is of interest include pulsed-neutron experiments, accelerator-driven subcritical systems, and fast burst reactors, where a steady-state assumption used in neutron transport is invalid for characterizing the time-dependent flux. To obtain the α-eigenvalue spectrum, the transition rate matrix method (TRMM) tallies transition rates describing neutron behavior in a discretized position-direction-energy phase space using Monte Carlo. Interpretation of the resulting Markov process transition rate matrix as the operator in the adjoint α-eigenvalue problem provides an avenue for determining a large finite set of α eigenvalues and eigenfunctions of a nuclear system. Results from the TRMM are verified using analytic solutions, time-dependent Monte Carlo simulations, and modal expansion from diffusion theory. For simplified infinite-medium and one-dimensional geometries, the TRMM accurately calculates eigenvalues, eigenfunctions, and eigenfunction expansion solutions. Applications and comparisons to measurements are made for the small fast burst reactor CALIBAN and the Fort St. Vrain high-temperature gas-cooled reactor. For large three-dimensional geometries, discretization of the large position-energy-direction phase space limits the accuracy of eigenfunction expansion solutions using the TRMM, but it can still generatemore » a fair estimate of the fundamental eigenvalue and eigenfunction. These results show that the TRMM generates an accurate estimate of a large number of α eigenvalues. This is not possible with existing Monte Carlo–based methods.« less

Authors:
ORCiD logo [1];  [2];  [2];  [3]
  1. Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States)
  2. Univ. of Michigan, Ann Arbor, MI (United States)
  3. Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
Publication Date:
Research Org.:
Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States)
Sponsoring Org.:
USDOE
OSTI Identifier:
1492187
DOE Contract Number:  
AC05-00OR22725
Resource Type:
Journal Article
Journal Name:
Nuclear Science and Engineering
Additional Journal Information:
Journal Volume: 192; Journal Issue: 2; Journal ID: ISSN 0029-5639
Publisher:
American Nuclear Society - Taylor & Francis
Country of Publication:
United States
Language:
English
Subject:
73 NUCLEAR PHYSICS AND RADIATION PHYSICS; Alpha eigenvalue; time-dependent transport; transition rate matrix; Markov process

Citation Formats

Betzler, Benjamin R., Kiedrowski, Brian C., Martin, William R., and Brown, Forrest B. Calculating Alpha Eigenvalues and Eigenfunctions with a Markov Transition Rate Matrix Monte Carlo Method. United States: N. p., 2018. Web. doi:10.1080/00295639.2018.1497397.
Betzler, Benjamin R., Kiedrowski, Brian C., Martin, William R., & Brown, Forrest B. Calculating Alpha Eigenvalues and Eigenfunctions with a Markov Transition Rate Matrix Monte Carlo Method. United States. doi:10.1080/00295639.2018.1497397.
Betzler, Benjamin R., Kiedrowski, Brian C., Martin, William R., and Brown, Forrest B. Fri . "Calculating Alpha Eigenvalues and Eigenfunctions with a Markov Transition Rate Matrix Monte Carlo Method". United States. doi:10.1080/00295639.2018.1497397.
@article{osti_1492187,
title = {Calculating Alpha Eigenvalues and Eigenfunctions with a Markov Transition Rate Matrix Monte Carlo Method},
author = {Betzler, Benjamin R. and Kiedrowski, Brian C. and Martin, William R. and Brown, Forrest B.},
abstractNote = {For a nuclear system in which the entire α-eigenvalue spectrum is known, eigenfunction expansion yields the time-dependent flux response to any arbitrary source. Applications in which this response is of interest include pulsed-neutron experiments, accelerator-driven subcritical systems, and fast burst reactors, where a steady-state assumption used in neutron transport is invalid for characterizing the time-dependent flux. To obtain the α-eigenvalue spectrum, the transition rate matrix method (TRMM) tallies transition rates describing neutron behavior in a discretized position-direction-energy phase space using Monte Carlo. Interpretation of the resulting Markov process transition rate matrix as the operator in the adjoint α-eigenvalue problem provides an avenue for determining a large finite set of α eigenvalues and eigenfunctions of a nuclear system. Results from the TRMM are verified using analytic solutions, time-dependent Monte Carlo simulations, and modal expansion from diffusion theory. For simplified infinite-medium and one-dimensional geometries, the TRMM accurately calculates eigenvalues, eigenfunctions, and eigenfunction expansion solutions. Applications and comparisons to measurements are made for the small fast burst reactor CALIBAN and the Fort St. Vrain high-temperature gas-cooled reactor. For large three-dimensional geometries, discretization of the large position-energy-direction phase space limits the accuracy of eigenfunction expansion solutions using the TRMM, but it can still generate a fair estimate of the fundamental eigenvalue and eigenfunction. These results show that the TRMM generates an accurate estimate of a large number of α eigenvalues. This is not possible with existing Monte Carlo–based methods.},
doi = {10.1080/00295639.2018.1497397},
journal = {Nuclear Science and Engineering},
issn = {0029-5639},
number = 2,
volume = 192,
place = {United States},
year = {2018},
month = {9}
}

Works referenced in this record:

Numerical Solution of the Time-Dependent Multigroup Diffusion Equations
journal, February 1968

  • Andrews, J. Barclay; Hansen, K. F.
  • Nuclear Science and Engineering, Vol. 31, Issue 2
  • DOI: 10.13182/NSE68-A18242

Time–eigenvalue calculations in multi-region Cartesian geometry using Green’s functions
journal, June 2005


Iterative computation of time-eigenvalues of the neutron transport equation
journal, November 2003


Diffusion theory methods for spatial kinetics calculations
journal, January 1996


On the spectrum of the linear transport operator
journal, November 1974

  • Larsen, Edward W.; Zweifel, Paul F.
  • Journal of Mathematical Physics, Vol. 15, Issue 11
  • DOI: 10.1063/1.1666570

Three-dimensional transport calculation of multiple alpha modes in subcritical systems
journal, December 2012


Time-Dependent, One-Speed Transport via Generalized Functions
journal, March 2001


An asymptotic expansion in the theory of neutron transport
journal, May 1958

  • Jörgens, Konrad
  • Communications on Pure and Applied Mathematics, Vol. 11, Issue 2
  • DOI: 10.1002/cpa.3160110206

the slab geometry
journal, March 1956


ENDF/B-VII.0: Next Generation Evaluated Nuclear Data Library for Nuclear Science and Technology
journal, December 2006


Calculating α Eigenvalues of One-Dimensional Media with Monte Carlo
journal, July 2014

  • Betzler, B. R.; Martin, W. R.; Kiedrowski, B. C.
  • Journal of Computational and Theoretical Transport, Vol. 43, Issue 1-7
  • DOI: 10.1080/00411450.2014.909851

Asymptotic solutions of the transport equation for thermal neutrons
journal, April 1963


Two Rossi- α Techniques for Measuring the Effective Delayed Neutron Fraction
journal, February 1993

  • Spriggs, Gregory D.
  • Nuclear Science and Engineering, Vol. 113, Issue 2
  • DOI: 10.13182/NSE93-2

A simple scheme for the direct evaluation of time-eigenvalues of neutron transport equation
journal, January 2003


Pulsed neutron source measurements in the subcritical ADS experiment YALINA-Booster
journal, December 2008


Iterative method for obtaining the prompt and delayed alpha-modes of the diffusion equation
journal, September 2011


Higher Order αMode Eigenvalue Calculation by Monte Carlo Power Iteration
journal, January 2011

  • Yamamoto, Toshihiro
  • Progress in Nuclear Science and Technology, Vol. 2, Issue 0
  • DOI: 10.15669/pnst.2.826

Monte Carlo simulation and benchmarking of pulsed neutron experiments in variable buckling Beo systems
journal, May 2009


On the spectrum of an unsymmetric operator arising in the transport theory of neutrons
journal, May 1955

  • Lehner, Joseph; Wing, G. Milton
  • Communications on Pure and Applied Mathematics, Vol. 8, Issue 2
  • DOI: 10.1002/cpa.3160080202

Calculating infinite-medium α-eigenvalue spectra with Monte Carlo using a transition rate matrix method
journal, December 2015


On the Relation between Decay Constants and Critical Parameters in Monoenergetic Neutron Transport
journal, March 1983

  • Dahl, E. B.; Protopopescu, V.; Sjöstrand, N. G.
  • Nuclear Science and Engineering, Vol. 83, Issue 3
  • DOI: 10.13182/NSE83-A17570