Abstract
In this paper, we propose a multi-level finite element method for the transmission eigenvalue problem of anisotropic media. The problem is non-standard and non-self-adjoint with important applications in inverse scattering theory. We employ a suitable finite element method to discretize the problem. The resulting generalized matrix eigenvalue problem is large, sparse and non-Hermitian. To compute the smallest real transmission eigenvalue, which is usually an interior eigenvalue, we devise a multi-level method using Arnoldi iteration. At the coarsest mesh, the eigenvalue is obtained using Arnoldi iteration with an adaptive searching technique. This value is used as the initial guess for Arnoldi iteration at the next mesh level. This procedure is then repeated until the finest mesh level. Numerical examples are presented to show the viability of the method.
Original language | English |
---|---|
Pages (from-to) | 422-435 |
Number of pages | 14 |
Journal | Journal of Computational Physics |
Volume | 255 |
DOIs | |
Publication status | Published - 15 Dec 2013 |
Externally published | Yes |
Keywords
- Anisotropic media
- Arnoldi iteration
- Finite element method
- Transmission eigenvalues