A multi-group extended linear discontinuous method for fixed-source discrete ordinates problems in slab geometry

Abstract

At present, neutron density calculation in non-multiplying media is relevant in many areas of engineering and science. In this paper, we propose the Extended Linear Discontinuous (ELD) method in multi-group discrete ordinates formulation, originally formulated for one-energy group fixed-source problems with isotropic scattering source in slab geometry. The proposed auxiliary equations are uncoupled on angular directions and combine the linear discontinuous approximation of the finite element method and the quasi-analytical general solution of the spectral nodal method. Thus, we can implement an efficient and simple algorithm using the conventional source iteration scheme for the sweeping equations. Numerical results for benchmark problems are presented to illustrate the accuracy and computational performance of the ELD method. The work shows that the main advantages of the proposed method are that the numerical scheme is stable for coarse-meshes, and its numerical results are more accurate than those generated by the Diamond Difference (DD) and Linear Discontinuous (LD) methods.