Evaluation of dominant time-eigenvalues of neutron transport equation by Meyer’s sub-space iterations

Abstract

The aim of this paper is to explore the use of Meyer’s sub-space iteration (SSI) method for the evaluation of dominant prompt time-eigenvalues of the neutron transport equation. The integro-differential form of the transport equation is considered. The SSI method is known to be an efficient technique to find the dominant eigenvalues of a non-symmetric matrix. It has been earlier used for eigenvalue problems in neutron diffusion theory. However, it does not seem to be tried in the transport theory case. Here, the use of SSI has been tested in transport theory for some 1-D mono-energetic homogeneous and heterogeneous benchmark problems. The space variable is discretised by finite differencing while neutron directions are discretised by discrete ordinates (S n -) method. The SSI method needs frequent multiplication of the relevant matrix operator with vectors. As known from earlier works in this area, this can be achieved in terms of external source calculations for which a 1-D programme was developed and used. With the availability of more versatile S n -method codes, it may perhaps be possible to extend use of SSI to more realistic cases.