Abstract:Piecewise divergence-free nonconforming virtual elements are designed for Stokes problem in any dimensions. After introducing a local energy projector based on the Stokes problem and the stabilization, a divergence-free nonconforming virtual element method is proposed for Stokes problem. A detailed and rigorous error analysis is presented for the discrete method. An important property in the analysis is that the local energy projector commutes with the divergence operator. With the help of a divergence-free interpolation operator onto a generalized Raviart-Thomas element space, a pressure-robust nonconforming virtual element method is developed by simply modifying the right hand side of the previous discretization. A reduced virtual element method is also discussed. Numerical results are provided to verify the theoretical convergence.
What problem does this paper attempt to address?
This paper aims to solve the problems of numerical simulation accuracy and stability in the Stokes problem. Specifically, the paper proposes a new piecewise divergence - free non - conforming virtual element method for solving the Stokes equations in arbitrary dimensions. This method ensures the strictly divergence - free property of the discrete velocity field by designing local energy projection operators and stabilization terms, which is crucial for the simulation of incompressible fluids. In addition, the paper also develops a pressure - robust non - conforming virtual element method by modifying the right - hand side of the discretization method. These methods not only improve the accuracy of the numerical solution but also ensure the stability of long - time simulations.
### Main contributions of the paper:
1. **Design of Divergence - Free Non - conforming Virtual Elements**: The paper designs a piecewise divergence - free non - conforming virtual element for the numerical solution of the Stokes problem. These virtual elements are divergence - free on each element, thus ensuring mass conservation.
2. **Local Energy Projection Operator**: A local energy projection operator $\Pi_K^k$ based on the Stokes problem is introduced. This operator commutes with the divergence operator, that is, $\text{div}(\Pi_K^k w)=\text{div}(w)$. This makes the projected velocity field still divergence - free on each element.
3. **Stabilization Term**: A stabilization term $S_K(w, v)$ is defined. This stabilization term involves all degrees of freedom of the velocity, except those corresponding to $G_{k - 2}(K):=\nabla P_{k - 1}(K)$. Through this stabilization term, the stability of the discrete method is ensured.
4. **Error Analysis**: A detailed error analysis of the proposed discrete method is carried out, proving the well - posedness and optimal error estimates of the discrete method. Specifically, it is proved that the discrete method satisfies the discrete inf - sup condition and the optimal error estimate is obtained.
5. **Pressure Robustness**: By modifying the right - hand side of the discretization method, a pressure - robust non - conforming virtual element method is developed. This method is more robust when dealing with the external force field $f$ and the accuracy of the velocity field will not be affected by the change of the pressure field.
### Summary of Mathematical Formulas:
- **Stokes Equations**:
\[
\begin{cases}
-\text{div}(\nu\epsilon(u))-\nabla p = f & \text{in }\Omega,\\
\text{div}u = 0 & \text{in }\Omega,\\
u = 0 & \text{on }\partial\Omega,
\end{cases}
\]
where $u$ is the velocity field, $p$ is the pressure, $\epsilon(u):=(\nabla u+(\nabla u)^T)/2$ is the symmetric gradient of the velocity field, $f$ is the external force field, and $\nu$ is the viscosity.
- **Local Energy Projection Operator**:
\[
\begin{cases}
(\epsilon(\Pi_K^k w),\epsilon(v))_K+(\text{div}v, P_K w)_K = (\epsilon(w),\epsilon(v))_K & \forall v\in P_k(K;\mathbb{R}^d),\\
\text{div}(\Pi_K^k w)=Q_K^{k - 1}(\text{div}w) & \forall q\in P_{k - 1}(K),\\
Q_K^0(\text{curl}\Pi_K^k w)=Q_K^0(\text{curl}w),\\
Q_K^0(\Pi_K^k w)=Q_K^0 w.
\end{cases}
\]
- **Stabilization Term**:
\[
S_K(w, v):=h_K^{-2}(Q_K^{G_{k - 2}} w,