DPD enables mesoscopic MRI simulation of slow flow

We present a novel method to simulate magnetic resonance imaging (MRI) for the assessment of slow flow at Reynolds number \(Re \approx 0.02\). We couple Bloch equations with dissipative particle dynamics (DPD) to study the effect of flow dynamics at the mesoscopic level on acquired MR images. The Bloch equations are used to propagate the evolution of the magnetization of particles while their trajectories are being computed simultaneously based on DPD interaction forces. The magnetic resonance assessment of fluid velocities is performed using a phase-contrast MRI technique, implemented by a spin echo single-sided bipolar gradient sequence. The computational cost for simulating the fluid flow is successfully reduced by an efficient implementation of a vectorized isochromat algorithm. We demonstrate successful simulation of laminar flow, flow with diffusion effects, and flow around an obstacle. The method can be used to simulate convective and diffusive flow MRI experiments at the mesoscopic level.