Abstract:In this paper we present a framework for the construction and implementation of general virtual element spaces based on projections built from constrained least squares problems. Building on the triples used for finite element spaces, we introduce the concept of a VEM tuple which encodes the necessary building blocks to construct these projections. Using this approach, a wide range of virtual element spaces can be defined. We discuss $H^k$-conforming spaces for $k=1,2$ as well as divergence and curl free spaces. This general framework has the advantage of being easily integrated into any existing finite element package and we demonstrate this within the open source software package DUNE.
What problem does this paper attempt to address?
### What problems does this paper attempt to solve?
This paper aims to solve a series of problems in the Virtual Element Method (VEM), especially the challenges encountered in constructing and implementing general virtual element spaces. Specifically, the goals of the paper include:
1. **Constructing a general framework**: The author proposes a general framework based on the projection construction of the constrained least - squares problem for implementing general virtual element spaces. This framework extends the traditional Finite Element Method (FEM), enabling the handling of more complex polygonal meshes and spaces with additional structures, such as arbitrary global regularity.
2. **Defining VEM tuples**: To describe how to construct these projection operators, the author introduces the concept of "VEM tuples". VEM tuples encode all the necessary components required to construct these projections. In this way, a wide range of virtual element spaces can be defined, including \(H^k\) - conforming spaces (\(k = 1,2\)), divergence - free and curl - free spaces, etc.
3. **Easy integration with existing software packages**: An important advantage of this framework is that it can be easily integrated into existing finite - element software packages. The author shows how this can be achieved in the open - source software package Dune, and users can switch seamlessly between finite - element and virtual - element discretizations.
4. **Solving nonlinear and high - order problems**: Traditional methods have limitations when dealing with nonlinear problems or partial differential equations with variable coefficients. This framework overcomes these limitations by introducing the constrained least - squares problem to calculate certain \(L^2\) projections. Moreover, this framework can be extended to fourth - order problems and three - dimensional problems.
5. **Improving flexibility and applicability**: By using VEM tuples, spaces with different properties can be flexibly constructed, such as spaces suitable for acoustic vibration problems, Stokes flow problems, etc.
### Formula summary
- **Projection operators**:
- The value projection \(\Pi_F^0\) and the gradient projection \(\Pi_E^1\) are defined by solving the constrained least - squares problem.
- The gradient projection is defined as follows:
\[
\int_E\Pi_E^1v_h\cdot q\,dx = -\int_E\Pi_E^0v_h\nabla\cdot q\,dx+\sum_{s\in S_h(E)}\int_s\Pi_s^0v_h(n\cdot q)\,ds
\]
where \(q\in B_E^1\) and \(n\) is the outer normal vector of the edge \(s\).
- **Constrained least - squares problem for value projection**:
- For the discrete function \(v_h\), the value projection \(\Pi_F^0\) is defined as the solution of the following constrained least - squares problem:
\[
\min\sum_{\lambda\in\Lambda_F}(\lambda(\Pi_F^0v_h - v_h))^2,\quad\text{subject to}\quad C(\Pi_F^0v_h - v_h)=0\quad\forall C\in C_F^0
\]
Through these formulas and methods, the author provides a general and flexible framework that can effectively construct and implement various virtual element spaces, thus solving many of the limitations in traditional methods.