Finite Elements with weighted bases for the fractional Laplacian

Abstract

This work presents a numerical study of the Dirichlet problem for the fractional Laplacian (-)s with s∈(0,1) using Finite Element methods with non-standard bases. Classical approaches based on piece-wise linear basis yield h 1 2 convergence rates in the Sobolev-Slobodeckij norm Hs due to the limited boundary regularity of the solution u(x), which behaves like dist(x,Rd )s, where h is the diameter of the mesh elements. To overcome this limitation, we propose a novel Finite Element basis of the form δs × (piece-wise linear functions), where δ is any suitably smooth approximation of dist(x,Rd ). This exploits the improved regularity of u/δs, achieving higher convergence rates. Under standard smoothness assumptions the method attains an order h2-s on quasi-uniform meshes, improving the rates with the piece-wise linear basis. We provide a rigorous theoretical error analysis with explicit rates and validate it through numerical experiments.

0

Turn this paper into a full lesson

ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…