On solutions to the Pn equations for thermal radiative transfer

  • Authors:

  • Ryan G. McClarren;James Paul Holloway;Thomas A. Brunner

  • Affiliations:

  • Computational Physics and Methods Group, Los Alamos National Laboratory, P.O. Box 1663, MS D413 Los Alamos, NM 87545, United States;Department of Nuclear Engineering and Radiological Sciences, College of Engineering, University of Michigan, 2355 Bonisteel Boulevard, Ann Arbor, MI 48109 2104, United States;Sandia National Laboratories, P.O. Box 5800, MS 1186, Albuquerque, NM 87185 1186, United States

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2008

Quantified Score

Hi-index 31.48

Visualization

Abstract

We present results for the spherical harmonics (P"n) method for solving problems of time-dependent thermal radiative transport. We prove a theorem that demonstrates that in the streaming limit, the spatially and temporally continuous P"n equations will allow negative energy densities for any finite order of n. We also develop an implicit numerical method for solving the P"n equations to explore the impact of the theorem. The numerical method uses a high-resolution Riemann solver to produce an upwinded discretization. We employ a quasi-linear approach to integrate the nonlinearites added to make the scheme non-oscillatory. We use the backward Euler method for time integration and treat the material interaction terms fully nonlinearly. Reflecting boundary conditions for the P"n equations are presented and we show how to implement this boundary condition using ghost cells. The implicit method was able to produce robust results to thermal transport problems in one and two dimensions. The numerical method is used to analyze the accuracy of various P"n expansion orders on several problems. In two-dimensional problems the numerical P"n solutions contained negative radiation energy densities as predicted by our theorem. The numerical results showed that the material temperature also became negative, a result outside the scope of the theorem. Our numerical method can handle these negative values, but they would cause problems in a radiation-hydrodynamics calculation.