Abstract:This paper introduces a nonconforming virtual element method for general second-order elliptic problems with variable coefficients on domains with curved boundaries and curved internal interfaces. We prove arbitrary order optimal convergence in the energy and $L^2$ norms, confirmed by numerical experiments on a set of polygonal meshes. The accuracy of the numerical approximation provided by the method is shown to be comparable with the theoretical analysis.
What problem does this paper attempt to address?
This paper aims to solve the numerical solution problem of general second - order elliptic problems on regions with curved boundaries. Specifically, the author proposes a non - conforming virtual element method (NVEM) to deal with general second - order elliptic problems with curved boundaries or internal interfaces. This method allows for high - precision numerical approximation on domains with curved boundaries and can maintain the optimal convergence of any order.
### Main Problem Description
The problem discussed in the paper can be formalized as the following Dirichlet problem:
\[
\begin{cases}
- \nabla \cdot (a \nabla u)+\nabla \cdot (b u)+c u = f & \text{in } \Omega, \\
u = g & \text{on } \partial \Omega,
\end{cases}
\]
where:
- \(\Omega\) is a region with a curved boundary in a \(d\)-dimensional space (\(d = 2,3\)).
- \(a=(a_{i,j}(x))\in[L^{\infty}(\Omega)]^{d\times d}\) is a symmetric matrix - valued function, and it is assumed to be strongly elliptic, that is, there exists a constant \(a_0>0\) such that \(a_0\xi^{T}\xi\leqslant\xi^{T}a\xi,\forall\xi\in\mathbb{R}^{d}\).
- \(b = (b_i(x))\in[L^{\infty}(\Omega)]^{d}\) is a vector - valued function.
- \(c = c(x)\in L^{\infty}(\Omega)\) is a scalar function.
- \(f\in L^{2}(\Omega)\) and \(g\in H^{1 / 2}(\partial \Omega)\) are the source term and the boundary condition respectively.
### Method Overview
The method proposed in the paper mainly includes the following aspects:
1. **Non - conforming virtual element space**: Define a non - conforming virtual element space \(V_h(K)\), where \(K\) is a polygonal element with a curved boundary.
2. **Projection operators**: Introduce the \(L^{2}\) projection operator and the Ritz - Galerkin operator to handle integrals and derivatives on the curved boundary.
3. **Discrete bilinear forms**: Define the discrete bilinear forms \(a_h(u_h,v_h)\), \(b_h(u_h,v_h)\) and \(c_h(u_h,v_h)\), and globally combine them into \(B_h(u_h,v_h)\).
4. **Stability analysis**: Prove the stability and convergence of the proposed projection and interpolation operators.
5. **Numerical experiments**: Verify the effectiveness of the method and the accuracy of the theoretical analysis through numerical experiments.
### Main Contributions
- Propose a non - conforming virtual element method suitable for general second - order elliptic problems on domains with curved boundaries.
- Prove the optimal convergence of any order of this method in the energy norm and the \(L^{2}\) norm.
- Verify the accuracy and effectiveness of the method through numerical experiments.
In short, the main goal of this paper is to develop an efficient and accurate numerical method for solving second - order elliptic problems with complex geometries.