Revisiting the calculation of non-fundamental flux modes in a Candu reactor

The determination of non-fundamental modes of the diffusion equation is required for computing Candu reactor power distribution from analysis of in-core detector readings. The legacy method for computing these modes is the Hotelling deflation technique based on bi-harmonic decontamination. However, the Hotelling technique becomes unstable as the number of modes increases or as their eigenvalues become closer. Effective and fast alternatives are provided with Implicit Arnoldi Restarted Methods (IRAM). Among them, we investigated the Krylov-Schur method available in the SLEPc library, and we are proposing a custom implementation of the augmented block Householder Arnoldi (ABHA) method, similar to the Open Source implementation of Prof. James Baglama. In our work, the ABHA method is applied to the neutron diffusion equation, discretized with the mesh-centered finite difference method. (author)