Low-Memory, Discrete Ordinates, Discontinuous Galerkin Methods for Radiative Transport
Abstract
The discrete ordinates discontinuous Galerkin (SN-DG) method is a well-established and practical approach for solving the radiative transport equation. In this paper, we study a low-memory variation of the upwind SN-DG method. The proposed method uses a smaller finite element space that is constructed by coupling spatial unknowns across collocation angles, thereby yielding an approximation with fewer degrees of freedom than the standard method. Like the original SN-DG method, the low-memory variation still preserves the asymptotic diffusion limit and maintains the characteristic structure needed for mesh sweeping algorithms. While we observe second-order convergence in the scattering dominated, diffusive regime, the low-memory method is in general only first-order accurate. To address this issue, we use upwind reconstruction to recover second-order accuracy. Finally, for both methods, numerical procedures based on upwind sweeps are proposed to reduce the system dimension in the underlying Krylov solver strategy.
- Authors:
-
- The Ohio State Univ., Columbus, OH (United States)
- Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States)
- Publication Date:
- Research Org.:
- Oak Ridge National Laboratory (ORNL), Oak Ridge, TN (United States); Oak Ridge Associated Univ., Oak Ridge, TN (United States)
- Sponsoring Org.:
- USDOE Office of Science (SC), Advanced Scientific Computing Research (ASCR); National Science Foundation (NSF)
- OSTI Identifier:
- 1902799
- Grant/Contract Number:
- AC05-00OR22725; SC0014664
- Resource Type:
- Accepted Manuscript
- Journal Name:
- SIAM Journal on Scientific Computing
- Additional Journal Information:
- Journal Volume: 42; Journal Issue: 4; Journal ID: ISSN 1064-8275
- Publisher:
- Society for Industrial and Applied Mathematics (SIAM)
- Country of Publication:
- United States
- Language:
- English
- Subject:
- 97 MATHEMATICS AND COMPUTING; radiative transport; discrete ordinates; discontinuous Galerkin; diffusion limit
Citation Formats
Sun, Zheng, and Hauck, Cory D. Low-Memory, Discrete Ordinates, Discontinuous Galerkin Methods for Radiative Transport. United States: N. p., 2020.
Web. doi:10.1137/19m1271956.
Sun, Zheng, & Hauck, Cory D. Low-Memory, Discrete Ordinates, Discontinuous Galerkin Methods for Radiative Transport. United States. https://doi.org/10.1137/19m1271956
Sun, Zheng, and Hauck, Cory D. Thu .
"Low-Memory, Discrete Ordinates, Discontinuous Galerkin Methods for Radiative Transport". United States. https://doi.org/10.1137/19m1271956. https://www.osti.gov/servlets/purl/1902799.
@article{osti_1902799,
title = {Low-Memory, Discrete Ordinates, Discontinuous Galerkin Methods for Radiative Transport},
author = {Sun, Zheng and Hauck, Cory D.},
abstractNote = {The discrete ordinates discontinuous Galerkin (SN-DG) method is a well-established and practical approach for solving the radiative transport equation. In this paper, we study a low-memory variation of the upwind SN-DG method. The proposed method uses a smaller finite element space that is constructed by coupling spatial unknowns across collocation angles, thereby yielding an approximation with fewer degrees of freedom than the standard method. Like the original SN-DG method, the low-memory variation still preserves the asymptotic diffusion limit and maintains the characteristic structure needed for mesh sweeping algorithms. While we observe second-order convergence in the scattering dominated, diffusive regime, the low-memory method is in general only first-order accurate. To address this issue, we use upwind reconstruction to recover second-order accuracy. Finally, for both methods, numerical procedures based on upwind sweeps are proposed to reduce the system dimension in the underlying Krylov solver strategy.},
doi = {10.1137/19m1271956},
journal = {SIAM Journal on Scientific Computing},
number = 4,
volume = 42,
place = {United States},
year = {Thu Jul 02 00:00:00 EDT 2020},
month = {Thu Jul 02 00:00:00 EDT 2020}
}
Works referenced in this record:
Positive and Asymptotic Preserving Approximation of the Radiation Transport Equation
journal, January 2020
- Guermond, Jean-Luc; Popov, Bojan; Ragusa, Jean
- SIAM Journal on Numerical Analysis, Vol. 58, Issue 1
Geometric Correction for Diffusive Expansion of Steady Neutron Transport Equation
journal, February 2015
- Wu, Lei; Guo, Yan
- Communications in Mathematical Physics, Vol. 336, Issue 3
On Locking and Robustness in the Finite Element Method
journal, October 1992
- Babuška, Ivo; Suri, Manil
- SIAM Journal on Numerical Analysis, Vol. 29, Issue 5
Runge-Kutta Discontinuous Galerkin Methods for Convection-Dominated Problems
journal, September 2001
- Cockburn, Bernardo; Shu, Chi-Wang
- Journal of Scientific Computing, Vol. 16, Issue 3, p. 173-261
Asymptotic Analysis of Upwind Discontinuous Galerkin Approximation of the Radiative Transport Equation in the Diffusive Limit
journal, January 2010
- Guermond, Jean-Luc; Kanschat, Guido
- SIAM Journal on Numerical Analysis, Vol. 48, Issue 1
Uniform asymptotic expansions in transport theory with small mean free paths, and the diffusion approximation
journal, April 1975
- Habetler, G. J.; Matkowsky, B. J.
- Journal of Mathematical Physics, Vol. 16, Issue 4
Unified Analysis of Discontinuous Galerkin Methods for Elliptic Problems
journal, January 2002
- Arnold, Douglas N.; Brezzi, Franco; Cockburn, Bernardo
- SIAM Journal on Numerical Analysis, Vol. 39, Issue 5
Boundary layers and homogenization of transport processes
journal, January 1979
- Bensoussan, Alain; Lions, Jacques-L.; Papanicolaou, George C.
- Publications of the Research Institute for Mathematical Sciences, Vol. 15, Issue 1
Numerical Schemes for Hyperbolic Conservation Laws with Stiff Relaxation Terms
journal, July 1996
- Jin, Shi; Levermore, C. David
- Journal of Computational Physics, Vol. 126, Issue 2
A robust S-DG-approximation for radiation transport in optically thick and diffusive regimes
journal, February 2012
- Ragusa, J. C.; Guermond, J. -L.; Kanschat, G.
- Journal of Computational Physics, Vol. 231, Issue 4
Asymptotic solutions of numerical transport problems in optically thick, diffusive regimes
journal, April 1987
- Larsen, Edward W.; Morel, J. E.; Miller, Warren F.
- Journal of Computational Physics, Vol. 69, Issue 2
Discontinuous Finite Element Transport Solutions in Thick Diffusive Problems
journal, March 2001
- Adams, Marvin L.
- Nuclear Science and Engineering, Vol. 137, Issue 3
Discrete-Ordinate Discontinuous Galerkin Methods for Solving the Radiative Transfer Equation
journal, January 2010
- Han, Weimin; Huang, Jianguo; Eichholz, Joseph A.
- SIAM Journal on Scientific Computing, Vol. 32, Issue 2
Asymptotic solutions of numerical transport problems in optically thick, diffusive regimes II
journal, July 1989
- larsen, Edward W.; Morel, J. E.
- Journal of Computational Physics, Vol. 83, Issue 1
Efficient Asymptotic-Preserving (AP) Schemes For Some Multiscale Kinetic Equations
journal, January 1999
- Jin, Shi
- SIAM Journal on Scientific Computing, Vol. 21, Issue 2
A Mixed Variational Framework for the Radiative Transfer Equation
journal, March 2012
- Egger, Herbert; Schlottbom, Matthias
- Mathematical Models and Methods in Applied Sciences, Vol. 22, Issue 03
Diffusion approximation and computation of the critical size
journal, February 1984
- Bardos, C.; Santos, R.; Sentis, R.
- Transactions of the American Mathematical Society, Vol. 284, Issue 2
On Locking and Robustness in the Finite Element Method
journal, October 1992
- Babuška, Ivo; Suri, Manil
- SIAM Journal on Numerical Analysis, Vol. 29, Issue 5
Boundary layers and homogenization of transport processes
journal, January 1979
- Bensoussan, Alain; Lions, Jacques-L.; Papanicolaou, George C.
- Publications of the Research Institute for Mathematical Sciences, Vol. 15, Issue 1