On the bilateral preconditioning for an L2-type all-at-once system arising from time-space fractional Bloch-Torrey equations
arXiv:2109.06510 · doi:10.1016/j.camwa.2023.08.014
Abstract
Time-space fractional Bloch-Torrey equations (TSFBTEs) are developed by some researchers to investigate the relationship between diffusion and fractional-order dynamics. In this paper, we first propose a second-order implicit difference scheme for TSFBTEs by employing the recently proposed L2-type formula [A.~A.~Alikhanov, C.~Huang, Appl.~Math.~Comput.~(2021) 126545]. Then, we prove the stability and the convergence of the proposed scheme. Based on such a numerical scheme, an L2-type all-at-once system is derived. In order to solve this system in a parallel-in-time pattern, a bilateral preconditioning technique is designed to accelerate the convergence of Krylov subspace solvers according to the special structure of the coefficient matrix of the system. We theoretically show that the condition number of the preconditioned matrix is uniformly bounded by a constant for the time fractional order . Numerical results are reported to show the efficiency of our method.
arXiv comment25 pages, 6 tables, 4 figures
References in corpus (8)
- A discrete Grönwall inequality with application to numerical schemes for subdiffusion problems
- A preconditioning technique for an all-at-once system from Volterra subdiffusion equations with graded time steps
- A preconditioning technique for all-at-once system from the nonlinear tempered fractional diffusion equation
- A fast implicit difference scheme for solving the generalized time-space fractional diffusion equations with variable coefficients
- Fast implicit difference schemes for time-space fractional diffusion equations with the integral fractional Laplacian
- A parallel-in-time two-sided preconditioning for all-at-once system from a non-local evolutionary equation with weakly singular kernel
- An efficient second-order energy stable BDF scheme for the space fractional Cahn-Hilliard equation
- Numerical algorithm for the space-time fractional Fokker-Planck system with two internal states