Abstract:In this paper, a new and efficient numerical algorithm by using weak Galerkin (WG) finite element methods is proposed for a type of fourth order problem arising from fluorescence tomography(FT). Fluorescence tomography is an emerging, in vivo non-invasive 3-D imaging technique which reconstructs images that characterize the distribution of molecules that are 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 in the usual $L^2$ norm are established, yielding optimal order of convergence for all the WG finite element algorithms except the one corresponding to the lowest order (i.e., piecewise quadratic elements). Some numerical experiments are presented to illustrate the efficiency and accuracy of the numerical scheme.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to develop a new and efficient numerical algorithm for solving a class of fourth - order problems originating from Fluorescence Tomography (FT). Specifically, a Weak Galerkin Finite Element Method (WG - FEM) is proposed in the paper to deal with fourth - order partial differential equations with Dirichlet and Neumann boundary conditions.
### Problem Background
Fluorescence Tomography is an emerging, non - invasive three - dimensional imaging technique. It captures molecule - specific information by using highly specific fluorescent probes and non - ionizing near - infrared (NIR) radiation. Compared with other medical imaging modalities such as X - rays or strong magnetic fields, Fluorescence Tomography is considered a potentially less - harmful medical imaging modality. Its main goal is to reconstruct the distribution of fluorescently - labeled target molecules from boundary measurements, so it has important application prospects in early cancer detection and drug monitoring, etc.
### Mathematical Model
The specific mathematical model studied in the paper is:
\[
(-\nabla\cdot(\kappa\nabla)+\mu)^2u = f\quad\text{in }\Omega,
\]
\[
u=\xi\quad\text{on }\partial\Omega,
\]
\[
\kappa\nabla u\cdot n=\nu\quad\text{on }\partial\Omega,
\]
where:
- $\Omega$ is an open bounded region located in $\mathbb{R}^d$ ($d = 2,3$) with a Lipschitz continuous boundary $\partial\Omega$.
- $n$ is the unit outer normal vector of $\partial\Omega$.
- $\kappa$ is a symmetric positive - definite matrix - valued function.
- $\mu$ is a non - negative real - valued function.
- $f$, $\xi$, and $\nu$ are given functions or known values on the boundary.
### Solution
To solve the above fourth - order problem, the paper proposes the following several innovative points:
1. **Introducing the Weak Second - Order Elliptic Operator**: A weak second - order elliptic operator suitable for discontinuous functions and its discrete version are defined, which enables calculations on polygonal or polyhedral meshes.
2. **Stability Term**: A stability term that takes into account the jumps of the coefficient matrix is introduced to improve the stability and accuracy of the algorithm.
3. **Variational Form**: Based on a specific variational form (rather than the variational form of the standard biharmonic equation), the Weak Galerkin Finite Element Method is designed, thus avoiding the C1 continuity requirement in traditional methods.
### Main Contributions
- A new Weak Galerkin finite element algorithm applicable to fourth - order problems arising from Fluorescence Tomography is proposed.
- The optimal - order error estimates are derived, and good convergence results are obtained in both the H2 - equivalent norm and the conventional L2 norm.
- The effectiveness and accuracy of the algorithm are verified through numerical experiments.
In summary, this paper aims to provide an efficient and accurate numerical algorithm through the Weak Galerkin finite element method to solve the fourth - order partial differential equation problems in Fluorescence Tomography, thereby promoting the further development of this field.