Abstract:The numerical simulation of the diffusion MRI signal arising from complex tissue micro-structures is helpful for understanding and interpreting imaging data as well as for designing and optimizing MRI sequences. The discretization of the Bloch-Torrey equation by finite elements is a more recently developed approach for this purpose, in contrast to random walk simulations, which has a longer history. While finite elements discretization is more difficult to implement than random walk simulations, the approach benefits from a long history of theoretical and numerical developments by the mathematical and engineering communities. In particular, software packages for the automated solutions of partial differential equations using finite elements discretization, such as FEniCS, are undergoing active support and development. However, because diffusion MRI simulation is a relatively new application area, there is still a gap between the simulation needs of the MRI community and the available tools provided by finite elements software packages. In this paper, we address two potential difficulties in using FEniCS for diffusion MRI simulation. First, we simplified software installation by the use of FEniCS containers that are completely portable across multiple platforms. Second, we provide a portable simulation framework based on Python and whose code is open source. This simulation framework can be seamlessly integrated with cloud computing resources such as Google Colaboratory notebooks working on a web browser or with Google Cloud Platform with MPI parallelization. We show examples illustrating the accuracy, the computational times, and parallel computing capabilities. The framework contributes to reproducible science and open-source software in computational diffusion MRI with the hope that it will help to speed up method developments and stimulate research collaborations.

Diffusion across semi-permeable barriers: spectral properties, efficient computation, and applications

Semigroup-theoretic approach to diffusion in thin layers separated by semi-permeable membranes

Surveying Diffusion in Complex Geometries. An Essay

Spectral properties of the Bloch-Torrey operator in three dimensions

Parallel diffusion operator for magnetized plasmas with improved spectral fidelity

Portable simulation framework for diffusion MRI

Local Time Statistics and Permeable Barrier Crossing: from Poisson to Birth-Death Diffusion Equations

A probabilistic model of diffusion through a semi-permeable barrier

Numerical Methods for Linear Diffusion Equations in the Presence of an Interface

Surface-to-volume ratio with oscillating gradients

Statistical algorithms for low-frequency diffusion data: A PDE approach

Diffusion maps for embedded manifolds with boundary with applications to PDEs

Nature of barriers determine first passage times in heterogeneous media

Spectral measures and iterative bounds for effective diffusivity of steady and space-time periodic flows

Mathematical Models of Diffusion in Physiology

Numerical Recovery of the Diffusion Coefficient in Diffusion Equations from Terminal Measurement

On the Mathematical foundations of Diffusion Monte Carlo

A probabilistic representation of the solution to a 1D evolution equation in a medium with negative index

Homogenization of locally stationary diffusions with possibly degenerate diffusion matrix