Francesca Gardini

Francesca Gardini, Gianmarco Manzini, Giuseppe Vacca

We analyse the nonconforming Virtual Element Method (VEM) for the approximation of elliptic eigenvalue problems. The nonconforming VEM allow to treat in the same formulation the two- and three-dimensional case.We present two possible formulations of the discrete problem, derived respectively by the nonstabilized and stabilized approximation of the L^2-inner product, and we study the convergence properties of the corresponding discrete eigenvalue problem. The proposed schemes provide a correct approximation of the spectrum, in particular we prove optimal-order error estimates for the eigenfunctions and the usual double order of convergence of the eigenvalues. Finally we show a large set of numerical tests supporting the theoretical results, including a comparison with the conforming Virtual Element choice.

O. Certik, F. Gardini, G. Manzini et al.

We discuss the $p$- and the $hp$-versions of the virtual element method for the approximation of eigenpairs of elliptic operators with a potential term on polygonal meshes. An application of this model is provided by the Schrödinger equation with a pseudo-potential term. We present in details the analysis of the p-version of the method, proving exponential convergence in the case of analytic eigenfunctions. The theoretical results are supplied with a wide set of experiments. We also show numerically that, in the case of eigenfunctions with finite Sobolev regularity, an exponential approximation of the eigenvalues in terms of the cubic root of the number of degrees of freedom can be obtained by employing $hp$-refinements. Importantly, the geometric flexibility of polygonal meshes is exploited in the construction of the $hp$-spaces.

Ondrej Certik, Francesca Gardini, Gianmarco Manzini et al.

We extend the conforming virtual element method to the numerical resolution of eigenvalue problems with potential terms on a polytopal mesh. An important application is that of the Schrodinger equation with a pseudopotential term. This model is a fundamental element in the numerical resolution of more complex problems from the Density Functional Theory. The VEM is based on the construction of the discrete bilinear forms of the variational formulation through certain polynomial projection operators that are directly computable from the degrees of freedom. The method shows a great flexibility with respect to the meshes and provide a correct spectral approximation with optimal convergence rates. This point is discussed from both the theoretical and the numerical viewpoint. The performance of the method is numerically investigated by solving the Quantum Harmonic Oscillator problem with the harmonic potential and a singular eigenvalue problem with zero potential for the first eigenvalues.

Francesca Gardini, Giuseppe Vacca

We introduce the Virtual Element Method (VEM) for elliptic eigenvalue problems. The main result of the paper states that VEM provides an optimal order approximation of the eigenmodes. A wide set of numerical tests confirm the theoretical analysis.

Daniele Boffi, Francesca Gardini, Lucia Gastaldi

The VEM approximation of eigenvalue problems usually involves the appropriate tuning of stabilization parameters, unless self-stabilizing or stabilization-free VEM are used. In this paper we prove that for elliptic self-adjoint eigenvalue problems the stabilization of the mass matrix is not necessary when lower order standard VEM spaces are adopted. Numerical evidence shows that also for higher order schemes the same result is true on various mesh sequences.

Dibyendu Adak, Daniele Boffi, Francesca Gardini et al.

In this work, we analyze the conforming and $C^0$-non-conforming Virtual Element Method for a fourth-order Steklov eigenvalue problem on a generally shaped, possibly nonconvex, polygonal domain. By employing an {\it enriching } operator, we derive the convergence analysis using the discrete $H^2$ seminorm, and the $H^1$ and $L^2$ norms. We use the Babuška--Osborn spectral theory \cite{BO} to prove that the numerical scheme approximates the spectrum without introducing any spurious eigenvalue. Moreover, we derive the optimal order of convergence for eigenfunctions and double order for eigenvalues. We assess the performance of the method on several numerical tests using different families of polygonal meshes.

Daniele Boffi, Dietmar Gallistl, Francesca Gardini et al.

It is shown that the h-adaptive mixed finite element method for the discretization of eigenvalue clusters of the Laplace operator produces optimal convergence rates in terms of nonlinear approximation classes. The results are valid for the typical mixed spaces of Raviart-Thomas or Brezzi-Douglas-Marini type with arbitrary fixed polynomial degree in two and three space dimensions.

Daniele Boffi, Ricardo G. Durán, Francesca Gardini et al.

In this paper we study an a posteriori error indicator introduced in E. Dari, R.G. Duran, C. Padra, Appl. Numer. Math., 2012, for the approximation of the Laplace eigenvalue problem with Crouzeix-Raviart non-conforming finite elements. In particular, we show that the estimator is robust also in presence of eigenvalues of multiplicity greater than one. Some numerical examples confirm the theory and illustrate the convergence of an adaptive algorithm when dealing with multiple eigenvalues.