Finite difference diagonalization to simulate nuclear magnetic resonance diffusion experiments in porous media

A finite difference approach for computing Laplacian eigenvalues and eigenvectors in discrete porous media is derived and used to approximately solve the Bloch-Torrey equations. Neumann, Dirichlet, and Robin boundary conditions are considered and applications to simulate nuclear magnetic resonance diffusion experiments are shown. The method is illustrated with MATLAB examples and computational tests in one and two dimensions and the extension to three dimensions is outlined. (c) 2015 Wiley Periodicals, Inc. Concepts Magn Reson Part A 44A: 160-180, 2015.