NOKIN1D: one-dimensional neutron kinetics based on a nodal collocation method

¹ Institute for the Industrial, Radiophysical and Environmental Safety Universitat Politècnica de València (UPV). Camí de Vera s/n, 46021 Valencia, Spain
² Department of Applied Mathematics Universitat Politècnica de València (UPV). Camí de Vera s/n, 46021 Valencia, Spain
³ SEA Propulsion S. L. Av. Atenas, 75. 28290 Las Rozas (Madrid), Spain.
⁴ Iberdrola Ingeniería y Construción S.A.U. Av. Manoteras, 20. 28050 Madrid, Spain.

^* Corresponding Author, gverdu@iqn.upv.es

Abstract

The TRAC-BF1 one-dimensional kinetic model is a formulation of the neutron diffusion equation in the two energy groups’ approximation, based on the analytical nodal method (ANM). The advantage compared with a zero-dimensional kinetic model is that the axial power profile may vary with time due to thermal-hydraulic parameter changes and/or actions of the control systems but at has the disadvantages that in unusual situations it fails to converge.

The nodal collocation method developed for the neutron diffusion equation and applied to the kinetics resolution of TRAC-BF1 thermal-hydraulics, is an adaptation of the traditional collocation methods for the discretization of partial differential equations, based on the development of the solution as a linear combination of analytical functions. It has chosen to use a nodal collocation method based on a development of Legendre polynomials of neutron fluxes in each cell.

The qualification is carried out by the analysis of the turbine trip transient from the NEA benchmark in Peach Bottom NPP using both the original 1D kinetics implemented in TRAC-BF1 and the 1D nodal collocation method.