On solutions to the P{sub n} equations for thermal radiative transfer
- 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)
We present results for the spherical harmonics (P{sub 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{sub n} equations will allow negative energy densities for any finite order of n. We also develop an implicit numerical method for solving the P{sub 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{sub 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{sub n} expansion orders on several problems. In two-dimensional problems the numerical P{sub 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.
- OSTI ID:
- 21028308
- Journal Information:
- Journal of Computational Physics, Vol. 227, Issue 5; Other Information: DOI: 10.1016/j.jcp.2007.11.027; PII: S0021-9991(07)00515-3; Copyright (c) 2007 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved; Country of input: International Atomic Energy Agency (IAEA); ISSN 0021-9991
- Country of Publication:
- United States
- Language:
- English
Similar Records
Implicit Riemann solvers for the Pn equations.
Upwind differencing scheme for the equations of ideal magnetohydrodynamics