Open AccessNumerical modeling of the radiative transfer in a turbid medium using the synthetic iteration
Author(s) -
Vladimir P. Budak,
Gennady A. Kaloshin,
Oleg V. Shagalov,
Victor S. Zheltov
Publication year - 2015
Publication title -
optics express
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.394
H-Index - 271
ISSN - 1094-4087
DOI - 10.1364/oe.23.00a829
Subject(s) - discretization , radiative transfer , convergence (economics) , mathematics , numerical analysis , metric (unit) , integral equation , iterative method , algorithm , optics , computer science , mathematical analysis , physics , operations management , economics , economic growth
In this paper we propose the fast, but the accurate algorithm for numerical modeling of light fields in the turbid media slab. For the numerical solution of the radiative transfer equation (RTE) it is required its discretization based on the elimination of the solution anisotropic part and the replacement of the scattering integral by a finite sum. The solution regular part is determined numerically. A good choice of the method of the solution anisotropic part elimination determines the high convergence of the algorithm in the mean square metric. The method of synthetic iterations can be used to improve the convergence in the uniform metric. A significant increase in the solution accuracy with the use of synthetic iterations allows applying the two-stream approximation for the regular part determination. This approach permits to generalize the proposed method in the case of an arbitrary 3D geometry of the medium.