Abstract:We present a finite element scheme for fractional diffusion problems with varying diffusivity and fractional order. We consider a symmetric integral form of these nonlocal equations defined on general geometries and in arbitrary bounded domains. A number of challenges are encountered when discretizing these equations. The first comes from the heterogeneous kernel singularity in the fractional integral operator. The second comes from the dense discrete operator with its quadratic growth in memory footprint and arithmetic operations. An additional challenge comes from the need to handle volume conditions-the generalization of classical local boundary conditions to the nonlocal setting. Satisfying these conditions requires that the effect of the whole domain, including both the interior and exterior regions, can be computed on every interior point in the discretization. Performed directly, this would result in quadratic complexity. To address these challenges, we propose a strategy that decomposes the stiffness matrix into three components. The first is a sparse matrix that handles the singular near-field separately and is computed by adapting singular quadrature techniques available for the homogeneous case to the case of spatially variable order. The second component handles the remaining smooth part of the near-field as well as the far field and is approximated by a hierarchical $\mathcal{H}^{2}$ matrix that maintains linear complexity in storage and operations. The third component handles the effect of the global mesh at every node and is written as a weighted mass matrix whose density is computed by a fast-multipole type method. The resulting algorithm has therefore overall linear space and time complexity. Analysis of the consistency of the stiffness matrix is provided and numerical experiments are conducted to illustrate the convergence and performance of the proposed algorithm.

An IRK‐QCD Scheme for the Space Fractional Diffusion Equation and Block‐Lower‐Triangle Preconditioners for the Corresponding Linear Systems

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

An Implicit Numerical Method for the Riemann–Liouville Distributed-Order Space Fractional Diffusion Equation

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

Fourth order quasi-compact difference schemes for (tempered) space fractional diffusion equations

Differential quadrature method for space-fractional diffusion equations on 2D irregular domains

A fourth-order compact solver for fractional-in-time fourth-order diffusion equations

Finite element discretizations for variable-order fractional diffusion problems

A COMPACT DIFFERENCE SCHEME FOR FOURTH-ORDER FRACTIONAL SUB-DIFFUSION EQUATIONS WITH NEUMANN BOUNDARY CONDITIONS

A numerical study on solving a fractional time-space diffusion equation via the finite difference method

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

Fast implicit difference schemes for time-space fractional diffusion equations with the integral fractional Laplacian

Efficient Numerical Schemes for Fractional Sub-Diffusion Equation with the Spatially Variable Coefficient

An Efficient QSC Approximation of Variable-Order Time-Fractional Mobile-Immobile Diffusion Equations with Variably Diffusive Coefficients

High order numerical methods based on quadratic spline collocation method and averaged L1 scheme for the variable-order time fractional mobile/immobile diffusion equation

A Temporal Second-Order Difference Scheme for Variable-Order-Time Fractional-Sub-Diffusion Equations of the Fourth Order

Strang-type Preconditioners for Solving Fractional Diffusion Equations by Boundary Value Methods

Efficient Variable Steps BDF2 Scheme for the Two-Dimensional Space Fractional Cahn-Hilliard Model

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

A Fourth-Order Compact ADI Scheme for Two-Dimensional Riesz Space Fractional Nonlinear Reaction-Diffusion Equation

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