Comparison of hemodynamics of intracranial aneurysms between MR fluid dynamics using 3D cine phase-contrast MRI and MR-based computational fluid dynamics

Characterization of Maximum Wall Shear Stress Points in Unruptured Cerebral Aneurysms Using Four-dimensional Flow Magnetic Resonance Imaging

BACKGROUND

Maximum wall shear stress (maxWSS) points of unruptured cerebral aneurysms (UCAs) may cause wall remodeling leading to rupture. We characterized maxWSS points and their inherent intra-aneurysmal flow structures in a sizable cohort of saccular UCAs using four-dimensional (4D) flow magnetic resonance imaging (MRI).

METHODS

After contrast administration, 50 saccular UCAs were subjected to 4D flow MRI using a 1.5 T MRI scanner. Post-processing of obtained data was performed using commercially available software. The maxWSS points and maxWSS values were evaluated. The maxWSS values were statistically compared between aneurysm groups.

RESULTS

The maxWSS point was located on the aneurysm apex in 9 (18.0%), body in 2 (4.0%), and neck in 39 (78.0%) UCAs. The inherent intra-aneurysmal flow structure of the maxWSS point was an inflow zone in 34 (68.0%) UCAs, an inflow jet in 8 (16.0%), and an impingement zone in 8 (16.0%). The maxWSS point on the neck had significantly higher maxWSS values than those points on the other wall areas (P = 0.008). The maxWSS values of the maxWSS points on the apex and on the impingement zone were not significantly different compared with those of the other maxWSS points.

CONCLUSION

The maxWSS points existed preferentially on the aneurysmal neck adjacent to the inflow zone with higher maxWSS values. The maxWSS points existed occasionally on the aneurysmal apex adjacent to the impingement zone. 4D flow MRI may be helpful to discriminate saccular UCAs with higher-risk maxWSS points that can cause wall remodeling leading to rupture.

Fluid dynamic analysis in predicting the recanalization of intracranial aneurysms after coil embolization - A study of spatiotemporal characteristics

Purpose

Hemodynamics play a key role in the management of cerebral aneurysm recanalization after coil embolization; however, the most reliable hemodynamic parameter remains unknown. Previous studies have explored the use of both spatiotemporally averaged and maximal definitions for hemodynamic parameters, based on computational fluid dynamics (CFD) analysis, to build predictive models for aneurysmal recanalization. In this study, we aimed to assess the influence of different spatiotemporal characteristics of hemodynamic parameters on predictive performance.

Methods

Hemodynamics were simulated using CFD for 66 cerebral aneurysms from 65 patients. We evaluated 14 types of spatiotemporal definitions for two hemodynamic parameters in the pre-coiling model and five in virtual post-coiling model (VM) created by cutting the aneurysm from the pre-coiling model. A total of 91 spatiotemporal hemodynamic features were derived and utilized to develop univariate predictor (UP) and multivariate logistic regression (LR) models. The model's performance was assessed using two metrics: the area under the receiver operating characteristic curve (AUROC) and the area under the precision-recall curve (AUPRC).

Results

Different spatiotemporal hemodynamic features exhibited a wide range of AUROC values ranging from 0.224 to 0.747, with 22 feature pairs showing a significant difference in AUROC value (P-value <0.05), despite being derived from the same hemodynamic parameter. PDave,q1 was identified as the strongest UP with AUROC/AUPRC values of 0.747/0.385, yielding sensitivity and specificity value of 0.889 and 0.614 at the optimal cut-off value, respectively. The LR model further improved the prediction performance, having AUROC/AUPRC values of 0.890/0.903. At the optimal cut-off value, the LR model achieved a specificity of 0.877, sensitivity of 0.719, outperforming the UP model.

Conclusion

Our research indicated that the characteristics of hemodynamic parameters in terms of space and time had a significant impact on the development of predictive model. Our findings suggest that LR model based on spatiotemporal hemodynamic features could be clinically useful in predicting recanalization after coil embolization in patients, without the need for invasive procedures.

Four-dimensional Flow MRI Assessment of Portal Hemodynamics and Hepatic Regeneration after Portal Vein Embolization

Background Percutaneous transhepatic portal vein (PV) embolization (PVE) is a standard preoperative procedure for advanced biliary cancer when the future liver remnant (FLR) is insufficient, yet the effect of this procedure on portal hemodynamics is still unclear. Purpose To assess whether four-dimensional (4D) MRI flowmetry can be used to estimate FLR volume and to identify the optimal time for this measurement. Materials and Methods This prospective single-center study enrolled consecutive adult patients with biliary cancer who underwent percutaneous transhepatic PVE for the right liver between June 2020 and November 2022. Portal hemodynamics were assessed using 4D flow MRI before PVE and within 1 day (0-day group) or 3-4 days (3-day group) after PVE. FLR volume was measured using CT before PVE and after PVE but before surgery. Blood flow changes were analyzed with the Wilcoxon signed rank test, and correlations with Spearman rank correlation. Results The 0-day group included 24 participants (median age, 72 years [IQR, 69-77 years]; 17 male participants), and the 3-day group included 13 participants (median age, 71 years [IQR, 68-78 years]; eight male participants). Both groups showed increased left PV (LPV) flow rate after PVE (0-day group: from median 3.72 mL/sec [IQR, 2.83-4.55 mL/sec] to 9.48 mL/sec [IQR, 8.12-10.7 mL/sec], P < .001; 3-day group: from median 3.65 mL/sec [IQR, 2.14-3.79 mL/sec] to 8.16 mL/sec [IQR, 6.82-8.98 mL/sec], P < .001). LPV flow change correlated with FLR volume change relative to the number of days from PVE to presurgery CT only in the 3-day group (ρ = 0.62, P = .02; 0-day group, P = .11). The output of the regression equation for estimating presurgery FLR volume correlated with CT-measured volume (ρ = 0.78; P = .002). Conclusion Four-dimensional flow MRI demonstrated increased blood flow in residual portal branches 3-4 days after PVE, offering insights for estimating presurgery FLR volume. Published under a CC BY 4.0 license. Supplemental material is available for this article. See also the editorial by Roldán-Alzate and Oechtering in this issue.

On the identification of hypoxic regions in subject-specific cerebral vasculature by combined CFD/MRI

A long-time exposure to lack of oxygen (hypoxia) in some regions of the cerebrovascular system is believed to be one of the causes of cerebral neurological diseases. In the present study, we show how a combination of magnetic resonance imaging (MRI) and computational fluid dynamics (CFD) can provide a non-invasive alternative for studying blood flow and transport of oxygen within the cerebral vasculature. We perform computer simulations of oxygen mass transfer in the subject-specific geometry of the circle of Willis. The computational domain and boundary conditions are based on four-dimensional (4D)-flow MRI measurements. Two different oxygen mass transfer models are considered: passive (where oxygen is treated as a dilute chemical species in plasma) and active (where oxygen is bonded to haemoglobin) models. We show that neglecting haemoglobin transport results in a significant underestimation of the arterial wall mass transfer of oxygen. We identified the hypoxic regions along the arterial walls by introducing the critical thresholds that are obtained by comparison of the estimated range of Damköhler number (Da ⊂ 〈9; 57〉) with the local Sherwood number. Finally, we recommend additional validations of the combined MRI/CFD approach proposed here for larger groups of subject- or patient-specific brain vasculature systems.

Distribution and regional variation of wall shear stress in the curved middle cerebral artery using four-dimensional flow magnetic resonance imaging

BACKGROUND

To investigate the distribution and regional variation of wall shear stress (WSS) in the curved middle cerebral artery (MCA) in healthy individuals using four-dimensional (4D) flow magnetic resonance imaging (MRI).

METHODS

A total of 44 healthy participants (18 males; mean ages: 27.16±5.69 years) were included in this cross-sectional study. The WSS parameters of mean, minimum, and maximum values, the coefficient of variation of time-averaged WSS (TAWSSCV), and the maximum values of the oscillatory shear index (OSI) were calculated and compared in the curved proximal (M1) segments. Three cross-sectional planes were selected: the location perpendicular to the beginning of the long axis of the curved M1 segment of the MCA (proximal section), the most curved M1 location (curved M1 section), and the location before the insular (M2) segment bifurcation (distal section). The WSS and OSI parameters of the proximal, curved, and distal sections of the curved M1 segment were compared, including the inner and outer curvatures of the curved M1 section.

RESULTS

Of the curved M1 segments, the curved M1 section had significantly lower minimum TAWSS values than the proximal (P=0.031) and distal sections (P=0.002), and the curved M1 section had significantly higher maximum OSI values than the distal section (P=0.001). The TAWSSCV values at the curved M1 section were significantly higher than the proximal (P=0.001) and distal sections (P<0.001). At the curved M1 section, the inner curvature showed a significantly lower minimum TAWSS (P=0.013) and higher maximum OSI values (P=0.002) than the outer curvature.

CONCLUSIONS

There are distribution variation of WSS and OSI parameters at the curved M1 section of the curved MCA, and the inner curvature of the curved M1 section has the lowest WSS and highest OSI distribution. The local hemodynamic features of the curved MCA may be related to the predilection for atherosclerotic plaque development.

Minimum wall shear stress points and their underlying intra-aneurysmal flow structures of unruptured cerebral aneurysms on 4D flow MRI

BACKGROUND AND PURPOSE

Minimum wall shear stress (Min-WSS) points may be associated with wall instability of unruptured cerebral aneurysms. We aimed to investigate the relationship between the locations of Min-WSS points and their underlying intra-aneurysmal flow structure patterns in unruptured cerebral aneurysms using four-dimensional (4D) flow magnetic resonance imaging (MRI).

MATERIALS AND METHODS

Min-WSS points and the intra-aneurysmal flow structure patterns were identified in 50 unruptured aneurysms by 4D flow MRI.

RESULTS

The Min-WSS points were located around a vortex core tip in 31 (62.0%) aneurysms and on an intra-bleb vortex center in 7 (14.0%). Sixteen (32.0%) aneurysms had the Min-WSS points on the aneurysmal apex, and in 24 (48.0%) were on the neck. The Min-WSS values of aneurysms with the Min-WSS points on an intra-bleb flow were significantly lower than those of the other groups (P = 0.030). Aneurysms with the Min-WSS points on the neck had significantly higher Min-WSS values than the other aneurysms (P = 0.008).

CONCLUSIONS

The location of the Min-WSS point was corresponding to the vortex core or center in 76% of all aneurysms. The underlying intra-aneurysmal flow structure and location of the Min-WSS point affect the Min-WSS value. Further studies are needed to characterize Min-WSS points to identify aneurysms with a higher risk of wall instability.

Quality Control for 4D Flow MR Imaging

In recent years, 4D flow MRI has become increasingly important in clinical applications for the blood vessels in the whole body, heart, and cerebrospinal fluid. 4D flow MRI has advantages over 2D cine phase-contrast (PC) MRI in that any targeted area of interest can be analyzed post-hoc, but there are some factors to be considered, such as ensuring measurement accuracy, a long imaging time and post-processing complexity, and interobserver variability.Due to the partial volume phenomenon caused by low spatial and temporal resolutions, the accuracy of flow measurement in 4D flow MRI is reduced. For spatial resolution, it is recommended to include at least four voxels in the vessel of interest, and if possible, six voxels. In large vessels such as the aorta, large voxels can be secured and SNR can be maintained, but in small cerebral vessels, SNR is reduced, resulting in reduced accuracy. A temporal resolution of less than 40 ms is recommended. The velocity-to-noise ratio (VNR) of low-velocity blood flow is low, resulting in poor measurement accuracy. The use of dual velocity encoding (VENC) or multi-VENC is recommended to avoid velocity wrap around and to increase VNR. In order to maintain sufficient spatio-temporal resolution, a longer imaging time is required, leading to potential patient movement during examination and a corresponding decrease in measurement accuracy.For the clinical application of new technologies, including various acceleration techniques, in vitro and in vivo accuracy verification based on existing accuracy-validated 2D cine PC MRI and 4D flow MRI, as well as accuracy verification on the conservation of mass' principle, should be performed, and intraobserver repeatability, interobserver reproducibility, and test-retest reproducibility should be checked.

Does the Pulsatile Non-uniform Flow Matter in MR Flowmetry?

3D cine phase-contrast (4D flow) MRI is a sequence with great potential for non-invasive time-resolved 3D flowmetry at arbitrary vessel sections in various blood vessels. However, it is not widely known that the flowmetry with 4D flow MRI is vulnerable to pulsatile and non-uniform flow. Due to the limited spatial and temporal resolutions, averaging within the 3D voxel is occurring during the flowmetry. A simple solution is to avoid setting the measurement plane in the area where non-uniform flow is dominant, which is possible with an aid of streamline depictions generated by computational fluid dynamics (CFD) or 4D flow MRI data. Unlike 4D flow MRI, flowmetry in CFD simulation can use higher spatial and temporal resolution depending on computer performance; therefore, it is robust to fluctuating non-uniform flow. However, the performance of CFD simulations might be limited due to inlet conditions with low temporal resolution. Difficulty applying complex blood flow such as reflection flow from periphery may also limit accurate simulation. Caution should be taken when comparing the result of CFD simulation to that of 4D flow measurement.

4D flow MRI hemodynamic biomarkers for cerebrovascular diseases

Alterations in cerebral blood flow are common in several neurological diseases among the elderly including stroke, cerebral small vessel disease, vascular dementia, and Alzheimer's disease. 4D flow magnetic resonance imaging (MRI) is a relatively new technique to investigate cerebrovascular disease, and makes it possible to obtain time-resolved blood flow measurements of the entire cerebral arterial venous vasculature and can be used to derive a repertoire of hemodynamic biomarkers indicative of cerebrovascular health. The information that can be obtained from one single 4D flow MRI scan allows both the investigation of aberrant flow patterns at a focal location in the vasculature as well as estimations of brain-wide disturbances in blood flow. Such focal and global hemodynamic biomarkers show the potential of being sensitive to impending cerebrovascular disease and disease progression and can also become useful during planning and follow-up of interventions aiming to restore a normal cerebral circulation. Here, we describe 4D flow MRI approaches for analyzing the cerebral vasculature. We then survey key hemodynamic biomarkers that can be reliably assessed using the technique. Finally, we highlight cerebrovascular diseases where one or multiple hemodynamic biomarkers are of central interest.

Hemodynamic Study of a Patient-Specific Intracranial Aneurysm: Comparative Assessment of Tomographic PIV, Stereoscopic PIV, In Vivo MRI and Computational Fluid Dynamics

Introduction

Wall shear stress (WSS) is associated with the growth and rupture of an intracranial aneurysm. To reveal their underlying connections, many image-based computational fluid dynamics (CFD) studies have been conducted. However, the methodological validations using both in vivo medical imaging and in vitro optical flow measurements were rarely accompanied in such studies.

Methods

In the present study, we performed a comparative assessment on the hemodynamics of a patient-specific intracranial saccular aneurysm using in vivo 4D Flow MRI, in silico CFD, in vitro stereoscopic and tomographic particle imaging velocimetry (Stereo-PIV and Tomo-PIV) techniques. PIV experiments and CFD were conducted under steady state corresponding to the peak systole of 4D Flow MRI.

Results

The results showed that all modalities provided similar flow features and overall surface distribution of WSS. However, a large variation in the absolute WSS values was found. 4D Flow MRI estimated a 2- to 4-fold lower peak WSS (3.99 Pa) and a 1.6- to 2-fold lower mean WSS (0.94 Pa) than Tomo-PIV, Stereo-PIV, and CFD. Bland-Altman plots of WSS showed that the differences between PIV-/CFD-based WSS and 4D Flow MRI-based WSS increase with higher WSS magnitude. Such proportional trend was absent in the Bland-Altman comparison of velocity where the resolutions of PIV and CFD datasets were matched to 4D Flow MRI. We also found that because of superior resolution in the out-of-plane direction, WSS estimation by Tomo-PIV was higher than Stereo-PIV.

Conclusions

Our results indicated that the differences in spatial resolution could be the main contributor to the discrepancies between each modality. The findings of this study suggest that with current techniques, care should be taken when using absolute WSS values to perform a quantitative risk analysis of aneurysm rupture.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13239-021-00583-2.

Improving Blood Flow Visualization of Recirculation Regions at Carotid Bulb in 4D Flow MRI Using Semi-Automatic Segmentation with ITK-SNAP

Assessment of carotid bulb hemodynamics using four-dimensional (4D) flow magnetic resonance imaging (MRI) requires accurate segmentation of recirculation regions that is frequently hampered by limited resolution. This study aims to improve the accuracy of 4D flow MRI carotid bulb segmentation and subsequent recirculation regions analysis. Time-of-flight (TOF) MRI and 4D flow MRI were performed on bilateral carotid artery bifurcations in seven healthy volunteers. TOF-MRI data was segmented into 3D geometry for computational fluid dynamics (CFD) simulations. ITK-SNAP segmentation software was included in the workflow for the semi-automatic generation of 4D flow MRI angiographic data. This study compared the velocities calculated at the carotid bifurcations and the 3D blood flow visualization at the carotid bulbs obtained by 4D flow MRI and CFD. By applying ITK-SNAP segmentation software, an obvious improvement in the 4D flow MRI visualization of the recirculation regions was observed. The 4D flow MRI images of the recirculation flow characteristics of the carotid artery bulbs coincided with the CFD. A reasonable agreement was found in terms of velocity calculated at the carotid bifurcation between CFD and 4D flow MRI. However, the dispersion of velocity data points relative to the local errors of measurement in 4D flow MRI remains. Our proposed strategy showed the feasibility of improving recirculation regions segmentation and the potential for reliable blood flow visualization in 4D flow MRI. However, quantitative analysis of recirculation regions in 4D flow MRI with ITK-SNAP should be enhanced for use in clinical situations.

Comparison of hemodynamic stress in healthy vessels after parent artery occlusion and flow diverter stent treatment for internal carotid artery aneurysm

OBJECTIVE

De novo aneurysms generally develop in healthy vessels after parent artery occlusion for large internal carotid artery (ICA) aneurysm, possibly owing to increased hemodynamic stress in the remaining vessels. In recent years, there has been a shift toward flow diverter stent treatment. However, there is a lack of direct evidence and data that prove this change in hemodynamic stress in healthy vessels after parent artery occlusion and flow diverter stent treatment. The authors compared hemodynamic stress in healthy-side vessels before and after parent artery occlusion and flow diverter treatments.

METHODS

The authors included patients who underwent 3D cine phase-contrast MRI before and after large ICA aneurysm treatment. Spatially and temporally averaged volume flow rates and spatially averaged systolic wall shear stress (WSS) in healthy-side ICA distal to the posterior communicating artery (C1 segment according to Fisher's classification) were measured before and after parent artery occlusion and flow diverter treatments.

RESULTS

Seventeen patients were included (5 patients in the parent artery occlusion group and 12 in the flow diverter group). At 1-2 months after treatment, median volume flow rate in healthy-side ICA increased from 5.36 ml/sec to 6.28 ml/sec (total increase 117%, p = 0.04) in the parent artery occlusion group and from 4.65 ml/sec to 4.93 ml/sec (total increase 106%, p = 0.02) in the flow diverter group. In the parent artery occlusion group, median WSS in the C1 segment of the healthy-side ICA increased from 3.91 Pa to 5.61 Pa (total increase 143%, p = 0.08); however, no significant increase was observed in the flow diverter group (4.29 Pa to 4.57 Pa [total increase 107%, p = 0.21]).

CONCLUSIONS

Postoperatively, volume flow rate and WSS in the C1 segment of the healthy-side ICA significantly increased in the parent artery occlusion group. Therefore, the parent artery occlusion group was more prone to de novo aneurysm than the flow diverter group.

Human upper-airway respiratory airflow: In vivo comparison of computational fluid dynamics simulations and hyperpolarized 129Xe phase contrast MRI velocimetry

Computational fluid dynamics (CFD) simulations of respiratory airflow have the potential to change the clinical assessment of regional airway function in health and disease, in pulmonary medicine and otolaryngology. For example, in diseases where multiple sites of airway obstruction occur, such as obstructive sleep apnea (OSA), CFD simulations can identify which sites of obstruction contribute most to airway resistance and may therefore be candidate sites for airway surgery. The main barrier to clinical uptake of respiratory CFD to date has been the difficulty in validating CFD results against a clinical gold standard. Invasive instrumentation of the upper airway to measure respiratory airflow velocity or pressure can disrupt the airflow and alter the subject's natural breathing patterns. Therefore, in this study, we instead propose phase contrast (PC) velocimetry magnetic resonance imaging (MRI) of inhaled hyperpolarized 129Xe gas as a non-invasive reference to which airflow velocities calculated via CFD can be compared. To that end, we performed subject-specific CFD simulations in airway models derived from 1H MRI, and using respiratory flowrate measurements acquired synchronously with MRI. Airflow velocity vectors calculated by CFD simulations were then qualitatively and quantitatively compared to velocity maps derived from PC velocimetry MRI of inhaled hyperpolarized 129Xe gas. The results show both techniques produce similar spatial distributions of high velocity regions in the anterior-posterior and foot-head directions, indicating good qualitative agreement. Statistically significant correlations and low Bland-Altman bias between the local velocity values produced by the two techniques indicates quantitative agreement. This preliminary in vivo comparison of respiratory airway CFD and PC MRI of hyperpolarized 129Xe gas demonstrates the feasibility of PC MRI as a technique to validate respiratory CFD and forms the basis for further comprehensive validation studies. This study is therefore a first step in the pathway towards clinical adoption of respiratory CFD.

Parametric Sequential Method for MRI-Based Wall Shear Stress Quantification

Wall shear stress (WSS) has been suggested as a potential biomarker in various cardiovascular diseases and it can be estimated from phase-contrast Magnetic Resonance Imaging (PC-MRI) velocity measurements. We present a parametric sequential method for MRI-based WSS quantification consisting of a geometry identification and a subsequent approximation of the velocity field. This work focuses on its validation, investigating well controlled high-resolution in vitro measurements of turbulent stationary flows and physiological pulsatile flows in phantoms. Initial tests for in vivo 2D PC-MRI data of the ascending aorta of three volunteers demonstrate basic applicability of the method to in vivo.

Comparison of Hemodynamic Visualization in Cerebral Arteries: Can Magnetic Resonance Imaging Replace Computational Fluid Dynamics?

A multimodality approach was applied using four-dimensional flow magnetic resonance imaging (4D flow MRI), time-of-flight magnetic resonance angiography (TOF-MRA) signal intensity gradient (SIG), and computational fluid dynamics (CFD) to investigate the 3D blood flow characteristics and wall shear stress (WSS) of the cerebral arteries. TOF-MRA and 4D flow MRI were performed on the major cerebral arteries in 16 healthy volunteers (mean age 34.7 ± 7.6 years). The flow rate measured with 4D flow MRI in the internal carotid artery, middle cerebral artery, and anterior cerebral artery were 3.8, 2.5, and 1.2 mL/s, respectively. The 3D blood flow pattern obtained through CFD and 4D flow MRI on the cerebral arteries showed reasonable consensus. CFD delivered much greater resolution than 4D flow MRI. TOF-MRA SIG and CFD WSS of the major cerebral arteries showed reasonable consensus with the locations where the WSS was relatively high. However, the visualizations were very different between TOF-MRA SIG and CFD WSS at the internal carotid artery bifurcations, the anterior cerebral arteries, and the anterior communicating arteries. 4D flow MRI, TOF-MRA SIG, and CFD are complementary methods that can provide additional insight into the hemodynamics of the human cerebral artery.

Inflow Hemodynamics of Intracranial Aneurysms: A Comparison of Computational Fluid Dynamics and 4D Flow Magnetic Resonance Imaging

PURPOSE

Although the inflow hemodynamics of cerebral aneurysms are key factors in their rupture and recurrence after endovascular treatments, the most available method for inflow hemodynamics evaluation remains unestablished. We compared the efficacy of inflow hemodynamics evaluation using computational fluid dynamics (CFD) analysis and that using four-dimensional (4D) flow magnetic resonance imaging (MRI).

METHODS

In 23 unruptured cerebral aneurysms, the inflow hemodynamics was evaluated using both CFD and 4D flow MRI. The evaluated parameters included visually classified inflow jet patterns, the inflow rate ratio (the ratio of the inflow rate at the aneurysmal orifice to the flow rate in the proximal parent artery), and the velocity ratio (the ratio of the inflow velocity to the velocity in the proximal parent artery). The Shapiro-Wilk test was used to assess the normality of variable data, and logarithmic transformation was performed for variables with non-normal distributions. Data analysis was performed using Pearson correlation analyses and the chi-square test.

RESULTS

There was a significant correlation between inflow jet patterns evaluated by CFD and 4D flow MRI (p = 0.008). Moreover, there was a strong correlation between the inflow rate ratios evaluated by CFD and 4D flow MRI (r = 0.801; p <0.001). Furthermore, there was a moderate correlation between the velocity ratios measured by CFD and 4D flow MRI (r = 0.559; p = 0.008).

CONCLUSION

Inflow hemodynamics evaluated by CFD analysis and 4D flow MRI showed good correlations in inflow jet pattern, inflow rate ratio, and velocity ratio.

A Transparent Vessel-on-a-Chip Device for Hemodynamic Analysis and Early Diagnosis of Intracranial Aneurysms by CFD and PC-MRI

Hemodynamics plays a critical role in early diagnosis and investigating the growth mechanism of intracranial aneurysms (IAs), which usually induce hemorrhagic stroke, serious neurological diseases, and even death. We developed a transparent blood vessel-on-a-chip (VOC) device for magnetic resonance imaging (MRI) to provide characteristic flow fields of early IAs as the reference for early diagnosis. This VOC device takes advantage of the transparent property to clearly exhibit the internal structure and identify the needless air bubbles in the biomimetic fluid experiment, which significantly affects the MRI image quality. Furthermore, the device was miniaturized and easily assembled with arbitrary direction using a 3D-printed scaffold in a radiofrequency coil. Computational fluid dynamics (CFD) simulations of the flow field were greatly consistent with those data from MRI. Both internal flow and wall shear stress (WSS) exhibited very low levels during the IA growth, thus leading to the growth and rupture of IAs. PC-MRI images can also provide a reasonable basis for the early diagnosis of IAs. Therefore, we believed that this proposed VOC-based MR imaging technique has great potential for early diagnostic of intracranial aneurysms.

Assessing the Risk of Intracranial Aneurysm Rupture Using Morphological and Hemodynamic Biomarkers Evaluated from Magnetic Resonance Fluid Dynamics and Computational Fluid Dynamics

PURPOSE

Evaluate in vivo hemodynamic and morphological biomarkers of intracranial aneurysms, using magnetic resonance fluid dynamics (MRFD) and MR-based patient specific computational fluid dynamics (CFD) in order to assess the risk of rupture.

METHODS

Forty-eight intracranial aneurysms (10 ruptured, 38 unruptured) were scrutinized for six morphological and 10 hemodynamic biomarkers. Morphological biomarkers were calculated based on 3D time-of-flight magnetic resonance angiography (3D TOF MRA) in MRFD analysis. Hemodynamic biomarkers were assessed using both MRFD and CFD analyses. MRFD was performed using 3D TOF MRA and 3D cine phase-contrast magnetic resonance imaging (3D cine PC MRI). CFD was performed utilizing patient specific inflow-outflow boundary conditions derived from 3D cine PC MRI. Univariate analysis was carried out to identify statistically significant biomarkers for aneurysm rupture and receiver operating characteristic (ROC) analysis was performed for the significant biomarkers. Binary logistic regression was performed to identify independent predictive biomarkers.

RESULTS

Morphological biomarker analysis revealed that aneurysm size [P = 0.021], volume [P = 0.035] and size ratio [P = 0.039] were statistically significantly different between the two groups. In hemodynamic biomarker analysis, MRFD results indicated that ruptured aneurysms had higher oscillatory shear index (OSI) [OSI.max, P = 0.037] and higher relative residence time (RRT) [RRT.ave, P = 0.035] compared with unruptured aneurysms. Correspondingly CFD analysis demonstrated significant differences for both average and maximum OSI [OSI.ave, P = 0.008; OSI.max, P = 0.01] and maximum RRT [RRT.max, P = 0.045]. ROC analysis revealed AUC values greater than 0.7 for all significant biomarkers. Aneurysm volume [AUC, 0.718; 95% CI, 0.491-0.946] and average OSI obtained from CFD [AUC, 0.774; 95% CI, 0.586-0.961] were retained in the respective logistic regression models.

CONCLUSION

Both morphological and hemodynamic biomarkers have significant influence on intracranial aneurysm rupture. Aneurysm size, volume, size ratio, OSI and RRT could be potential biomarkers to assess aneurysm rupture risk.

From 2D to 4D Phase-Contrast MRI in the Neurovascular System: Will It Be a Quantum Jump or a Fancy Decoration?

Considering the crosstalk between the flow and vessel wall, hemodynamic assessment of the neurovascular system may offer a well-integrated solution for both diagnosis and management by adding prognostic significance to the standard CT/MR angiography. 4D flow MRI or time-resolved 3D velocity-encoded phase-contrast MRI has long been promising for the hemodynamic evaluation of the great vessels, but challenged in clinical studies for assessing intracranial vessels with small diameter due to long scan times and low spatiotemporal resolution. Current accelerated MRI techniques, including parallel imaging with compressed sensing and radial k-space undersampling acquisitions, have decreased scan times dramatically while preserving spatial resolution. 4D flow MRI visualized and measured 3D complex flow of neurovascular diseases such as aneurysm, arteriovenous shunts, and atherosclerotic stenosis using parameters including flow volume, velocity vector, pressure gradients, and wall shear stress. In addition to the noninvasiveness of the phase contrast technique and retrospective flow measurement through the wanted windows of the analysis plane, 4D flow MRI has shown several advantages over Doppler ultrasound or computational fluid dynamics. The evaluation of the flow status and vessel wall can be performed simultaneously in the same imaging modality. This article is an overview of the recent advances in neurovascular 4D flow MRI techniques and their potential clinical applications in neurovascular disease. LEVEL OF EVIDENCE: 5 TECHNICAL EFFICACY STAGE: 3.

Four-Dimensional Magnetic Resonance Imaging Assessment of Intracranial Aneurysms: A State-of-the-Art Review

Treatment of unruptured intracranial aneurysms can reduce the risk of subarachnoid hemorrhage and its associated morbidity and mortality. However, current methods to predict the risk of rupture and optimize treatment strategies for intracranial aneurysms are limited. Assessment of intra-aneurysmal flow using 4-dimensional magnetic resonance imaging (4D MRI) is a novel tool that could be used to guide therapy. A systematic search of the literature was performed to provide a state-of-the-art review on 4D MRI assessment of unruptured intracranial aneurysms. A total of 18 studies were available for review. Eccentric flow on 4D MRI is associated with a greater aspect ratio and peak wall shear stress (WSS). WSS, vorticity, and peak velocity are greater in saccular than fusiform aneurysms. Unstable aneurysms are associated with greater WSS, peak wall stress, and flow jet angle and may exhibit wall enhancement. In comparison to computational fluid dynamics (CFD), 4D MRI has a lower spatial resolution and reports lower WSS and velocity magnitudes, but these parameters equalize when spatial resolution is matched. 4D MRI demonstrates the intra-aneurysmal hemodynamic changes associated with flow diversion, including significantly decreased flow velocity. Thus, 4D MRI is a novel, noninvasive imaging tool used for the evaluation of hemodynamics within intracranial aneurysms. Hemodynamic indices derived from 4D MRI appear to correlate well with the simulated (CFD) values and may be used to measure the success of endovascular therapies and risk factors for aneurysm growth and rupture.

Abnormal Flow Dynamics Result in Low Wall Shear Stress and High Oscillatory Shear Index in Abdominal Aortic Dilatation: Initial in vivo Assessment with 4D-flow MRI

PURPOSE

To characterize the non-laminar flow dynamics and resultant decreased wall shear stress (WSS) and high oscillatory shear index (OSI) of the infrarenal abdominal aortic dilatation, cardiac phase-resolved 3D phase-contrast MRI (4D-flow MRI) was performed.

METHODS

The prospective single-arm study was approved by the Institutional Review Board and included 18 subjects (median 67.5 years) with the dilated infrarenal aorta (median diameter 35 mm). 4D-flow MRI was conducted on a 1.5T MRI system. On 3D streamline images, laminar and non-laminar (i.e., vortex or helical) flow patterns were visually assessed both for the dilated aorta and for the undilated upstream aorta. Cardiac phase-resolved flow velocities, WSS and OSI, were also measured for the dilated aorta and the upstream undilated aorta.

RESULTS

Non-laminar flow represented by vortex or helical flow was more frequent and overt in the dilated aorta than in the undilated upstream aorta (P < 0.0156) with a very good interobserver agreement (weighted kappa: 0.82-1.0). The WSS was lower, and the OSI was higher on the dilated aortic wall compared with the proximal undilated segments. In mid-systole, mean spatially-averaged WSS was 0.20 ± 0.016 Pa for the dilated aorta vs. 0.68 ± 0.071 Pa for undilated upstream aorta (P < 0.0001), and OSI on the dilated aortic wall was 0.093 ± 0.010 vs. 0.041 ± 0.0089 (P = 0.013). The maximum values and the amplitudes of the WSS at the dilated aorta were inversely proportional to the ratio of dilated/undilated aortic diameter (r = -0.694, P = 0.0014).

CONCLUSION

4D-flow can characterize abnormal non-laminar flow dynamics within the dilated aorta in vivo. The wall of the infrarenal aortic dilatation is continuously and increasingly affected by atherogenic stimuli due to the flow disturbances represented by vortex or helical flow, which is reflected by lower WSS and higher OSI.

Haemodynamics in a patient-specific intracranial aneurysm according to experimental and numerical approaches: A comparison of PIV, CFD and PC-MRI

BACKGROUND

The haemodynamics determined by different approaches for studying fluid dynamics - i.e. computational fluid dynamics (CFD), particle image velocimetry (PIV), and phase-contrast magnetic resonance imaging (PC-MRI) - have rarely been thoroughly compared; nor have the factors that affect accuracy and precision in each method. As each method has its own advantages and limitations, this knowledge is important for future studies to be able to achieve valid analyses of fluid flows.

OBJECTIVE

To gauge the capacity of these methods for analysing aneurysmal flows, we compared the haemodynamic behaviours determined by each method within a patient-specific aneurysm model.

METHODS

An in vitro silicone aneurysm model was fabricated for PIV and PC-MRI, and an in silico aneurysm model with the same geometry was reconstructed for CFD. With the same fluid model prepared numerically and physically, CFD, PIV and PC-MRI were performed to study aneurysmal haemodynamics.

RESULTS

2D velocity vectors and magnitudes show good agreement between PIV and CFD, and 3D flow patterns show good similarity between PC-MRI and CFD.

CONCLUSIONS

These findings give confidence to future haemodynamic studies using CFD technology. For the first time, the morphological inconsistency between the PCMRI model and others is found to affect the measurement of local flow patterns.

4D Flow with MRI

Magnetic resonance imaging (MRI) has become an important tool for the clinical evaluation of patients with cardiac and vascular diseases. Since its introduction in the late 1980s, quantitative flow imaging with MRI has become a routine part of standard-of-care cardiothoracic and vascular MRI for the assessment of pathological changes in blood flow in patients with cardiovascular disease. More recently, time-resolved flow imaging with velocity encoding along all three flow directions and three-dimensional (3D) anatomic coverage (4D flow MRI) has been developed and applied to enable comprehensive 3D visualization and quantification of hemodynamics throughout the human circulatory system. This article provides an overview of the use of 4D flow applications in different cardiac and vascular regions in the human circulatory system, with a focus on using 4D flow MRI in cardiothoracic and cerebrovascular diseases.

Effect of Ambient Temperature Changes on Blood Flow in Anterior Cerebral Artery of Patients with Skull Prosthesis

BACKGROUND

In many cases, an injury to the head leads to the replacement of a part of the skull with materials such as titanium and polyether ether ketone.

METHODS

Three-dimensional heads model of 15 healthy individuals and 13 patients were prepared. The models were simulated using thermal fluid structure interaction analysis to evaluate the effects of cold (5°C) and hot (55°C) temperatures of the skull on the conditions of blood flow in the anterior cerebral artery.

RESULTS

The results showed negligible changes (<3%) in wall shear stress (WSS) vessel and von Mises stress between the healthy individuals and patients both at 25°C and 55°C. However, at 5°C, the values of these 2 parameters in the patients were 2.1 and 2.5 times those in healthy individuals, respectively. The value of WSS in healthy individuals and the patients in cold temperature was 1.2 and 2.9 times those at normal temperature. The corresponding values for von Mises stress were 1.1 and 2.2, respectively. Accordingly, the stress changes between cold and hot ambient temperatures were found to be negligible in all samples.

CONCLUSIONS

The changes in stress were significant only for the patients when exposed to cold ambient temperature, and only in patients, exposure to a cold ambient temperature significantly increased the risks of vascular aneurysm and damage to the brain tissue surrounding the blood vessels. These risks were found to be negligible for both healthy individuals and patients when exposed to hot ambient temperature and also for healthy individuals exposed to cold ambient temperature.

Four-Dimensional Flow Magnetic Resonance Imaging for Assessment of Velocity Magnitudes and Flow Patterns in The Human Carotid Artery Bifurcation: Comparison with Computational Fluid Dynamics

Purpose: Knowledge of the hemodynamics in the vascular system is important to understand and treat vascular pathology. The present study aimed to evaluate the hemodynamics in the human carotid artery bifurcation measured by four-dimensional (4D) flow magnetic resonance imaging (MRI) as compared to computational fluid dynamics (CFD). Methods: This protocol used MRI data of 12 healthy volunteers for the 3D vascular models and 4D flow MRI measurements for the boundary conditions in CFD simulation. We compared the velocities measured at the carotid bifurcation and the 3D velocity streamlines of the carotid arteries obtained by these two methods. Results: There was a good agreement for both maximum and minimum velocity values between the 2 methods for velocity magnitude at the bifurcation plane. However, on the 3D blood flow visualization, secondary flows, and recirculation regions are of poorer quality when visualized through the 4D flow MRI. Conclusion: 4D flow MRI and CFD show reasonable agreement in demonstrated velocity magnitudes at the carotid artery bifurcation. However, the visualization of blood flow at the recirculation regions and the assessment of secondary flow characteristics should be enhanced for the use of 4D flow MRI in clinical situations.

Accuracy of the Flow Velocity and Three-directional Velocity Profile Measured with Three-dimensional Cine Phase-contrast MR Imaging: Verification on Scanners from Different Manufacturers

Purpose:

The accuracy of flow velocity and three-directional velocity components are important for the precise visualization of hemodynamics by 3D cine phase-contrast MRI (3D cine PC MRI, also referred to as 4D-flow). The aim of this study was to verify the accuracy of these measurements of prototype or commercially available 3D cine PC MRI obtained by three different manufactures’ MR scanners.

Methods:

The verification of the accuracy of flow velocity in 3D cine PC MRI was performed by circulating blood mimicking fluid through a straight-tube phantom in a slanting position, such that the three-directional velocity components were simultaneously measurable, using three 3T MR scanners from different manufacturers. The data obtained were processed by phase correction, and the velocity and three-directional velocity components in the center of the tube on the central cross section of a slab were calculated. The velocity profile in each three directions and the composite velocity profiles were compared with the calculated reference values, using the Hagen–Poiseuille equation. In addition, velocity profiles and the spatially time-averaged velocity perpendicular to the tube were compared with the theoretical values and measured values by a flowmeter, respectively.

Results:

An underestimation of the maximum velocity in the center of the tube and an overestimation of the velocity near the tube wall due to partial volume effects were observed in all three scanners. A roughening and flattening of profiles in the center of the tube were observed in one scanner, due, presumably, to the low signal-to-noise ratio. However, the spatially time-averaged velocities corresponded well with the measured values by the flowmeter in all three scanners.

Conclusion:

In this study, we have demonstrated that the accuracy of flow velocity and three-directional velocity components in 3D cine PC MRI was satisfactory in all three MR scanners.

Assessment of boundary conditions for CFD simulation in human carotid artery

Computational fluid dynamics (CFD) is an increasingly used method for investigation of hemodynamic parameters and their alterations under pathological conditions, which are important indicators for diagnosis of cardiovascular disease. In hemodynamic simulation models, the employment of appropriate boundary conditions (BCs) determines the computational accuracy of the CFD simulation in comparison with pressure and velocity measurements. In this study, we have first assessed the influence of inlet boundary conditions on hemodynamic CFD simulations. We selected two typical patients suspected of carotid artery disease, with mild stenosis and severe stenosis. Both patients underwent digital subtraction angiography (DSA), magnetic resonance angiography, and the invasive pressure guide wire measured pressure profile. We have performed computational experiments to (1) study the hemodynamic simulation outcomes of distributions of wall shear stress, pressure, pressure gradient and (2) determine the differences in hemodynamic performances caused by inlet BCs derived from DSA and Womersley analytical solution. Our study has found that the difference is related to the severity of the stenosis; the greater the stenosis, the more the difference ensues. Further, in our study, the two typical subjects with invasively measured pressure profile and thirty subjects with ultrasound Doppler velocimeter (UDV) measurement served as the criteria to evaluate the hemodynamic outcomes of wall shear stress, pressure, pressure gradient and velocity due to different outlet BCs based on the Windkessel model, structured-tree model, and fully developed flow model. According to the pressure profiles, the fully developed model appeared to have more fluctuations compared with the other two models. The Windkessel model had more singularities before convergence. The three outlet BCs models also showed good correlation with the UDV measurement, while the Windkessel model appeared to be slightly better ([Formula: see text]). The structured-tree model was seen to have the best performance in terms of available computational cost and accuracy. The results of our numerical simulation and the good correlation with the computed pressure and velocity with their measurements have highlighted the effectiveness of CFD simulation in patient-specific human carotid artery with suspected stenosis.

A comparison of 4D flow MRI-derived wall shear stress with computational fluid dynamics methods for intracranial aneurysms and carotid bifurcations - A review

BACKGROUND

4D flow MRI is a relatively quick method for obtaining wall shear stress (WSS) in vivo, a hemodynamic parameter which has shown promise in risk stratification for rupture of cerebrovascular diseases such as intracranial aneurysms and atherosclerotic plaques. The accuracy of such measurements is still largely unknown.

OBJECTIVE

To quantify the accuracy of 4D flow MRI-derived wall shear stress values for intracranial aneurysms and carotid bifurcations.

METHOD

We performed a review of all original research articles which compared the magnitudes of WSS derived from 4D flow MRI with corresponding values derived from computational fluid dynamics (CFD) within both intracranial aneurysms and carotid bifurcations.

RESULT

For intracranial aneurysms and carotid bifurcations, 4D flow MRI-derived WSS estimations are generally lower in magnitude compared to WSS derived by CFD methods. These differences are more pronounced in regions of higher WSS. However, the relative distributions of WSS derived from both methods are reasonably similar.

CONCLUSION

Pooled analysis suggests that WSS magnitudes obtained by 4D flow MRI are underestimated, while the relative distribution is reasonably accurate, the latter being an important factor for determining the natural history of intracranial aneurysms and other cerebrovascular diseases. 4D flow MRI shows enormous potential in providing new risk stratification parameters which could have significant impact on individualized treatment decisions and improved patient outcomes.

Transitional hemodynamics in intracranial aneurysms - Comparative velocity investigations with high resolution lattice Boltzmann simulations, normal resolution ANSYS simulations, and MR imaging

PURPOSE

Blood flow in intracranial aneurysms has, until recently, been considered to be disturbed but still laminar. Recent high resolution computational studies have demonstrated, in some situations, however, that the flow may exhibit high frequency fluctuations that resemble weakly turbulent or transitional flow. Due to numerous assumptions required for simplification in computational fluid dynamics (CFD) studies, the occurrence of these events, in vivo, remains unsettled. The detection of these fluctuations in aneurysmal blood flow, i.e., hemodynamics by CFD, poses additional challenges as such phenomena cannot be captured in clinical data acquisition with magnetic resonance (MR) due to inadequate temporal and spatial resolutions. The authors' purpose was to address this issue by comparing results from highly resolved simulations, conventional resolution laminar simulations, and MR measurements, identify the differences, and identify their causes.

METHODS

Two aneurysms in the basilar artery, one with disturbed yet laminar flow and the other with transitional flow, were chosen. One set of highly resolved direct numerical simulations using the lattice Boltzmann method (LBM) and another with adequate resolutions under laminar flow assumption were conducted using a commercially available ANSYS Fluent solver. The velocity fields obtained from simulation results were qualitatively and statistically compared against each other and with MR acquisition.

RESULTS

Results from LBM, ANSYS Fluent, and MR agree well qualitatively and quantitatively for one of the aneurysms with laminar flow in which fluctuations were <80 Hz. The comparisons for the second aneurysm with high fluctuations of > ∼ 600 Hz showed vivid differences between LBM, ANSYS Fluent, and magnetic resonance imaging. After ensemble averaging and down-sampling to coarser space and time scales, these differences became minimal.

CONCLUSIONS

A combination of MR derived data and CFD can be helpful in estimating the hemodynamic environment of intracranial aneurysms. Adequately resolved CFD would suffice gross assessment of hemodynamics, potentially in a clinical setting, and highly resolved CFD could be helpful in a detailed and retrospective understanding of the physiological mechanisms.

Influence of Spatial Resolution in Three-dimensional Cine Phase Contrast Magnetic Resonance Imaging on the Accuracy of Hemodynamic Analysis

Introduction:

We aim to elucidate the effect of spatial resolution of three-dimensional cine phase contrast magnetic resonance (3D cine PC MR) imaging on the accuracy of the blood flow analysis, and examine the optimal setting for spatial resolution using flow phantoms.

Materials and Methods:

The flow phantom has five types of acrylic pipes that represent human blood vessels (inner diameters: 15, 12, 9, 6, and 3 mm). The pipes were fixed with 1% agarose containing 0.025 mol/L gadolinium contrast agent. A blood-mimicking fluid with human blood property values was circulated through the pipes at a steady flow. Magnetic resonance (MR) images (three-directional phase images with speed information and magnitude images for information of shape) were acquired using the 3-Tesla MR system and receiving coil. Temporal changes in spatially-averaged velocity and maximum velocity were calculated using hemodynamic analysis software. We calculated the error rates of the flow velocities based on the volume flow rates measured with a flowmeter and examined measurement accuracy.

Results:

When the acrylic pipe was the size of the thoracicoabdominal or cervical artery and the ratio of pixel size for the pipe was set at 30% or lower, spatially-averaged velocity measurements were highly accurate. When the pixel size ratio was set at 10% or lower, maximum velocity could be measured with high accuracy. It was difficult to accurately measure maximum velocity of the 3-mm pipe, which was the size of an intracranial major artery, but the error for spatially-averaged velocity was 20% or less.

Conclusions:

Flow velocity measurement accuracy of 3D cine PC MR imaging for pipes with inner sizes equivalent to vessels in the cervical and thoracicoabdominal arteries is good. The flow velocity accuracy for the pipe with a 3-mm-diameter that is equivalent to major intracranial arteries is poor for maximum velocity, but it is relatively good for spatially-averaged velocity.

4D Flow MRI in Neuroradiology: Techniques and Applications

Assessment of the intracranial flow is important for the understanding and management of cerebral vascular diseases. From brain aneurysms and arteriovenous malformations lesions to intracranial and cervical stenosis, the appraisal of the blood flow can be crucial and influence positively on patients' management. The determination of the intracranial hemodynamics and the collateral pattern seems to play to a major role in the management of these lesions. 4D flow magnetic resonance imaging is a noninvasive phase contrast derived method that has been developed and applied in neurovascular diseases. It has a great potential if followed by further technical improvements and comprehensive and systematic clinical studies.

Patient-specific assessment of hemodynamics by computational fluid dynamics in patients with bicuspid aortopathy

OBJECTIVE

Hemodynamics related to eccentric blood flow may factor into the development of bicuspid aortic valve aortopathy. We investigated wall shear stress distribution by means of magnetic resonance imaging-based computational fluid dynamics in patients with a bicuspid aortic valve.

METHODS

Included were 12 patients with a bicuspid aortic valve (aortic stenosis, n = 11; root enlargement, n = 1). Three patients with a normal tricuspid aortic valve (arch aneurysm, n = 1; descending aortic aneurysm, n = 2) were included for comparison. The thoracic aorta geometry was reconstructed by means of 3-dimensional computed tomography angiography, and the bicuspid aortic valve orifice was modeled. Flow rates at the sinotubular junction and 3 aortic branches were measured at various time points by cine phase-contrast magnetic resonance imaging to define boundary conditions for computational fluid dynamics, and the flow was simulated.

RESULTS

Bicuspid aortic valve cusp configurations were type 0 lateral (n = 4), type 0 anterior-posterior (n = 2), type 1 L-R (n = 4), and type 1 R-N (n = 2). Abnormal aortic helical flow was seen in the ascending aorta and transverse arch in all patients with bicuspid aortic valves and was right handed in 11 patients (91%). No such flow was seen in the patients with tricuspid aortic valves. The patients with bicuspid aortic valves were likely to have jet flow/wall impingement against the greater curvature of the proximal ascending aorta, resulting in remarkably increased wall shear stress around the impingement area.

CONCLUSIONS

Computational fluid dynamics simulation is useful for precise evaluation of hemodynamics related to bicuspid aortic valve aortopathy. Such evaluation will advance our understanding of the disease pathophysiology and may facilitate molecular biological investigation.

Characterizing saccular aortic arch aneurysms from the geometry-flow dynamics relationship

OBJECTIVE

Low wall shear stress (WSS) has been reported to be associated with accelerated atherosclerosis, aneurysm growth, or rupture. We evaluated the geometry of aortic arch aneurysms and their relationship with WSS by using the 4-dimensional flow magnetic resonance imaging to better characterize the saccular aneurysms.

METHODS

We analyzed the geometry in 100 patients using multiplanar reconstruction of computed tomography. We evaluated WSS and vortex flow using 4-dimensional flow magnetic resonance imaging in 16 of them, which were compared with 8 age-matched control subjects and eight healthy young volunteers.

RESULTS

Eighty-two patients had a saccular aneurysm, and 18 had a fusiform aneurysm. External diameter/aneurysm length ratio and sac depth/neck width ratio of the fusiform aneurysms were constant at 0.76 ± 0.18 and 0.23 ± 0.09, whereas those of saccular aneurysms, especially those involving the outer curvature, were higher and more variable. Vortex flow was always present in the aneurysms, resulting in low WSS. When the sac depth/neck width ratio was less than 0.8, peak WSS correlated inversely with luminal diameter even in the saccular aneurysms. When this ratio exceeded 0.8, which was the case only with the saccular aneurysms, such correlation no longer existed and WSS was invariably low.

CONCLUSIONS

Fusiform aneurysms elongate as they dilate, and WSS is lower as the diameter is larger. Saccular aneurysms dilate without proportionate elongation, and they, especially those occupying the inner curvature, have higher and variable sac depth/neck width ratio. When this ratio exceeds 0.8, WSS is low regardless of diameter, which may explain their malignant clinical behavior.

Overcoming spatio-temporal limitations using dynamically scaled in vitro PC-MRI - A flow field comparison to true-scale computer simulations of idealized, stented and patient-specific left main bifurcations

The majority of patients with angina or heart failure have coronary artery disease. Left main bifurcations are particularly susceptible to pathological narrowing. Flow is a major factor of atheroma development, but limitations in imaging technology such as spatio-temporal resolution, signal-to-noise ratio (SNRv), and imaging artefacts prevent in vivo investigations. Computational fluid dynamics (CFD) modelling is a common numerical approach to study flow, but it requires a cautious and rigorous application for meaningful results. Left main bifurcation angles of 40°, 80° and 110° were found to represent the spread of an atlas based 100 computed tomography angiograms. Three left mains with these bifurcation angles were reconstructed with 1) idealized, 2) stented, and 3) patient-specific geometry. These were then approximately 7× scaled-up and 3D printing as large phantoms. Their flow was reproduced using a blood-analogous, dynamically scaled steady flow circuit, enabling in vitro phase-contrast magnetic resonance (PC-MRI) measurements. After threshold segmentation the image data was registered to true-scale CFD of the same coronary geometry using a coherent point drift algorithm, yielding a small covariance error (σ² <;5.8×10^-4). Natural-neighbour interpolation of the CFD data onto the PC-MRI grid enabled direct flow field comparison, showing very good agreement in magnitude (error 2-12%) and directional changes (r² 0.87-0.91), and stent induced flow alternations were measureable for the first time. PC-MRI over-estimated velocities close to the wall, possibly due to partial voluming. Bifurcation shape determined the development of slow flow regions, which created lower SNRv regions and increased discrepancies. These can likely be minimised in future by testing different similarity parameters to reduce acquisition error and improve correlation further. It was demonstrated that in vitro large phantom acquisition correlates to true-scale coronary flow simulations when dynamically scaled, and thus can overcome current PC-MRI's spatio-temporal limitations. This novel method enables experimental assessment of stent induced flow alternations, and in future may elevate CFD coronary flow simulations by providing sophisticated boundary conditions, and enable investigations of stenosis phantoms.

Non-invasive evaluation of pulmonary arterial blood flow and wall shear stress in pulmonary arterial hypertension with 3D phase contrast magnetic resonance imaging

Background

Recently, time-resolved 3D phase contrast magnetic resonance imaging (4D-flow) allows flow dynamics in patients with pulmonary arterial hypertension to be measured. Abnormal flow dynamics, such as vortex blood flow pattern in the pulmonary artery (PA), may reflect progression of pulmonary arterial hypertension (PAH). Some reports suggested that abnormal blood flow parameters including wall shear stress (WSS) could be markers of PAH. However, it was not fully assessed clinical usefulness of these variables. We aimed to assess whether these flow dynamic parameters, such as vortex formation time (VFT) and WSS, were associated with right ventricular (RV) function.

Results

Fifteen subjects, nine with PAH and six healthy volunteers, underwent 4D-flow. Differences of Blood flow patterns, blood flow velocities and WSS between PAH patients and healthy volunteers were evaluated. We also assessed the association between VFT, WSS and RV function in PAH patients. Both vortex blood flow patterns and early systolic retrograde flow in the main PA were observed in all patients with PAH. The PA flow velocities and WSS in patients with PAH were lower than those in healthy volunteers, but that blood flow volumes in the MPA, RPA and LPA and SV in the MPA were broadly comparable between the groups. The mean VFT was 35.0 ± 16.6 % of the cardiac cycle. The VFT significantly correlated with RV ejection fraction, RV end systolic volume, and RV end systolic volume index (RVEF = 75.1 + (−85.7)·VFT, p = 0.003, RVESV = 12.4 + 181.8·VFT, p = 0.037 and RVESVI = 10.6 + 114.8·VFT, p = 0.038, respectively) in PAH patients, whereas WSS did not correlate with RV function.

Conclusions

We confirmed that abnormal blood flow dynamics, including the vortex formation and the early onset of retrograde flow, low WSS in the PA were characteristics of PAH. The VFT may be associated with right ventricular dysfunction, whereas WSS was not. Our results suggest that 4D-flow is an effective means of detecting right heart failure as well as diagnosing PAH.

Clinical trial registration URL: https://upload.umin.ac.jp/cgi-open-bin/ctr/ctr.cgi. Unique identifier: UMIN000011128

Electronic supplementary material

The online version of this article (doi:10.1186/s40064-016-2755-7) contains supplementary material, which is available to authorized users.

Advanced flow MRI: emerging techniques and applications

Magnetic resonance imaging (MRI) techniques provide non-invasive and non-ionising methods for the highly accurate anatomical depiction of the heart and vessels throughout the cardiac cycle. In addition, the intrinsic sensitivity of MRI to motion offers the unique ability to acquire spatially registered blood flow simultaneously with the morphological data, within a single measurement. In clinical routine, flow MRI is typically accomplished using methods that resolve two spatial dimensions in individual planes and encode the time-resolved velocity in one principal direction, typically oriented perpendicular to the two-dimensional (2D) section. This review describes recently developed advanced MRI flow techniques, which allow for more comprehensive evaluation of blood flow characteristics, such as real-time flow imaging, 2D multiple-venc phase contrast MRI, four-dimensional (4D) flow MRI, quantification of complex haemodynamic properties, and highly accelerated flow imaging. Emerging techniques and novel applications are explored. In addition, applications of these new techniques for the improved evaluation of cardiovascular (aorta, pulmonary arteries, congenital heart disease, atrial fibrillation, coronary arteries) as well as cerebrovascular disease (intra-cranial arteries and veins) are presented.

Inflow Jet Patterns of Unruptured Cerebral Aneurysms Based on the Flow Velocity in the Parent Artery: Evaluation Using 4D Flow MRI

BACKGROUND AND PURPOSE

Inflow jet characteristics may be related to aneurysmal bleb formation and rupture. We investigated the visualization threshold on the basis of the flow velocity in the parent artery to classify the inflow jet patterns observed on 4D flow MR imaging.

MATERIALS AND METHODS

Fifty-seven unruptured aneurysms (24 bifurcation and 33 sidewall aneurysms) were subjected to 4D flow MR imaging to visualize inflow streamline bundles whose velocity exceeded visualization thresholds corresponding to 60%, 75%, and 90% of the maximum flow velocity in the parent artery. The shape of the streamline bundle was determined visually, and the inflow jet patterns were classified as concentrated, diffuse, neck-limited, and unvisualized.

RESULTS

At the 75% threshold, bifurcation aneurysms exhibited a concentrated inflow jet pattern at the highest rate. At this threshold, the inflow jets were concentrated in 13 aneurysms (group C, 22.8%), diffuse in 18 (group D, 31.6%), neck-limited in 11 (group N, 19.3%), and unvisualized in 15 (group U, 26.3%). In 16 (28.1%) of the 57 aneurysms, the inflow jet pattern was different at various thresholds. Most inflow parameters, including the maximum inflow velocity and rate, the inflow velocity ratio, and the inflow rate ratio, were significantly higher in groups C and D than in groups N and U.

CONCLUSIONS

The inflow jet pattern may depend on the threshold applied to visualize the inflow streamlines on 4D flow MR imaging. For the classification of the inflow jet patterns on 4D flow MR imaging, the 75% threshold may be optimal among the 3 thresholds corresponding to 60%, 75%, and 90% of the maximum flow velocity in the parent artery.

Relation between wall shear stress and carotid artery wall thickening MRI versus CFD

Wall shear stress (WSS), a parameter associated with endothelial function, is calculated by computational fluid dynamics (CFD) or phase-contrast (PC) MRI measurements. Although CFD is common in WSS (WSSCFD) calculations, PC-MRI-based WSS (WSSMRI) is more favorable in population studies; since it is straightforward and less time consuming. However, it is not clear if WSSMRI and WSSCFD show similar associations with vascular pathology. Our aim was to test the associations between wall thickness (WT) of the carotid arteries and WSSMRI and WSSCFD. The subjects (n=14) with an asymptomatic carotid plaque who underwent MRI scans two times within 4 years of time were selected from the Rotterdam Study. We compared WSSCFD and WSSMRI at baseline and follow-up. Baseline WSSMRI and WSSCFD values were divided into 3 categories representing low, medium and high WSS tertiles. WT of each tertile was compared by a one-way ANOVA test. The WSSMRI and WSSCFD were 0.50±0.13Pa and 0.73±0.25Pa at baseline. Although WSSMRI was underestimated, a significant regression was found between WSSMRI and WSSCFD (r(2)=0.71). No significant difference was found between baseline and follow-up WSS by CFD and MRI-based calculations. The WT at baseline was 1.36±0.16mm and did not change over time. The WT was 1.55±0.21mm in low, 1.33±0.20mm in medium and 1.21±0.21mm in the high WSSMRI tertiles. Similarly, the WT was 1.49±0.21mm in low, 1.33±0.20mm in medium and 1.26±0.21mm in high WSSCFD tertiles. We found that WSSMRI and WSSCFD were inversely related with WT. WSSMRI and WSSCFD patterns were similar although MRI-based calculations underestimated WSS.

Velocity Measurement in Carotid Artery: Quantitative Comparison of Time-Resolved 3D Phase-Contrast MRI and Image-based Computational Fluid Dynamics

Background:

Understanding hemodynamic environment in vessels is important for realizing the mechanisms leading to vascular pathologies.

Objectives:

Three-dimensional velocity vector field in carotid bifurcation is visualized using TR 3D phase-contrast magnetic resonance imaging (TR 3D PC MRI) and computational fluid dynamics (CFD). This study aimed to present a qualitative and quantitative comparison of the velocity vector field obtained by each technique.

Subjects and Methods:

MR imaging was performed on a 30-year old male normal subject. TR 3D PC MRI was performed on a 3 T scanner to measure velocity in carotid bifurcation. 3D anatomical model for CFD was created using images obtained from time-of-flight MR angiography. Velocity vector field in carotid bifurcation was predicted using CFD and PC MRI techniques. A statistical analysis was performed to assess the agreement between the two methods.

Results:

Although the main flow patterns were the same for the both techniques, CFD showed a greater resolution in mapping the secondary and circulating flows. Overall root mean square (RMS) errors for all the corresponding data points in PC MRI and CFD were 14.27% in peak systole and 12.91% in end diastole relative to maximum velocity measured at each cardiac phase. Bland-Altman plots showed a very good agreement between the two techniques. However, this study was not aimed to validate any of methods, instead, the consistency was assessed to accentuate the similarities and differences between Time-resolved PC MRI and CFD.

Conclusion:

Both techniques provided quantitatively consistent results of in vivo velocity vector fields in right internal carotid artery (RCA). PC MRI represented a good estimation of main flow patterns inside the vasculature, which seems to be acceptable for clinical use. However, limitations of each technique should be considered while interpreting results.

Reproducibility and interobserver variability of systolic blood flow velocity and 3D wall shear stress derived from 4D flow MRI in the healthy aorta

PURPOSE

To investigate the reproducibility and interobserver variability of 3D aortic velocity vector fields and wall shear stress (WSS) averaged over five systolic timeframes derived from noncontrast 4D flow magnetic resonance imaging (MRI).

MATERIALS AND METHODS

Fourteen controls underwent test-retest 4D flow MRI examinations separated by 16 ± 3 days (resolution = 3.0-3.6 × 2.3-2.6 × 2.5-2.7 mm(3) ; TE/TR/FA = 2.5/4.9 msec/7°; Venc = 150 cm/s). Two observers segmented the aorta, and WSS was calculated for both series of scans and both segmentations. Test-retest and interobserver velocity and WSS vectors were compared on a voxel-by-voxel basis in the aorta and on a regional basis by subdividing the aortas in six segments.

RESULTS

Test-retest: voxel-by-voxel Bland-Altman analysis revealed small differences (-0.03/-0.02 m/s/Pa), limits of agreement (LOA) of 0.25 m/s/0.29 Pa, and coefficients of variation (CV) of 20% for velocity/WSS. Voxel-by-voxel orthogonal regression analysis showed moderate agreement (slope: 1.14/1.16, intraclass correlation coefficient [ICC]: 0.76/0.67 for velocity/WSS). The regional analysis revealed a CV of 9%/8% and ICC of 0.9/0.9 for velocity/WSS. Interobserver: voxel-by-voxel difference for WSS was 0, LOA: 0.17/0.19 Pa, CV: 12/13%, slope: 1.01/1.09, ICC: 0.87/0.85 for test/retest. The CV/ICC for WSS in the regional analysis was 4%/1.0 for test and 3%/1.0 for retest.

CONCLUSION

Systolic velocity and WSS derived from 4D flow MRI are reproducible between consecutive visits, with low interobserver variability in healthy volunteers.

The predictive value of magnetic resonance imaging in evaluating intracranial arteriovenous malformation obliteration after stereotactic radiosurgery

OBJECT

The current gold standard for diagnosing arteriovenous malformation (AVM) and assessing its obliteration after stereotactic radiosurgery (SRS) is digital subtraction angiography (DSA). Recently, MRI and MR angiography (MRA) have become increasingly popular imaging modalities for the follow-up of patients with an AVM because of their convenient setup and noninvasiveness. In this study, the authors assessed the sensitivity and specificity of MRI/MRA in evaluating AVM nidus obliteration as assessed by DSA.

METHODS

The authors study a consecutive series of 136 patients who underwent SRS between January 2000 and December 2012 and who underwent regular clinical examinations, several MRI studies, and at least 1 post-SRS DSA follow- up evaluation at the University of Virginia. The average follow-up time was 47.3 months (range 10.1–165.2 months). Two blinded observers were enrolled to interpret the results of MRI/MRA compared with those of DSA. The sensitivity, specificity, positive predictive value, and negative predictive value for the obliteration of AVM were reported.

RESULTS

On the basis of DSA, 73 patients (53.7%) achieved final angiographic obliteration in a median of 28.8 months. The sensitivity (the probability of finding obliteration on MRI/MRA among those for whom complete obliteration was shown on DSA) was 84.9% for one observer (Observer 1) and 76.7% for the other (Observer 2). The specificity was 88.9% and 95.2%, respectively. The false-negative interpretations were significantly related to the presence of draining veins, perinidal edema on T2-weighted images, and the interval between the MRI/MRA and DSA studies.

CONCLUSIONS

MRI/MRA predicted AVM obliteration after SRS in most patients and can be used in their follow-up. However, because the specificity of MRI/MRA is not perfect, DSA should still be performed to confirm AVM nidus obliteration after SRS.

4D flow imaging: current status to future clinical applications

4D flow MRI permits a comprehensive in-vivo assessment of three-directional blood flow within 3-dimensional vascular structures throughout the cardiac cycle. Given the large coverage permitted from a 4D flow acquisition, the distribution of vessel wall and flow parameters along an entire vessel of interest can thus be derived from a single measurement without being dependent on multiple predefined 2D acquisitions. In addition to qualitative 3D visualizations of complex cardiac and vascular flow patterns, quantitative flow analysis can be performed and is complemented by the ability to compute sophisticated hemodynamic parameters, such as wall shear stress or 3D pressure difference maps. These metrics can provide information previously unavailable with conventional modalities regarding the impact of cardiovascular disease or therapy on global and regional changes in hemodynamics. This review provides an introduction to the methodological aspects of 4D flow MRI to assess vascular hemodynamics and describes its potential for the assessment and understanding of altered hemodynamics in the presence of cardiovascular disease.

A practical introduction to the hemodynamic analysis of the cardiovascular system with 4D Flow MRI

The 4D Flow MRI technique provides a three-dimensional representation of blood flow over time, making it possible to evaluate the hemodynamics of the cardiovascular system both qualitatively and quantitatively. In this article, we describe the application of the 4D Flow technique in a 3T scanner; in addition to the technical parameters, we discuss the advantages and limitations of the technique and its possible clinical applications. We used 4D Flow MRI to study different body areas (chest, abdomen, neck, and head) in 10 volunteers. We obtained 3D representations of the patterns of flow and quantitative hemodynamic measurements. The technique makes it possible to evaluate the pattern of blood flow in large and midsize vessels without the need for exogenous contrast agents.

Flow visualization of recurrent aneurysms after coil embolization by 3D phase-contrast MRI

BACKGROUND

Flow patterns in cerebral aneurysms are clinically important. Information on inflow patterns into aneurysms is especially helpful in preventing a recurrence after coil embolization. Computational fluid dynamics (CFD) simulations of patient-specific cerebral aneurysms are feasible and provide information on flow patterns. However, flow visualization by CFD simulations is challenging for recurrent aneurysms after coil embolization because coils make it difficult to obtain precise geometry of the recurrent aneurysms. In this study, we assessed the feasibility of flow visualization of recurrent aneurysms using 3D phase-contrast magnetic resonance imaging (PC-MRI).

METHOD

Time-of-flight magnetic resonance angiography and 3D PC-MRI were performed in eight cases of recurrent aneurysms after coil embolization. We attempted to visualize flow inside the aneurysms using data of 3D PC-MRI and evaluated the visualization. Additionally, CFD simulations were performed in a single case.

RESULTS

Inflow into aneurysms was visualized in all eight cases (100%). Flow patterns inside aneurysms were visualized in six cases (75%), and these were associated with a large size of recurrent aneurysms (mean size, 10.3 mm for visualized cases vs. 4.8 mm for unvisualized cases; p = 0.046, Mann-Whitney test). Flow patterns were similar between PC-MRI and CFD simulations. PC-MRI was faster and easier for observing inflow patterns than CFD simulations.

CONCLUSIONS

This is the first study to demonstrate that flow visualization of recurrent aneurysms by 3D PC-MRI is feasible. This technique may be more practical and easier than CFD simulations, and may provide clinically helpful information.