Finite element methods for eigenvalue problems have been studied extensively [1, 2, 9]. In this paper, we propose a new finite element approach for the Dirichlet eigenvalue problem. The problem is formulated as an eigenvalue problem of a holomorphic Fredholm operator function . Using Lagrange finite elements, the convergence is proved by the abstract approximation theory for holomorphic operator functions [7, 8].
The new approach has the following characteristics: 1) it provides a new finite element methodology which is different than the classic Babuška-Osborn theory; 2) it can be applied to a large class of nonlinear eigenvalue problems ; and 3) combined with the spectral indicator method [5, 6, 4], it can be parallelized to compute many eigenvalues effectively.
The rest of the paper is arranged as follows. In Section 2, preliminaries of holomorphic Fredholm operator functions and the associated abstract approximation are presented. In Section 3, we reformulate the Dirichlet eigenvalue problem as the eigenvalue problem of a holomorphic Fredholm operator function of index zero. The linear Lagrange finite element is used for discretization and the convergence is proved using the abstract approximation results of Karma [7, 8]. In Section 4, the spectral indicator method [5, 6, 4] is applied to compute the eigenvalues of the unit square.
We present some preliminaries on the eigenvalue approximation theory of holomorphic Fredholm operator functions following [7, 8]. Let be complex Banach spaces, be compact. Denote by the space of bounded linear operators and the set of holomorphic Fredholm operator functions of index zero . Assume that . We consider the problem of finding , such that
The resolvent set and the spectrum of are defined as
Throughout the paper, we assume that . Then the spectrum has no cluster points in , and every is an eigenvalue .
To approximate the eigenvalues of , we consider a sequence of operator functions . Assume the following properties hold.
There exist Banach spaces , and linear bounded mappings with the property
is equibounded on , i.e., there exists a constant such that
approximates for every , i.e.,
is regular for every , i.e., if and is compact in the sense of Karma , then
3 Finite Element Approximation
Let be a bounded Lipschitz domain. The Dirichlet eigenvalue problem is to find and such that
The associated source problem is, given , to find such that
For , the weak formulation of (3.4) is finding such that
Due to the wellposeness of (3.5) (see, e.g., ), there exists a linear compact solution operator such that . The Dirichelt eigenvalue problem (3.3) is equivalent to the operator eigenvalue problem: .
Assume that . Define a nonlinear operator function by
where is the identity operator. Clearly, is a Dirichlet eigenvalue of (3.3) if and only if is an eigenvalue of .
Let be a compact set. Then is a holomorphic Fredholm operator function of index zero.
Proof. Since is compact and is the identity operator, is a Fredholm operator of index zero. Clearly, is holomorphic in .
In the rest of the paper, stands for and is a generic constant. Let be a regular triangular mesh for with mesh size and is the linear Lagrange element space associated with . Then the discrete formulation of (3.5) is to find , such that
where is the -projection operator. Obviously, is bounded and as , for any . Let . The well-posedness of (3.7) implies that .
Let , be the solution operator of (3.7). Define an operator function
The error estimate (3.8) indicates that. This implies that
which is due to for all .
There exists small enough such that, for every compact set
Proof. Let be sufficiently small and . Then
Since is compact and , (3.11) holds.
Let . Then
Proof. From (3.8), we have that
Now we are ready to present the main convergence theorem.
Let . Assume that is small enough. Then there exists such that as . For any sequence the following estimate hold
where is the maximum rank of eigenvectors.
is the maximum rank of eigenvectors.
Proof. Let be a sequence of sufficiently small positive numbers with as and , and . Clearly, (b1) holds with , , and . (b2) and (b3) hold due to Lemma 3.2 and 3.3.
To verify (b4), assume that , with and
for some . We estimate as follows by considering and separately.
If , then exists and is bounded. Let . We have
Recalling from (3.10) it holds
Using (3.14) and Lemma 3.3 we have as .
Assume that . Let
be the finite dimensional generalized eigenspace of. We denote by the projection from to , by the inverse of from to . Due (3.14), we have that
Since is closed, . Let and . Since
by Lemma 3.3, similar to (3.15), we deduce that
On the other hand, since is finite dimensional, there is a subsequence and such that as . Therefore we have
Under the same conditions of the above theorem, it is possible to obtain error estimates for generalized eigenvectors , which is not included in this paper due to simplicity.
4 Numerical Results
Let be the unit square . The smallest eigenvalue is and its rank . We use the linear Lagrange element on a series of uniformly refined meshes for discretization. The spectral indicator method [5, 6, 4] is employed to compute the smallest eigenvalue of (3.9). The results are shown in Table 1, which confirms the second order convergence.
The research of B. Gong is supported partially by China Postdoctoral Science Foundation Grant 2019M650460. The research of J. Sun is supported partially by MTU REF. The research of Z. Zhang is supported partially by the National Natural Science Foundation of China grants NSFC 11871092, NSAF U1930402, and NSF 11926356.
-  I. Babuška and J. Osborn, Eigenvalue Problems, Handbook of Numerical Analysis, Vol. II, Elsevier Science Publishers, North-Holland, 1991.
-  D. Boffi, Finite element approximation of eigenvalue problems. Acta Numer. 19 (2010), 1-120.
-  I. Gohberg and J. Leiterer, Holomorphic operator functions of one variable and applications. 192. Birkhäuser Verlag, Basel, 2009.
-  B. Gong, J. Sun, T. Turner, C. Zheng, Finite element approximation of the nonlinear transmission eigenvalue problem for anisotropic media. preprint, 2019.
-  R. Huang, A. Struthers, J. Sun, and R. Zhang, Recursive integral method for transmission eigenvalues. J. Comput. Phys. 327 (2016), 830-840.
-  R. Huang, J. Sun and C. Yang, Recursive integral method with Cayley transformation. Numer. Linear Algebra Appl., 25 (2018), no. 6, e2199.
-  O. Karma, Approximation in eigenvalue problems for holomorphic Fredholm operator functions. I. Numer. Funct. Anal. Optim. 17 (1996), no. 3-4, 365-387.
-  O. Karma, Approximation in eigenvalue problems for holomorphic Fredholm operator functions. II. (Convergence rate). Numer. Funct. Anal. Optim. 17 (1996), no. 3-4, 389-408.
-  J. Sun and A. Zhou, Finite element methods for eigenvalue problem. CRC Press, Taylor Francis Group, Boca Raton, London, New York, 2016.