High-resolution MRI velocimetry compared with numerical simulations

Numerical investigation of patient-specific thoracic aortic aneurysms and comparison with normal subject via computational fluid dynamics (CFD)

Vascular hemodynamics play an important role in cardiovascular diseases. This work aimed to investigate the effects of an increase in ascending aortic diameter (AAD) on hemodynamics throughout a cardiac cycle for real patients. In this study, two scans of thoracic aortic aneurysm (TAA) subject with different AADs (42.94 mm and 48.01 mm) and a scan of a normal subject (19.81 mm) were analyzed to assess the effects of hemodynamics on the progression of TAA with the same flow rate. Real-patient aortic geometries were scanned by computed tomography angiography (CTA), and steady and pulsatile flow conditions were used to simulate real patient aortic geometries. Aortic arches were obtained from routine clinical scans. Computational fluid dynamics (CFD) simulations were performed with in vivo boundary conditions, and 3D Navier-Stokes equations were solved by a UDF (user-defined function) code defining a real cardiac cycle of one patient using Fourier series (FS). Wall shear stress (WSS) and pressure distributions were presented from normal subject to TAA cases. The results show that during the peak systolic phase pressure load increased by 18.56% from normal subject to TAA case 1 and by 23.8% from normal subject to TAA case 2 in the aneurysm region. It is concluded that although overall WSS increased in aneurysm cases but was low in dilatation areas. As a result, abnormal changes in WSS and higher pressure load may lead to rupture and risk of further dilatation. CFD simulations were highly effective to guide clinical predictions and assess the progress of aneurysm regions in case of early surgical intervention. Graphical abstract.

Accuracy of phase-contrast velocimetry in systems with skewed intravoxel velocity distributions

Phase contrast velocimetry (PCV) has been widely used to investigate flow properties in numerous systems. Several authors have reported errors in velocity measurements and have speculated on the sources, which have ranged from eddy current effects to acceleration artefacts. An often overlooked assumption in the theory of PCV, which may not be met in complex or unsteady flows, is that the intravoxel displacement distributions (propagators) are symmetric. Here, the effect of the higher moments of the displacement distribution (variance, skewness and kurtosis) on the accuracy of PCV is investigated experimentally and theoretically. Phase and propagator measurements are performed on tailored intravoxel distributions, achieved using a simple phantom combined with a single large voxel. Asymmetric distributions (Skewness ≠ 0) are shown to generate important phase measurement errors that lead to significant velocimetry errors. Simulations of the phase of the spin vector sum, based on experimentally measured propagators, are shown to quantitatively reproduce the relationship between measured phase and experimental parameters. These allow relating the observed velocimetry errors to a discrepancy between the average phase of intravoxel spins considered in PCV theory and the vector phase actually measured by a PFG experiment. A theoretical expression is derived for PCV velocimetry errors as a function of the moments of the displacement distribution. Positively skewed distributions result in an underestimation of the true mean velocity, while negatively skewed distributions result in an overestimation. The magnitude of these errors is shown to increase with the variance and decrease with the kurtosis of the intravoxel displacement distribution.

Measurement with microscopic MRI and simulation of flow in different aneurysm models

PURPOSE

The impact and the development of aneurysms depend to a significant degree on the exchange of liquid between the regular vessel and the pathological extension. A better understanding of this process will lead to improved prediction capabilities. The aim of the current study was to investigate fluid-exchange in aneurysm models of different complexities by combining microscopic magnetic resonance measurements with numerical simulations. In order to evaluate the accuracy and applicability of these methods, the fluid-exchange process between the unaltered vessel lumen and the aneurysm phantoms was analyzed quantitatively using high spatial resolution.

METHODS

Magnetic resonance flow imaging was used to visualize fluid-exchange in two different models produced with a 3D printer. One model of an aneurysm was based on histological findings. The flow distribution in the different models was measured on a microscopic scale using time of flight magnetic resonance imaging. The whole experiment was simulated using fast graphics processing unit-based numerical simulations. The obtained simulation results were compared qualitatively and quantitatively with the magnetic resonance imaging measurements, taking into account flow and spin-lattice relaxation.

RESULTS

The results of both presented methods compared well for the used aneurysm models and the chosen flow distributions. The results from the fluid-exchange analysis showed comparable characteristics concerning measurement and simulation. Similar symmetry behavior was observed. Based on these results, the amount of fluid-exchange was calculated. Depending on the geometry of the models, 7% to 45% of the liquid was exchanged per second.

CONCLUSIONS

The result of the numerical simulations coincides well with the experimentally determined velocity field. The rate of fluid-exchange between vessel and aneurysm was well-predicted. Hence, the results obtained by simulation could be validated by the experiment. The observed deviations can be caused by the noise in the measurement and by the limited resolution of the simulation. The resulting differences are small enough to allow reliable predictions of the flow distribution in vessels with stents and for pulsed blood flow.