Development of two-dimensional, multigroup neutron diffusion computer code based on GFEM with unstructured triangle elements

Abstract: Various methods for solving the forward/adjoint equation in hexagonal and rectangular geometries are known in the literatures. In this paper, the solution of multigroup forward/adjoint equation using Finite Element Method (FEM) for hexagonal and rectangular reactor cores is reported. The spatial discretization of equations is based on Galerkin FEM (GFEM) using unstructured triangle elements. Calculations are performed for both linear and quadratic approximations of the shape function; based on which results are compared. Using power iteration method for the forward and adjoint calculations, the forward and adjoint fluxes with the corresponding eigenvalues are obtained. The results are then benchmarked against the valid results for IAEA-2D, BIBLIS-2D and IAEA-PWR benchmark problems. Convergence rate of GFEM in linear and quadratic approximations of the shape function are calculated and results are quantitatively compared. A sensitivity analysis of the calculations to the number and arrangement of elements has been performed. [Copyright &y& Elsevier]