High-order numerical methods for the Riesz space fractional advection-dispersion equations

03/31/2020
by   Libo Feng, et al.
0

In this paper, we propose high-order numerical methods for the Riesz space fractional advection-dispersion equations (RSFADE) on a finite domain. The RSFADE is obtained from the standard advection-dispersion equation by replacing the first-order and second-order space derivative with the Riesz fractional derivatives of order α∈(0,1) and β∈(1,2], respectively. Firstly, we utilize the weighted and shifted Grünwald difference operators to approximate the Riesz fractional derivative and present the finite difference method for the RSFADE. Specifically, we discuss the Crank-Nicolson scheme and solve it in matrix form. Secondly, we prove that the scheme is unconditionally stable and convergent with the accuracy of O(τ^2+h^2). Thirdly, we use the Richardson extrapolation method (REM) to improve the convergence order which can be O(τ^4+h^4). Finally, some numerical examples are given to show the effectiveness of the numerical method, and the results are excellent with the theoretical analysis.

READ FULL TEXT

Please sign up or login with your details

Forgot password? Click here to reset