Axially Quadratic Discontinuous Finite Elements for Neutron Transport Reactor Simulations

摘要

This work serves to indicate the increase in computational efficiency of the axially quadratic LDQZ spatial discretization method in PDT with respect to the fully linear LD method. The LDQZ method utilizes the set of LD basis functions and appends a fifth axially quadratic basis function, which reduces both the difficulty of implementation and the computational resources needed to build the linear system matrices for each spatial cell. An increase in efficiency by a factor of about 100 is seen with the use of LDQZ over LD for the 3D C5G7 UO_2 pin with the figure of merit dependent on reactivity differences. However, the LDQZ method may result in higher efficiencies if used with geometries with higher axial heterogeneity.