1 / 20

A finite element approach for modeling Diffusion equation

A finite element approach for modeling Diffusion equation. Subha Srinivasan 10/30/09. Forward Problem Definition.

chuong
Download Presentation

A finite element approach for modeling Diffusion equation

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. A finite element approach for modeling Diffusion equation SubhaSrinivasan 10/30/09

  2. Forward Problem Definition • Given a distribution of light sources on the boundary of an object and a distribution of tissue parameters within , to find the resulting measurement set on

  3. Light Propagation in a 3-D Breast Model using BEM

  4. Light Propagation in a 3-D Breast Model using BEM

  5. Inverse Problem Definition • Given a distribution of light sources and a distribution of measurements on the boundary , to derive the distribution of tissue parameters within

  6. Diffusion equation in frequency domain • is the isotropic fluence, is the Diffusion coefficient, is the absorption coefficient and is the isotropic source • is the reduced scattering coefficient

  7. Solutions to Diffusion equation • Analytical solutions exist in simple geometries • Finite difference methods (FDM) use approximations for differentiation and integration. Works well for 2D problems with regular boundaries parallel to coordinate axis, cumbersome for regions with curved or irregular boundaries • Finite element methods (FEM) can be easily applied to complicated and inhomogeneous domains and boundaries. Versatile and computationally feasible (compared to Monte Carlo methods)

  8. Using FEM for Modeling • Main concept: divide a volume/area into elements and build behavior in entire area by characterizing each element (Mosaic) • Uses integral formulation to generate a set of equations • Uses continuous piecewise smooth functions for approximating unknown quantities

  9. Basis Functions For a set of basis functions, we can choose anything. For simplicity here, shown are piecewise linear “hat functions”. Our solution will be a linear combination of these functions. φ3 φ1 φ2 x1=0 x2=L/2 x3=L

  10. Derivation of FEM formulation for Diffusion Equation • The approximate solution is: • And for flux: • Galerkin formulation gives the weighted residual to equal zero: • Galerkin weak form: • Green’s identity: • Substituting:

  11. Matrix form of FEM Model • Discretizing parameters: • Overall: • Matrix form: • For A,B detailed, refer to Paulsen et al, Med Phys, 1995 • Need to apply BCs

  12. Boundary Conditions • Type III BC, Robin type • α incorporates reflection at the boundary due to refractive index change

  13. Source Modeling • Point source: contribution of source to element in which it falls • Gaussian source: modeled with known FWHM • Distributed source model: • Hybrid monte-carlo model: Monte carlo model close to source & diffusion model away from source Paulsen et al, Med Phys, 1995

  14. Forward Model:

  15. Forward Model for Homogeneous Domain: Multiple Sources

  16. Forward Model with Inclusion

  17. Boundary Measurements

  18. Hybrid Source Model Ashley Laughney summer project, 2007

  19. Plots near the source Ashley Laughney summer project, 2007

  20. References • Arridge et al, Med Phys, 20(2), 1993 • Schweiger et al, Med Phys, 22(11), 1995 • Paulsen et al, Med Phys, 22(6), 1995 • Wang et al, JOSA(A), 10(8), 1993

More Related