No abstract

Abstract. We introduce the nonconforming Virtual Element Method (VEM) for the approximation of second order elliptic problems. We present the construction of the new element in two and three dimensions, highlighting the main differences with the conforming VEM and the classical nonconforming finite element methods. We provide the error analysis and establish the equivalence with a family of mimetic finite difference methods.

In this paper we analyze a discontinuous finite element method recently introduced by Bassi and Rebay for the approximation of elliptic problems. Stability and error estimates in various norms are proven.

We present in a unified framework new conforming and nonconforming Virtual Element Methods (VEM) for general second order elliptic problems in two and three dimensions. The differential operator is split into its symmetric and non-symmetric parts and conditions for stability and accuracy on their discrete counterparts are established. These conditions are shown to lead to optimal H 1 -and L 2 -error estimates, confirmed by numerical experiments on a set of polygonal meshes. The accuracy of the numerical approximation provided by the two methods is shown to be comparable.

Generalized barycentric coordinates such as Wachspress and mean value coordinates have been used in polygonal and polyhedral finite element methods. Recently, mimetic finite difference schemes were cast within a variational framework, and a consistent and stable finite element method on arbitrary polygonal meshes was devised. The method was coined as the virtual element method (VEM), since it did not require the explicit construction of basis functions. This advance provides a more in-depth understanding of mimetic schemes, and also endows polygonal-based Galerkin methods with greater flexibility than three-node and four-node finite element methods. In the VEM, a projection operator is used to realize the decomposition of the stiffness matrix into two terms: a consistent matrix that is known, and a stability matrix that must be positive semi-definite and which is only required to scale like the consistent matrix. In this paper, we first present an overview of previous developments on conforming polygonal and polyhedral finite elements, and then appeal to the exact decomposition in the VEM to obtain a robust and efficient generalized barycentric coordinate-based Galerkin method on polygonal and polyhedral elements. The consistent matrix of the VEM is adopted, and numerical quadrature with generalized barycentric coordinates is used to compute the stability matrix. This facilitates post-processing of field variables and visualization in the VEM, and on the other hand, provides a means to exactly satisfy the patch test with efficient numerical integration in polygonal and polyhedral finite elements. We present numerical examples that demonstrate the sound accuracy and performance of the proposed method. For Poisson problems in ℝ2 and ℝ3, we establish that linearly complete generalized barycentric interpolants deliver optimal rates of convergence in the L2-norm and the H1-seminorm.

In this paper we address the numerical approximation of linear fourth-order elliptic problems on polygonal meshes. In particular, we present a novel nonconforming virtual element discretization of arbitrary order of accuracy for biharmonic problems. The approximation space is made of possibly discontinuous functions, thus giving rise to the fully nonconforming virtual element method. We derive optimal error estimates in a suitable (broken) energy norm and present numerical results to assess the validity of the theoretical estimates.

A posteriori error estimation and adaptivity are very useful in the context of the virtual element and mimetic discretization methods due to the flexibility of the meshes to which these numerical schemes can be applied. Nevertheless, developing error estimators for virtual and mimetic methods is not a straightforward task due to the lack of knowledge of the basis functions. In the new virtual element setting, we develop a residual based a posteriori error estimator for the Poisson problem with (piecewise) constant coefficients, that is proven to be reliable and efficient. We moreover show the numerical performance of the proposed estimator when it is combined with an adaptive strategy for the mesh refinement.

SUMMARYIn this paper, we establish the connections between the virtual element method (VEM) and the hourglass control techniques that have been developed since the early 1980s to stabilize underintegrated C 0 Lagrange finite element methods. In the virtual element method, the bilinear form is decomposed into two parts: a consistent term that reproduces a given polynomial space and a correction term that provides stability. The essential ingredients of C 0 -continuous virtual element methods on polygonal and polyhedral meshes are described, which reveals that the variational approach adopted in the VEM affords a generalized and robust means to stabilize underintegrated finite elements. We focus on the heat conduction (Poisson) equation, and present a virtual element approach for the isoparametric four-node quadrilateral and the eight-node hexahedral element. In addition, we show quantitative * Correspondence to: N. Sukumar,