Two kinds of numerical algorithms for ultra-slow diffusion equations

MaRDI QIDQ6434496zbMATH OpenFDO

Authors Min Cai, Changpin Li

Publication date 27 April 2023

Mathematics Subject Classification ID

Fractional partial differential equations (35R11) Finite difference methods for initial value and initial-boundary value problems involving PDEs (65M06)

Abstract: In this article, two kinds of numerical algorithms are derived for the ultra-slow (or superslow) diffusion equation in one and two space dimensions, where the ultra-slow diffusion is characterized by the Caputo-Hadamard fractional derivative of order $alphain(0,1)$. To describe the spatial interaction, the Riesz fractional derivative and the fractional Laplacian are used in one and two space dimensions, respectively. The Caputo-Hadamard derivative is discretized by two typical approximate formulae, i.e., L2-1${}_{sigma}$ and L1-2 methods. The spatial fractional derivatives are discretized by the 2-nd order finite difference methods. When L2-1${}_{sigma}$ discretization is used, the derived numerical scheme is unconditionally stable with error estimate $mathcalO(a{u}^2+h^2)$ for all $alphain(0,1)$, in which $au$ and $h$ are temporal and spatial stepsizes, respectively. When L1-2 discretization is used, the derived numerical scheme is stable with error estimate $mathcalO(a{u}^{3-alpha}+h^2)$ for $alphain(0,0.3738)$. The illustrative examples displayed are in line with the theoretical analysis.

This page was built for publication: Two kinds of numerical algorithms for ultra-slow diffusion equations

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6434496)