Abstract:In this paper, we develop a numerical scheme for the space-time fractional parabolic equation, i.e., an equation involving a fractional time derivative and a fractional spatial operator. Both the initial value problem and the non-homogeneous forcing problem (with zero initial data) are considered. The solution operator $E(t)$ for the initial value problem can be written as a Dunford-Taylor integral involving the Mittag-Leffler function $e_{\alpha,1}$ and the resolvent of the underlying (non-fractional) spatial operator over an appropriate integration path in the complex plane. Here $\alpha$ denotes the order of the fractional time derivative. The solution for the non-homogeneous problem can be written as a convolution involving an operator $W(t)$ and the forcing function $F(t)$. We develop and analyze semi-discrete methods based on finite element approximation to the underlying (non-fractional) spatial operator in terms of analogous Dunford-Taylor integrals applied to the discrete operator. The space error is of optimal order up to a logarithm of $1/h$. The fully discrete method for the initial value problem is developed from the semi-discrete approximation by applying an exponentially convergent sinc quadrature technique to approximate the Dunford-Taylor integral of the discrete operator and is free of any time stepping. To approximate the convolution appearing in the semi-discrete approximation to the non-homogeneous problem, we apply a pseudo midpoint quadrature. This involves the average of $W_h(s)$, (the semi-discrete approximation to $W(s)$) over the quadrature interval. This average can also be written as a Dunford-Taylor integral. We first analyze the error between this quadrature and the semi-discrete approximation. To develop a fully discrete method, we then introduce sinc quadrature approximations to the Dunford-Taylor integrals for computing the averages.

Approximation and application of the Riesz-Caputo fractional derivative of variable order with fixed memory

Caputo-Hadamard fractional derivatives of variable order

Numerical Algorithms for Approximation of Fractional Integrals and Derivatives Based on Quintic Spline Interpolation

Three Approximations of Numerical Solution's by Finite Element Method for Resolving Space-Time Partial Differential Equations Involving Fractional Derivative's Order

A fast method for variable-order Caputo fractional derivative with applications to time-fractional diffusion equations

Numerical simulation of fractional-order reaction–diffusion equations with the Riesz and Caputo derivatives

New recursive approximations for variable-order fractional operators with applications

Numerical Approximation of Space-time Fractional Parabolic Equations

On new approximations of Caputo–Prabhakar fractional derivative and their application to reaction–diffusion problems with variable coefficients

A simple and fast finite difference method for the integral fractional Laplacian of variable order

Numerical approximation for a nonlinear variable-order fractional differential equation via an integral equation method

Fractional integration and differentiation

On the variable-order fractional derivatives with respect to another function

A cubic spline interpolation based numerical method for fractional differential equations

Asymptotic expansions and approximations for the Caputo derivative

Mathematical analysis and efficient finite element approximation for variable‐order time‐fractional reaction–diffusion equation with nonsingular kernel

Numerical Approximations for Fractional Differential Equations

A Second Order Approximation for the Caputo Fractional Derivative

An Innovative Fixed-Pole Numerical Approximation for Fractional Order Systems.

Numerical approximations for the tempered fractional Laplacian: Error analysis and applications

Fast second-order evaluation for variable-order Caputo fractional derivative with applications to fractional sub-diffusion equations