Numerical Method to Modify the Fractional-Order Diffusion Equation

Time or space or time-space fractional-order diffusion equations (FODEs) are widely used to describe anomalous diffusion processes in many physical and biological systems. In recent years, many authors have proposed different numerical methods to solve the modified fractional-order diffusion equations, and some achievements have been obtained. However, to our knowledge of the literature, up to date, all the proposed numerical methods to modify FODE have achieved at most a second-order time accuracy. In this study, we focus mainly on the numerical methods based on numerical integration in order to modify the fractional-order diffusion equation: 1+1/12δx2pjk−1+1/12δx2pjk−1=μα∑l=0kλαlδx2pjk−l+ μβ∑l=0kλβlδx2pjk−l+τ/21+1/12δx2fjk−1+fjk, k=1,2,…,K;j=1,2,…,j−1,pj0=ωj,j=0,1,…,J, p0k=φtk,pjk=ψtk,k=0,1,…,K,fjl=fxj,tl,ωj=ωxj. Accordingly, numerical methods can be built to modify FODE with second-order time accuracy and fourth-order spatial accuracy in ∂px,t/∂t=∂1−α/∂t1−α+B∂1−β/∂t1−β∂2px,t/∂x2+fx,t,0