A static memory sparse spectral method for time-fractional PDEs

Timon S. Gutleb,José A. Carrillo

DOI: https://doi.org/10.1016/j.jcp.2023.112522

2023-10-11

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.

Numerical Analysis

What problem does this paper attempt to address?

The main problem that this paper attempts to solve is to develop an efficient and accurate numerical method for solving time - fractional partial differential equations (time - fractional PDEs). Specifically, this method aims to reduce the excessive memory requirements in traditional methods due to non - local fractional - order derivative operators. By combining the multivariable sparse spectral method and the fast Gaussian quadrature node calculation method, as well as the recursive non - classical method to handle the Caputo fractional - order derivative, this method can provide accurate numerical solutions without storing the historical values of functions.

### Main contributions of the paper
1. **Reduction of memory requirements**: A static - memory sparse spectral method is proposed, which can significantly reduce the memory requirements when solving time - fractional PDEs.
2. **Efficient numerical calculation**: The multivariable sparse spectral method and the fast Gaussian quadrature node calculation method are combined to improve the computational efficiency.
3. **Strong versatility**: This method has minimal theoretical overhead in any domain where an orthogonal polynomial basis can be used and is applicable to multiple types of PDEs.
4. **Error analysis and convergence strategy**: A detailed error analysis is provided, and a reasonable convergence strategy is proposed.
5. **Application prospects**: Special attention is paid to the propagation model of sound waves in viscoelastic media, which has important applications in medical ultrasonic treatment devices.

### Method overview
- **Caputo fractional - order derivative**: The Caputo fractional - order derivative used in the paper is defined as:
\[
\frac{\partial^\alpha f(t)}{\partial t^\alpha}=D_C^\alpha f(t)=\frac{1}{\Gamma(\lceil\alpha\rceil - \alpha)}\int_0^t(t - s)^{\lceil\alpha\rceil - \alpha - 1}f^{(\lceil\alpha\rceil)}(s)\,ds
\]
- **Sparse spectral method**: Functions are expanded using an orthogonal polynomial basis, so that linear operators (such as differentiation and integration) have a banded matrix form.
- **Recursive non - classical method**: By introducing an auxiliary function \(\psi_j(t)\) and using a recursive relationship, the storage of the historical values of functions is avoided, thus achieving a "history - free" calculation method.

### Numerical experiments
The paper verifies the effectiveness and accuracy of the method through multiple numerical experiments, including the solution of time - fractional PDEs on intervals, triangles and disks, and presents the error analysis results.

### Future research directions
The paper also discusses future research directions, especially how to apply this method to more complex multi - domain models, such as three - dimensional models of the human skull and brain.

### Conclusion
The method proposed in this paper is not only innovative in theory but also performs well in practical applications, especially in scenarios requiring high efficiency and low memory consumption, such as the simulation of medical ultrasonic treatment devices.