Affiliation:
1. Department of Mechanical Engineering, University of California, Berkeley, CA 94720
Abstract
This paper introduces a powerful but simple methodology for solving the general equation of radiative transfer for scattering and/or absorbing one-dimensional systems. Existing methods, usually designed to handle specific boundary and energy equilibrium conditions, either provide crude estimates or involve intricate mathematical analysis coupled with numerical techniques. In contrast, the present scheme, which uses a discrete-ordinate technique to reduce the integro-differential equation to a system of ordinary differential equations, utilizes readily available software routines to solve the resulting set of coupled first-order ordinary differential equations as a two-point boundary value problem. The advantage of this approach is that the user is freed from having to understand complicated mathematical analysis and perform extensive computer programming. Additionally, the software used is state of the art, which is less prone to numerical instabilities and inaccuracies. Any degree of scattering anisotropy and albedo can be incorporated along with different conditions of energy equilibrium or specified temperature distributions and boundary conditions. Examples are presented where the radiative transfer is computed by using different quadratures such as Gaussian, Lobatto, Fiveland, Chebyshev, and Newton-Cotes. Comparison with benchmark cases shows that in a highly forward scattering medium Gaussian quadrature provides the most accurate and stable solutions.
Subject
Mechanical Engineering,Mechanics of Materials,Condensed Matter Physics,General Materials Science
Cited by
69 articles.
订阅此论文施引文献
订阅此论文施引文献,注册后可以免费订阅5篇论文的施引文献,订阅后可以查看论文全部施引文献