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.

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.

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.