Abstract:In Gander et al. [2004] we proposed a new non-conforming domain decomposition paradigm, the New Interface Cement Equilibrated Mortar (NICEM) method, based on Schwarz type methods that allows for the use of Robin interface conditions on non-conforming grids. The error analysis was done for P1 finite elements, in 2D and 3D. In this paper, we provide new numerical analysis results that allow to extend this error analysis in 2D for piecewise polynomials of higher order and also prove the convergence of the iterative algorithm in all these cases.
What problem does this paper attempt to address?
### Problems the paper attempts to solve
This paper aims to extend and improve a new non - conforming grid decomposition method (New Interface Cement Equilibrated Mortar, NICEM method), especially to conduct numerical analysis on high - order piecewise polynomials in two - dimensional space and prove the convergence of the iterative algorithm. Specifically:
1. **Extension of the NICEM method**:
- The NICEM method is a new non - conforming grid decomposition method based on the Schwarz - type method, which allows the use of Robin interface conditions on non - conforming grids.
- Previous work has carried out error analysis on P1 finite elements, and this paper extends it to higher - order piecewise polynomials (Pp finite elements, \( p \geq 1 \)).
2. **Results of numerical analysis**:
- New numerical analysis results are provided to deal with higher - order piecewise polynomials in two - dimensional space.
- The convergence of the iterative algorithm in these cases is proved.
3. **Error estimate analysis**:
- The previous error estimate analysis is extended to be applicable to high - order piecewise polynomials in two - dimensional space.
4. **Theoretical and practical verification**:
- After theoretical analysis, the correctness of the theoretical estimates is verified through numerical simulation.
### Representation of mathematical formulas
- **Problem definition**:
Consider the following partial differential equation problem:
\[
(I - \Delta)u = f \quad \text{in } \Omega
\]
\[
u = 0 \quad \text{on } \partial\Omega
\]
where \( f \in L^2(\Omega) \).
- **Variational form**:
The variational form of the problem is: find \( u \in H_0^1(\Omega) \) such that
\[
\int_\Omega (\nabla u \cdot \nabla v + uv) \, dx = \int_\Omega fv \, dx, \quad \forall v \in H_0^1(\Omega)
\]
- **Interface condition**:
The Robin interface condition is:
\[
\frac{\partial u_k}{\partial n_k} + \alpha u_k = -\frac{\partial u_\ell}{\partial n_\ell} + \alpha u_\ell \quad \text{on } \Gamma_{k,\ell}
\]
- **Discretization space**:
Define the discrete space \( V_h \) as:
\[
V_h = \left\{ (u_h, p_h) \in \left( \prod_{k = 1}^K X_k^h \right) \times \left( \prod_{k = 1}^K \tilde{W}_k^h \right), \int_{\Gamma_{k,\ell}} (p_{h,k} + \alpha u_{h,k}) \psi_{h,k,\ell} = \int_{\Gamma_{k,\ell}} (-p_{h,\ell} + \alpha u_{h,\ell}) \psi_{h,k,\ell}, \forall \psi_{h,k,\ell} \in \tilde{W}_{k,\ell}^h \right\}
\]
- **Discrete problem**:
Find \( (u_h, p_h) \in V_h \) such that
\[
\sum_{k = 1}^K \int_{\Omega_k} (\nabla u_{h,k} \cdot \nabla v_{h,k} + u_{h,k} v_{h,k})