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

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.