A note on parallel preconditioning for the all-at-once solution of Riesz fractional diffusion equations

The $p$-step backwards difference formula (BDF) for solving the system of ODEs can result in a kind of all-at-once linear systems, which are solved via the parallel-in-time preconditioned Krylov subspace solvers (see McDonald, Pestana, and Wathen [SIAM J. Sci. Comput., 40(2) (2018): A1012-A1033] and Lin and Ng [arXiv:2002.01108, 17 pages]. However, these studies ignored that the $p$-step BDF ($p\geq 2$) is not selfstarting, when they are exploited to solve time-dependent PDEs. In this note, we focus on the 2-step BDF which is often superior to the trapezoidal rule for solving the Riesz fractional diffusion equations, but its resultant all-at-once discretized system is a block triangular Toeplitz system with a low-rank perturbation. Meanwhile, we first give an estimation of the condition number of the all-at-once systems and then adapt the previous work to construct two block circulant (BC) preconditioners. Both the invertibility of these two BC preconditioners and the eigenvalue distributions of preconditioned matrices are discussed in details. The efficient implementation of these BC preconditioners is also presented especially for handling the computation of dense structured Jacobi matrices. Finally, numerical experiments involving both the one- and two-dimensional Riesz fractional diffusion equations are reported to support our theoretical findings.Comment: 18 pages. 2 figures. 6 Table. Tech. Rep.: Institute of Mathematics, Southwestern University of Finance and Economics. Revised-1: refine/shorten the contexts and correct some typos; Revised-2: correct some reference

A Finite Element Method for the Multiterm Time-Space Riesz Fractional Advection-Diffusion Equations in Finite Domain

We present an effective finite element method (FEM) for the multiterm time-space Riesz fractional advection-diffusion equations (MT-TS-RFADEs). We obtain the weak formulation of MT-TS-RFADEs and prove the existence and uniqueness of weak solution by the Lax-Milgram theorem. For multiterm time discretization, we use the Diethelm fractional backward finite difference method based on quadrature. For spatial discretization, we show the details of an FEM for such MT-TS-RFADEs. Then, stability and convergence of such numerical method are proved, and some numerical examples are given to match well with the main conclusions

Numerical methods for solving space fractional partial differential equations by using Hadamard finite-part integral approach

We introduce a novel numerical method for solving two-sided space fractional partial differential equation in two dimensional case. The approximation of the space fractional Riemann-Liouville derivative is based on the approximation of the Hadamard finite-part integral which has the convergence order $O(h^{3- \alpha})$, where $h$ is the space step size and $\alpha\in (1, 2)$ is the order of Riemann-Liouville fractional derivative. Based on this scheme, we introduce a shifted finite difference method for solving space fractional partial differential equation. We obtained the error estimates with the convergence orders $O(\tau +h^{3-\alpha}+ h^{\beta})$, where $\tau$ is the time step size and $\beta >0$ is a parameter which measures the smoothness of the fractional derivatives of the solution of the equation. Unlike the numerical methods for solving space fractional partial differential equation constructed by using the standard shifted Gr\"unwald-Letnikov formula or higher order Lubich'e methods which require the solution of the equation satisfies the homogeneous Dirichlet boundary condition in order to get the first order convergence, the numerical method for solving space fractional partial differential equation constructed by using Hadamard finite-part integral approach does not require the solution of the equation satisfies the Dirichlet homogeneous boundary condition. Numerical results show that the experimentally determined convergence order obtained by using the Hadamard finite-part integral approach for solving space fractional partial differential equation with non-homogeneous Dirichlet boundary conditions is indeed higher than the convergence order obtained by using the numerical methods constructed with the standard shifted Gr\"unwald-Letnikov formula or Lubich's higer order approximation schemes

Nonlinear Evolution Equations: Analysis and Numerics

The workshop was devoted to the analytical and numerical investigation of nonlinear evolution equations. The main aim was to stimulate a closer interaction between experts in analytical and numerical methods for areas such as wave and Schrödinger equations or the Navier–Stokes equations and fluid dynamics