Summary: | A numerical method for solving fractional partial differential equations (fPDEs) of the diffusion and reaction–diffusion type, subject to Dirichlet boundary data, in three dimensions is developed. Such fPDEs may describe fluid flows through porous media better than classical diffusion equations. This is a new, fractional version of the Alternating Direction Implicit (ADI) method, where the source term is balanced, in that its effect is split in the three space directions, and it may be relevant, especially in the case of anisotropy. The method is unconditionally stable, second-order in space, and third-order in time. A strategy is devised in order to improve its speed of convergence by means of an extrapolation method that is coupled to the <i>PageRank</i> algorithm. Some numerical examples are given.
|