Finite difference method for inhomogeneous fractional Dirichlet problem

01/27/2021
by   Jing Sun, et al.
0

We make the split of the integral fractional Laplacian as (-Δ)^s u=(-Δ)(-Δ)^s-1u, where s∈(0,1/2)∪(1/2,1). Based on this splitting, we respectively discretize the one- and two-dimensional integral fractional Laplacian with the inhomogeneous Dirichlet boundary condition and give the corresponding truncation errors with the help of the interpolation estimate. Moreover, the suitable corrections are proposed to guarantee the convergence in solving the inhomogeneous fractional Dirichlet problem and an 𝒪(h^1+α-2s) convergence rate is obtained when the solution u∈ C^1,α(Ω̅^δ_n), where n is the dimension of the space, α∈(max(0,2s-1),1], δ is a fixed positive constant, and h denotes mesh size. Finally, the performed numerical experiments confirm the theoretical results.

READ FULL TEXT

page 23

page 24

research
10/06/2021

Besov regularity for the Dirichlet integral fractional Laplacian in Lipschitz domains

We prove Besov regularity estimates for the solution of the Dirichlet pr...
research
09/27/2022

Two-level error estimation for the integral fractional Laplacian

For the singular integral definition of the fractional Laplacian, we con...
research
08/13/2019

An Auxiliary Space Preconditioner for Fractional Laplacian of Negative Order

Coupled multiphysics problems often give rise to interface conditions na...
research
09/22/2020

A unified meshfree pseudospectral method for solving both classical and fractional PDEs

In this paper, we propose a meshfree method based on the Gaussian radial...
research
07/26/2023

A grid-overlay finite difference method for the fractional Laplacian on arbitrary bounded domains

A grid-overlay finite difference method is proposed for the numerical ap...
research
01/26/2022

Schlömilch integrals and probability distributions on the simplex

The Schlömilch integral, a generalization of the Dirichlet integral on t...

Please sign up or login with your details

Forgot password? Click here to reset