Asymptotic convergence of the angular discretization error in the scalar flux computed from the particle transport equation with the method of discrete ordinates

Hu, Xiaoyu; Azmy, Yousry Y.

doi:10.1016/j.anucene.2019.107199

Asymptotic convergence of the angular discretization error in the scalar flux computed from the particle transport equation with the method of discrete ordinates

Journal Article · Mon Nov 18 23:00:00 EST 2019 · Annals of Nuclear Energy (Oxford)

DOI:https://doi.org/10.1016/j.anucene.2019.107199· OSTI ID:1801122

Hu, Xiaoyu ^[1]; Azmy, Yousry Y. ^[2]

North Carolina State Univ., Raleigh, NC (United States); OSTI
North Carolina State Univ., Raleigh, NC (United States)

The asymptotic convergence of the angular discretization error in the scalar flux solution of the particle transport equation computed with the Discrete Ordinates (S_N) method with increasing quadrature order is examined. Five selected angular quadrature types are considered: Level Symmetric (LS), Legendre-Chebyshev Quadrangular (LCQ), Legendre-Chebyshev Triangular(LCT), Quadruple Range (QR) and Quadruple Range Spence-type (QRS) quadrature sets. We relate the S_N angular discretization error to the quadrature error, and split the total flux into the uncollided flux and the fully collided flux, and then we verify the uncollided scalar flux error and the fully collided scalar flux error separately. After developing the theoretical basis for the relationship between true solution regularity and quadrature rule error, we employ a two-dimensional problem to verify the observed and theoretical convergence orders for the region-averaged uncollided and fully collided scalar flux errors. Numerical results show that the angular discretization errors in the region-averaged scalar flux obtained by different quadrature types converge with different rates, that are commensurate with the regularity order of the exact angular flux within the quadrature-designed integration interval. The angular discretization error in the uncollided region-averaged scalar flux obtained by LC class quadratures converges linearly, and the error obtained by QR class quadrature sets converges quadratically. The angular discretization error in the fully collided region-averaged scalar flux converges linearly for LC class quadratures, and faster-than-second order for QR class quadrature sets.

Research Organization:: North Carolina State Univ., Raleigh, NC (United States)

Sponsoring Organization:: USDOE National Nuclear Security Administration (NNSA)

Grant/Contract Number:: NA0002576

OSTI ID:: 1801122

Alternate ID(s):: OSTI ID: 2325393

Journal Information:: Annals of Nuclear Energy (Oxford), Journal Name: Annals of Nuclear Energy (Oxford) Vol. 138; ISSN 0306-4549

Publisher:: ElsevierCopyright Statement

Country of Publication:: United States

Language:: English

References (13)

Arbitrarily high order characteristic methods for solving the neutron transport equation Azmy, Yousry Y. Annals of Nuclear Energy, Vol. 19, Issue 10-12 https://doi.org/10.1016/0306-4549(92)90004-U	journal	October 1992
The generation of arbitrary order, non-classical, Gauss-type quadrature for transport applications Spence, Peter J. Journal of Computational Physics, Vol. 296 https://doi.org/10.1016/j.jcp.2015.04.038	journal	September 2015
Recent Studies on the Asymptotic Convergence of the Spatial Discretization for Two-Dimensional Discrete Ordinates Solutions Barichello, L. B.; Tres, A.; Picoloto, C. B. Journal of Computational and Theoretical Transport, Vol. 45, Issue 4 https://doi.org/10.1080/23324309.2016.1171242	journal	May 2016
Angular Quadratures for Improved Transport Computations Abu-Shumays, I. K. Transport Theory and Statistical Physics, Vol. 30, Issue 2-3 https://doi.org/10.1081/TT-100105367	journal	June 2001
Pointwise Convergence of the Three-Dimensional Discrete Ordinate Method Madsen, N. K. SIAM Journal on Numerical Analysis, Vol. 8, Issue 2 https://doi.org/10.1137/0708027	journal	June 1971
Asymptotic Error Bounds for a Class of Interpolatory Quadratures Kütz, Martin SIAM Journal on Numerical Analysis, Vol. 21, Issue 1 https://doi.org/10.1137/0721011	journal	February 1984
Numerical Evaluation of Multiple Integrals Haber, Seymour SIAM Review, Vol. 12, Issue 4 https://doi.org/10.1137/1012102	journal	October 1970
Analysis of Ray-Effect Mitigation Techniques Morel, J. E.; Wareing, T. A.; Lowrie, R. B. Nuclear Science and Engineering, Vol. 144, Issue 1 https://doi.org/10.13182/NSE01-48	journal	May 2003
Error Comparison of Diamond Difference, Nodal, and Characteristic Methods for Solving Multidimensional Transport Problems with the Discrete Ordinates Approximation Duo, J. I.; Azmy, Y. Y. Nuclear Science and Engineering, Vol. 156, Issue 2 https://doi.org/10.13182/NSE05-91	journal	June 2007
A Quantitative Theory of Angular Truncation Errors in Three-Dimensional S _N Calculations Larsen, Edward W.; Wollaber, Allan B. Nuclear Science and Engineering, Vol. 160, Issue 3 https://doi.org/10.13182/NSE160-267	journal	November 2008
Ray-Effect Mitigation Methods for Two-Dimensional Neutron Transport Theory Miller, W. F.; Reed, Wm. H. Nuclear Science and Engineering, Vol. 62, Issue 3 https://doi.org/10.13182/NSE62-391	journal	March 1977
Compatible Product Angular Quadrature for Neutron Transport in x-y Geometry Abu-Shumays, I. K. Nuclear Science and Engineering, Vol. 64, Issue 2 https://doi.org/10.13182/NSE64-299	journal	October 1977
Spatial Convergence Properties of the Diamond Difference Method in x,y Geometry Larsen, Edward W. Nuclear Science and Engineering, Vol. 80, Issue 4 https://doi.org/10.13182/NSE82-A18980	journal	April 1982

Similar Records

Discrete Ordinate Quadrature Selection for Reactor-based Eigenvalue Problems

Conference · Mon Dec 31 23:00:00 EST 2012 · OSTI ID:1083720

Discrete ordinate quadrature selection for reactor-based Eigenvalue problems

Conference · Mon Jul 01 00:00:00 EDT 2013 · OSTI ID:22212741

Recent Studies on the Asymptotic Convergence of the Spatial Discretization for Two-Dimensional Discrete Ordinates Solutions

Journal Article · Wed Jul 13 00:00:00 EDT 2016 · Journal of Computational and Theoretical Transport · OSTI ID:1437435

Related Subjects

42 ENGINEERING
angular discretization error
discrete ordinates method
numerical quadrature
transport equation

Asymptotic convergence of the angular discretization error in the scalar flux computed from the particle transport equation with the method of discrete ordinates

Citation Formats

References (13)

Similar Records

Related Subjects