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 of the derivatives beyond Taylor expansion

Stability of Multiquadric Quasi-Interpolation to Approximate High Order Derivatives

On approximate implicit Taylor methods for ordinary differential equations

Implicit Summation by Parts Operators for Finite Difference Approximations of First and Second Derivatives

Enhancing Interpolation and Approximation Error Estimates Using a Novel Taylor-like Formula

Approximation to the K-Th Derivatives by Multiquadric Quasi-Interpolation Method

Taylor-Series Expansion Based Numerical Methods: A Primer, Performance Benchmarking and New Approaches for Problems with Non-smooth Solutions

Stability and Error Analysis for a Second-Order Fast Approximation of the Local and Nonlocal Diffusion Equations on the Real Line

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

A set of compact finite-difference approximations to first and second derivatives, related to the extended Numerov method of Chawla on nonuniform grids

First Derivative Approximations and Applications

On obtaining convergence order of a fourth and sixth order method of Hueso et al. without using Taylor series expansion

A Local Fractional Taylor Expansion and Its Computation for Insufficiently Smooth Functions

-analysis of H3N3-2 -based difference method for fractional hyperbolic equations

Numerical Approximation of Space-time Fractional Parabolic Equations

An Efficient Numerical Method for Solving a Class of Nonlinear Fractional Differential Equations and Error Estimates

Convergence of summation-by-parts finite difference methods for the wave equation

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

New higher-order implict method for approximating solutions of the initial value problems

On full history recursive multilevel Picard approximations and numerical approximations for high-dimensional nonlinear parabolic partial differential equations and high-dimensional nonlinear backward stochastic differential equations

Summation by Parts Operators for Finite Difference Approximations of Second-Derivatives with Variable Coefficients