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 finite difference scheme for the approximation of the third initial boundary value parabolic problem

Compact finite-difference scheme for some Sobolev type equations with Dirichlet boundary conditions

The Unconditional Stable and Convergent Difference Methods with Intrinsic Parallelism for Quasilinear Parabolic Systems

Finite difference method for boundary value problem for nonlinear elliptic equation with nonlocal conditions

A finite difference scheme for two-dimensional singularly perturbed convection-diffusion problem with discontinuous source term

An iterative numerical method for an inverse source problem for a multidimensional nonlinear parabolic equation

Construction Of Difference Schemes For Nonlinear Singular Perturbed Equations By Approximation Of Coefficients

On an Iterative Method of Solving Direct and Inverse Problems for Parabolic Equations

Numerical Approximation of Space-time Fractional Parabolic Equations

A nine-point finite difference scheme for one-dimensional wave equation

Finite Difference Method for a Semilinear Parabolic Equations with Variable Coefficients

Advanced numerical scheme and its convergence analysis for a class of two-point singular boundary value problems

An efficient fourth-order convergent scheme based on half-step spline function for two-point mixed boundary value problems

Three-Step Approximation Methods from Continuous and Discrete Perspective for Quasi-Variational Inequalities

Penalty method for the Navier-Stokes-Fourier system with Dirichlet boundary conditions: convergence and error estimates

Linear regularized finite difference scheme for the quasilinear subdiffusion equation

Finite Element Approximation of Hamilton-Jacobi-Bellman equations with nonlinear mixed boundary conditions

High accuracy two-level compact implicit method in exponential form for 2D fourth order quasi-linear parabolic equations

Convergence of a numerical scheme for a nonlinear oblique derivative boundary value problem

A convergent finite difference-quadrature scheme for the porous medium equation with nonlocal pressure

A uniformly convergent difference method for singularly perturbed parabolic partial differential equations with large delay and integral boundary condition