Quantitative analysis of PC MRI velocity maps: pulsatile flow in cylindrical vessels

The accuracy of MR phase contrast (PC) velocity mapping, and the subsequent derivation of wall shear stress (WSS) values, has been quantitatively assessed. Using a retrospectively gated PC gradient-echo technique, the temporal-spatial velocity fields were measured for pulsatile flow in a rigid cylindrical vessel. The experimental data were compared with values derived from the Womersley solution of the Navier-Stokes equations. For a sinusoidal waveform, the overall root-mean-square (rms) difference between the measured and analytical velocities corresponded to 13% of the peak fluid velocity. The WSS derived from the data displayed a 14% rms difference with the analytical model. As an example of a more complicated flow, a triangular saw-tooth waveform was deconstructed into its Fourier components. Velocity maps and the WSS were calculated by the superposition of the individual solutions, weighted by the Fourier series coefficient, for each harmonic. The velocity and experimentally derived WSS agreed with the analytical results (4% and 12% rms difference, respectively). Evaluation of the analytical models allowed an estimate of the inherent accuracy in the measurement of velocity maps and WSS values.