In this paper, an innovative and effective numerical algorithm by the use of weak Galerkin (WG) finite element methods is proposed for a type of fourth order problem arising from fluorescence tomography. Fluorescence tomography is emerging as an in vivo non-invasive 3D imaging technique reconstructing images that characterize the distribution of molecules tagged by fluorophores. Weak second order elliptic operator and its discrete version are introduced for a class of discontinuous functions defined on a finite element partition of the domain consisting of general polygons or polyhedra. An error estimate of optimal order is derived in an Hκ2-equivalent norm for the WG finite element solutions. Error estimates of optimal order except the lowest order finite element in the usual L2 norm are established for the WG finite element approximations. Numerical tests are presented to demonstrate the accuracy and efficiency of the theory established for the WG numerical algorithm.
- Finite element methods
- Fluorescence tomography
- Fourth order problem
- Polygonal or polyhedral meshes
- Weak Galerkin
- Weak second order elliptic operator