skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: On the eigenstructure of spherical harmonic equations for radiative transport

ORCiD logo;
Publication Date:
Sponsoring Org.:
OSTI Identifier:
Grant/Contract Number:
Resource Type:
Journal Article: Publisher's Accepted Manuscript
Journal Name:
Computers and Mathematics with Applications (Oxford)
Additional Journal Information:
Journal Name: Computers and Mathematics with Applications (Oxford); Journal Volume: 72; Journal Issue: 2; Related Information: CHORUS Timestamp: 2017-10-06 15:59:44; Journal ID: ISSN 0898-1221
Country of Publication:
United Kingdom

Citation Formats

Garrett, C. Kristopher, and Hauck, Cory D. On the eigenstructure of spherical harmonic equations for radiative transport. United Kingdom: N. p., 2016. Web. doi:10.1016/j.camwa.2015.05.030.
Garrett, C. Kristopher, & Hauck, Cory D. On the eigenstructure of spherical harmonic equations for radiative transport. United Kingdom. doi:10.1016/j.camwa.2015.05.030.
Garrett, C. Kristopher, and Hauck, Cory D. 2016. "On the eigenstructure of spherical harmonic equations for radiative transport". United Kingdom. doi:10.1016/j.camwa.2015.05.030.
title = {On the eigenstructure of spherical harmonic equations for radiative transport},
author = {Garrett, C. Kristopher and Hauck, Cory D.},
abstractNote = {},
doi = {10.1016/j.camwa.2015.05.030},
journal = {Computers and Mathematics with Applications (Oxford)},
number = 2,
volume = 72,
place = {United Kingdom},
year = 2016,
month = 7

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record at 10.1016/j.camwa.2015.05.030

Citation Metrics:
Cited by: 1work
Citation information provided by
Web of Science

Save / Share:
  • Highlights: {yields} We developed a variable order global basis scheme to solve light transport in 3D. {yields} Based on finite elements, the method can be applied to a wide class of geometries. {yields} It is computationally cheap when compared to the fixed order scheme. {yields} Comparisons with local basis method and other models demonstrate its accuracy. {yields} Addresses problems encountered n modeling of light transport in human brain. - Abstract: We propose the P{sub N} approximation based on a finite element framework for solving the radiative transport equation with optical tomography as the primary application area. The key idea ismore » to employ a variable order spherical harmonic expansion for angular discretization based on the proximity to the source and the local scattering coefficient. The proposed scheme is shown to be computationally efficient compared to employing homogeneously high orders of expansion everywhere in the domain. In addition the numerical method is shown to accurately describe the void regions encountered in the forward modeling of real-life specimens such as infant brains. The accuracy of the method is demonstrated over three model problems where the P{sub N} approximation is compared against Monte Carlo simulations and other state-of-the-art methods.« less
  • A new two-dimensional monochromatic method that computes the transfer of solar or thermal radiation through atmospheres with arbitrary optical properties is described. The model discretizes the radiative transfer equation by expanding the angular part of the radiance field in a spherical harmonic series and representing the spatial part with a discrete grid. The resulting sparse coupled system of equations is solved iteratively with the conjugate gradient method. A Monte Carlo model is used for extensive verification of outgoing flux and radiance values from both smooth and highly variable (multifractal) media. The spherical harmonic expansion naturally allows for different levels ofmore » approximation, but tests show that the 2D equivalent of the two-stream approximation is poor at approximating variations in the outgoing flux. The model developed here is shown to be highly efficient so that media with tens of thousands of grid points can be computed in minutes. The large improvement in efficiency will permit quick, accurate radiative transfer calculations of realistic cloud fields and improve our understanding of the effect of inhomogeneity on radiative transfer in cloudy atmospheres. 21 refs., 7 figs., 3 tabs.« less
  • We develop a time-dependent, multi-group, multi-dimensional relativistic radiative transfer code, which is required to numerically investigate radiation from relativistic fluids that are involved in, e.g., gamma-ray bursts and active galactic nuclei. The code is based on the spherical harmonic discrete ordinate method (SHDOM) which evaluates a source function including anisotropic scattering in spherical harmonics and implicitly solves the static radiative transfer equation with ray tracing in discrete ordinates. We implement treatments of time dependence, multi-frequency bins, Lorentz transformation, and elastic Thomson and inelastic Compton scattering to the publicly available SHDOM code. Our code adopts a mixed-frame approach; the source functionmore » is evaluated in the comoving frame, whereas the radiative transfer equation is solved in the laboratory frame. This implementation is validated using various test problems and comparisons with the results from a relativistic Monte Carlo code. These validations confirm that the code correctly calculates the intensity and its evolution in the computational domain. The code enables us to obtain an Eddington tensor that relates the first and third moments of intensity (energy density and radiation pressure) and is frequently used as a closure relation in radiation hydrodynamics calculations.« less