Multigrid Waveform Relaxation for the Time-Fractional Heat Equation

From MaRDI portal
Publication:5348253

DOI10.1137/16M1090193zbMATH Open1371.65103arXiv1608.05400OpenAlexW2963027558MaRDI QIDQ5348253FDOQ5348253


Authors: Carmen Rodrigo, Francisco José Gaspar Edit this on Wikidata


Publication date: 15 August 2017

Published in: SIAM Journal on Scientific Computing (Search for Journal in Brave)

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(NMlog(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.


Full work available at URL: https://arxiv.org/abs/1608.05400




Recommendations




Cites Work


Cited In (21)

Uses Software





This page was built for publication: Multigrid Waveform Relaxation for the Time-Fractional Heat Equation

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q5348253)