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

Title: RAYLEIGH QUOTIENT ITERATION IN 3D, DETERMINISTIC NEUTRON TRANSPORT

Abstract

Today’s “grand challenge” neutron transport problems require 3-D meshes with billions of cells, hundreds of energy groups, and accurate quadratures and scattering expansions. Leadership-class computers provide platforms on which high-fidelity fluxes can be calculated. However, appropriate methods are needed that can use these machines effectively. Such methods must be able to to use hundreds of thousands of cores and have good convergence properties. Rayleigh quotient iteration (RQI) is an eigenvalue solver that has been added to the SN code Denovo to address convergence. Rayleigh quotient iteration is an optimal shifted inverse iteration method that should converge in fewer iterations than the more common power method and other shifted inverse iteration methods for many problems of interest. Denovo’s RQI uses a new multigroup Krylov solver for the fixed source solutions inside every iteration that allows parallelization in energy in addition to space and angle. This Krylov solver has been shown to scale successfully to 200,000 cores: for example one test problem scaled from 69,120 cores to 190,080 cores with 98% efficiency. This paper shows that RQI works for some small problems. However, the Krylov method upon which it relies does not always converge because RQI creates ill-conditioned systems. This result leadsmore » to the conclusion that preconditioning is needed to allow this method to be applicable to a wider variety of problems.« less

Authors:
 [1];  [2];  [2];  [3]
  1. Univ. of Wisconsin, Madison, WI (United States). Dept. of Nuclear Engineering and Engineering Physics
  2. Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States). Radiation Transport Group
  3. Univ. of Wisconsin, Madison, WI (United States)
Publication Date:
Research Org.:
Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States). Oak Ridge Leadership Computing Facility (OLCF)
Sponsoring Org.:
USDOE Office of Science (SC)
OSTI Identifier:
1567684
Resource Type:
Conference
Resource Relation:
Conference: PHYSOR 2012 Advances in Reactor Physics Linking Research, Industry, and Education, Knoxville, Tennessee, USA, April 15-20, 2012
Country of Publication:
United States
Language:
English

Citation Formats

Slaybaugh, R. N., Evans, T. M., Davidson, G. G., and Wilson, Paul P.H.. RAYLEIGH QUOTIENT ITERATION IN 3D, DETERMINISTIC NEUTRON TRANSPORT. United States: N. p., 2012. Web. doi:10.6084/m9.figshare.1504109.
Slaybaugh, R. N., Evans, T. M., Davidson, G. G., & Wilson, Paul P.H.. RAYLEIGH QUOTIENT ITERATION IN 3D, DETERMINISTIC NEUTRON TRANSPORT. United States. https://doi.org/10.6084/m9.figshare.1504109
Slaybaugh, R. N., Evans, T. M., Davidson, G. G., and Wilson, Paul P.H.. 2012. "RAYLEIGH QUOTIENT ITERATION IN 3D, DETERMINISTIC NEUTRON TRANSPORT". United States. https://doi.org/10.6084/m9.figshare.1504109.
@article{osti_1567684,
title = {RAYLEIGH QUOTIENT ITERATION IN 3D, DETERMINISTIC NEUTRON TRANSPORT},
author = {Slaybaugh, R. N. and Evans, T. M. and Davidson, G. G. and Wilson, Paul P.H.},
abstractNote = {Today’s “grand challenge” neutron transport problems require 3-D meshes with billions of cells, hundreds of energy groups, and accurate quadratures and scattering expansions. Leadership-class computers provide platforms on which high-fidelity fluxes can be calculated. However, appropriate methods are needed that can use these machines effectively. Such methods must be able to to use hundreds of thousands of cores and have good convergence properties. Rayleigh quotient iteration (RQI) is an eigenvalue solver that has been added to the SN code Denovo to address convergence. Rayleigh quotient iteration is an optimal shifted inverse iteration method that should converge in fewer iterations than the more common power method and other shifted inverse iteration methods for many problems of interest. Denovo’s RQI uses a new multigroup Krylov solver for the fixed source solutions inside every iteration that allows parallelization in energy in addition to space and angle. This Krylov solver has been shown to scale successfully to 200,000 cores: for example one test problem scaled from 69,120 cores to 190,080 cores with 98% efficiency. This paper shows that RQI works for some small problems. However, the Krylov method upon which it relies does not always converge because RQI creates ill-conditioned systems. This result leads to the conclusion that preconditioning is needed to allow this method to be applicable to a wider variety of problems.},
doi = {10.6084/m9.figshare.1504109},
url = {https://www.osti.gov/biblio/1567684}, journal = {},
number = ,
volume = ,
place = {United States},
year = {2012},
month = {1}
}

Conference:
Other availability
Please see Document Availability for additional information on obtaining the full-text document. Library patrons may search WorldCat to identify libraries that hold this conference proceeding.

Save / Share: