Angular quadrature perturbations in radiative transfer theory

Abstract

A new numerical method is presented for solving the general equation of radiative transfer. The approximation, which replaces the integral term over angle in the transfer equation by a quadrature sum, is studied; an estimate of the error involved is obtained and this error, which may be thought of as a further source or sink of photons (depending upon the sign), can then be used to evaluate a corection to the radiation field originally determined. This process may then be continued as a perturbation series. The method is found to give a final solution, when starting from the Eddington approximation, at least as accurate as that obtained using variable Eddington factors. Furthermore, the technique involves very little extra computing over that required using the Eddington approximation, and may be trivially generalized to any radiative transfer problem. It can also be used in conjunction with any of the existing methods for solving the equation of transfer. Examples are given in the context of spectral line formation in slab geometry.