Development of a finite volume inter-cell polynomial expansion method for the neutron diffusion equation

Abstract

Heterogeneous nuclear reactors require numerical methods to solve the neutron diffusion equation (NDE) to obtain the neutron flux distribution inside them, by discretizing the heterogeneous geometry in a set of homogeneous regions. This discretization requires additional equations at the inner faces of two adjacent cells: neutron flux and current continuity, which imply an excess of equations. The finite volume method (FVM) is suitable to be applied to NDE, because it can be easily applied to any mesh and it is typically used in the transport equations due to the conservation of the transported quantity within the volume. However, the gradient and face-averaged values in the FVM are typically calculated as a function of the cell-averaged values of adjacent cells. So, if the materials of the adjacent cells are different, the neutron current condition could not be accomplished. Therefore, a polynomial expansion of the neutron flux is developed in each cell for assuring the accomplishment of the flux and current continuity and calculating both analytically. In this polynomial expansion, the polynomial terms for each cell were assigned previously and the constant coefficients are determined by solving the eigenvalue problem with SLEPc. A sensitivity analysis for determining the best set of polynomial terms is performed.