A monotone finite volume method for time fractional Fokker-Planck equations

We develop a monotone finite volume method for the time fractional Fokker-Planck equations and theoretically prove its unconditional stability. We show that the convergence rate of this method is of order 1 in the space and if the space grid becomes suffciently fine, the convergence rate can be improved to order 2. Numerical results are given to support our theoretical findings. One characteristic of our method is that it has monotone property such that it keeps the nonnegativity of some physical variables such as density, concentration, etc.