DOI QR코드

DOI QR Code

Discontinuous finite-element quadrature sets based on icosahedron for the discrete ordinates method

  • Dai, Ni (School of Nuclear Science and Engineering, North China Electric Power University) ;
  • Zhang, Bin (School of Nuclear Science and Engineering, North China Electric Power University) ;
  • Chen, Yixue (School of Nuclear Science and Engineering, North China Electric Power University)
  • Received : 2019.08.06
  • Accepted : 2019.11.20
  • Published : 2020.06.25

Abstract

The discrete ordinates method (SN) is one of the major shielding calculation method, which is suitable for solving deep-penetration transport problems. Our objective is to explore the available quadrature sets and to improve the accuracy in shielding problems involving strong anisotropy. The linear discontinuous finite-element (LDFE) quadrature sets based on the icosahedron (in short, ICLDFE quadrature sets) are developed by defining projected points on the surfaces of the icosahedron. Weights are then introduced in the integration of the discontinuous finite-element basis functions in the relevant angular regions. The multivariate secant method is used to optimize the discrete directions and their corresponding weights. The numerical integration of polynomials in the direction cosines and the Kobayashi benchmark are used to analyze and verify the properties of these new quadrature sets. Results show that the ICLDFE quadrature sets can exactly integrate the zero-order and first-order of the spherical harmonic functions over one-twentieth of the spherical surface. As for the Kobayashi benchmark problem, the maximum relative error between the fifth-order ICLDFE quadrature sets and references is only -0.55%. The ICLDFE quadrature sets provide better integration precision of the spherical harmonic functions in local discrete angle domains and higher accuracy for simple shielding problems.

Keywords

References

  1. D.G. Cacuci, Handbook of Nuclear Engineering, Springer, New York, 2010.
  2. E.B. Podgorsak, Radiation physics for medical physicists, Med. Phys. 33 (1) (2005), 249-249.
  3. D. Mihalas, B.W. Mihalas, Foundations of Radiation Hydrodynamics, Oxford University Press, New York, 1984.
  4. E.E. Lewis, W.F. Miller, Computational Methods of Neutron Transport, American Nuclear Society, Inc. La Grange Park, Illinois, 1993.
  5. B.G. Carlson, Lee, et al., Mechanical Quadrature and the Transport Equation, Report LAMS-2573, Los Alamos Scientific Laboratory, 1961.
  6. W.F. Walters, Use of the Chebyshev-Legendre Quadrature Set in Discrete Ordinates Codes, Report LA-UR-87-3621, Los Alamos Scientific Laboratory, 1987.
  7. C.P. Thurgood, A. Pollard, H.A. Becker, The $T_N$ quadrature set for the discrete ordinates method, J. Heat Transf. 117 (4) (1995) 1068-1070. https://doi.org/10.1115/1.2836285
  8. G. Longoni, A. Haghighat, Development and Application of the Regional Angular Refinement Technique and its Application to Non-conventional Problems, 2002, PHYSOR, Seoul, 2002.
  9. J.J. Jarell, M.L. Adams, Discrete Ordinates Quadrature Sets Based on Linear Discontinuous Finite Elements, Mathematics and Computational Methods Applied to Nuclear Science and Engineering (M&C 2011), Rio de Janeiro, 2011.
  10. C.D. Ahrens, Highly efficient, exact quadratures for three-dimensional discrete ordinates transport calculations, Nucl. Sci. Eng. 170 (1) (2012) 98-101. https://doi.org/10.13182/NSE10-69TN
  11. K.E. Atkinson, W. Han, I. Ebrary, Spherical Harmonics and Approximations on the Unit Sphere : an Introduction, Springer Verlag, 2012.
  12. C.Y. Lau, M.L. Adams, Discrete Ordinates Quadratures Based on Linear and Quadratic Discontinuous Finite Elements over Spherical Quadrilaterals, LaGrange Park: 2015 ANS Joint International Conference on Mathematics and Computation, M&C 2015), 2015.
  13. J. Jarrell, An Adaptive Angular Discretization Method for Neutral-Particle Transport in Three-Dimensional Geometries, Dissertations & Theses-Gradworks, 2012.
  14. K. Kobayashi, N. Sugimura, Y. Nagaya, 3D radiation transport benchmark problems and results for simple geometries with void region, Prog. Nucl. Energy 39 (2) (2001) 119-144. https://doi.org/10.1016/S0149-1970(01)00007-5
  15. E.W. Larsen, Spatial convergence properties of the diamond difference method in x, y geometry, Nucl. Sci. Eng. 80 (4) (1982) 710-713. https://doi.org/10.13182/NSE82-A18980
  16. Y.X. Chen, B. Zhang, L. Zhang, et al., ARES: a parallel discrete ordinates transport code for radiation shielding applications and reactor physics analysis, Sci. Technol.Nucl. Install. 2017 (2017) 1-11.