RESUMO
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.
RESUMO
We present an analysis of a number of different approximations for the plane transmittance Tp and diffuse attenuation coefficient Kd of a semi-infinite, unbounded, plane-parallel, and optically homogeneous layer. The maximally wide optical conditions (from the full absorption to the full scattering and from the fully forward to the fully backward scattering) were considered. The approximations were analyzed from the point of view of their physical limitations and closeness to the numerical solution of the radiative transfer equation for the plane transmittance. The main criterion for inclusion of the models for analysis was the possibility of practical use, i.e., approximations were well parameterized and included only easily measured or estimated parameters. A detailed analysis of errors for different Tp and Kd models showed that the two-stream radiative transfer Ben-David model yields the best results over all optical conditions and depths. However, the quasi-single-scattering and polynomial Gordon's approximations proved to be the best for the depths close to zero.
