skip to main content

DOE PAGESDOE PAGES

Title: An assessment of coupling algorithms for nuclear reactor core physics simulations

This paper evaluates the performance of multiphysics coupling algorithms applied to a light water nuclear reactor core simulation. The simulation couples the k-eigenvalue form of the neutron transport equation with heat conduction and subchannel flow equations. We compare Picard iteration (block Gauss–Seidel) to Anderson acceleration and multiple variants of preconditioned Jacobian-free Newton–Krylov (JFNK). The performance of the methods are evaluated over a range of energy group structures and core power levels. A novel physics-based approximation to a Jacobian-vector product has been developed to mitigate the impact of expensive on-line cross section processing steps. Furthermore, numerical simulations demonstrating the efficiency of JFNK and Anderson acceleration relative to standard Picard iteration are performed on a 3D model of a nuclear fuel assembly. Both criticality (k-eigenvalue) and critical boron search problems are considered.
Authors:
 [1] ;  [1] ;  [1] ;  [2] ;  [3] ;  [3] ;  [1] ;  [1]
  1. Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States)
  2. Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)
  3. North Carolina State Univ., Raleigh, NC (United States)
Publication Date:
Report Number(s):
SAND-2016-7579J
Journal ID: ISSN 0021-9991; NT0304000; NEAF343
Grant/Contract Number:
AC05-00OR22725; AC04-94AL85000
Type:
Accepted Manuscript
Journal Name:
Journal of Computational Physics
Additional Journal Information:
Journal Volume: 311; Journal Issue: C; Journal ID: ISSN 0021-9991
Publisher:
Elsevier
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); Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)
Sponsoring Org:
USDOE Office of Nuclear Energy (NE)
Country of Publication:
United States
Language:
English
Subject:
22 GENERAL STUDIES OF NUCLEAR REACTORS; Multiphysics; Jacobian-free Newton–Krylov; Anderson acceleration; Nuclear reactor analysis; 73 NUCLEAR PHYSICS AND RADIATION PHYSICS; 97 MATHEMATICS AND COMPUTING; multiphysics; nuclear reactor analysis
OSTI Identifier:
1247927
Alternate Identifier(s):
OSTI ID: 1291340; OSTI ID: 1347462

Hamilton, Steven, Berrill, Mark, Clarno, Kevin, Pawlowski, Roger, Toth, Alex, Kelley, C. T., Evans, Thomas, and Philip, Bobby. An assessment of coupling algorithms for nuclear reactor core physics simulations. United States: N. p., Web. doi:10.1016/j.jcp.2016.02.012.
Hamilton, Steven, Berrill, Mark, Clarno, Kevin, Pawlowski, Roger, Toth, Alex, Kelley, C. T., Evans, Thomas, & Philip, Bobby. An assessment of coupling algorithms for nuclear reactor core physics simulations. United States. doi:10.1016/j.jcp.2016.02.012.
Hamilton, Steven, Berrill, Mark, Clarno, Kevin, Pawlowski, Roger, Toth, Alex, Kelley, C. T., Evans, Thomas, and Philip, Bobby. 2016. "An assessment of coupling algorithms for nuclear reactor core physics simulations". United States. doi:10.1016/j.jcp.2016.02.012. https://www.osti.gov/servlets/purl/1247927.
@article{osti_1247927,
title = {An assessment of coupling algorithms for nuclear reactor core physics simulations},
author = {Hamilton, Steven and Berrill, Mark and Clarno, Kevin and Pawlowski, Roger and Toth, Alex and Kelley, C. T. and Evans, Thomas and Philip, Bobby},
abstractNote = {This paper evaluates the performance of multiphysics coupling algorithms applied to a light water nuclear reactor core simulation. The simulation couples the k-eigenvalue form of the neutron transport equation with heat conduction and subchannel flow equations. We compare Picard iteration (block Gauss–Seidel) to Anderson acceleration and multiple variants of preconditioned Jacobian-free Newton–Krylov (JFNK). The performance of the methods are evaluated over a range of energy group structures and core power levels. A novel physics-based approximation to a Jacobian-vector product has been developed to mitigate the impact of expensive on-line cross section processing steps. Furthermore, numerical simulations demonstrating the efficiency of JFNK and Anderson acceleration relative to standard Picard iteration are performed on a 3D model of a nuclear fuel assembly. Both criticality (k-eigenvalue) and critical boron search problems are considered.},
doi = {10.1016/j.jcp.2016.02.012},
journal = {Journal of Computational Physics},
number = C,
volume = 311,
place = {United States},
year = {2016},
month = {4}
}