Numerical analysis and physical simulations for the time fractional radial diffusion equation

We do the numerical analysis and simulations for the time fractional radial diffusion equation used to describe the anomalous subdiffusive transport processes on the symmetric diffusive field. Based on rewriting the equation in a new form, we first present two kinds of implicit finite difference schemes for numerically solving the equation. Then we strictly establish the stability and convergence results. We prove that the two schemes are both unconditionally stable and second order convergent with respect to the maximum norm. Some numerical results are presented to confirm the rates of convergence and the robustness of the numerical schemes. Finally, we do the physical simulations. Some interesting physical phenomena are revealed; we verify that the long time asymptotic survival probability ∝t^−α, but independent of the dimension, where α is the anomalous diffusion exponent.