Numerical Solution of Fractional Order Anomalous Subdiffusion Problems Using Radial Kernels and Transform

By coupling of radial kernels and localized Laplace transform, a numerical scheme for the approximation of time fractional anomalous subdiffusion problems is presented. The fractional order operators are well suited to handle by Laplace transform and radial kernels are also built for high dimensions. The numerical computations of inverse Laplace transform are carried out by contour integration technique. The computation can be done in parallel and no time sensitivity is involved in approximating the time fractional operator as contrary to finite differences. The proposed numerical scheme is stable and accurate.