A ForthandBack Implicit ΛIteration
Abstract
In many radiative transfer (RT) problems, the sources contain a scattering term that couples all the specific RT equations, one for each frequency and direction, so that solving the problem means solving the system formed by these equations. Each of them is a firstorder linear differential equation with its own initial condition assigned at a different point of the medium, which makes the solution of the system extraordinarily difficult.
One simple way to achieve a solution is with the socalled Λiteration: sources and sinks given as a first approximation > computation of the specific intensities from their own RT equations > computation of the scattering terms > recomputation of the sources and sinks. This scheme is straightforward, but unfortunately in practice its convergence rate is too slow to be of value in the case of optically thick systems.
The aim of this paper is to show that a forthandback approach (the natural approach to describing sequentially the two intensities propagating along the two directions of a straight line), together with an implicit representation of the source function in the computation of the intensities within the above iterative scheme, can dramatically accelerate the convergence of the iterative process while retaining the straightforwardness of ordinary Λiteration.
 Publication:

The Astrophysical Journal
 Pub Date:
 October 1997
 DOI:
 10.1086/304626
 Bibcode:
 1997ApJ...487..735A
 Keywords:

 Methods: Numerical;
 Radiative Transfer