LAUSR

The finite-difference time-domain (FDTD) method, which is often used to solve Maxwell’s equations, has been proposed to solve the optical diffusion equations to simulate light propagation in turbid media such as biological tissue. Finite-difference methods can be numerically unstable and calculation can diverge unless a stability condition is satisfied. Until now, the stability condition for FDTD of the optical diffusion equations has not been derived. In this paper, we derive a stability condition for FDTD of the optical diffusion equations using the von Neumann stability analysis and verify the condition numerically.