Calculation of the eigenfunctions and corresponding eigenvalues of the two-group diffusion equation in heterogeneous systems
Paper in proceedings, 2007

This paper deals with the extension of an existing code developed at the Department of Nuclear Engineering, Chalmers University of Technology, to the calculation of the higher order eigenfunctions and corresponding eigenvalues of the two-group diffusion equations. More specifically, this code, originally meant for calculating the two-group static neutron flux and its adjoint for any two-dimensional heterogeneous system, is modified in order to calculate any higher mode of both the neutron flux and its adjoint. A one-dimensional two-region system is first considered for benchmark purposes, where the reference solution can be determined semi-analytically. The numerical solution is computed using a modified power iteration method. The agreement between the numerical solution and the reference solution is found to be excellent, both regarding the eigenfunctions and the eigenvalues, and both for the forward and adjoint problems. The modified power iteration method is then applied to the calculation of the different modes of the neutron flux and its adjoint in a two-dimensional heterogeneous representation of the Swedish Forsmark-1 BWR core corresponding to the 1996/1997 channel instability event.

diffusion theory

numerical methods

eigenvalues

eigenfunctions

Author

Christophe Demaziere

Chalmers, Applied Physics, Nuclear Engineering

Carl Sunde

Chalmers, Applied Physics, Nuclear Engineering

Proc. Joint Int. Topl. Mtg. Mathematics & Computation and Supercomputing in Nuclear Applications (M & C + SNA 2007)

Subject Categories

Other Engineering and Technologies not elsewhere specified

Areas of Advance

Energy

ISBN

0-89448-059-6

More information

Created

10/8/2017