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 sparse representation related with fractional heat equations

Sparse representations of approximation to identity via time–space fractional heat equations

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

A sparse approximation for fractional Fourier transform

A frame approach for equations involving the fractional Laplacian

The numerical solution of fractional integral equations via orthogonal polynomials in fractional powers

Sparse Representations of Solutions to a class of Random Boundary Value Problems

Sparse Series Solutions of Random Boundary and Initial Value Problems

A static memory sparse spectral method for time-fractional PDEs

Fourier pseudospectral methods for the spatial variable-order fractional wave equations

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

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

Discretization of Fractional Fully Nonlinear Equations by Powers of Discrete Laplacians

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

Numerical methods for forward fractional Feynman–Kac equation

Numerical Approximation of Space-time Fractional Parabolic Equations

A stabilized semi-implicit Fourier spectral method for nonlinear space-fractional reaction-diffusion equations

Numerical solutions to time-fractional stochastic partial differential equations

A weak solution for the fractional N-Laplacian flow