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

Title: Stability and accuracy of 3D neutron transport simulations using the 2D/1D method in MPACT

Abstract

We derived a consistent “2D/1D” neutron transport method from the 3D Boltzmann transport equation, to calculate fuel-pin-resolved neutron fluxes for realistic full-core Pressurized Water Reactor (PWR) problems. The 2D/1D method employs the Method of Characteristics to discretize the radial variables and a lower order transport solution to discretize the axial variable. Our paper describes the theory of the 2D/1D method and its implementation in the MPACT code, which has become the whole-core deterministic neutron transport solver for the Consortium for Advanced Simulations of Light Water Reactors (CASL) core simulator VERA-CS. We also performed several applications on both leadership-class and industry-class computing clusters. Results are presented for whole-core solutions of the Watts Bar Nuclear Power Station Unit 1 and compared to both continuous-energy Monte Carlo results and plant data.

Authors:
ORCiD logo [1];  [1]; ORCiD logo [2];  [2]; ORCiD logo [2]; ORCiD logo [2];  [2];  [2];  [1]
  1. Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States)
  2. Univ. of Michigan, Ann Arbor, MI (United States). Dept. of Nuclear Engineering and Radiological Sciences
Publication Date:
Research Org.:
Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States). Oak Ridge Leadership Computing Facility (OLCF); Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States). Consortium for Advanced Simulation of LWRs (CASL)
Sponsoring Org.:
USDOE Office of Science (SC)
OSTI Identifier:
1326469
Alternate Identifier(s):
OSTI ID: 1359305
Grant/Contract Number:  
AC05-00OR22725
Resource Type:
Accepted Manuscript
Journal Name:
Journal of Computational Physics
Additional Journal Information:
Journal Volume: 326; Journal Issue: C; Journal ID: ISSN 0021-9991
Publisher:
Elsevier
Country of Publication:
United States
Language:
English
Subject:
97 MATHEMATICS AND COMPUTING; 71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; 2D/1D method; MOC; MPACT

Citation Formats

Collins, Benjamin, Stimpson, Shane, Kelley, Blake W., Young, Mitchell T. H., Kochunas, Brendan, Graham, Aaron, Larsen, Edward W., Downar, Thomas, and Godfrey, Andrew. Stability and accuracy of 3D neutron transport simulations using the 2D/1D method in MPACT. United States: N. p., 2016. Web. doi:10.1016/j.jcp.2016.08.022.
Collins, Benjamin, Stimpson, Shane, Kelley, Blake W., Young, Mitchell T. H., Kochunas, Brendan, Graham, Aaron, Larsen, Edward W., Downar, Thomas, & Godfrey, Andrew. Stability and accuracy of 3D neutron transport simulations using the 2D/1D method in MPACT. United States. https://doi.org/10.1016/j.jcp.2016.08.022
Collins, Benjamin, Stimpson, Shane, Kelley, Blake W., Young, Mitchell T. H., Kochunas, Brendan, Graham, Aaron, Larsen, Edward W., Downar, Thomas, and Godfrey, Andrew. Thu . "Stability and accuracy of 3D neutron transport simulations using the 2D/1D method in MPACT". United States. https://doi.org/10.1016/j.jcp.2016.08.022. https://www.osti.gov/servlets/purl/1326469.
@article{osti_1326469,
title = {Stability and accuracy of 3D neutron transport simulations using the 2D/1D method in MPACT},
author = {Collins, Benjamin and Stimpson, Shane and Kelley, Blake W. and Young, Mitchell T. H. and Kochunas, Brendan and Graham, Aaron and Larsen, Edward W. and Downar, Thomas and Godfrey, Andrew},
abstractNote = {We derived a consistent “2D/1D” neutron transport method from the 3D Boltzmann transport equation, to calculate fuel-pin-resolved neutron fluxes for realistic full-core Pressurized Water Reactor (PWR) problems. The 2D/1D method employs the Method of Characteristics to discretize the radial variables and a lower order transport solution to discretize the axial variable. Our paper describes the theory of the 2D/1D method and its implementation in the MPACT code, which has become the whole-core deterministic neutron transport solver for the Consortium for Advanced Simulations of Light Water Reactors (CASL) core simulator VERA-CS. We also performed several applications on both leadership-class and industry-class computing clusters. Results are presented for whole-core solutions of the Watts Bar Nuclear Power Station Unit 1 and compared to both continuous-energy Monte Carlo results and plant data.},
doi = {10.1016/j.jcp.2016.08.022},
journal = {Journal of Computational Physics},
number = C,
volume = 326,
place = {United States},
year = {Thu Aug 25 00:00:00 EDT 2016},
month = {Thu Aug 25 00:00:00 EDT 2016}
}

Journal Article:

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

Save / Share:

Works referenced in this record:

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

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

Apollo2 year 2010
journal, October 2010

  • Sanchez, Richard; Zmijarevi, Igor; Coste-Delclaux, M.
  • Nuclear Engineering and Technology, Vol. 42, Issue 5
  • DOI: 10.5516/NET.2010.42.5.474

Assembly homogenization techniques for light water reactor analysis
journal, January 1986


A new characteristics algorithm for 3D transport calculations
journal, January 2003


Solution of the C5G7MOX benchmark three-dimensional extension problems by the DeCART direct whole core calculation code
journal, July 2006


Practical numerical reactor employing direct whole core neutron transport and subchannel thermal/hydraulic solvers
journal, December 2013


High-Fidelity Light Water Reactor Analysis with the Numerical Nuclear Reactor
journal, March 2007

  • Weber, David P.; Sofu, Tanju; Yang, Won Sik
  • Nuclear Science and Engineering, Vol. 155, Issue 3
  • DOI: 10.13182/NSE07-A2672

A consistent 2D/1D approximation to the 3D neutron transport equation
journal, December 2015


Derivation of Optimum Polar Angle Quadrature Set for the Method of Characteristics Based on Approximation Error for the Bickley Function
journal, February 2007

  • Yamamoto, Akio; Tabuchi, Masato; Sugimura, Naoki
  • Journal of Nuclear Science and Technology, Vol. 44, Issue 2
  • DOI: 10.1080/18811248.2007.9711266

Diffusion theory methods for spatial kinetics calculations
journal, January 1996


Monte Carlo Criticality Methods and Analysis Capabilities in SCALE
journal, May 2011

  • Goluoglu, Sedat; Petrie, Lester M.; Dunn, Michael E.
  • Nuclear Technology, Vol. 174, Issue 2
  • DOI: 10.13182/NT10-124

Works referencing / citing this record:

Consistent pCMFD Acceleration Schemes of the Three-Dimensional Transport Code PROTEUS-MOC
journal, February 2019


The Subray Method of Characteristics
journal, January 2019

  • Graham, Aaron M.; Collins, Benjamin S.; Downar, Thomas J.
  • Nuclear Science and Engineering, Vol. 193, Issue 6
  • DOI: 10.1080/00295639.2018.1550988

A 2D/1D Algorithm for Effective Cross-Section Generation in Fast Reactor Neutronic Transport Calculations
journal, July 2018


The RAPID Fission Matrix Approach to Reactor Core Criticality Calculations
journal, July 2018

  • Walters, William J.; Roskoff, Nathan J.; Haghighat, Alireza
  • Nuclear Science and Engineering, Vol. 192, Issue 1
  • DOI: 10.1080/00295639.2018.1497395

Improved Accuracy in the 2-D/1-D Method with Anisotropic Transverse Leakage and Cross-Section Homogenization
journal, September 2018


Polar Parity for Efficient Evaluation of Anisotropic Transverse Leakage in the 2D/1D Transport Method
journal, July 2019


Quasi-Diffusion Method with 3-D Cross Sections for TREAT Core Analysis
journal, December 2019