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, nonselfadjoint, 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. This scheme is based on Lagrange finite elements and is one of the less expensive methods in terms of the amount of degrees of freedom. Due to the nonselfadjointness, the discretization of elastic transmission eigenvalue problem leads to a non-classical mixed method which does not fit into the framework of classical theoretical analysis. In stead, we obtain the convergence analysis based on the spectral approximation theory of compact operators. Numerical examples are presented to verify the theory. Both real and complex eigenvalues can be obtained.