open access

https://hdl.handle.net/2115/67483

Renormalization of the highly forward-peaked phase function using the double exponential formula for radiative transfer

英語

Double-exponential formula

2016-11

Springer

Journal of Mathematical Chemistry

54(10)

2048-2061

Numerical calculation of photon migration in biological tissue using the radiative transfer equation (RTE) has attracted great interests in biomedical optics and imaging. Because biological tissue is a highly forward-peaked scattering medium, renormalization of the phase function in numerical calculation of the RTE is crucial. This paper proposes a simple approach of renormalizing the phase function by the double exponential formula, which was heuristically modified from the original one. Firstly, the validity of the proposed approach was tested by comparing numerical results for an average cosine of the polar scattering angle calculated by the proposed approach with those by the conventional approach in highly forward-peaked scattering. The results show that calculation of the average cosine converged faster using the proposed approach than using the conventional one as a total number of discrete angular directions increases. Next, the accuracy of the numerical solutions of the RTE using the proposed approach was examined by comparing the numerical solutions with the analytical solutions of the RTE in a homogeneous medium with highly forward-peaked scattering. It was found that the proposed approach reduced the errors of the numerical solutions from those using the conventional one especially at a small value of the total number of the directions. This result suggests that the proposed approach has a possibility to improve the accuracy for the numerical results of the RTE in the highly scattering medium.

学術雑誌論文

Accepted Manuscript

The final publication is available at Springer via http://dx.doi.org/10.1007/s10910-016-0670-3

0259-9791