Nuclear Science and Engineering / Volume 155 / Number 2 / February 2007 / Pages 168-178
Technical Paper / Mathematics and Computation, Supercomputing, Reactor Physics and Nuclear and Biological Applications
A lumped linear-discontinuous spatial finite element discretization of the Sn equations in r-z geometry on triangular meshes is derived and computationally tested. An asymptotic analysis indicates that the scheme preserves the thick diffusion limit and behaves well with unresolved boundary layers. Computational results are presented that indicate the scheme is second-order accurate in the transport regime and that confirm the main predictions of the asymptotic diffusion-limit analysis.