A Posteriori Error Estimation for a Nodal Method in Neutron Transport Calculations

An a posteriori error analysis of the spatial approximation is developed for the one-dimensional Arbitrarily High Order Transport-Nodal method. The error estimator preserves the order of convergence of the method when the mesh size tends to zero with respect to the L{sup 2} norm. It is based on the difference between two discrete solutions that are available from the analysis. The proposed estimator is decomposed into error indicators to allow the quantification of local errors. Some test problems with isotropic scattering are solved to compare the behavior of the true error to that of the estimated error

Improving the Accuracy of High-Order Nodal Transport Methods

This paper outlines some recent advances towards improving the accuracy of neutron transport calculations using the Arbitrarily High Order Transport-Nodal (AHOT-N) Method. These advances consist of several contributions: (a) A formula for the spatial weights that allows for the polynomial order to be raised arbitrarily high without suffering adverse effects from round-off error; (b) A reconstruction technique for the angular flux, based upon a recursive formula, that reduces the pointwise error by one ordeq (c) An a posterior error indicator that estimates the true error and its distribution throughout the domain, so that it can be used for adaptively refining the approximation. Present results are mainly for ID, extension to 2D-3D is in progress