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 integral discretization scheme on a graded mesh for a fractional differential equation with integral boundary conditions

A Compact Difference Scheme on Graded Meshes for the Nonlinear Fractional Integro-differential Equation with Non-smooth Solutions

An efficient numerical method for a time-fractional diffusion equation

Fast compact finite difference schemes on graded meshes for fourth-order multi-term fractional sub-diffusion equations with the first Dirichlet boundary conditions

Finite Difference Schemes for the Variable Coefficients Single and Multi-Term Time-Fractional Diffusion Equations with Non-Smooth Solutions on Graded and Uniform Meshes

Numerical approximation for a nonlinear variable-order fractional differential equation via an integral equation method

Local analysis of L1-finite difference method on graded meshes for multi-term two-dimensional time-fractional initial-boundary value problem with Neumann boundary conditions

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

High order approximation to Caputo derivative on graded mesh and time-fractional diffusion equation for non-smooth solutions

A second-order finite difference scheme for nonlinear tempered fractional integrodifferential equations in three dimensions

Finite Difference Scheme on Graded Meshes to the Time-Fractional Neutron Diffusion Equation with Non-Smooth Solutions

Error analysis for a fractional-derivative parabolic problem on quasi-graded meshes using barrier functions

A finite difference scheme for fractional sub-diffusion equations on an unbounded domain using artificial boundary conditions

Two Finite Difference Schemes for Multi-Dimensional Fractional Wave Equations with Weakly Singular Solutions

Finite element method combined with time graded meshes for the time-fractional coupled Burgers' equations

Finite element discretizations for variable-order fractional diffusion problems

Successively Iterative Method for Fractional Differential Equations with Integral Boundary Conditions

Some Grönwall Inequalities for a Class of Discretizations of Time Fractional Equations on Nonuniform Meshes

Difference Approximation for 2D Time-Fractional Integro-Differential Equation with Given Initial and Boundary Conditions

Numerical Methods for Caputo–Hadamard Fractional Differential Equations with Graded and Non-Uniform Meshes

A Discretization of Caputo Derivatives with Application to Time Fractional SDEs and Gradient Flows