DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: 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:
ORCiD logo [1]; ORCiD logo [2]
  1. The Ohio State Univ., Columbus, OH (United States)
  2. 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
  • DOI: 10.1137/19M1260785

Geometric Correction for Diffusive Expansion of Steady Neutron Transport Equation
journal, February 2015


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
  • DOI: 10.1137/0729075

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
  • DOI: 10.1023/A:1012873910884

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
  • DOI: 10.1137/090746938

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
  • DOI: 10.1063/1.522618

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
  • DOI: 10.1137/S0036142901384162

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
  • DOI: 10.2977/prims/1195188427

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
  • DOI: 10.1006/jcph.1996.0149

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
  • DOI: 10.1016/j.jcp.2011.11.017

Asymptotic solutions of numerical transport problems in optically thick, diffusive regimes
journal, April 1987


Discontinuous Finite Element Transport Solutions in Thick Diffusive Problems
journal, March 2001

  • Adams, Marvin L.
  • Nuclear Science and Engineering, Vol. 137, Issue 3
  • DOI: 10.13182/NSE00-41

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
  • DOI: 10.1137/090767340

Asymptotic solutions of numerical transport problems in optically thick, diffusive regimes II
journal, July 1989


Efficient Asymptotic-Preserving (AP) Schemes For Some Multiscale Kinetic Equations
journal, January 1999


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
  • DOI: 10.1142/S021820251150014X

Diffusion approximation and computation of the critical size
journal, February 1984


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
  • DOI: 10.1137/0729075

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
  • DOI: 10.2977/prims/1195188427