Abstract
We present a method for handling nonscattering regions within diffusing domains. The method develops from an iterative radiosity-diffusion approach using Green's functions that was computationally slow. Here we present an improved implementation using a finite element method (FEM) that is direct. The fundamental idea is to introduce extra equations into the standard diffusion FEM to represent nondiffusive light propagation across a nonscattering region. By appropriate mesh node ordering the computational time is not much greater than for diffusion alone. We compare results from this method with those from a discrete ordinate transport code, and with Monte Carlo calculations. The agreement is very good, and, in addition, our scheme allows us to easily model time-dependent and frequency domain problems.
Original language | English |
---|---|
Pages (from-to) | 252-264 |
Number of pages | 13 |
Journal | Medical Physics |
Volume | 27 |
Issue number | 1 |
DOIs | |
Publication status | Published - Jan 2000 |
Keywords
- Diffusion
- Finite element method
- Light propagation
- Transport equation
- Voids
ASJC Scopus subject areas
- Biophysics
- Radiology Nuclear Medicine and imaging