American Nuclear Society

Home / Store / Journals / Electronic Articles / Nuclear Science and Engineering / Volume 147 / Number 1 / Pages 26-42

Krylov Subspace Iterations for Deterministic k-Eigenvalue Calculations

James S. Warsa, Todd A. Wareing, Jim E. Morel, John M. McGhee, Richard B. Lehoucq

Nuclear Science and Engineering / Volume 147 / Number 1 / Pages 26-42

May 2004


Member Price:$27.00
Member Savings:$3.00

The Implicitly Restarted Arnoldi Method (IRAM), a Krylov subspace iterative method, applied to k-eigenvalue calculations for criticality problems in deterministic transport codes is discussed. A computationally efficient alternative to the power iteration method that is typically used for such problems, the IRAM not only finds the largest eigenvalue but also several additional higher order eigenvectors with little extra computational cost. Implementation requires only modest changes to existing power iteration coding present in an SN transport program. Numerical results are presented for three-dimensional SN transport on unstructured tetrahedral meshes to compare the IRAM results with those computed using the traditional, unaccelerated power iteration method. The results indicate that the IRAM can be an efficient and powerful technique, especially for problems with dominance ratios approaching unity.

Questions or comments about the site? Contact the ANS Webmaster.