Abstract:We study the asymptotic error between the finite element solutions of nonlocal models with a bounded interaction neighborhood and the exact solution of the limiting local model. The limit corresponds to the case when the horizon parameter, the radius of the spherical nonlocal interaction neighborhood of the nonlocal model, and the mesh size simultaneously approach zero. Two important cases are discussed: one involving the original nonlocal models and the other for nonlocal models with polygonal approximations of the nonlocal interaction neighborhood. Results of numerical experiments are also reported to substantiate the theoretical studies.
What problem does this paper attempt to address?
### What problems does this paper attempt to solve?
This paper aims to study the error estimates in solving non - local models by the finite element method, especially when the interaction neighborhood of the non - local model is an exact or approximate spherical region. Specifically, the author focuses on the asymptotic error between the finite element solution of the non - local model and the exact solution of the corresponding local model.
#### Main research contents:
1. **Relationship between non - local and local models**:
- The paper explores how the solution of the non - local model approaches the solution of the local model when the horizon parameter \(\delta\) in the non - local model and the mesh size \(h\) both approach zero.
- The local model refers to the classical partial differential equation (PDE), while the non - local model involves integral operators which describe the interactions within a certain distance.
2. **Research on two different situations**:
- **Original non - local model**: Consider the non - local model with an exact spherical interaction neighborhood.
- **Non - local model with polygon - approximate interaction neighborhood**: Consider the case of using polygons to approximate the spherical interaction neighborhood.
3. **Numerical experiments**:
- The paper also reports the results of numerical experiments to verify the correctness of the theoretical analysis.
#### Specific mathematical expressions:
- The non - local operator \(L_\delta\) is defined as:
\[
L_\delta u_\delta(x)=2\int_{\mathbb{R}^d}(u_\delta(y) - u_\delta(x))\gamma_\delta(x, y)\,dy\quad\forall x\in\Omega,
\]
where \(\gamma_\delta(x, y)\) is a non - negative symmetric kernel function.
- The non - local Dirichlet volume - constrained diffusion problem is defined as:
\[
\begin{cases}
-L_\delta u_\delta(x)=f_\delta(x)&\text{on }\Omega,\\
u_\delta(x)=g_\delta(x)&\text{on }\Omega_c^\delta,
\end{cases}
\]
where \(\Omega_c^\delta = \{y\in\mathbb{R}^d\setminus\Omega:\text{dist}(y, \partial\Omega)<\delta\}\) represents the interaction domain.
#### Research motivation:
- When the horizon parameter \(\delta\) approaches zero, the non - local effect disappears and the solution of the non - local model should converge to the solution of the classical local PDE. This property is very important for practical applications (such as multi - scale modeling and benchmark testing).
- An important problem in numerical analysis is to study the extent to which various discrete approximation methods maintain this limit property at the discrete level.
Through the above research, the paper provides a theoretical basis for the numerical solution of non - local problems and guidance for parameter selection in practical applications.