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.

Estimation of a heat source in a parabolic equation with nonlocal Wentzell boundary condition using a spectral technique

Inverse time-dependent source problem for the heat equation with a nonlocal Wentzell-Neumann boundary condition

Identifying a Time-Dependent Heat Source with Nonlocal Boundary and Overdetermination Conditions by the Variational Iteration Method

Method of virtual sources using on-surface radiation conditions for the Helmholtz equation

A Generalized Quasi-Boundary Value Method for Recovering a Source in a Fractional Diffusion-Wave Equation

Determination of space–time-dependent heat source in a parabolic inverse problem via the Ritz–Galerkin technique

Combined finite difference and spectral methods for the numerical solution of hyperbolic equation with an integral condition

Numerical and mathematical analysis of nonlocal singular Emden–Fowler type BVPs by improved Taylor-wavelet method

Application of two-sided approximations method to solution of first boundary value problem for one-dimensional nonlinear heat conductivity equation

Convergence of numerical schemes for the solution of partial integro-differential equations used in heat transfer

Determination of an unknown heat source term from boundary data

Gaussian beam approach for the boundary value problem of high frequency Helmholtz equation

A Mollification Regularization Method for Identifying the Time-Dependent Heat Source Problem

Approximate Solution of the Nonlinear Heat Conduction Equation in a Semi-Infinite Domain

Identifying an Unknown Source Term in A Heat Equation

Error Estimation of a Class of Stable Spectral Approximation to the Cahn-Hilliard Equation

Solution to a Two-Dimensional Nonlinear Parabolic Heat Equation Subject to a Boundary Condition Specified on a Moving Manifold

An energy-based summation-by-parts finite difference method for the wave equation in second order form

Application of the time discontinuous Galerkin finite element method to heat wave simulation

Numerical Approximation of Space-time Fractional Parabolic Equations

A Quasi‐boundary Value Regularization Method for Determining the Heat Source