Abstract:We propose an implementation of linear finite element method for nonlocal diffusion problem in 2D space. In the implementation, we reduce the integral from 4D to 2D which would simplify the computation significantly.
What problem does this paper attempt to address?
This paper attempts to solve the computational complexity problems encountered in the numerical implementation of non - local diffusion models in two - dimensional space using the linear finite element method. Specifically, the author proposes a method to simplify the four - dimensional integral into a two - dimensional integral, thereby significantly simplifying the computational process.
### Background and Problem Description of the Paper
The non - local diffusion model is an important class of mathematical models and is widely used in fields such as physics and materials science. For such models, finite - element discretization usually has many excellent properties, such as asymptotic preservation and a good convergence rate. However, directly applying the finite - element method to solve non - local diffusion problems is difficult, mainly because it is necessary to calculate high - dimensional integrals (usually 4D integrals), which are very time - consuming and difficult to handle in practical applications.
### Main Contributions of the Paper
1. **Simplifying the Integral Dimension**: The author proposes a method to simplify the originally required 4D integral into a 2D integral. This simplification greatly reduces the computational complexity, making the numerical solution of non - local diffusion problems more feasible.
2. **Specific Models and Equations**: In this paper, taking the Poisson equation with Neumann boundary conditions as an example, the non - local corresponding Poisson equation form is given:
\[
\frac{1}{\delta^2}\int_{\Omega}R_{\delta}(x,y)(u(x) - u(y))dy+\int_{\Omega}\bar{R}_{\delta}(x,y)u(y)dy=\int_{\Omega}\bar{R}_{\delta}(x,y)f(y)dy + 2\int_{\partial\Omega}\bar{R}_{\delta}(x,y)g(y)d\tau_y
\]
where the kernel functions \(R_{\delta}(x,y)\) and \(\bar{R}_{\delta}(x,y)\) are defined as follows:
\[
R_{\delta}(x,y)=C_{\delta}R\left(\frac{\|x - y\|^2}{4\delta^2}\right),\quad\bar{R}_{\delta}(x,y)=C_{\delta}\bar{R}\left(\frac{\|x - y\|^2}{4\delta^2}\right)
\]
where \(C_{\delta}=\alpha^2\delta^{- 2}\), and \(\bar{R}(r)=\int_{r}^{+\infty}R(s)ds\), \(\bar{\bar{R}}(r)=\int_{r}^{+\infty}\bar{R}(s)ds\).
3. **Discretization Method**: In order to perform discretization, the author decomposes the computational domain \(\Omega\) into triangular meshes and constructs linear basis functions \(\phi_{i,k}\). Through these basis functions, the author derives the discretized form of the non - local diffusion model and details how to further simplify the integral calculation through techniques such as integration by parts.
4. **Calculating the Stiffness Matrix**: The paper discusses in detail how to calculate the coefficients in the stiffness matrix, especially how to transform complex 4D integrals into more tractable 2D integrals. This includes the treatment of zero - order terms and source terms, as well as the calculation of boundary terms.
### Conclusion
Through the above methods, the author successfully simplifies the 4D integral of the non - local diffusion model into a 2D integral, thereby significantly reducing the computational complexity. This method is not only applicable to non - local diffusion models with Neumann boundary conditions, but can also be extended to models with Dirichlet boundary conditions and other higher - dimensional problems.
In summary, the main contribution of this paper is to provide an effective numerical method, making the solution of non - local diffusion problems more efficient and feasible.