A multigrid solver for subdiffusion equations

Mokhtari, Reza; Ramezani, Mohadeseh; Haase, Gundolf

doi:10.22034/cmde.2025.65520.3017

A multigrid solver for subdiffusion equations

Document Type : Research Paper

Authors

¹ 1. Department of Mathematical Sciences, Isfahan University of Technology, Isfahan 84156-83111, Iran. 2. Department for Mathematics and Scientific Computing, University of Graz, Graz, Austria.

² Department of Mathematical Sciences, Isfahan University of Technology, Isfahan 84156-83111, Iran.

10.22034/cmde.2025.65520.3017

Abstract

In this paper, we investigate the S3-FD method for solving time-fractional diffusion equations in both one-dimensional (1D) and two-dimensional (2D) spatial domains, achieving high-order temporal accuracy. We leverage the S3 formula, which has a temporal accuracy of \(4 - \alpha\), to approximate the Caputo fractional derivative of order \(\alpha \in (0,1)\), and we employ the finite difference approach for spatial discretization. We develop a fully discrete scheme for both uniform and non-uniform spatial meshes. Our analysis begins with the 1D subdiffusion problem, where we employ the cyclic reduction method alongside OpenMP-based parallel programming to reduce computational costs. Leveraging this groundwork, we extend our technique to the 2D subdiffusion problem using a multigrid method and domain decomposition strategy paired with MPI programming. This innovative method yields an impressive temporal convergence order of \(\mathcal{O}(\Delta t^{4-\alpha})\). The performance and efficiency of the proposed S3-FD algorithm are demonstrated through numerical experiments, highlighting its potential for large-scale fractional diffusion problems.

Keywords

Main Subjects

Fractional partial differential equations

References

[1] A. A. Alikhanov, A new difference scheme for the time fractional diffusion equation, J. Comput. Phys., 280 (2015), 424–438.
[2] C. C. Douglas, G. Haase, and U. Langer, A tutorial on elliptic PDE solvers and their parallelization, Society for Industrial and Applied Mathematics, 2003.
[3] D. Gan, G. F. Zhang, and Z. Z. Liang, An efficient preconditioner for linear systems arising from high-order accurate schemes of time fractional diffusion equations, J. Appl. Math. Comput., 70 (2024), 5129–5151.
[4] G. H. Gao, Z. Z. Sun, and H. W. Zhang, A new fractional numerical differentiation formula to approximate the Caputo fractional derivative and its applications, J. Comput. Phys., 259 (2014), 33–50.
[5] X. Hu, C. Rodrigo, and F. J. Gaspar, Using hierarchical matrices in the solution of the time-fractional heat equation by multigrid waveform relaxation, J. Comput. Phys., 416 (2020), 109540.
[6] Y. Jiang and H. Bai, Multigrid methods for time fractional conservation laws, Numer. Algorithms, 97 (2024), 1301–1322.
[7] Y. Jiang, H. Chen, T. Sun, and C. Huang, Efficient L1-ADI finite difference method for the two-dimensional nonlinear time-fractional diffusion equation, Appl. Math. Comput., 471 (2024), 128609.
[8] A. Khibiev, A. Alikhanov, and C. Huang, A second-order difference scheme for generalized time-fractional diffusion equation with smooth solutions, Comput. Methods Appl. Math., 24 (2024), 101–117.
[9] Y. Li, L. Zikatanov, and C. Zuo, A reduced conjugate gradient basis method for fractional diffusion, SIAM J. Sci. Comput., 46 (2024), S68–S87.
[10] T. A. M. Langlands and B. I. Henry, The accuracy and stability of an implicit solution method for the fractional diffusion equation, J. Comput. Phys., 205 (2005), 719–736.
[11] K. Oldham and J. Spanier, The fractional calculus: Theory and applications of differentiation and integration to arbitrary order, Elsevier, 1974.
[12] J. Liu, H. Fu, and J. Zhang, A QSC method for fractional subdiffusion equations with fractional boundary conditions and its application in parameters identification, Math. Comput. Simul., 174 (2020), 153–174.
[13] S. Maji and S. Natesan, Adaptive-grid technique for the numerical solution of a class of fractional boundary-value- problems, Comput. Methods Differ. Equ., 12 (2024), 338–349.
[14] R. Mokhtari and F. Mostajeran, A high order formula to approximate the Caputo fractional derivative, Commun. Appl. Math. Comput., 2 (2020), 1–29.
[15] R. Mokhtari, M. Ramezani, and G. Haase, Stability and convergence analyses of the FDM based on some L-type formulae for solving the subdiffusion equation, Numer. Math. Theory Methods Appl., 14 (2021), 1–27.
[16] K. Pan, H. W. Sun, Y. Xu, and Y. Xu, An efficient multigrid solver for two-dimensional spatial fractional diffusion equations with variable coefficients, Appl. Math. Comput., 402 (2021), 126091.
[17] M. Ramezani and R. Mokhtari, A novel high-order finite-difference method for the time-fractional diffusion equation with smooth/nonsmooth solutions, Bull. Iran. Math. Soc., 48 (2022), 3987–4013.
[18] M. Ramezani, R. Mokhtari, and G. Haase, Some high order formulae for approximating Caputo fractional derivatives, Appl. Numer. Math., 153 (2020), 300–318.
[19] M. Ramezani, R. Mokhtari, and G. Haase, Analysis of stability and convergence for L-type formulas combined with a spatial finite element method for solving subdiffusion problems, Electron. Trans. Numer. Anal., (2022), 568–584.
[20] M. Ramezani, R. Mokhtari, and Y. Yan, Correction of a high-Order numerical method for approximating time- fractional wave equation, J. Sci. Comput., 100 (2024), 71.
[21] A. Singh, S. Kumar, and H. Ramos, An efficient computational method based on exponential B-splines for a class of fractional sub-diffusion equations, Comput. Methods Differ. Equ., 12 (2024), 719–740.
[22] M. Stynes, E. O’Riordan, and J. L. Gracia, Error analysis of a finite difference method on graded meshes for a time-fractional diffusion equation, SIAM J. Numer. Anal., 55(2) (2017), 1057–1079.
[23] S. P. Tang and Y. M. Huang, A fast preconditioning iterative method for solving the discretized second-order space-fractional advection–diffusion equations, J. Comput. Appl. Math., 438 (2024), 115513.
[24] J. Tan and J. Liu, An efficient numerical solver for anisotropic subdiffusion problems, J. Comput. Appl. Math., 364 (2020), 112318.
[25] Y. Wang, N. An, and C. Huang, Unconditional optimal error bounds of the fast nonuniform Alikhanov scheme for a nonlinear time-fractional biharmonic equation, J. Appl. Math. Comput., 70 (2024), 4053–4071.
[26] Y. Xu, S. L. Lei, and H. W. Sun, An efficient multigrid method with preconditioned smoother for two-dimensional anisotropic space-fractional diffusion equations, Comput. Math. Appl., 124 (2022), 218–226.
[27] Y. Yan, M. Khan, and N. J. Ford, An analysis of the modified L1 scheme for time-fractional partial differential equations with nonsmooth data, SIAM J. Numer. Anal., 56 (2018), 210–227.
[28] S. Yeganeh, R. Mokhtari, and J. S. Hesthaven, A local discontinuous Galerkin method for two-dimensional time fractional diffusion equations, Commun. Appl. Math. Comput., 2 (2020), 689–709.
[29] Y. Zhao, Y. Zhang, F. Liu, I. Turner, Y. Tang, and V. Anh, Convergence and superconvergence of a fully-discrete scheme for multi-term time fractional diffusion equations, Comput. Math. Appl., 73 (2017), 1087–1099.
[30] Z. Zhou, S. Zhang, and W. Li, The splitting characteristic finite difference domain decomposition scheme for solving time-fractional MIM nonlinear advection–diffusion equations, J. Sci. Comput., 100 (2024), 49.