Pei-Yong Zhu

Yong-Liang Zhao, Pei-Yong Zhu, Xian-Ming Gu et al.

A block lower triangular Toeplitz system arising from time-space fractional diffusion equation is discussed. For efficient solutions of such the linear system, the preconditioned biconjugate gradient stabilized method and flexible general minimal residual method are exploited. The main contribution of this paper has two aspects: (i) A block bi-diagonal Toeplitz preconditioner is developed for the block lower triangular Toeplitz system, whose storage is of $\mathcal{O}(N)$ with $N$ being the spatial grid number; (ii) A new skew-circulant preconditioner is designed to fast calculate the inverse of the block bi-diagonal Toeplitz preconditioner multiplying a vector. Numerical experiments are given to demonstrate the efficiency of our preconditioners.

Yong-Liang Zhao, Pei-Yong Zhu, Xian-Ming Gu et al.

An implicit finite difference scheme based on the $L2$-$1_σ$ formula is presented for a class of one-dimensional time fractional reaction-diffusion equations with variable coefficients and time drift term. The unconditional stability and convergence of this scheme are proved rigorously by the discrete energy method, and the optimal convergence order in the $L_2$-norm is $\mathcal{O}(τ^2 + h^2)$ with time step $τ$ and mesh size $h$. Then, the same measure is exploited to solve the two-dimensional case of this problem and a rigorous theoretical analysis of the stability and convergence is carried out. Several numerical simulations are provided to show the efficiency and accuracy of our proposed schemes and in the last numerical experiment of this work, three preconditioned iterative methods are employed for solving the linear system of the two-dimensional case.

Yong-Liang Zhao, Pei-Yong Zhu, Xian-Ming Gu et al.

A kind of spatial fractional diffusion equations in this paper are studied. Firstly, an L1 formula is employed for the spatial discretization of the equations. Then, a second order scheme is derived based on the resulting semi-discrete ordinary differential system by using the implicit integration factor method, which is a class of efficient semi-implicit temporal scheme. Numerical results show that the proposed scheme is accurate even for the discontinuous coefficients.