Numerical approximation of the integral fractional Laplacian

We propose a new nonconforming finite element algorithm to approximate the solution to the elliptic problem involving the fractional Laplacian. We first derive an integral representation of the bilinear form corresponding to the variational problem. The numerical approximation of the action of the c...

Ausführliche Beschreibung

Gespeichert in:
Bibliographische Detailangaben
Veröffentlicht in:Numerische Mathematik Jg. 142; H. 2; S. 235 - 278
Hauptverfasser: Bonito, Andrea, Lei, Wenyu, Pasciak, Joseph E.
Format: Journal Article
Sprache:Englisch
Veröffentlicht: Berlin/Heidelberg Springer Berlin Heidelberg 01.06.2019
Springer Nature B.V
Schlagworte:
Algorithms
Approximation
Error analysis
Finite element method
Integrals
Mathematical and Computational Engineering
Mathematical and Computational Physics
Mathematical Methods in Physics
Mathematics
Mathematics and Statistics
Numerical Analysis
Numerical and Computational Physics
Parameters
Partial differential equations
Representations
Simulation
Stiffness matrix
Theoretical
65N15
65N12
65R20
35S15
65N30
ISSN:0029-599X, 0945-3245
Online-Zugang:Volltext
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
Beschreibung
Zusammenfassung:We propose a new nonconforming finite element algorithm to approximate the solution to the elliptic problem involving the fractional Laplacian. We first derive an integral representation of the bilinear form corresponding to the variational problem. The numerical approximation of the action of the corresponding stiffness matrix consists of three steps: (1) apply a sinc quadrature scheme to approximate the integral representation by a finite sum where each term involves the solution of an elliptic partial differential equation defined on the entire space, (2) truncate each elliptic problem to a bounded domain, (3) use the finite element method for the space approximation on each truncated domain. The consistency error analysis for the three steps is discussed together with the numerical implementation of the entire algorithm. The results of computations are given illustrating the error behavior in terms of the mesh size of the physical domain, the domain truncation parameter and the quadrature spacing parameter.
Bibliographie:ObjectType-Article-1
SourceType-Scholarly Journals-1
ObjectType-Feature-2
content type line 14
ISSN:0029-599X
0945-3245
DOI:10.1007/s00211-019-01025-x