A Multilevel Correction Method for Steklov Eigenvalue Problem by Nonconforming Finite Element Methods

In this paper, a multilevel correction scheme is proposed to solve the Steklov eigenvalue problem by nonconforming finite element methods. With this new scheme, the accuracy of eigenpair approximations can be improved after each correction step which only needs to solve a source problem on finer finite element space and an Steklov eigenvalue problem on the coarsest finite element space. This correction scheme can increase the overall efficiency of solving eigenvalue problems by the nonconforming finite element method. Furthermore, as same as the direct eigenvalue solving by nonconforming finite element methods, this multilevel correction method can also produce the lower-bound approximations of the eigenvalues.     

Adaptive finite element methods for computing band gaps in photonic crystals

In this paper we propose and analyse adaptive finite element methods for computing the band structure of 2D periodic photonic crystals. The problem can be reduced to the computation of the discrete spectra of each member of a family of periodic Hermitian eigenvalue problems on a unit cell, parametrised by a two-dimensional parameter - the quasimomentum. These eigenvalue problems involve non-coercive elliptic operators with generally discontinuous coefficients and are solved by adaptive finite elements. We propose an error estimator of residual type and show it is reliable and efficient for each eigenvalue problem in the family. In particular we prove that if the error estimator converges to zero then the distance of the computed eigenfunction from the true eigenspace also converges to zero and the computed eigenvalue converges to a true eigenvalue with double the rate. We also prove that if the distance of a computed sequence of approximate eigenfunctions from the true eigenspace approaches zero, then so must the error estimator. The results hold for eigenvalues of any multiplicity. We illustrate the benefits of the resulting adaptive method in practice, both for fully periodic structures and also for the computation of eigenvalues in the band gap of structures with defect, using the supercell method.     

Extrapolation of Mixed Finite Element Approximations for the Maxwell Eigenvalue Problem

In this paper,a general method to derive asymptotic error expansion formulas for the mixed finite element approximations of the Maxwell eigenvalue problem is established.Abstract lemmas for the error of the eigenvalue approximations are obtained.Based on the asymptotic error expansion formulas,the Richardson extrapolation method is employed to improve the accuracy of the approximations for the eigenvalues of the Maxwell system from θ(h2) to θ(h4) when applying the lowest order Nédé1ec mixed finite element and a nonconforming mixed finite element.To our best knowledge,this is the first superconvergence result of the Maxwell eigenvalue problem by the extrapolation of the mixed finite element approximation.Numerical experiments are provided to demonstrate the theoretical results.     

Extremal Eigenvalues of the Sturm-Liouville Problems with Discontinuous Coefficients

In this paper, an extremal eigenvalue problem to the Sturm-Liouville equations with discontinuous coefficients and volume constraint is investigated. Liouville transformation is applied to change the problem into an equivalent minimization problem. Finite element method is proposed and the convergence for the finite element solution is established. A monotonic decreasing algorithm is presented to solve the extremal eigenvalue problem. A global convergence for the algorithm in the continuous case is proved. A few numerical results are given to depict the efficiency of the method.     

Wilson元特征值下逼近准确特征值

该文讨论矩形域上Laplace算子特征值问题有限元近似.证明了Wilson非协调有限元特征值下逼近准确特征值,从而解决了有限元法中长期存在的一个猜想.     

Approximation of an Eigenvalue Problem Associated with the Stokes Problem by the Stream Function-Vorticity-Pressure Method

By means of eigenvalue error expansion and integral expansion techniques, we propose and analyze the stream function-vorticity-pressure method for the eigenvalue problem associated with the Stokes equations on the unit square. We obtain an optimal order of convergence for eigenvalues and eigenfuctions. Furthermore, for the bilinear finite element space, we derive asymptotic expansions of the eigenvalue error, an efficient extrapolation and an a posteriori error estimate for the eigenvalue. Finally, numerical experiments are reported. The first author was supported by China Postdoctoral Sciences Foundation.     

On a variational approximation method for a class of elliptic eigenvalue problems in composite structures

We consider a second-order elliptic eigenvalue problem on a convex polygonal domain, divided in nonoverlapping subdomains. The conormal derivative of the unknown function is continuous on the interfaces, while the function itself is discontinuous. We present a general finite element method to obtain a numerical solution of the eigenvalue problem, starting from a nonstandard formally equivalent variational formulation in an abstract setting in product Hilbert spaces. We use standard Lagrange finite element spaces on the subdomains. Moreover, the bilinear forms are approximated by suitable numerical quadrature formulas. We obtain error estimates for both the eigenfunctions and the eigenvalues, allowing for the case of multiple exact eigenvalues, by a pure variational method.

     

带间断系数椭圆问题的P1非协调四边形元的加性 Schwarz 方法

本文讨论了带间断系数的二阶椭网问题的P1非协调四边形元的加性Schwarz方法.通过分析加性Schwaxz预处理后系统的特征值分布,我们证明了除少数小特征值外,其余所有特征值都有正的关于间断系数和网格尺寸拟一致的上下界.数值试验验证了我们的结论.     

Eigenvalues of the discretized Navier-Stokes equation with application to the detection of Hopf bifurcations

This paper is concerned with a numerical approach to the problem of finding the leftmost eigenvalues of large sparse nonsymmetric generalised eigenvalue problems which arise in stability studies of incompressible fluid flow problems. The matrices have a special block structure that is typical of mixed finite element discretizations for such problems. The numerical approach is an extension of the hybrid technique introduced by Saad [22] and utilizes the idea of preconditioning the eigenvalue problem before applying Arnoldi's method. Two preconditioners, one a modified Cayley transform, the other a Chebyshev polynomial transform, are compared in numerical experiments on a double diffusive convection problem and the Cayley transform proves superior. The Cayley transform is then used to provide numerical results for the finite Taylor problem.     

THE LOWER APPROXIMATION OF EIGENVALUE BY LUMPED MASS FINITE ELEMENT METHOD

In the present paper, we investigate properties of lumped mass finite element method (LFEM hereinafter) eigenvalues of elliptic problems. We propose an equivalent formulation of LFEM and prove that LFEM eigenvalues are smaller than the standard finite element method (SFEM hereinafter) eigenvalues. It is shown, for model eigenvalue problems with uniform meshes, that LFEM eigenvalues are not greater than exact solutions and that they are increasing functions of the number of elements of the triangulation, and numerical examples show that this result equally holds for general problems.     

Computing the lower and upper bounds of Laplace eigenvalue problem: by combining conforming and nonconforming finite element methods

We introduce some ways to compute the lower and upper bounds of the Laplace eigenvalue problem.By using the special nonconforming finite elements,i.e.,enriched Crouzeix-Raviart element and extended Q1ro t,we get the lower bound of the eigenvalue.Additionally,we use conforming finite elements to do the postprocessing to get the upper bound of the eigenvalue,which only needs to solve the corresponding source problems and a small eigenvalue problem if higher order postprocessing method is implemented.Thus,we can obtain the lower and upper bounds of the eigenvalues simultaneously by solving eigenvalue problem only once.Some numerical results are also presented to demonstrate our theoretical analysis.     

A multilevel preconditioner for the C-R FEM for elliptic problems with discontinuous coefficients

In this paper, we propose a multilevel preconditioner for the Crouzeix-Raviart finite element approximation of second-order elliptic partial differential equations with discontinuous coefficients. Since the finite element spaces are nonnested, weighted intergrid transfer operators, which are stable under the weighted L2 norm, are introduced to exchange information between different meshes. By analyzing the eigenvalue distribution of the preconditioned system, we prove that except a few small eigenvalues, all the other eigenvalues are bounded below and above nearly uniformly with respect to the jump and the mesh size. As a result, we get that the convergence rate of the preconditioned conjugate gradient method is quasi-uniform with respect to the jump and the mesh size. Numerical experiments are presented to confirm our theoretical analysis.     

Boundary element approximation for Maxwell's eigenvalue problem

We introduce a new method for computing eigenvalues of the Maxwell operator with boundary finite elements. On bounded domains with piecewise constant material coefficients, the Maxwell solution for fixed wave number can be represented by boundary integrals, which allows to reduce the eigenvalue problem to a nonlinear problem for determining the wave number along with boundary and interface traces. A Galerkin discretization yields a smooth nonlinear matrix eigenvalue problem that is solved by Newton's method or, alternatively, the contour integral method. Several numerical results including an application to the band structure computation of a photonic crystal illustrate the efficiency of this approach. Copyright © 2013 John Wiley & Sons, Ltd.     

The extrapolation of numerical eigenvalues by finite elements for differential operators

This paper discusses the extrapolation of numerical eigenvalues by finite elements for differential operators and obtains the following new results: (a) By extending a theorem of eigenvalue error estimate, which was established by Osborn, a new expansion of eigenvalue error is obtained. Many achievements, which are about the asymptotic expansions of finite element methods of differential operator eigenvalue problems, are brought into the framework of functional analysis. (b) The Richardson extrapolation of nonconforming finite elements for multiple eigenvalues and splitting extrapolation of finite elements based on domain decomposition of non-selfadjoint differential operators for multiple eigenvalues are achieved. In addition, numerical examples are provided to support the theoretical analysis.     

Lower bounds for eigenvalues and postprocessing by an integral type nonconforming FEM

In this paper, we analyze some approximation properties of a nonconforming piecewise linear finite element with integral degrees of freedom. A nonconforming finite element method (FEM) is applied to second-order eigenvalue problems (EVPs).We prove that the eigenvalues computed by means of this element are smaller than the exact ones if the mesh size is small enough. The case when an EVP is defined on a nonconvex domain is considered. A superconvergent rate is established to a second-order elliptic problem by the introduction of nonstandard interpolated elements based on the integral type linear element. A simple postprocessing method applied to second-order EVPs is also proposed and analyzed. Finally, computational aspects are discussed and numerical examples are presented.     

Distribution of eigenvalues for the discontinuous boundary-value problem with functional-manypoint conditions

In this study, we investigate the boundary-value problem with eigenvalue parameter generated by the differential equation with discontinuous coefficients and boundary conditions which contains not only endpoints of the considered interval, but also a point of discontinuity, a finite number of internal points and abstract linear functionals. So our problem is not a pure boundary-value one. We single out a class of linear functionals and find simple algebraic conditions on the coefficients which guarantee the existence of an infinite number of eigenvalues. Also, the asymptotic formulas for the eigenvalues are found. The results obtained in this paper are new, even in the case of boundary conditions either without internal points or without linear functionals.     

Boundary layer analysis for the eigenvalues of hanging rod

The motion of a naturally straight inextensible flexible elastic hanging rod is formulated and then linearized about the straight solution. To solve this equation by separation of variables, an eigenvalue problem is derived. When the stiffness of the rod is small, the eigenvalue equation is a singular perturbation problem. This paper is devoted to solving this eigenvalue problem by boundary layer analysis when the stiffness is suitably small, especially on the analytic approximate solutions of the first several eigenvalues and eigenfunctions. The first three eigenvalues are also compared with the numerical results computed by a finite difference method. The excellent agreement shows the efficiency of the boundary layer analysis.     

A mixed method of subspace iteration for Dirichlet eigenvalue problems

A full multigrid scheme was used in computing some eigenvalues of the Laplace eigenvalue problem with the Dirichlet boundary condition. We get a system of algebraic equations with an aid of finite difference method and apply subspace iteration method to the system to compute first some eigenvalues. The result shows that this is very effective in calculating some eigenvalues of this problem.     

Divergence-free finite elements for the numerical solution of a hydroelastic vibration problem

In this paper, we analyze a divergence-free finite element method to solve a fluid–structure interaction spectral problem in the three-dimensional case. The unknowns of the resulting formulation are the fluid and solid displacements and the fluid pressure on the interface separating both media. The resulting mixed eigenvalue problem is approximated by using appropriate basis of the divergence-free lowest order Raviart–Thomas elements for the fluid, piecewise linear elements for the solid and piecewise constant elements for the interface pressure. It is proved that eigenvalues and eigenfunctions are correctly approximated and some numerical results are reported in order to assess the performance of the method.     

DERIVATIVES OF EIGENPAIRS OF SYMMETRIC QUADRATIC EIGENVALUE PROBLEM

Derivatives of eigenvalues and eigenvectors with respect to parameters in symmetric quadratic eigenvalue problem are studied. The first and second order derivatives of eigenpairs are given. The derivatives are calculated in terms of the eigenvalues and eigenvectors of the quadratic eigenvalue problem, and the use of state space representation is avoided, hence the cost of computation is greatly reduced. The efficiency of the presented method is demonstrated by considering a spring-mass-damper system.