American Nuclear Society

Home / Publications / Journals / Nuclear Science and Engineering / Volume 189 / Number 3

An Accurate Globally Conservative Subdomain Discontinuous Least-Squares Scheme for Solving Neutron Transport Problems

Weixiong Zheng, Ryan G. McClarren, Jim E. Morel

Nuclear Science and Engineering / Volume 189 / Number 3 / March 2018 / Pages 259-271

Technical Paper /

Received:December 6, 2016
Accepted:November 16, 2017
Published:February 8, 2018

In this work, we present a subdomain discontinuous least-squares (SDLS) scheme for neutronics problems. Least-squares (LS) methods are known to be inaccurate for problems with sharp total cross-section interfaces. In addition, the LS scheme is known not to be globally conservative in heterogeneous problems. In problems where global conservation is important, e.g., k-eigenvalue problems, a conservative treatment must be applied. In this study, we propose an SDLS method that retains global conservation and, as a result, gives high accuracy on eigenvalue problems. Such a method resembles the LS formulation in each subdomain without a material interface and differs from LS in that an additional LS interface term appears for each interface. The scalar flux is continuous in each subdomain with the continuous finite element method while discontinuous on interfaces for every pair of contiguous subdomains. The SDLS numerical results are compared with those obtained from other numerical methods with test problems having material interfaces. High accuracy of scalar flux in fixed-source problems and in eigenvalue problems is demonstrated.

Questions or comments about the site? Contact the ANS Webmaster.