Numerical approximation of nonlinear fractional differential equations with subdiffusion and superdiffusion

https://doi.org/10.1016/j.camwa.2011.02.045Get rights and content
Under an Elsevier user license
open archive

Abstract

In this paper, we study the time–space fractional order (fractional for simplicity) nonlinear subdiffusion and superdiffusion equations, which can relate the matter flux vector to concentration gradient in the general sense, describing, for example, the phenomena of anomalous diffusion, fractional Brownian motion, and so on. The semi-discrete and fully discrete numerical approximations are both analyzed, where the Galerkin finite element method for the space Riemann–Liouville fractional derivative with order 1+β[1,2] and the finite difference scheme for the time Caputo derivative with order α(0,1) (for subdiffusion) and (1,2) (for superdiffusion) are analyzed, respectively. Results on the existence and uniqueness of the weak solutions, the numerical stability, and the error estimates are presented. Numerical examples are included to confirm the theoretical analysis. During our simulations, an interesting diffusion phenomenon of particles is observed, that is, on average, the diffusion velocity for 0<α<1 is slower than that for α=1, but the diffusion velocity for 1<α<2 is faster than that for α=1. For the spatial diffusion, we have a similar observation.

Keywords

Time–space fractional diffusion equation
Subdiffusion
Superdiffusion
Riemann–Liouville derivative
Caputo derivative
Difference method
Finite element method

Cited by (0)

This work was partially supported by the Natural Science Foundation of China under Grant No. 10872119, the Key Disciplines of Shanghai Municipality under Grant No. S30104, and the Innovation Foundation of Shanghai University under Grant No. SHUCX101070.