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.

ERROR ANALYSIS OF SOLUTION FOR GENERALIZED DENSITY EVOLUTION EQUATION VIA A FAMILY OF Δ SEQUENCES

SOLUTION OF GENERALIZED DENSITY EVOLUTION EQUATION VIA A FAMILY OF Δ SEQUENCES

Solution Of Generalized Density Evolution Equation Via A Family Of Delta Sequences

An Ensemble Evolution Numerical Method for Solving Generalized Density Evolution Equation

Improvement to the discretized initial condition of the generalized density evolution equation

A Difference-Wavelet Method for Solving Generalized Density Evolution Equation in Stochastic Structural Analysis

Density functions for the overdamped generalized Langevin equation and its Euler--Maruyama method: smoothness and convergence

A Numerical Approach with Residual Error Estimation for Eolution of High-order Linear Differential-difference Equations by Using Gegenbauer Polynomials

A modified Chebyshev collocation method for the generalized probability density evolution equation

The Principle of Preservation of Probability and the Generalized Density Evolution Equation

Error estimates of the Euler's method for stochastic differential equations with multiplicative noise via relative entropy

ADVANCES OF NUMERICAL METHODS FOR GENERAL DYNAMIC EQUATION DESCRIBING DISPERSED SYSTEM

A Unified Formalism of the GE-GDEE for Generic Continuous Responses and First-Passage Reliability Analysis of Multi-Dimensional Nonlinear Systems Subjected to Non-White-noise Excitations

Convergence analysis of exponential time differencing scheme for the nonlocal Cahn-Hilliard equation

A Priori and a Posteriori Error Analysis of the First Order Hyperbolic Equation by Using DG Method

Probabilistic evolution of the error of numerical method for linear stochastic differential equation

The Galerkin Analysis for the Random Periodic Solution of Semilinear Stochastic Evolution Equations

A Coupling Method of Difference with High Order Accuracy and Boundary Integral Equation for Evolutionary Equation and Its Error Estimates

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

A Reproducing Kernel Particle Method for Solving Generalized Probability Density Evolution Equation in Stochastic Dynamic Analysis

Approximation of SDEs: a stochastic sewing approach