TY - JOUR

T1 - A Finite Difference Method for Space Fractional Differential Equations with Variable Diffusivity Coefficient

AU - Mustapha, K. A.

AU - Furati, K. M.

AU - Knio, Omar

AU - Le Maître, O. P.

N1 - KAUST Repository Item: Exported on 2021-02-25
Acknowledged KAUST grant number(s): KAUST005

PY - 2020/5/29

Y1 - 2020/5/29

N2 - Anomalous diffusion is a phenomenon that cannot be modeled accurately by second-order diffusion equations, but is better described by fractional diffusion models. The nonlocal nature of the fractional diffusion operators makes substantially more difficult the mathematical analysis of these models and the establishment of suitable numerical schemes. This paper proposes and analyzes the first finite difference method for solving variable-coefficient one-dimensional (steady state) fractional differential equations (DEs) with two-sided fractional derivatives (FDs). The proposed scheme combines first-order forward and backward Euler methods for approximating the left-sided FD when the right-sided FD is approximated by two consecutive applications of the first-order backward Euler method. Our scheme reduces to the standard second-order central difference in the absence of FDs. The existence and uniqueness of the numerical solution are proved, and truncation errors of order h are demonstrated (h denotes the maximum space step size). The numerical tests illustrate the global O(h) accuracy, except for nonsmooth cases which, as expected, have deteriorated convergence rates.

AB - Anomalous diffusion is a phenomenon that cannot be modeled accurately by second-order diffusion equations, but is better described by fractional diffusion models. The nonlocal nature of the fractional diffusion operators makes substantially more difficult the mathematical analysis of these models and the establishment of suitable numerical schemes. This paper proposes and analyzes the first finite difference method for solving variable-coefficient one-dimensional (steady state) fractional differential equations (DEs) with two-sided fractional derivatives (FDs). The proposed scheme combines first-order forward and backward Euler methods for approximating the left-sided FD when the right-sided FD is approximated by two consecutive applications of the first-order backward Euler method. Our scheme reduces to the standard second-order central difference in the absence of FDs. The existence and uniqueness of the numerical solution are proved, and truncation errors of order h are demonstrated (h denotes the maximum space step size). The numerical tests illustrate the global O(h) accuracy, except for nonsmooth cases which, as expected, have deteriorated convergence rates.

UR - http://hdl.handle.net/10754/626452

UR - http://link.springer.com/10.1007/s42967-020-00066-6

U2 - 10.1007/s42967-020-00066-6

DO - 10.1007/s42967-020-00066-6

M3 - Article

VL - 2

SP - 671

EP - 688

JO - Communications on Applied Mathematics and Computation

JF - Communications on Applied Mathematics and Computation

SN - 2096-6385

IS - 4

ER -