Radiative transfer with delta-Eddington-type phase functions

Appl Math Comput. 2017 May 1:300:70-78. doi: 10.1016/j.amc.2016.12.001. Epub 2016 Dec 26.

Abstract

The radiative transfer equation (RTE) arises in a wide variety of applications, in particular, in biomedical imaging applications associated with the propagation of light through the biological tissue. However, highly forward-peaked scattering feature in a biological medium makes it very challenging to numerically solve the RTE problem accurately. One idea to overcome the difficulty associated with the highly forward-peaked scattering is through the use of a delta-Eddington phase function. This paper is devoted to an RTE framework with a family of delta-Eddington-type phase functions. Significance in biomedical imaging applications of the RTE with delta-Eddington-type phase functions are explained. Mathematical studies of the problems include solution existence, uniqueness, and continuous dependence on the problem data: the inflow boundary value, the source function, the absorption coefficient, and the scattering coefficient. Numerical results are presented to show that employing a delta-Eddington-type phase function with properly chosen parameters provides accurate simulation results for light propagation within highly forward-peaked scattering media.

Keywords: continuous dependence; existence; generalized delta-Eddington phase function; radiative transfer equation; uniqueness.