The second order weighted and shifted Grünwald difference (WSGD) operators are developed in [Tian, Zhou and Deng, Math. Comput., 84 (2015), pp. 1703–1727] to solve space fractional partial differential equations. Along this direction, we further design a new family of second order WSGD operators; by properly choosing the weighted parameters, they can be effectively used to discretize space (Riemann-Liouville) fractional derivatives. Based on the new second order WSGD operators, we derive a family of difference schemes for the space fractional advection diffusion equation. By von Neumann stability analysis, it is proved that the obtained schemes are unconditionally stable. Finally, extensive numerical experiments are performed to demonstrate the performance of the schemes and confirm the convergence orders.

In this paper, a new numerical algorithm for solving the time fractional Fokker-Planck equation is proposed. The analysis of local truncation error and the stability of this method are investigated. Theoretical analysis and numerical experiments show that the proposed method has higher order of accuracy for solving the time fractional Fokker-Planck equation.