A Lowest-Order Mixed Finite Element Method for the Elastic Transmission Eigenvalue Problem

Abstract

The goal of this paper is to develop numerical methods computing a few smallest elastic interior transmission eigenvalues, which are of practical importance in inverse elastic scattering theory. The problem is challenging since it is nonlinear, non-self-adjoint, and of fourth order. In this paper, we construct a lowest-order mixed finite element method which is close to the Ciarlet-Raviart mixed finite element method. The scheme is based on Lagrange finite element and is one of the less expensive methods in terms of the amount of degrees of freedom. Due to the non-self-adjointness, the discretization of elastic transmission eigenvalue problem leads to a non-classical mixed problem which does not fit into the framework of classical theoretical analysis. Instead, we obtain the convergence analysis based on the spectral approximation theory of compact operator. Numerical examples are presented to verify the theory. Both real and complex eigenvalues can be obtained.