skip to main content


This content will become publicly available on April 1, 2017

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

Here we evaluate 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 was developed to mitigate the impact of expensive on-line cross section processing steps. 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. Finally, both criticality (k-eigenvalue) and critical boron search problems are considered.
 [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:
OSTI Identifier:
Grant/Contract Number:
AC05-00OR22725; AC04-94AL85000
Accepted Manuscript
Journal Name:
Journal of Computational Physics
Additional Journal Information:
Journal Volume: 311; Journal Issue: C; Journal ID: ISSN 0021-9991
Research Org:
Oak Ridge National Laboratory (ORNL), Oak Ridge, TN (United States). Oak Ridge Leadership Computing Facility (OLCF); Consortium for Advanced Simulation of LWRs (CASL)
Sponsoring Org:
Country of Publication:
United States
22 GENERAL STUDIES OF NUCLEAR REACTORS Multiphysics; Jacobian-free Newton–Krylov; Anderson acceleration; Nuclear reactor analysis