EPJ Web Conf.
Volume 247, 2021PHYSOR2020 – International Conference on Physics of Reactors: Transition to a Scalable Nuclear Future
|Number of page(s)||8|
|Published online||22 February 2021|
DEVELOPMENT AND TESTING OF THE MODIFIED SIMPSON’S RULE FOR DISCRETE ORDINATES TRANSPORT APPLICATIONS
Department of Nuclear Engineering, North Carolina State University 2500 Stinson Drive, Campus Box 7909, Raleigh, NC 27695-7909, USA
Published online: 22 February 2021
A new angular quadrature type termed Modified Simpson Trapezoidal (MST) is developed based on the conventional Simpson’s 1/3 rule where the angular pattern over polar levels has a trapezoid shape. An adaptive coefficient correction scheme is developed to enable our new quadrature to integrate the angular flux over subintervals separated by the interior jump irregularities. A two-dimensional test problem is employed to verify the angular discretization error in the uncollided SN scalar flux computed with our new quadrature sets, as well as conventional angular quadrature types. Numerical results show that the MST quadrature error in the point-wise scalar flux converges with second order against increasing number of discrete angles, while the error obtained with other conventional quadrature types converges slower than first order depending on the regularity of the exact point-wise uncollided angular flux. In order to reduce the number of discrete points needed, a variant of the MST quadrature, namely MSTP30, is developed by using the Quadruple Range  polar quadrature with fixed 30 polar angles and applying the MST quadrature to the azimuthal dependence in each polar level. The angular discretization error in the point-wise SN scalar flux obtained with MSTP30 sets converges with fourth order because the polar discretization error is sufficiently reduced that MSTP30 behaves like a one-dimensional quadrature. Furthermore, because MSTP30 computes the integral over subintervals that keep the true solution’s irregularity at the boundaries, this fourth order convergence rate is unaffected by such inevitable irregularities.
Key words: Modified Simpson’s Rule / Discrete Ordinates / Quadrature Error / Uncollided Point-wise Flux
© The Authors, published by EDP Sciences, 2021
This is an Open Access article distributed under the terms of the Creative Commons Attribution License 4.0, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.