Analytical model of optical fluence inside multiple cylindrical inhomogeneities embedded in an otherwise homogeneous turbid medium for quantitative photoacoustic imaging.

We present an analytical model of optical fluence for multiple cylindrical inhomogeneities embedded in an otherwise homogeneous turbid medium. The model is based on the diffusion equation and represents the optical fluence distribution inside and outside inhomogeneities as a series of modified Bessel functions. We take into account the interplay between cylindrical inhomogeneities by introducing new boundary conditions on the surface of inhomogeneities. The model is compared with the numerical solution of the diffusion equation with NIRFAST software. The fluences inside the inhomogeneities obtained by the two methods are in close agreement. This permits the use of the model as a forward model for quantitative photoacoustic imaging.