Abstract:Finite difference method as a popular numerical method has been widely used to solve fractional diffusion equations. In the general spatial error analyses, an assumption $u\in C^{4}(\bar{\Omega})$ is needed to preserve $\mathcal{O}(h^{2})$ convergence when using central finite difference scheme to solve fractional sub-diffusion equation with Laplace operator, but this assumption is somewhat strong, where $u$ is the exact solution and $h$ is the mesh size. In this paper, a novel analysis technique is proposed to show that the spatial convergence rate can reach $\mathcal{O}(h^{\min(\sigma+\frac{1}{2}-\epsilon,2)})$ in both $l^{2}$-norm and $l^{\infty}$-norm in one-dimensional domain when the initial value and source term are both in $\hat{H}^{\sigma}(\Omega)$ but without any regularity assumption on the exact solution, where $\sigma\geq 0$ and $\epsilon>0$ being arbitrarily small. After making slight modifications on the scheme, acting on the initial value and source term, the spatial convergence rate can be improved to $\mathcal{O}(h^{2})$ in $l^{2}$-norm and $\mathcal{O}(h^{\min(\sigma+\frac{3}{2}-\epsilon,2)})$ in $l^{\infty}$-norm. It's worth mentioning that our spatial error analysis is applicable to high dimensional cube domain by using the properties of tensor product. Moreover, two kinds of averaged schemes are provided to approximate the Riemann--Liouville fractional derivative, and $\mathcal{O}(\tau^{2})$ convergence is obtained for all $\alpha\in(0,1)$. Finally, some numerical experiments verify the effectiveness of the built theory.

Localized Bases for Finite-Dimensional Homogenization Approximations with Nonseparated Scales and High Contrast.

Localized bases for finite dimensional homogenization approximations with non-separated scales and high-contrast

Polyharmonic homogenization, rough polyharmonic splines and sparse super-localization

Metric-Based Upscaling

Hybrid Localized Spectral Decomposition for multiscale problems

Localized Harmonic Characteristic Basis Functions for Multiscale Finite Element Methods

Local Convergence of Inexact Methods under the Hölder Condition

Adaptive Local (AL) Basis for Elliptic Problems with $L^\infty$-Coefficients

Homogenization of Nondivergence-Form Elliptic Equations with Discontinuous Coefficients and Finite Element Approximation of the Homogenized Problem

N A ] 1 1 M ay 2 00 5 Metric based upscaling

Homogenization of Parabolic Equations with a Continuum of Space and Time Scales

Local Hardy spaces associated with inhomogeneous higher order elliptic operators

Hierarchical Super-Localized Orthogonal Decomposition Method

Computational multiscale methods for nondivergence-form elliptic partial differential equations

ITERATED NUMERICAL HOMOGENIZATION FOR MULTISCALE ELLIPTIC EQUATIONS WITH MONOTONE NONLINEARITY

Super-localized orthogonal decomposition for convection-dominated diffusion problems

Homogenization of Fundamental Solutions for Parabolic Operators Involving Non-Self-similar Scales

Localized subspace iteration methods for elliptic multiscale problems

Compactness and Stable Regularity in Multiscale Homogenization

Sharp error estimates for spatial-temporal finite difference approximations to fractional sub-diffusion equation without regularity assumption on the exact solution