Thermoacoustic tomography forward modeling with the spectral element method.


Journal Article

A thermoacoustic tomography (TAT) forward solver based on the spectral element method (SEM) with perfectly matched layer absorbing boundary condition has been developed. The TAT forward solver is intended to model acoustically inhomogeneous media with high accuracy in the frequency domain. The high-order basis functions used in the SEM are Gauss-Lobatto-Legendre (GLL) polynomials. Due to the orthogonality of the GLL basis functions and GLL quadrature integration, the mass matrix is diagonal and the stiffness matrix is sparse. Thus, the proposed method greatly reduces the memory requirement and computational time in comparison with the conventional finite element method (FEM). Numerical results show that the high-order SEM is able to achieve the same accuracy as the FEM but with a much smaller number of unknowns. Therefore, the TAT forward solver based on SEM is able to simulate a large-scale and realistic TAT problem.

Full Text

Duke Authors

Cited Authors

  • Lim, KH; Lee, J-H; Liu, QH

Published Date

  • January 2008

Published In

Volume / Issue

  • 35 / 1

Start / End Page

  • 4 - 12

PubMed ID

  • 18293555

Pubmed Central ID

  • 18293555

Electronic International Standard Serial Number (EISSN)

  • 2473-4209

International Standard Serial Number (ISSN)

  • 0094-2405

Digital Object Identifier (DOI)

  • 10.1118/1.2805478


  • eng