Stabilization-free virtual element methods based on finite element interpolation
Abstract
In this paper, we introduce a new framework for designing stabilization-free virtual element methods (VEMs) based on an finite element interpolation-based strategy, where we can simultaneously eliminate the stabilization terms in the discretizations of diffusion and reaction terms. The core idea is to construct a computable, polynomial-preserving, and norm-equivalent interpolation operator from the virtual element space to a (local) finite element space. Leveraging the properties of this operator, we design two types of stabilization-free schemes. The first scheme requires the interpolation to preserve the polynomial consistency related to the bilinear forms, thereby maintaining both consistency and stability as in the standard VEM. The second scheme relaxes this consistency requirement. While it may not satisfy the standard polynomial consistency, the second scheme retains optimal convergence with simpler construction, fewer degrees of freedom and, more importantly, applicable to more complex problems such as those involving nonlinearities or variable coefficients. We construct concrete interpolation operators for both conforming and nonconforming virtual elements in two and three dimensions. These operators are then employed to realize stabilization-free schemes for conforming and nonconforming VEMs. Numerical experiments confirm the optimal convergence rates of the proposed methods. The presented framework can be extended to design stabilization-free schemes for other polytopal discretization methods, such as the hybrid high-order method and the weak Galerkin method.
Related
- C1 Virtual Element Methods for the Optimal Control of Oseen Equations with Stream-Function Formulation
- The pressure-robust virtual element methods with reconstruction for the Stokes equation
- A two-level preconditioner for discontinuous virtual element methods with irregular subdomains
- A quasi-interpolation method for finite element spaces
- Space-time finite element interpolation of nonsmooth functions satisfying boundary conditions
BibTeX
@misc{zhao2026stabilizationfree,
title = {{Stabilization-free virtual element methods based on finite element interpolation}},
author = {Jikun Zhao and Wenhao Zhu and Bei Zhang and Shipeng Mao},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2606.01614},
eprint = {2606.01614},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2606.01614},
}
RIS
TY - UNPB
TI - Stabilization-free virtual element methods based on finite element interpolation
AU - Jikun Zhao
AU - Wenhao Zhu
AU - Bei Zhang
AU - Shipeng Mao
PY - 2026
DA - 2026/10/04
PB - arXiv
DO - 10.48550/arxiv.2606.01614
UR - https://arxiv.org/abs/2606.01614
AB - In this paper, we introduce a new framework for designing stabilization-free virtual element methods (VEMs) based on an finite element interpolation-based strategy, where we can simultaneously eliminate the stabilization terms in the discretizations of diffusion and reaction terms. The core idea is to construct a computable, polynomial-preserving, and norm-equivalent interpolation operator from the virtual element space to a (local) finite element space. Leveraging the properties of this operator, we design two types of stabilization-free schemes. The first scheme requires the interpolation to preserve the polynomial consistency related to the bilinear forms, thereby maintaining both consistency and stability as in the standard VEM. The second scheme relaxes this consistency requirement. While it may not satisfy the standard polynomial consistency, the second scheme retains optimal convergence with simpler construction, fewer degrees of freedom and, more importantly, applicable to more complex problems such as those involving nonlinearities or variable coefficients. We construct concrete interpolation operators for both conforming and nonconforming virtual elements in two and three dimensions. These operators are then employed to realize stabilization-free schemes for conforming and nonconforming VEMs. Numerical experiments confirm the optimal convergence rates of the proposed methods. The presented framework can be extended to design stabilization-free schemes for other polytopal discretization methods, such as the hybrid high-order method and the weak Galerkin method.
ER -