Radiative transfer equation modeling by streamline diffusion modified continuous Galerkin method

Abstract

Optical tomography has a wide range of biomedical applications. Accurate prediction of photon transport in media is critical, as it directly affects the accuracy of the reconstructions. The radiative transfer equation (RTE) is the most accurate deterministic forward model, yet it has not been widely employed in practice due to the challenges in robust and efficient numerical implementations in high dimensions. Herein, we propose a method that combines the discrete ordinate method (DOM) with a streamline diffusion modified continuous Galerkin method to numerically solve RTE. Additionally, a phase function normalization technique was employed to dramatically reduce the instability of the DOM with fewer discrete angular points. To illustrate the accuracy and robustness of our method, the computed solutions to RTE were compared with Monte Carlo (MC) simulations when two types of sources (ideal pencil beam and Gaussian beam) and multiple optical properties were tested. Results show that with standard optical properties of human tissue, photon densities obtained using RTE are, on average, around 5% of those predicted by MC simulations in the entire/deeper region. These results suggest that this implementation of the finite element method-RTE is an accurate forward model for optical tomography in human tissues.

Accurate and efficient computation of the 3D radiative transfer equation in highly forward-peaked scattering media using a renormalization approach

“…The LSE and EO sets have limitations on N Ω under the condition that w l is positive, i.e., maximum values of N Ω are 360 and 288 for the LSE and for the EO sets, respectively. On the other hand, it has been reported in (Gregersen and York, 2005;Sanchez, 2012;Long et al, 2016) that the Lebedev quadrature set can provide the same accuracy as the LSE set and at the same time reduce computation loads to two-thirds of that by the LSE set. The Lebedev quadrature set is determined based on the spherical harmonics expansions to satisfy the invariance under Fig.…”

Single-cell computational analysis of light harvesting in a flat-panel photo-bioreactor

“…5a , which is an integro-differential equation, in MATLAB. The integral part, i.e., the in-scattering term, was solved using the Lebedev Quadrature technique [ 36 , 37 ]. Figure 3 compares solutions of the RTE at different algae concentrations and wavelengths of irradiated light.…”

Fluorescence Molecular Tomography Reconstruction of Small Targets Using Stacked Auto-Encoder Neural Networks

“…FMT problem, describing the photon migration in highly scattering media, is non-linear and can be linearized by using the coupled diffusion equations (DEs) with the Robintype boundary condition as a lower-order approximation. For time-domain FMT with point excitation sources, the coupled diffuse equations provide a description of the forward problem [18], [19].…”