Solution of the fixed source neutron diffusion equation by using the pseudo-harmonics method

Lima, Zelmo Rodrigues de; Silva, F.C. da; Alvim, A.C.M.

doi:10.1016/j.anucene.2003.08.008

Cited by 5 publications

(1 citation statement)

References 7 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…One of the first methods used for solving this eigenvalue problem was the finite difference method (FDM), however as it requires very mesh points to model accurately, it becomes a bit slow from the computational point of view [1]. Over time other methodologies have been implemented for solving the neutron diffusion eigenvalue problem, as for example: Finite Element Method [2], [Pseudo-Harmonics Expansion Method [1], Finite Volume Method [3], Lagrange Polynomial Algorithms (LAP) using the same approach in FDM [4], Taylor Series Expansion Method [5], Isogeometric Analysis [6], Integral Transform Techniques [7] and Generalized Integral Transform Technique (GITT) combined with Laplace Transform given by [8] and Nodal Method [9].…”

Section: Introductionmentioning

confidence: 99%

Solution of the Multigroup Neutron Diffusion Eigenvalue Problem in Slab Geometry by Modified Power Method

Zanette¹,

Petersen²,

Tavares³

2021

Braz. J. Rad. Sci.

View full text Add to dashboard Cite

We describe in this work the application of the modified power method for solve the multigroup neutron diffusion eigenvalue problem in slab geometry considering two-dimensions for nuclear reactor global calculations. It is well known that criticality calculations can often be best approached by solving eigenvalue problems. The criticality in nuclear reactors physics plays a relevant role since establishes the ratio between the numbers of neutrons generated in successive fission reactions. In order to solve the eigenvalue problem, a modified power method is used to obtain the dominant eigenvalue (effective multiplication factor) and its corresponding eigenfunction (scalar neutron flux), which is non-negative in every domain, that is, physically relevant. The innovation of this work is solving the neutron diffusion equation in analytical form for each new iteration of the power method. For solve this problem we propose to apply the Finite Fourier Sine Transform on one of the variables obtaining a transformed problem which is resolved by well-established methods for ordinary differential equations. The inverse Fourier Transform is used to reconstruct the solution for the original problem. It is known that the power method is an iterative source method in which is updated by the neutron flux expression of previous iteration. Thus, for each new iteration, the neutron flux expression becomes larger and more complex due to analytical solution what makes us propose that it be reconstructed through a polynomial interpolation. The methodology is implemented to solve a homogeneous problem and the results are compared with works presents in the literature.

show abstract