A multigrid solver for subdiffusion equations

Mokhtari, Reza; Ramezani, Mohadeseh; Haase, Gundolf

doi:10.22034/cmde.2025.65520.3017

فهرست نشریات دارای اعتبار وزارت علوم، تحقیقات و فناوری

تعداد نشریات	45
تعداد شماره‌ها	1,531
تعداد مقالات	18,710
تعداد مشاهده مقاله	61,574,108
تعداد دریافت فایل اصل مقاله	22,291,060

	A multigrid solver for subdiffusion equations
Computational Methods for Differential Equations
مقاله 26، دوره 14، شماره 3، مهر 2026، صفحه 1426-1441 اصل مقاله (471.15 K)
نوع مقاله: Research Paper
شناسه دیجیتال (DOI): 10.22034/cmde.2025.65520.3017
نویسندگان
Reza Mokhtari^* ¹؛ Mohadeseh Ramezani²؛ Gundolf Haase³
¹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.
³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.
چکیده
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.
کلیدواژه‌ها
S3 formula؛ Subdiffusion equation؛ Multigrid method؛ Domain decomposition؛ Parallel processing

مراجع
[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.
آمار تعداد مشاهده مقاله: 278 تعداد دریافت فایل اصل مقاله: 256

سامانه مدیریت نشریات علمی. قدرت گرفته از سیناوب

پیوندهای مفید

آمار

A multigrid solver for subdiffusion equations