A response matrix method for the refined Analytic Function Expansion Nodal (AFEN) method in the two-dimensional hexagonal geometry and its numerical performance

In order to improve calculational efficiency of the CAPP code in the analysis of the hexagonal reactor core, we have tried to implement a refined AFEN method with transverse gradient basis functions and interface flux moments in the hexagonal geometry. The numerical scheme for the refined AFEN method adopted here is the response matrix method that uses the interface partial currents as nodal unknowns instead of the interface fluxes used in the original AFEN method. Since the response matrix method is single-node based, it has good properties such as good calculational efficiency and parallel computing affinity. Because a refined AFEN method equivalent nonlinear FDM response matrix method tried first could not provide a numerically stable solution, a direct formulation of the refined AFEN response matrix were developed. To show the numerical performance of this response matrix method against the original AFEN method, the numerical error analyses were performed for several benchmark problems including the VVER-440 LWR benchmark problem and the MHTGR-350 HTGR benchmark problem. The results showed a more than three times speedup in computing time for the LWR and HTGR benchmark problems due to good convergence and excellent calculational efficiency of the refined AFEN response matrix method.