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

•An efficient solver for the time-fractional heat equation discretized on non- uniform temporal meshes is proposed.•The solver is a space-time multigrid method based on the waveform relaxation technique.•Hierarchical matrix approximations are considered to maintain optimal complexity of the algorith...

Celý popis

Uloženo v:
Podrobná bibliografie
Vydáno v:Journal of computational physics Ročník 416; s. 109540
Hlavní autoři: Hu, Xiaozhe, Rodrigo, Carmen, Gaspar, Francisco J.
Médium: Journal Article
Jazyk:angličtina
Vydáno: Cambridge Elsevier Inc 01.09.2020
Elsevier Science Ltd
Témata:
Algorithms
Computational physics
Differential equations
Discretization
Graded meshes
Hierarchical matrices
Multigrid waveform relaxation
Operators (mathematics)
Semi-algebraic mode analysis
Thermodynamics
Time-fractional heat equation
Waveforms
Multigrid waveform relaxation
Hierarchical matrices
Graded meshes
Time-fractional heat equation
Semi-algebraic mode analysis
ISSN:0021-9991, 1090-2716
On-line přístup:Získat plný text
Tagy: Přidat tag
Žádné tagy, Buďte první, kdo vytvoří štítek k tomuto záznamu!
Popis
Shrnutí:•An efficient solver for the time-fractional heat equation discretized on non- uniform temporal meshes is proposed.•The solver is a space-time multigrid method based on the waveform relaxation technique.•Hierarchical matrix approximations are considered to maintain optimal complexity of the algorithm.•A semi-algebraic analysis supports the numerical findings. This work deals with the efficient numerical solution of the time-fractional heat equation discretized on non-uniform temporal meshes. Non-uniform grids are essential to capture the singularities of “typical” solutions of time-fractional problems. We propose an efficient space-time multigrid method based on the waveform relaxation technique, which accounts for the nonlocal character of the fractional differential operator. To maintain an optimal complexity, which can be obtained for the case of uniform grids, we approximate the coefficient matrix corresponding to the temporal discretization by its hierarchical matrix (H-matrix) representation. In particular, the proposed method has a computational cost of O(kNMlog⁡(M)), where M is the number of time steps, N is the number of spatial grid points, and k is a parameter which controls the accuracy of the H-matrix approximation. The efficiency and the good convergence of the algorithm, which can be theoretically justified by a semi-algebraic mode analysis, are demonstrated through numerical experiments in both one- and two-dimensional spaces.
Bibliografie:ObjectType-Article-1
SourceType-Scholarly Journals-1
ObjectType-Feature-2
content type line 14
ISSN:0021-9991
1090-2716
DOI:10.1016/j.jcp.2020.109540