3D calculation of the lambda eigenvalues and eigenmodes of the two-group neutron diffusion equation by coarse-mesh nodal methods

Abstract This paper shows the development and roots of the lambda eigenvector and eigenmodes calculation by coarse-mesh finite difference nodal methods. In addition, this paper shows an inter-comparison of the eigenvalues and power profiles obtained by different 3D nodal methods with two neutron energy groups. The methods compared are: the nodal collocation method (Verdu et al., 1998, 1993; Hebert, 1987) with different orders of the Legendre expansions, the modified coarse-mesh nodal method explained in this paper, and the method implemented in the PARCS code by Wysocki et al. (2014, 2015). In this paper we have developed a program NODAL-LAMBDA that uses a two-group modified coarse-mesh finite difference method with albedo boundary conditions. Some of the approximations performed originally by Borressen (1971) have been discarded and more exact expressions have been used. We compare for instance the eigenvalues and power profiles obtained with Borressen original approach of 1.5 group and with 2 groups. Also, some improvements in the albedo boundary conditions suggested by (Turney 1975; Chung et al., 1981 ) have been incorporated to the code as an option. The goal is to obtain the eigenvalues and the sub-criticalities ( k i − k 0 ) of the harmonic modes that can be excited during an instability event in a fast way and with an acceptable precision.

[1]  Gumersindo Verdú,et al.  Resolution of the Generalized Eigenvalue Problem in the Neutron Diffusion Equation Discretized by the Finite Volume Method , 2014 .

[2]  Toshikazu Takeda,et al.  Neutronic model for modal multichannel analysis of out-of-phase instability in boiling water reactor cores , 1997 .

[3]  Gene H. Golub,et al.  Matrix computations , 1983 .

[4]  Danny C. Sorensen,et al.  Implicit Application of Polynomial Filters in a k-Step Arnoldi Method , 1992, SIAM J. Matrix Anal. Appl..

[5]  The Laplace Transform Method for the Albedo Boundary Conditions in Neutron Diffusion Eigenvalue Problems , 2010 .

[6]  Gumersindo Verdú,et al.  Coupling of density wave oscillations in parallel channels with high order modal kinetics: application to BWR out of phase oscillations , 2000 .

[7]  Y. Saad,et al.  Numerical Methods for Large Eigenvalue Problems , 2011 .

[8]  A. Manera,et al.  A Physical mechanism for rotating lines of symmetry in BWR OIJT-of-phase limit cycle oscillations , 2015 .

[9]  Peter Arbenz,et al.  Lecture Notes on Solving Large Scale Eigenvalue Problems , 2012 .

[10]  Gumersindo Verdú,et al.  Solution of the Lambda modes problem of a nuclear power reactor using an h–p finite element method , 2014 .

[11]  Paritosh Pandey,et al.  NEMSQR: A 3-D multi group diffusion theory code based on nodal expansion method for square geometry , 2014 .

[12]  Petr N. Vabishchevich,et al.  Spectral properties of dynamic processes in a nuclear reactor , 2017 .

[13]  Christian Constanda,et al.  Integral Methods in Science and Engineering, Volume 2 , 2010 .

[14]  Alain Hébert,et al.  Development of the nodal collocation method for solving the neutron diffusion equation , 1987 .

[15]  A. Manera,et al.  TRACE/PARCS analysis of out-of-phase power oscillations with a rotating line of symmetry , 2014 .

[16]  W. Arnoldi.  The principle of minimized iterations in the solution of the matrix eigenvalue problem , 1951 .

[17]  W. B. Terney.  Albedo-Adjusted Fast-Neutron Diffusion Coefficients in Reactor Reflectors , 1975 .

[18]  Gumersindo Verdú,et al.  3D λ-modes of the neutron-diffusion equation , 1994 .

[19]  Gumersindo Verdú,et al.  A nodal modal method for the neutron diffusion equation. Application to BWR instabilities analysis , 2002 .

[20]  M. Döring,et al.  Subspace iteration for nonsymmetric eigenvalue problems applied to the λ-eigenvalue problem , 1993 .

[21]  R. Herrero.  Métodos modales para el estudio de inestabilidades en reactores nucleares bwr , 2002 .

[22]  S. Børresen,et al.  A Simplified, Coarse-Mesh, Three-Dimensional Diffusion Scheme for Calculating the Gross Power Distribution in a Boiling Water Reactor , 1971 .

[23]  Nam Zin Cho,et al.  FUNDAMENTALS AND RECENT DEVELOPMENTS OF REACTOR PHYSICS METHODS , 2005 .

[24]  A. F. Henry,et al.  The replacement of reflectors by Albedo-type boundary conditions. , 1975 .

[25]  A. Hébert,et al.  A USER GUIDE FOR TRIVAC VERSION4 , 2014 .

[26]  Gumersindo Verdú,et al.  The implicit restarted Arnoldi method, an efficient alternative to solve the neutron diffusion equation , 1999 .