A static memory sparse spectral method for time-fractional PDEs

Abstract

We introduce a method which provides accurate numerical solutions to fractional-in-time partial differential equations posed on $[0,T]\times \Omega$ with $\Omega \subset {\mathbb{R}}^d$ without the excessive memory requirements associated with the nonlocal fractional derivative operator. Our approach combines recent advances in the development and utilization of multivariate sparse spectral methods as well as fast methods for the computation of Gauss quadrature nodes with recursive non-classical methods for the Caputo fractional derivative of general fractional order $\alpha >0$. An attractive feature of the method is that it has minimal theoretical overhead when using it on any domain $\Omega$ on which an orthogonal polynomial basis is already available. We discuss the memory requirements of the method, present several numerical experiments demonstrating the method's performance in solving time-fractional PDEs on intervals, triangles and disks and derive error bounds which suggest sensible convergence strategies. As an important model problem for this approach we consider a type of wave equation with time-fractional dampening related to acoustic waves in viscoelastic media with applications in the physics of medical ultrasound and outline future research steps required to use such methods for the reverse problem of image reconstruction from sensor data.

arXiv comment27 pages, 13 figures

References in corpus (9)

• Ten equivalent definitions of the fractional Laplace operator
• Modified Szabo's wave equation models for lossy media obeying frequency power law
• The automatic solution of partial differential equations using a global spectral method
• Fractional derivative models for atmospheric dispersion of pollutants
• Tensor calculus in polar coordinates using Jacobi polynomials
• A contour method for time-fractional PDEs and an application to fractional viscoelastic beam equations
• Computing Equilibrium Measures with Power Law Kernels
• Computation of Power Law Equilibrium Measures on Balls of Arbitrary Dimension
• A fast sparse spectral method for nonlinear integro-differential Volterra equations with general kernels

Cited by in corpus (1)

• Modelling power-law ultrasound absorption using a time-fractional, static memory, Fourier pseudo-spectral method