An optimal preconditioner for high-order scheme arising from multi-dimensional Riesz space fractional diffusion equations with variable coefficients

Abstract

In this paper, we propose an efficient method for solving multi-dimensional Riesz space fractional diffusion equations with variable coefficients. The Crank-Nicolson (CN) method is used for temporal discretization, while the fourth-order fractional centered difference (4FCD) method is employed for spatial discretization. Using a novel technique, we show that the CN-4FCD scheme for the multi-dimensional case is unconditionally stable and convergent, achieving second-order accuracy in time and fourth-order accuracy in space with respect to the discrete L2-norm. Moreover, leveraging the symmetric multi-level Toeplitz-like structure of the coefficient matrix in the discrete linear systems, we enhance the computational efficiency of the proposed scheme with a sine transform-based preconditioner, ensuring a mesh-size-independent convergence rate for the conjugate gradient method. Finally, two numerical examples validate the theoretical analysis and demonstrate the superior performance of the proposed preconditioner compared to existing methods.

0

Turn this paper into a full lesson

ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…