Abstract:This paper presents and analyzes a parallelizable iterative procedure based on domain decomposition for primal-dual weak Galerkin (PDWG) finite element methods applied to the Poisson equation. The existence and uniqueness of the PDWG solution are established. Optimal order of error estimates are derived in both a discrete norm and the $L^2$ norm. The convergence analysis is conducted for domain decompositions into individual elements associated with the PDWG methods, which can be extended to larger subdomains without any difficulty.
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is to develop and analyze a parallel iterative algorithm based on domain decomposition for the primal - dual weak Galerkin (PDWG) finite - element method of solving the Poisson equation. Specifically, the author focuses on the following points:
1. **Establishing existence and uniqueness**: Proved the existence and uniqueness of the PDWG solution.
2. **Error estimation**: Derived the optimal - order error estimates in the discrete norm and \(L^2\) norm.
3. **Convergence analysis**: Conducted a convergence analysis on the domain decomposition based on a single element and pointed out that this method can be extended to larger sub - domains without difficulty.
4. **Parallelization**: Designed a parallel iterative algorithm, enabling the PDWG method to be efficiently implemented on large - scale parallel computers, with each sub - domain being assigned to a processor for calculation.
### Main contributions
- **Theoretical contribution**: Through strict mathematical derivations, proved the effectiveness and stability of the PDWG method in solving the Poisson equation.
- **Numerical contribution**: Proposed a parallel iterative algorithm, which improves computational efficiency and is especially suitable for large - scale computational problems.
- **Application prospects**: Provides a new and efficient numerical method for scientific computing and engineering applications, especially having significant advantages in dealing with large - scale partial differential equation problems.
### Technical details
- **Weak Galerkin method**: Utilizes weak derivatives and weak continuity to design numerical schemes, making the method flexible and accurate in solving various partial differential equations.
- **Domain decomposition**: Decomposes the computational domain into multiple sub - domains, and the problems on each sub - domain can be independently solved. Finally, the sub - domain solutions are combined through an iterative method.
- **Parallel computing**: By assigning sub - domains to different processors, realizes the parallelization of computational tasks, greatly improving computational efficiency.
### Formula examples
- The form of the Poisson equation:
\[
\Delta u = f \quad \text{in } \Omega, \quad u = g \quad \text{on } \Gamma
\]
- The weak form:
\[
(\Delta_w \sigma, \phi)_T = (\sigma_0, \Delta \phi)_T - \langle \sigma_b, \nabla \phi \cdot n \rangle_{\partial T} + \langle \sigma_n, \phi \rangle_{\partial T}
\]
- The discrete weak Laplacian operator:
\[
(\Delta_{w,r,T} \sigma, w)_T = (\sigma_0, \Delta w)_T - \langle \sigma_b, \nabla w \cdot n \rangle_{\partial T} + \langle \sigma_n, w \rangle_{\partial T}, \quad \forall w \in P_r(T)
\]
These technical details and formulas demonstrate the innovativeness and practicality of the paper in theoretical and practical applications.