Lourenço Beirão da Veiga,Franco Dassi,Daniele A. Di Pietro,Jérôme Droniou
Abstract:This paper contains two major contributions. First we derive, following the discrete de Rham (DDR) and Virtual Element (VEM) paradigms, pressure-robust methods for the Stokes equations that support arbitrary orders and polyhedral meshes. Unlike other methods presented in the literature, pressure-robustness is achieved here without resorting to an $\boldsymbol{H}({\rm div})$-conforming construction on a submesh, but rather projecting the volumetric force onto the discrete $\boldsymbol{H}({\bf curl})$ space. The cancellation of the pressure error contribution stems from key commutation properties of the underlying DDR and VEM complexes. The pressure-robust error estimates in $h^{k+1}$ (with $h$ denoting the meshsize and $k\ge 0$ the polynomial degree of the DDR or VEM complex) are proven theoretically and supported by a panel of three-dimensional numerical tests. The second major contribution of the paper is an in-depth study of the relations between the DDR and VEM approaches. We show, in particular, that a complex developed following one paradigm admits a reformulation in the other, and that couples of related DDR and VEM complexes satisfy commuting diagram properties with the degrees of freedom maps.
What problem does this paper attempt to address?
### What problems does this paper attempt to solve?
The main purpose of this paper is to develop and analyze pressure - robust discrete de Rham (DDR) and virtual element methods (VEM) applicable to any order and polyhedral meshes for solving the Stokes problem. Specifically, the article addresses the following two main issues:
1. **Develop new pressure - robust numerical methods**:
- By projecting the body force onto the discrete \( \mathbf{H}(\text{curl}) \) space, the author achieves pressure robustness that does not depend on the sub - grid \( \mathbf{H}(\text{div}) \)-conforming construction.
- This method takes advantage of the key commutativity properties of the DDR and VEM complexes, thereby eliminating the influence of the pressure error term.
2. **Study the relationship between DDR and VEM methods**:
- The article delves into the connection between DDR and VEM methods, showing that a complex developed based on one paradigm can be reformulated in another paradigm, and the relevant DDR and VEM complexes satisfy the commutative diagram property under the degree - of - freedom mapping.
### Markdown representation of formulas
To ensure the correctness and readability of formulas, the following are several key formulas represented in Markdown format:
- **Weak form of the Stokes problem**:
\[
\begin{cases}
\int_{\Omega} \text{curl} \, \mathbf{u} \cdot \text{curl} \, \mathbf{v} + \int_{\Omega} \nabla p \cdot \mathbf{v} = \int_{\Omega} \mathbf{f} \cdot \mathbf{v} & \forall \mathbf{v} \in \mathbf{H}(\text{curl}; \Omega), \\
\int_{\Omega} \nabla q \cdot \mathbf{u} = 0 & \forall q \in H^1(\Omega) \cap L^2_0(\Omega).
\end{cases}
\]
- **Definition of discrete spaces**:
\[
\mathbf{X}_h^{\text{curl}, k} = \left\{ \mathbf{v}_h = \left( (\mathbf{v}_{R,T}, \mathbf{v}_{c,T})_{T \in \mathcal{T}_h}, (\mathbf{v}_{R,F}, \mathbf{v}_{c,F})_{F \in \mathcal{F}_h}, (v_E)_{E \in \mathcal{E}_h} \right) :
\begin{aligned}
&\mathbf{v}_{R,T} \in R_{k - 1}(T), \mathbf{v}_{c,T} \in R_{c,k}(T) \text{ for all } T \in \mathcal{T}_h, \\
&\mathbf{v}_{R,F} \in R_{k - 1}(F), \mathbf{v}_{c,F} \in R_{c,k}(F) \text{ for all } F \in \mathcal{F}_h, \\
&v_E \in P_k(E) \text{ for all } E \in \mathcal{E}_h
\end{aligned}
\right\}.
\]
- **Discrete gradient operator**:
\[
\mathbf{G}_h^k : X_h^{\text{grad}, k} \to \mathbf{X}_h^{\text{curl}, k},
\]