A Numerical Simulation of the Diffusion of Near Infrared Light Through Brain Tissue

A numerical simulation of the transport of light energy in the near infrared region of the spectrum through human brain tissue is presented. This simulation models the use of near infrared spectroscopic techniques to quantify the levels of oxygen present in brain tissue. Successful application of the technique requires knowledge of the optical pathlength in the tissue, and it is the goal of this simulation to quantify the relationship of the optical pathlength and the oxygenation state of the tissue. Both implicit and explicit finite element schemes for unstructured grids are implemented and discussed. Several application simulations using three tissue grids of varying degrees of physiological accuracy are then conducted, and figures for the optical pathlength of light through the tissue at varying levels of oxygenation are computed. These results are then used to develop a quantitative relationship between the pathlength and the absorption parameter for the three tissue models which we explore.