In this article, two block-centered finite difference schemes are introduced and analyzed to solve the parabolic integro-differential equation arising in modeling non-Fickian flow in porous media. One ...