Abstract:We study a nonconforming virtual element method (VEM) for advection-diffusion-reaction problems with continuous interior penalty (CIP) stabilization. The design of the method is based on a standard variational formulation of the problem (no skew-symmetrization), and boundary conditions are imposed with a Nitsche technique. We use the enhanced version of VEM, with a ``DoFi-DoFi'' stabilization in the diffusion and reaction terms. We prove stability of the proposed method and derive $h$-version error estimates.
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is that for advection - diffusion - reaction problems in the case of convection - dominated, the standard Galerkin methods (such as the Finite Element Method, FEM) will lead to non - physical oscillation phenomena in the numerical solution. Specifically, when the diffusion term becomes very small relative to the convection and reaction terms, traditional discretization methods may produce unstable solutions unless an extremely fine grid is used.
To solve this problem, the authors propose a nonconforming virtual element method (VEM) and combine it with the Continuous Interior Penalty (CIP) stabilization technique to deal with these problems. This method aims to provide a quasi - robust method, that is, it can still maintain good numerical stability when the diffusion parameter is small, and can obtain error estimates independent of the diffusion parameter.
### Main contributions
1. **Nonconforming VEM**: Compared with previous work, this paper extends the application of nonconforming VEM, making it applicable not only to the H1 - conforming case but also to a wider range of nonconforming cases.
2. **CIP stabilization**: The CIP stabilization term is introduced to ensure the stability of the method in the case of convection - dominated.
3. **Nitsche method**: The symmetric version of the Nitsche method is used to impose Dirichlet boundary conditions, ensuring that the optimal L2 error estimate can still be maintained when the convection tends to zero.
4. **Theoretical analysis**: The stability of the proposed method is proved, and the h - version error estimate is derived.
### Mathematical model
Consider the steady - state advection - diffusion - reaction equation in the following form:
\[
\begin{cases}
- \varepsilon \Delta u+\beta \cdot \nabla u+\sigma u = f & \text{in } \Omega, \\
u = g & \text{on } \Gamma,
\end{cases}
\]
where $\Omega \subset \mathbb{R}^d$ is a polygonal domain, $\Gamma$ is its boundary, $\varepsilon$ is the diffusion coefficient, $\beta$ is the convection coefficient, $\sigma$ is the reaction coefficient, $f$ is the source term, and $g$ is the boundary condition.
### Discretization method
To discretize the above problem, the authors adopt a standard discretization method based on the variational principle and do not skew - symmetrize the discrete convection form. In addition, they use the Nitsche method to impose Dirichlet boundary conditions and introduce the CIP stabilization term to enhance the stability of the method.
Through these improvements, this method can effectively suppress the non - physical oscillation of the numerical solution in the case of convection - dominated, while maintaining high computational efficiency and accuracy.
### Summary
The main purpose of this paper is to develop, analyze, and numerically test a nonconforming virtual element method for advection - diffusion - reaction problems, combined with the CIP stabilization technique, in order to overcome the numerical instability and non - physical oscillation problems that may occur in traditional methods in the case of convection - dominated.