AbstractPurposeComputational simulation of phase‐contrast MRI (PC‐MRI) is an attractive way to physically interpret properties and errors in MRI‐reconstructed flow velocity fields. Recent studies have developed PC‐MRI simulators that solve the Bloch equation, with the magnetization transport being modeled using a Lagrangian approach. Because this method expresses the magnetization as spatial distribution of particles, influences of particle densities and their spatial uniformities on numerical accuracy are well known. This study developed an alternative method for PC‐MRI modeling using an Eulerian approach in which the magnetization is expressed as a spatially smooth continuous function.MethodsThe magnetization motion was described using the Bloch equation with an advection term and computed on a fixed grid using a finite difference method, and k‐space sampling was implemented using the spoiled gradient echo sequence. PC‐MRI scans of a fully developed flow in straight and stenosed cylinders were acquired to provide numerical examples.ResultsReconstructed flow in a straight cylinder showed excellent agreement with input velocity profiles and mean errors were less than 0.5% of the maximum velocity. Numerical cases of flow in a stenosed cylinder successfully demonstrated the velocity profiles, with displacement artifacts being dependent on scan parameters and intravoxel dephasing due to flow disturbances. These results were in good agreement with those obtained using the Lagrangian approach with a sufficient particle density.ConclusionThe feasibility of the Eulerian approach to PC‐MRI modeling was successfully demonstrated.