Abstract:In this paper, we develop a fast numerical method for solving the time-dependent Riesz space fractional diffusion equations with a nonlinear source term in the convex domain. An implicit finite difference method is employed to discretize the Riesz space fractional diffusion equations with a penalty term in a rectangular region by the volume-penalization approach. The stability and the convergence of the proposed method are studied. As the coefficient matrix is with the Toeplitz-like structure, the generalized minimum residual method with a preconditioner based on the sine transform is exploited to solve the discretized linear system, where the preconditioner is constructed in view of the combination of two approximate inverse ${\tau}$ matrices, which can be diagonalized by the sine transform. The spectrum of the preconditioned matrix is also investigated. Numerical experiments are carried out to demonstrate the efficiency of the proposed method.

What problem does this paper attempt to address?

The main problem that this paper attempts to solve is: the numerical calculation challenges encountered when solving the time - dependent Riesz - space fractional - order diffusion equations (RSFDEs) with nonlinear source terms in convex domains. Specifically, the article focuses on the following issues: 1. **Development of efficient numerical methods**: A fast numerical method based on the implicit finite - difference method is proposed to discretize the Riesz - space fractional - order diffusion equation, and a penalty term is introduced to handle the boundary conditions in convex domains. 2. **Special structure of the coefficient matrix**: Since the coefficient matrix has a Toeplitz - like structure, direct solution will lead to high computational cost and storage requirements. Therefore, the article explores how to use the generalized minimum residual method (GMRES) combined with a preconditioner based on the sine transform to accelerate the solution of the linear system. 3. **Stability and convergence analysis**: The stability and convergence of the proposed numerical method are studied to ensure its reliability and accuracy in practical applications. 4. **Construction and optimization of preconditioners**: In order to improve the convergence speed of the iterative method, the article proposes a method for constructing preconditioners based on the sine transform, and further optimizes the preconditioner by combining two approximate inverse τ matrices, enabling it to be calculated in O(N log N) time. 5. **Verification by numerical experiments**: The effectiveness of the proposed method is verified through a series of numerical experiments, demonstrating its superior performance in solving Riesz - space fractional - order diffusion equations in convex domains. In summary, this paper aims to provide an efficient and stable numerical method for solving Riesz - space fractional - order diffusion equations in convex domains, and reduce the computational complexity and improve the solution efficiency through reasonable preconditioner design.

A preconditioner based on sine transform for two-dimensional Riesz space factional diffusion equations in convex domains

Spectral analysis for preconditioning of multi-dimensional Riesz fractional diffusion equations

Sine Transform Based Preconditioning for an Inverse Source Problem of Time-Space Fractional Diffusion Equations

A -Preconditioner for Space Fractional Diffusion Equation with Non-separable Variable Coefficients

A $τ$-preconditioner for space fractional diffusion equation with non-separable variable coefficients

Generalized Finite Difference Method for a Class of Multidimensional Space-Fractional Diffusion Equations

A novel fourth-order scheme for two-dimensional Riesz space fractional nonlinear reaction-diffusion equations and its optimal preconditioned solver

On $τ$-preconditioners for a quasi-compact difference scheme to Riesz fractional diffusion equations with variable coefficients

Circulant-based approximate inverse preconditioners for a class of fractional diffusion equations

Sine-transform-based fast solvers for Riesz fractional nonlinear Schrödinger equations with attractive nonlinearities

An α -robust fast algorithm for distributed-order time–space fractional diffusion equation with weakly singular solution

Multilevel Tau preconditioners for symmetrized multilevel Toeplitz systems with applications to solving space fractional diffusion equations

Stiff-cut leap-frog scheme for fractional Laplacian diffusion equations

A fast algorithm for two-dimensional distributed-order time-space fractional diffusion equations

A fast finite volume method for spatial fractional diffusion equations on nonuniform meshes

Double fast algorithm for solving time-space fractional diffusion problems with spectral fractional Laplacian

A τ Matrix Based Approximate Inverse Preconditioning for Tempered Fractional Diffusion Equations

An efficient multigrid method with preconditioned smoother for two-dimensional anisotropic space-fractional diffusion equations

Two fast and unconditionally stable finite difference methods for Riesz fractional diffusion equations with variable coefficients

Accurate Boundary Treatment for Riesz Space Fractional Diffusion Equations

R. Chan's circulant-based approximate inverse preconditioning iterative method for solving second-order space fractional advection–dispersion equations with variable coefficients