Abstract:This paper addresses the non-homogeneous Neumann boundary condition on domains with curved boundaries. We consider the Raviart-Thomas element (RTk ) of degree $k \geq 1 $on triangular mesh. on a triangular mesh. A key feature of our boundary value correction method is the shift from the true boundary to a surrogate boundary. We present a high-order version of the method, achieving an $O(h^k+1/2)$ convergence in $L^2$-norm estimate for the velocity field and an $O(h^k )$ convergence in $H^1$-norm estimate for the pressure. Finally, numerical experiments validate our theoretical results.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to deal with the non - homogeneous Neumann boundary condition problem on the region with curved boundaries. Specifically, the paper focuses on how to improve the accuracy of the numerical solution by the boundary value correction method when using Raviart - Thomas elements (RTk) on triangular meshes.
### Problem Background
In science and engineering, many practical problems involve regions with curved boundaries. For a region Ω with a curved boundary, approximating the boundary Γ with a polygon will lead to geometric errors, which will affect the accuracy of high - order elements. To solve this problem, two main strategies are proposed in the existing literature:
1. **Isoparametric finite element method** and **isogeometric analysis**: Reduce geometric errors by improving geometric approximation without changing the variational form.
2. **Boundary value correction method**: Solve directly on the polygon - approximated domain and modify the variational form.
### Specific Contributions of the Paper
This paper proposes a new boundary value correction method, which is achieved by mapping the real boundary Γ to an alternative boundary Γh. The key of this method lies in using Taylor expansion to deal with the normal flux differences between Γ and Γh. In addition, the paper also conducts a strict theoretical analysis and proves that this method can achieve the following convergence rates:
- In the L2 norm, the convergence rate of the velocity field is O(h^{k + 1/2}).
- In the H1 norm, the convergence rate of the pressure field is O(h^k).
### Mathematical Formula Representation
To ensure the correctness and readability of the formulas, the following is the Markdown - format representation of the main formulas involved in the paper:
1. **Taylor Expansion** :
\[
v(M_h(x_h))=\sum_{j = 0}^m\frac{\delta_h(x_h)^j}{j!}\partial_\nu^j v(x_h)+R_m v(x_h),\quad\text{on }\Gamma_h,
\]
where \(\partial_\nu^j\) represents the \(j\)-th partial derivative in the direction of the unit vector \(\nu_h\), and \(R_m v(x_h)\) is the remainder term and satisfies \(|R_m v(x_h)| = o(\delta^m)\).
2. **Discrete Space Definition** :
\[
V_h=\{v_h\in H(\text{div},\Omega_h):v_h|_K\in RT_k(K),K\in T_h\},
\]
\[
Q_h=\{q_h\in L^2(\Omega_h):q_h|_K\in P_k(K),K\in T_h\}.
\]
3. **Weak - form Equation** :
\[
a_h(u_h,v_h)+b_{h1}(v_h,p_h)=\sum_{e\in E_b^h}\langle h^{-1}_K e_{gN},T_k v_h\cdot e_n\rangle_e,\quad\forall v_h\in V_h,
\]
\[
b_{h0}(u_h,q_h)=-(f_E,q_h)_{\Omega_h},\quad\forall q_h\in Q_0^h.
\]
### Conclusion
By introducing the boundary value correction method, the paper successfully improves the accuracy of the numerical solution of the non - homogeneous Neumann boundary condition problem on the region with curved boundaries. Numerical experiments verify the effectiveness of the theoretical results and show the potential of this method in practical applications.