A scalable domain decomposition method for FEM discretizations of nonlocal equations of integrable and fractional type

Manuel Klar,Giacomo Capodaglio,Marta D'Elia,Christian Glusa,Max Gunzburger,Christian Vollmann
2023-06-01
Abstract:Nonlocal models allow for the description of phenomena which cannot be captured by classical partial differential equations. The availability of efficient solvers is one of the main concerns for the use of nonlocal models in real world engineering applications. We present a domain decomposition solver that is inspired by substructuring methods for classical local equations. In numerical experiments involving finite element discretizations of scalar and vectorial nonlocal equations of integrable and fractional type, we observe improvements in solution time of up to 14.6x compared to commonly used solver strategies.
Numerical Analysis
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: in non - local models, how to efficiently solve non - local equations. Specifically, the paper proposes a domain decomposition solver based on the sub - structure method for non - local integral - type and fractional - order equations with finite - element discretization. This method aims to improve the efficiency of solving non - local models, especially in practical engineering applications. ### Paper Background Non - local models can describe phenomena that cannot be captured by classical partial differential equations, such as long - range interactions and jump discontinuities. However, an efficient solver is one of the key challenges in using non - local models for practical engineering applications. ### Main Contributions 1. **Domain Decomposition Method**: The paper proposes a domain decomposition solver inspired by the sub - structure method of classical local equations. 2. **Numerical Experiments**: Through numerical experiments, it shows a significant improvement in the solution time of this method, which can be up to 14.6 times faster than the commonly used solution strategies. 3. **Theoretical Framework**: It provides detailed steps and algorithm descriptions, which are applicable to non - local diffusion models and can be extended to other non - local models. ### Key Technologies - **Non - local Diffusion Model**: Considered the non - local diffusion setting, where the interaction distance between points does not exceed a given constant \( \delta \). - **Domain Decomposition Algorithm**: Proposed a non - local domain decomposition algorithm to solve the overlapping sub - domain problem in non - local models. - **Finite Element Method**: Combined with the finite element method, defined the discrete weak form of non - local models. - **FETI Method**: Generalized the classical FETI method to non - local models, providing a scalable multi - domain solution method. ### Challenges - **Non - locality**: The finite - element stiffness matrix in non - local models is less sparse than that in local PDE models, and the condition number depends not only on the mesh size but also on the horizon \( \delta \). - **Overlapping Sub - domains**: Sub - domains in non - local models not only overlap on the common boundary but also exchange information between internal points of adjacent sub - domains. ### Applications The application fields of non - local models are extensive, including multi - scale and multi - physical systems, underground fluid flow, corrosion and fracture propagation in solid mechanics, turbulence, phase change, image processing and stochastic processes, etc. ### Conclusion The method proposed in the paper performs well in solving non - local models, especially in improving the solution efficiency. Through detailed algorithm descriptions and numerical experiment verifications, it shows the effectiveness and practicality of this method.