Nuclear Science and Engineering / Volume 147 / Number 1 / May 2004 / Pages 26-42
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.