Variable-order fractional wave equation: Analysis, numerical approximation, and fast algorithm

TL;DR

This paper studies a variable-order fractional wave equation modeling diffusive wave propagation in viscoelastic media, providing analysis, a finite element scheme, and a fast algorithm to improve computational efficiency.

Contribution

It introduces a novel fast divide-and-conquer algorithm for the variable-order fractional wave equation, reducing computational complexity significantly.

Findings

The model is well-posed with high order regularity.

The finite element scheme achieves rigorous error estimates.

The fast algorithm reduces computational complexity from O(MN^2) to O(MN log^2 N).

Abstract

We investigate a local modification of a variable-order fractional wave equation, which describes the propagation of diffusive wave in viscoelastic media with evolving physical property. We incorporate an equivalent formulation to prove the well-posedness of the model as well as its high order regularity estimates. To accommodate the convolution term in the reformulated model, we adopt the Ritz-Volterra finite element projection and then derive the rigorous error estimate for the fully-discretized finite element scheme. To circumvent the high computational cost from the temporal integral term, we exploit the translational invariance of the discrete coefficients associated with the convolution structure and construct a fast divide-and-conquer algorithm which reduces the computational complexity from $O(MN^2)$ to $O(MN\log^2 N)$. Numerical experiments are provided to verify the theoretical results and to demonstrate the accuracy and efficiency of the proposed method.