Solving time-fractional differential equation via rational approximation

by   Ustim Khristenko, et al.

Fractional differential equations (FDEs) describe subdiffusion behavior of dynamical systems. Its non-local structure requires to take into account the whole evolution history during the time integration, which then possibly causes additional memory use to store the history, growing in time. An alternative to a quadrature of the history integral is to approximate the fractional kernel with the sum of exponentials, which is equivalent to consider the FDE solution as a sum of solutions to a system of ODEs. One possibility to construct this system is to approximate the Laplace spectrum of the fractional kernel with a rational function. In this paper, we use the adaptive Antoulas–Anderson (AAA) algorithm for the rational approximation of the kernel spectrum which yields only a small number of real valued poles. We propose a numerical scheme based on this idea and study its stability and convergence properties. Moreover, we apply the algorithm to a time-fractional Cahn-Hilliard problem.


Approximate solutions of one dimensional systems with fractional derivative

The fractional calculus is useful to model non-local phenomena. We const...

Splitting Schemes for Non-Stationary Problems with a Rational Approximation for Fractional Powers of the Operator

Problems of the numerical solution of the Cauchy problem for a first-ord...

Highly Accurate Global Padé Approximations of Generalized Mittag-Leffler Function and its Inverse

The two-parametric Mittag-Leffler function (MLF), E_α,β, is fundamental ...

Rational approximation preconditioners for multiphysics problems

We consider a class of mathematical models describing multiphysics pheno...

Error Estimates for a Linearized Fractional Crank-Nicolson FEM for Kirchhoff type Quasilinear Subdiffusion Equation with Memory

In this paper, we develop a linearized fractional Crank-Nicolson-Galerki...

Equivalence between a time-fractional and an integer-order gradient flow: The memory effect reflected in the energy

Time-fractional partial differential equations are nonlocal in time and ...

A fast and oblivious matrix compression algorithm for Volterra integral operators

The numerical solution of dynamical systems with memory requires the eff...