skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Conforming and nonconforming virtual element methods for elliptic problems

Abstract

Here we present, in a unified framework, new conforming and nonconforming virtual element methods for general second-order elliptic problems in two and three dimensions. The differential operator is split into its symmetric and nonsymmetric 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.

Authors:
 [1];  [2];  [1]
  1. Univ. of Leicester, Leicester (United Kingdom)
  2. Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
Publication Date:
Research Org.:
Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
Sponsoring Org.:
USDOE Laboratory Directed Research and Development (LDRD) Program
OSTI Identifier:
1331260
Report Number(s):
LA-UR-15-23951
Journal ID: ISSN 0272-4979
Grant/Contract Number:
AC52-06NA25396
Resource Type:
Journal Article: Accepted Manuscript
Journal Name:
IMA Journal of Numerical Analysis
Additional Journal Information:
Journal Name: IMA Journal of Numerical Analysis; Journal ID: ISSN 0272-4979
Publisher:
Oxford University Press/Institute of Mathematics and its Applications
Country of Publication:
United States
Language:
English
Subject:
97 MATHEMATICS AND COMPUTING; mathematics; elliptic problems; virtual element methods; polygonal and polyhedral meshes; convection-diffusion-reaction equations

Citation Formats

Cangiani, Andrea, Manzini, Gianmarco, and Sutton, Oliver J. Conforming and nonconforming virtual element methods for elliptic problems. United States: N. p., 2016. Web. doi:10.1093/imanum/drw036.
Cangiani, Andrea, Manzini, Gianmarco, & Sutton, Oliver J. Conforming and nonconforming virtual element methods for elliptic problems. United States. doi:10.1093/imanum/drw036.
Cangiani, Andrea, Manzini, Gianmarco, and Sutton, Oliver J. Wed . "Conforming and nonconforming virtual element methods for elliptic problems". United States. doi:10.1093/imanum/drw036. https://www.osti.gov/servlets/purl/1331260.
@article{osti_1331260,
title = {Conforming and nonconforming virtual element methods for elliptic problems},
author = {Cangiani, Andrea and Manzini, Gianmarco and Sutton, Oliver J.},
abstractNote = {Here we present, in a unified framework, new conforming and nonconforming virtual element methods for general second-order elliptic problems in two and three dimensions. The differential operator is split into its symmetric and nonsymmetric parts and conditions for stability and accuracy on their discrete counterparts are established. These conditions are shown to lead to optimal H1- and L2-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.},
doi = {10.1093/imanum/drw036},
journal = {IMA Journal of Numerical Analysis},
number = ,
volume = ,
place = {United States},
year = {Wed Aug 03 00:00:00 EDT 2016},
month = {Wed Aug 03 00:00:00 EDT 2016}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

Citation Metrics:
Cited by: 6works
Citation information provided by
Web of Science

Save / Share:
  • This thesis deals with the condition numbers and singular value distribution of the preconditioned operators B[sub h][sup [minus]1]A[sub h] and A[sub h]B[sub h][sup [minus]1], where A[sub h] and B[sub h] are nonconforming finite element discretizations of second-order elliptic operators A and B. It generalizes the works of Manteuffel and Parter, Goldstein, Manteuffel and Parter, as well as Parter and Wong. Three nonconforming finite element methods are considered. The first two methods, the penalty method and the method of nearly-zero boundary conditions, deal with Dirichlet boundary conditions on curved domains. It is shown that if the leading part of A ismore » a smooth function times the leading part of B and if the boundary condition of B is the same as that of A, then the L[sup 2]-condition number of A[sub h] B[sub h][sup [minus]1] is bounded independent of the mesh size h, and its L[sup 2]-singular values cluster and fill-in some bounded, estimable interval. If the boundary condition of the adjoint of B is the same as that of the adjoint of A, then the above conclusion holds for B[sub h][sup [minus]1]A[sub h] instead of A[sub h]B[sub h][sup [minus]1]. If B is self-adjoint and positive definite, then the B[sub h]-condition number of B[sub h][sup [minus]1]A[sub h] is bounded independent of h and the B[sub h]-singular values of B[sub h][sup [minus]1]A[sub h] cluster in some bounded, estimable interval. None of the results require full H[sup 2]-regularity nor uniform grids. The third method uses the P[sub 1] nonconforming finite elements, which consist of piecewise linear functions continuous only at mid-points of interelement boundaries. The generalized results of Parter and Wong are valid for this method with quasi-uniform grids and simpler boundary conditions. Through these three particular examples, one has learned what difficulties to anticipate when dealing with nonconforming finite element methods in general.« less
  • In this paper, we present the nonconforming virtual element method (VEM) for the numerical approximation of velocity and pressure in the steady Stokes problem. The pressure is approximated using discontinuous piecewise polynomials, while each component of the velocity is approximated using the nonconforming virtual element space. On each mesh element the local virtual space contains the space of polynomials of up to a given degree, plus suitable nonpolynomial functions. The virtual element functions are implicitly defined as the solution of local Poisson problems with polynomial Neumann boundary conditions. As typical in VEM approaches, the explicit evaluation of the non-polynomial functionsmore » is not required. This approach makes it possible to construct nonconforming (virtual) spaces for any polynomial degree regardless of the parity, for two- and three-dimensional problems, and for meshes with very general polygonal and polyhedral elements. We show that the nonconforming VEM is inf-sup stable and establish optimal a priori error estimates for the velocity and pressure approximations. Finally, numerical examples confirm the convergence analysis and the effectiveness of the method in providing high-order accurate approximations.« less
  • In this paper, we study a coupling of spectral methods and the [rho]-version finite element methods for elliptic boundary value problems containing singularities. The method of auxiliary mapping, which is a recent development to deal with domain singularities in the [rho]-version of the finite element method, is employed to remove the pollution effect caused by singularities. An iterative interfacial coupling between spectral methods and the [rho]-version of the finite element method is used and investigated numerically. The advantages of such an approach are demonstrated by the high accuracy of spectral methods for the smooth part of solutions and the flexibilitymore » of the [rho]-version of the finite element method for dealing with singularities and irregular domains. The efficiency of the coupling method is also evaluated by comparing results obtained by this method with those obtained by the full finite element algorithm. 26 refs., 9 figs., 6 tabs.« less
  • Engineering applications frequently require the numerical solution of elliptic boundary value problems in irregularly shaped domains. For smooth problems spectral element methods have proved very successful, since they can accommodate fairly complicated geometries while retaining a rapid rate of convergence. Geometric singularities, however, often give rise to singular solutions. The accuracy of the spectral element methods is then degraded, and they offer no apparent advantage over low-order finite element methods. In many cases, however, the singular structure of the solution is known, and its form may be exploited by the spectral element method. Among the various ways of doing somore » (through supplementary basis functions, eigenfunction expansions, and graded meshes), the method of auxiliary mapping proves to be particularly effective. For certain simple cases, the problem is transformed to a coordinate system in which the solution is analytic, and exponential convergence is recovered. Even when this is not possible, the singularity is usually much weaker after mapping, so that other treatments are more effective in the new coordinate system. In this paper, the authors study different ways of treating singularities, and in particular, the method of auxiliary mapping coupled with the use of supplementary basis functions. Error estimates are presented explaining why the combined approach is more effective, and these estimates are confirmed through a number of numerical experiments for the Laplace. Poisson, and Helmholtz equations. 28 refs., 13 figs., 5 tabs.« less