Implementation of the radiative transfer equation on block-structured grids for modeling fluorescence light propagation in tissue with arbitrary shape

We developed a method for solving the fluorescence equation of radiative transfer in the frequency domain on blockstructured grids. In this way fluorescence light propagation in arbitrarily shaped tissue can be modeled with high accuracy without compromising on the convergence speed of these codes. The block-structure grid generator is developed as a multi-purpose tool that can be used with many numerical schemes. We present results from numerical studies that show that it is possible to resolve curved boundaries with grids that maintain much of the intrinsic structure of Cartesian grids. The natural ordering of this grid allows for simplified algorithms. In simulation studies we found that we can reduce the error in boundary fluence by a factor of five by using a two-level block structured grid. The increase in computational cost is only two-fold. We compare benchmark solutions to results with various levels of refinement, boundary conditions, and different geometries.

[1]  Alexander D. Klose,et al.  Radiative transfer of luminescence light in biological tissue , 2009 .

[2]  Andreas H Hielscher,et al.  Optical tomographic imaging of small animals. , 2005, Current opinion in biotechnology.

[3]  William S. Russell,et al.  Polynomial interpolation schemes for internal derivative distributions on structured grids , 1995 .

[4]  A H Hielscher,et al.  Iterative reconstruction scheme for optical tomography based on the equation of radiative transfer. , 1999, Medical physics.

[5]  Andreas H. Hielscher,et al.  Three-dimensional optical tomography with the equation of radiative transfer , 2000, J. Electronic Imaging.

[6]  Andreas H. Hielscher,et al.  A PDE-constrained SQP algorithm for optical tomography based on the frequency-domain equation of radiative transfer , 2008 .

[7]  Guillaume Bal,et al.  Algorithm for solving the equation of radiative transfer in the frequency domain. , 2004, Optics letters.

[8]  L H Howell.  A Discrete Ordinates Algorithm for Radiation Transport Using Block-Structured Adaptive Mesh Refinement , 2003 .

[9]  G. C. Pomraning,et al.  Linear Transport Theory , 1967 .

[10]  Jürgen Beuthan,et al.  Optimal source-modulation frequencies for transport-theory-based optical tomography of small-tissue volumes. , 2008, Optics express.

[11]  Vasilis Ntziachristos,et al.  The inverse source problem based on the radiative transfer equation in optical molecular imaging , 2005 .

[12]  M. Berger,et al.  Adaptive mesh refinement for hyperbolic partial differential equations , 1982 .

[13]  Andreas H Hielscher,et al.  Frequency-domain sensitivity analysis for small imaging domains using the equation of radiative transfer. , 2007, Applied optics.

[14]  Alexander D Klose,et al.  Fluorescence tomography with simulated data based on the equation of radiative transfer. , 2003, Optics letters.

[15]  P. Colella,et al.  An Adaptive Mesh Refinement Algorithm for the Radiative Transport Equation , 1998 .