Uncollided flux techniques for arbitrary finite element meshes
Abstract
The uncollided angular flux can be difficult to compute accurately in discrete-ordinate radiation transport codes, especially in weakly-scattering configurations with localized sources. It has long been recognized that an analytical or semi-analytical treatment of the uncollided flux, coupled with a discrete-ordinate solution for the collided flux, can yield dramatic improvements in solution accuracy and computational efficiency. In this paper, we present such an algorithm for the semi-analytical calculation of the uncollided flux. This algorithm is unique in several aspects: (1) it applies to arbitrary polyhedral cells (and can be thus coupled with collided flux solvers that support arbitrary polyhedral meshes without the need for explicit tetrahedral re-meshing), (2) it provides accurate uncollided solutions near sources, (3) it is devised with parallel implementation in mind, and (4) it minimizes the total number of traced rays and maintains a reasonable ray density on each local subdomain. This paper provides a complete derivation of the algorithm and demonstrates its important features on a set of simple examples and a standard transport benchmark. Assessment of its parallel performance will be the subject of a subsequent paper.
- Authors:
-
- Texas A & M Univ., College Station, TX (United States)
- Publication Date:
- Research Org.:
- Texas A & M Univ., College Station, TX (United States)
- Sponsoring Org.:
- USDOE National Nuclear Security Administration (NNSA)
- OSTI Identifier:
- 1801088
- Alternate Identifier(s):
- OSTI ID: 1561374
- Grant/Contract Number:
- NA0002376
- Resource Type:
- Accepted Manuscript
- Journal Name:
- Journal of Computational Physics
- Additional Journal Information:
- Journal Volume: 398; Journal ID: ISSN 0021-9991
- Publisher:
- Elsevier
- Country of Publication:
- United States
- Language:
- English
- Subject:
- 42 ENGINEERING; uncollided flux; ray tracing; arbitrary finite elements; first-collision source; radiation transport
Citation Formats
Hanuš, Milan, Harbour, Logan H., Ragusa, Jean C., Adams, Michael P., and Adams, Marvin L. Uncollided flux techniques for arbitrary finite element meshes. United States: N. p., 2019.
Web. doi:10.1016/j.jcp.2019.07.046.
Hanuš, Milan, Harbour, Logan H., Ragusa, Jean C., Adams, Michael P., & Adams, Marvin L. Uncollided flux techniques for arbitrary finite element meshes. United States. https://doi.org/10.1016/j.jcp.2019.07.046
Hanuš, Milan, Harbour, Logan H., Ragusa, Jean C., Adams, Michael P., and Adams, Marvin L. Tue .
"Uncollided flux techniques for arbitrary finite element meshes". United States. https://doi.org/10.1016/j.jcp.2019.07.046. https://www.osti.gov/servlets/purl/1801088.
@article{osti_1801088,
title = {Uncollided flux techniques for arbitrary finite element meshes},
author = {Hanuš, Milan and Harbour, Logan H. and Ragusa, Jean C. and Adams, Michael P. and Adams, Marvin L.},
abstractNote = {The uncollided angular flux can be difficult to compute accurately in discrete-ordinate radiation transport codes, especially in weakly-scattering configurations with localized sources. It has long been recognized that an analytical or semi-analytical treatment of the uncollided flux, coupled with a discrete-ordinate solution for the collided flux, can yield dramatic improvements in solution accuracy and computational efficiency. In this paper, we present such an algorithm for the semi-analytical calculation of the uncollided flux. This algorithm is unique in several aspects: (1) it applies to arbitrary polyhedral cells (and can be thus coupled with collided flux solvers that support arbitrary polyhedral meshes without the need for explicit tetrahedral re-meshing), (2) it provides accurate uncollided solutions near sources, (3) it is devised with parallel implementation in mind, and (4) it minimizes the total number of traced rays and maintains a reasonable ray density on each local subdomain. This paper provides a complete derivation of the algorithm and demonstrates its important features on a set of simple examples and a standard transport benchmark. Assessment of its parallel performance will be the subject of a subsequent paper.},
doi = {10.1016/j.jcp.2019.07.046},
journal = {Journal of Computational Physics},
number = ,
volume = 398,
place = {United States},
year = {Tue Aug 06 00:00:00 EDT 2019},
month = {Tue Aug 06 00:00:00 EDT 2019}
}
Web of Science
Works referenced in this record:
The Solid Angle of a Plane Triangle
journal, February 1983
- Van Oosterom, A.; Strackee, J.
- IEEE Transactions on Biomedical Engineering, Vol. BME-30, Issue 2
A numerical algorithm for the construction of efficient quadrature rules in two and higher dimensions
journal, January 2010
- Xiao, Hong; Gimbutas, Zydrunas
- Computers & Mathematics with Applications, Vol. 59, Issue 2
Fast Ray-Tetrahedron Intersection Using Plucker Coordinates
journal, January 2003
- Platis, Nikos; Theoharis, Theoharis
- Journal of Graphics Tools, Vol. 8, Issue 4
FNSUNCL3: First Collision Source Code for TORT
journal, March 2000
- Kosako, Kazuaki; Konno, Chikara
- Journal of Nuclear Science and Technology, Vol. 37, Issue sup1
Fast, Minimum Storage Ray-Triangle Intersection
journal, January 1997
- Möller, Tomas; Trumbore, Ben
- Journal of Graphics Tools, Vol. 2, Issue 1
Gmsh: A 3-D finite element mesh generator with built-in pre- and post-processing facilities
journal, September 2009
- Geuzaine, Christophe; Remacle, Jean-François
- International Journal for Numerical Methods in Engineering, Vol. 79, Issue 11
Discontinuous finite element solution of the radiation diffusion equation on arbitrary polygonal meshes and locally adapted quadrilateral grids
journal, January 2015
- Ragusa, Jean C.
- Journal of Computational Physics, Vol. 280
AETIUS solutions for Kobayashi 3D benchmarks with the first collision source method on the volume source and unstructured tetrahedral mesh
journal, March 2018
- Kim, Jong Woon; Lee, Young-Ouk
- Annals of Nuclear Energy, Vol. 113
Quadratic serendipity discontinuous finite element discretization for S transport on arbitrary polygonal grids
journal, December 2018
- Hackemack, Michael W.; Ragusa, Jean C.
- Journal of Computational Physics, Vol. 374
Denovo: A New Three-Dimensional Parallel Discrete Ordinates Code in SCALE
journal, August 2010
- Evans, Thomas M.; Stafford, Alissa S.; Slaybaugh, Rachel N.
- Nuclear Technology, Vol. 171, Issue 2
TORT solutions with FNSUNCL3 for Kobayashi's 3D benchmarks
journal, January 2001
- Konno, Chikara
- Progress in Nuclear Energy, Vol. 39, Issue 2
TetGen, a Delaunay-Based Quality Tetrahedral Mesh Generator
journal, February 2015
- Si, Hang
- ACM Transactions on Mathematical Software, Vol. 41, Issue 2
3D radiation transport benchmark problems and results for simple geometries with void region
journal, January 2001
- Kobayashi, Keisuke; Sugimura, Naoki; Nagaya, Yasunobu
- Progress in Nuclear Energy, Vol. 39, Issue 2
Ray-Effect Mitigation Methods for Two-Dimensional Neutron Transport Theory
journal, March 1977
- Miller, W. F.; Reed, Wm. H.
- Nuclear Science and Engineering, Vol. 62, Issue 3
On the singular structure of the uncollided and first-collided components of the Green's function
journal, September 2000
- Sanchez, R.
- Annals of Nuclear Energy, Vol. 27, Issue 13
A robust segment/triangle intersection algorithm for interference tests. Efficiency study
journal, July 2010
- Jiménez, Juan J.; Segura, Rafael J.; Feito, Francisco R.
- Computational Geometry, Vol. 43, Issue 5