Abstract:In this work, we propose an efficient and robust multigrid method for solving the time-fractional heat equation. Due to the nonlocal property of fractional differential operators, numerical methods usually generate systems of equations for which the coefficient matrix is dense. Therefore, the design of efficient solvers for the numerical simulation of these problems is a difficult task. We develop a parallel-in-time multigrid algorithm based on the waveform relaxation approach, whose application to time-fractional problems seems very natural due to the fact that the fractional derivative at each spatial point depends on the values of the function at this point at all earlier times. Exploiting the Toeplitz-like structure of the coefficient matrix, the proposed multigrid waveform relaxation method has a computational cost of $O(N M \log(M))$ operations, where $M$ is the number of time steps and $N$ is the number of spatial grid points. A semi-algebraic mode analysis is also developed to theoretically confirm the good results obtained. Several numerical experiments, including examples with non-smooth solutions and a nonlinear problem with applications in porous media, are presented.

What problem does this paper attempt to address?

### What problem does this paper attempt to solve? This paper aims to solve the problem of efficient numerical solution of the time - fractional heat equation. Due to the non - local nature of the fractional - order differential operator, numerical methods usually generate systems of equations with dense coefficient matrices, which makes it a difficult problem to design efficient solvers. Specifically: 1. **Problem background**: - Fractional - order calculus has frequently emerged in fields such as fluid mechanics, signal processing, viscoelasticity, porous media flow, quantum mechanics, biology, medicine, physics, and engineering in recent years. - Fractional - order partial differential equations (PDEs), due to their non - local characteristics, lead to the coefficient matrices of the linear systems generated by numerical methods usually being dense. Traditional solution methods such as Gaussian elimination have a computational complexity as high as $O(n^3)$ and a storage cost of $O(n^2)$. 2. **Research objectives**: - Propose a parallel - time multigrid waveform relaxation method based on the waveform relaxation method for efficiently solving the time - fractional heat equation. - Utilize the Toeplitz - like structure of the coefficient matrix to make the proposed multigrid waveform relaxation method have a computational complexity of $O(NM\log(M))$, where $M$ is the number of time steps and $N$ is the number of spatial grid points. - Verify the good convergence of the proposed method through semi - algebraic mode analysis (SAMA) theory. 3. **Main contributions**: - Propose a time - parallel multigrid waveform relaxation method directly applied to the original discretized problem. - This method is applicable to nonlinear problems and shows good robustness and efficiency when dealing with time - fractional PDEs. 4. **Practical applications**: - The paper presents several numerical experiments, including examples of non - smooth solutions and nonlinear problems applied in porous media, which verify the effectiveness of the proposed method. In summary, the main purpose of this paper is to develop an efficient and robust multigrid method to solve the numerical solution problem of the time - fractional heat equation, especially making progress in dealing with the computational challenges brought by the non - locality caused by fractional - order derivatives and dense matrices.

Multigrid waveform relaxation for the time-fractional heat equation

Using hierarchical matrices in the solution of the time-fractional heat equation by multigrid waveform relaxation

Multigrid methods for time fractional conservation laws

Space-time waveform relaxation multigrid for Navier-Stokes

Fully Finite Element Adaptive Algebraic Multigrid Method for Time-Space Caputo-Riesz Fractional Diffusion Equations

An efficient multigrid solver for two-dimensional spatial fractional diffusion equations with variable coefficients

Efficient Multigrid Reduction-in-Time for Method-of-Lines Discretizations of Linear Advection

A multigrid-reduction-in-time solver with a new two-level convergence for unsteady fractional Laplacian problems

Multigrid for two-sided fractional differential equations discretized by finite volume elements on graded meshes

A Fast Multigrid Finite Element Method for the Time-Dependent Tempered Fractional Problem

A static memory sparse spectral method for time-fractional PDEs

Double fast algorithm for solving time-space fractional diffusion problems with spectral fractional Laplacian

Multigrid method for symmetric Toeplitz tridiagonal matrix: Convergence analysis and application to fractional Feynman-Kac equation

Fast time-stepping discontinuous Galerkin method for the subdiffusion equation

Fourier pseudospectral methods for the spatial variable-order fractional wave equations

A Space-Time Multigrid Method for Space-Time Finite Element Discretizations of Parabolic and Hyperbolic PDEs

Efficient computational hybrid method for the solution of 2D multi-term fractional order advection-diffusion equation

Multigrid reduction-in-time convergence for advection problems: A Fourier analysis perspective

Finite element discretizations for variable-order fractional diffusion problems

A fast finite volume method for spatial fractional diffusion equations on nonuniform meshes

Two-grid methods for nonlinear time fractional diffusion equations by L 1 -Galerkin FEM