L. Mascotto

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.

L. Mascotto, I. Perugia, A. Pichler

We introduce a novel virtual element method (VEM) for the two dimensional Helmholtz problem endowed with impedance boundary conditions. Local approximation spaces consist of Trefftz functions, i.e., functions belonging to the kernel of the Helmholtz operator. The global trial and test spaces are not fully discontinuous, but rather interelement continuity is imposed in a nonconforming fashion. Although their functions are only implicitly defined, as typical of the VEM framework, they contain discontinuous subspaces made of functions known in closed form and with good approximation properties (plane waves, in our case). We carry out an abstract error analysis of the method, and derive $h$-version error estimates. Moreover, we initiate its numerical investigation by presenting a first test, which demonstrates the theoretical convergence rates.

L. Mascotto, I. Perugia, A. Pichler

We discuss the implementation details and the numerical performance of the recently introduced nonconforming Trefftz virtual element method for the 2D Helmholtz problem. In particular, we present a strategy to significantly reduce the ill-conditioning of the original method; such a recipe is based on an automatic filtering of the basis functions edge by edge, and therefore allows for a notable reduction of the number of degrees of freedom. A widespread set of numerical experiments, including an application to acoustic scattering, the $h$-, $p$-, and $hp$-versions of the method, is presented. Moreover, a comparison with other Trefftz-based methods for the Helmholtz problem shows that this novel approach results in robust and effective performance.

L. Mascotto, A. Pichler

We extend the nonconforming Trefftz virtual element method introduced in arXiv:1805.05634 to the case of the fluid-fluid interface problem, that is, a Helmholtz problem with piecewise constant wave number. With respect to the original approach, we address two additional issues: firstly, we define the coupling of local approximation spaces with piecewise constant wave numbers, secondly, we enrich such local spaces with special functions capturing the physical behaviour of the solution to the target problem. As these two issues are directly related to an increase of the number of degrees of freedom, we use a reduction strategy inspired by arXiv:1807.11237, which allows to mitigate the growth of the dimension of the approximation space when considering $h$- and $p$-refinements. This renders the new method highly competitive in comparison to other Trefftz and quasi-Trefftz technologies tailored for the Helmholtz problem with piecewise constant wave number. A wide range of numerical experiments, including the $p$-version with quasi-uniform meshes and the $hp$-version with isotropic and anisotropic mesh refinements, is presented.

P. F. Antonietti, L. Mascotto, M. Verani

We present a multigrid algorithm for the solution of the linear systems of equations stemming from the $p-$version of the Virtual Element discretization of a two-dimensional Poisson problem. The sequence of coarse spaces are constructed decreasing progressively the polynomial approximation degree of the Virtual Element space, as in standard $p$-multigrid schemes. The construction of the interspace operators relies on auxiliary Virtual Element spaces, where it is possible to compute higher order polynomial projectors. We prove that the multigrid scheme is uniformly convergent, provided the number of smoothing steps is chosen sufficiently large. We also demonstrate that the resulting scheme provides a uniform preconditioner with respect to the number of degrees of freedom that can be employed to accelerate the convergence of classical Krylov-based iterative schemes. Numerical experiments validate the theoretical results.

L. Beirão da Veiga, A. Chernov, L. Mascotto et al.

In the present paper we initiate the study of $hp$ Virtual Elements. We focus on the case with uniform polynomial degree across the mesh and derive theoretical convergence estimates that are explicit both in the mesh size $h$ and in the polynomial degree $p$ in the case of finite Sobolev regularity. Exponential convergence is proved in the case of analytic solutions. The theoretical convergence results are validated in numerical experiments. Finally, an initial study on the possible choice of local basis functions is included.