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.

A New Spectral Method Using Nonstandard Singular Basis Functions for Time-Fractional Differential Equations

A Multi-Domain Spectral Method for Time-Fractional Differential Equations

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

Spectral approximation of a new class of stochastic fractional evolution equations

Generalized Jacobi Functions and Their Applications to Fractional Differential Equations.

A Novel Method for Linear Systems of Fractional Ordinary Differential Equations with Applications to Time-Fractional PDEs

A novel spectral method for the subdiffusion equation

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

Exponentially accurate spectral and spectral element methods for fractional ODEs

Numerical solution of fractional Fredholm integro-differential equations by spectral method with fractional basis functions

A static memory sparse spectral method for time-fractional PDEs

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

A pseudo-spectral scheme for variable order fractional stochastic Volterra integro-differential equations

Numerical Approximation of Space-time Fractional Parabolic Equations

Semi-implicit Galerkin–Legendre Spectral Schemes for Nonlinear Time-Space Fractional Diffusion–Reaction Equations with Smooth and Nonsmooth Solutions

Quasi-Linear (Stochastic) Partial Differential Equations with Time-Fractional Derivatives

Fourier-Gegenbauer Pseudospectral Method for Solving Time-Dependent One-Dimensional Fractional Partial Differential Equations with Variable Coefficients and Periodic Solutions

A search for a spectral technique to solve nonlinear fractional differential equations

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

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

Numerical approximations for the nonlinear time fractional reaction–diffusion equation