A Second-Order Accurate Numerical Approximation for Two-Sided Fractional Boundary Value Advection-Diffusion Problem

Subject Areas : Applied Mathematics

Elyas Shivanian ¹ , Hamid Reza Khodabandehlo ²

1 - Department of Mathematics, Imam Khomeini International University, Qazvin, 34149-16818, Iran
2 - Department of Mathematics, Imam Khomeini International University, Qazvin, 34149-16818, Iran

Received: 2016-01-24 Accepted : 2016-06-17 Published : 2019-05-01

Keywords: Numerical fractional PDE, Two-sided fractional partial differential equation, Shifted Gr"{u}nwald-Letnikov formula, Fractional diffusion, Crank-Nicolson method,

Abstract :

Fractional order partial differential equations are generalization of classical partialdifferential equations. Increasingly, these models are used in applications such as fluid flow, financeand others. In this paper we examine some practical numerical methods to solve a class of initial-boundary value fractional partial differential equations with variable coefficients on a finite domain. An approach based on the classical Crank-Nicolson method combined with spatial extrapolation is used to obtain temporally and spatially second-order accurate numerical estimates. Stability, consistency, and convergence of the method are examined. It is shown that the fractional Crank-Nicolson method based on the shifted Gr"{u}nwaldformula is unconditionally stable. Some numerical examples are presented and compared with the exact analytical solutionfor its order of convergence.Fractional order partial differential equations are generalization of classical partialdifferential equations. Increasingly, these models are used in applications such as fluid flow, financeand others. In this paper we examine some practical numerical methods to solve a class of initial-boundary value fractional partial differential equations with variable coefficients on a finite domain. An approach based on the classical Crank-Nicolson method combined with spatial extrapolation is used to obtain temporally and spatially second-order accurate numerical estimates. Stability, consistency, and convergence of the method are examined. It is shown that the fractional Crank-Nicolson method based on the shifted Gr"{u}nwaldformula is unconditionally stable. Some numerical examples are presented and compared with the exact analytical solutionfor its order of convergence.