Yidu Yang

Yidu Yang, Jiayu Han, Hai Bi

The transmission eigenvalue problem is an important and challenging topic arising in the inverse scattering theory. In this paper, for the Helmholtz transmission eigenvalue problem, we give a weak formulation which is a nonselfadjoint linear eigenvalue problem. Based on the weak formulation, we first discuss the non-conforming finite element approximation, and prove the error estimates of the discrete eigenvalues obtained by the Adini element, Morley-Zienkiewicz element, modified-Zienkiewicz element et. al. And we report some numerical examples to validate the efficiency of our approach for solving transmission eigenvalue problem.

Hai Bi, Hao Li, Yidu Yang

This paper proposes and analyzes an a posteriori error estimator for the finite element multi-scale discretization approximation of the Steklov eigenvalue problem. Based on the a posteriori error estimates, an adaptive algorithm of shifted inverse iteration type is designed. Finally, numerical experiments comparing the performances of three kinds of different adaptive algorithms are provided, which illustrate the efficiency of the adaptive algorithm proposed here.

Jiayu Han, Yidu Yang, Hai Bi

Numerical methods for the transmission eigenvalue problems are hot topics in recent years. Based on the work of Lin and Xie [Math. Comp., 84(2015), pp. 71-88], we build a multigrid method to solve the problems. With our method, we only need to solve a series of primal and dual eigenvalue problems on a coarse mesh and the associated boundary value problems on the finer and finer meshes. Theoretical analysis and numerical results show that our method is simple and easy to implement and is efficient for computing real and complex transmission eigenvalues.

Yidu Yang, Jiayu Han, Hai Bi

In this paper, using the linearization technique we write the Helmholtz transmission eigenvalue problem as an equivalent nonselfadjoint linear eigenvalue problem whose left-hand side term is a selfadjoint, continuous and coercive sesquilinear form. To solve the resulting nonselfadjoint eigenvalue problem, we give an $H^{2}$ conforming finite element discretization and establish a two grid discretization scheme. We present a complete error analysis for both discretization schemes, and theoretical analysis and numerical experiments show that the methods presented in this paper can efficiently compute real and complex transmission eigenvalues.

Yu Zhang, Hai Bi, Yidu Yang

We propose a multigrid correction scheme to solve a new Steklov eigenvalue problem in inverse scattering. With this scheme, solving an eigenvalue problem in a fine finite element space is reduced to solve a series of boundary value problems in fine finite element spaces and a series of eigenvalue problems in the coarsest finite element space. And the coefficient matrices associated with those linear systems are constructed to be symmetric and positive definite. We prove error estimates of eigenvalues and eigenfunctions. Numerical results coincide in theoretical analysis and indicate our scheme is highly efficient in solving the eigenvalue problem.

Hai Bi, Yu Zhang, Yidu Yang

In this paper, for a new Stekloff eigenvalue problem which is non-selfadjoint and not $H^1$-elliptic, we establish and analyze two kinds of two-grid discretization scheme and a local finite element scheme. We present the error estimates of approximations of two-grid discretizations. We also prove a local error estimate which is suitable for the case that the local refined region contains singular points lying on the boundary of domain. Numerical experiments are reported finally to show the efficiency of our schemes.

Jiayu Han, Yidu Yang

In this paper we develop an $H^m$-conforming ($m\ge1$) spectral element method on multi-dimensional domain associated with the partition into multi-dimensional rectangles. We construct a set of basis functions on the interval $[-1,1]$ that is made up of the generalized Jacobi polynomials (GJPs) and the nodal basis functions. So the basis functions on multi-dimensional rectangles consist of the tensorial product of the basis functions on the interval $[-1,1]$. Then we construct the spectral element interpolation operator and prove the associated interpolation error estimates. Finally we apply the $H^2$-conforming spectral element method to the Helmholtz transmission eigenvalues that is a hot topic in the field engineering and mathematics.

Liangkun Xu, Shixi Wang, Yidu Yang et al.

In this paper, we discuss a novel higher-order stabilization-free virtual element method for general second-order elliptic eigenvalue problems. Optimal a priori error estimates are derived for both the approximate eigenspace and eigenvalues. Numerical experiments are conducted on regular convex polygonal meshes, convex-concave polygonal meshes, and concave polygonal meshes. The numerical results validate the effectiveness of the proposed method.

Hai Bi, Yidu Yang, Yuanyuan Yu

In this paper we make a further discussion on the finite elements approximation for the Steklov eigenvalue problem on concave polygonal domain. We make full use of the regularity estimate and the characteristic of edge average interpolation operator of nonconforming Crouzeix-Raviart element, which is different from the existing proof argument, and prove a new and optimal error estimate in $\|\cdot\|_{0,\partialΩ}$ for the eigenfunction of linear conforming finite element and the nonconforming Crouzeix-Raviart element, which is an improvement of the current results. Finally, we present some numerical experiments to support the theoretical analysis.

Yao Luo, Yidu Yang

The Helmholtz transmission eigenvalue problem has received much concern in materials science, so it's significant to explore the efficient calculational method of the problem to mathematics and mechanics community. In this paper, based on a variational formulation proposed by Cakon, Monk and Sun, we introduce a mixed-element two-grid discretization and prove error estimates for this method theoretically. Some numerical results are presented to confirm the theoretical analysis and show that the method here is efficient.

Shixi Wang, Hai Bi, Yidu Yang

In this paper, we first discuss the optimal convergence of the adaptive finite element methods for non-self-adjoint eigenvalue problems. We present new theoretical error estimators and computable error estimators for multiple and clustered eigenvalues with the help of the error estimators of finite element solutions for the corresponding source problems, and prove the equivalence between these two estimators. We propose an adaptive algorithm for the eigenvalue cluster and demonstrate that it achieves the optimal convergence rate.We also provide numerical experiments to support our theoretical findings.

Yidu Yang, Yu Zhang, Hai Bi

In this paper, we use the non-conforming Crouzeix-Raviart element method to solve a Stekloff eigenvalue problem arising in inverse scattering. The weak formulation corresponding to this problem is non-selfadjoint and does not satisfy $H^{1}$-elliptic condition,and its Crouzeix-Raviart element discretization does not meet the Strang lemma condition. We use the standard duality techniques to prove an extension of Strang lemma. And we prove the convergence and error estimate of discrete eigenvalues and eigenfunctions using the spectral perturbation theory for compact operators. Finally, we present some numerical examples not only on uniform meshes but also in an adaptive refined meshes to show that the Crouzeix-Raviart method is efficient for computing real and complex eigenvalues as expected.

Hao Li, Yidu Yang

The interior penalty methods using $C^0$ Lagrange elements ($C^0$IPG) developed in the last decade for the fourth order problems are an interesting topic in academia at present. In this paper, we discuss the adaptive fashion of $C^0$IPG method for the Helmholtz transmission eigenvalue problem.We give the a posteriori error indicators for primal and dual eigenfunctions, and prove their reliability and efficiency. We also give the a posteriori error indicator for eigenvalues and design a $C^0$IPG adaptive algorithm. Numerical experiments show that this algorithm is efficient and can get the optimal convergence rate.

Yidu Yang, Hao Li, Hai Bi

In this paper, we prove that the Morley element eigenvalues approximate the exact ones from below on regular meshes, including adaptive local refined meshes, for the fourth-order elliptic eigenvalue problems with the clamped boundary condition in any dimension. And we implement the adaptive computation to obtain lower bounds of the Morley element eigenvalues for the vibration problem of clamped plate under tension.

Yidu Yang, Jiayu Han

Based on the work of Xu and Zhou [Math.Comput., 69(2000), pp.881-909], we establish new three-level and multilevel finite element discretizations by local defect-correction technique. Theoretical analysis and numerical experiments show that the schemes are simple and easy to carry out, and can be used to solve singular nonsymmetric eigenvalue problems efficiently. We also discuss the local error estimates of finite element approximations; it's a new feature here that the estimates apply to the local domains containing corner points.