Improving computation of cardiovascular relative pressure fields from velocity MRI

To evaluate a multigrid-based solver for the pressure Poisson equation (PPE) with Galerkin coarsening, which works directly on the specified domain, for the computation of relative pressure fields from velocity MRI data.

We compared the proposed structure-defined Poisson solver to other popular Poisson solvers working on unmodified rectangular and modified quasirectangular domains using synthetic and in vitro phantoms in which the mathematical solution of the pressure field is known, as well as on in vivo MRI velocity measurements of aortic blood flow dynamics.

All three PPE solvers gave accurate results for convex computational domains. Using a rectangular or quasirectangular domain on a more complicated domain, like a c-shape, revealed a systematic underestimation of the pressure amplitudes, while the proposed PPE solver, working directly on the specified domain, provided accurate estimates of the relative pressure fields.

Popular iterative approaches with quasirectangular computational domains can lead to significant systematic underestimation of the pressure amplitude. We suggest using a multigrid-based PPE solver with Galerkin coarsening, which works directly on the structure-defined computational domain. This solver provides accurate estimates of the relative pressure fields for both simple and complex geometries with additional significant improvements with respect to execution speed.