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.

The Spectral Collocation Method for Efficiently Solving PDEs with Fractional Laplacian.

A novel and simple spectral method for nonlocal PDEs with the fractional Laplacian

An Efficient and Accurate Numerical Method for the Spectral Fractional Laplacian Equation

Spectral Collocation Method for Fractional Differential/Integral Equations with Generalized Fractional Operator

Fractional-order dependent Radial basis functions meshless methods for the integral fractional Laplacian

A Highly Efficient Numerical Method for the Time-Fractional Diffusion Equation on Unbounded Domains

Jacobi spectral collocation method of space-fractional Navier-Stokes equations

Radial basis function collocation method for decoupled fractional Laplacian wave equations

Numerical Approximation of Space-time Fractional Parabolic Equations

An efficient spectral method for the fractional Schrödinger equation on the real line

Efficient and Accurate Spectral Method Using Generalized Jacobi Functions for Solving Riesz Fractional Differential Equations

Finite Difference Scheme and Finite Volume Scheme for Fractional Laplacian Operator and Some Applications

A PDE approach to fractional diffusion: a space-fractional wave equation

Double fast algorithm for solving time-space fractional diffusion problems with spectral fractional Laplacian

A sparse spectral method for fractional differential equations in one-spatial dimension

Collocation methods for integral fractional Laplacian and fractional PDEs based on radial basis functions

A static memory sparse spectral method for time-fractional PDEs

Exponentially accurate spectral and spectral element methods for fractional ODEs

Spectral solver for Cauchy problems in polar coordinates using discrete Hankel transforms

On a collocation method for the time-fractional convection-diffusion equation with variable coefficients

Spectral collocation methods for fractional multipantograph delay differential equations*