Numerical approximation of nonlinear fractional differential equations with subdiffusion and superdiffusion

详细信息    查看全文

• 作者：Changpin Li^a ; ^{lcp@shu.edu.cn"" rel=""nofollow} ; Zhengang Zhao^a ; ^{zgzhao999@yahoo.com.cn"" rel=""nofollow} ; YangQuan Chen^b ; ^{yangquan.chen@usu.edu"" rel=""nofollow}
• 关键词：Time– ; space fractional diffusion equation ; Subdiffusion ; Superdiffusion ; Riemann– ; Liouville derivative ; Caputo derivative ; Difference method ; Finite element method
• 刊名：Computers and Mathematics with Applications
• 出版年：2011
• 出版时间：August 2011
• 年：2011
• 卷：62
• 期：3
• 页码：855-875
• 全文大小：570 K

文摘

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.