Wall shear stress calculations based on 3D cine phase contrast MRI and computational fluid dynamics: a comparison study in healthy carotid arteries

Regional Aortic Wall Shear Stress Increases over Time in Patients with a Bicuspid Aortic Valve

BACKGROUND

Aortic wall shear stress (WSS) is a known predictor of ascending aortic growth in patients with a bicuspid aortic valve (BAV). The aim of this study was to study regional WSS and changes over time in BAV patients.

METHODS

BAV patients and age-matched healthy controls underwent 4D flow CMR. Regional, peak systolic ascending aortic WSS, aortic valve function, aortic stiffness measures and aortic dimensions were assessed. In BAV patients, 4D flow CMR was repeated after three years follow-up and both at baseline and follow-up computed tomography angiography (CTA) was acquired. Aortic growth (volume increase of ≥5%) was measured on CTA. Regional WSS differences within patients' aorta and WSS changes over time were analysed using linear mixed-effect models and were associated with clinical parameters.

RESULTS

Thirty BAV patients (aged 34 years [IQR 25-41]) were included in the follow-up analysis. Additionally, another 16 BAV patients and 32 healthy controls (aged 33 years [IQR 28-48]) were included for other regional analyses. Magnitude, axial, and circumferential WSS increased over time (all p<0.001) irrespective of aortic growth. The percentage of regions exposed to a magnitude WSS >95th percentile of healthy controls increased from 21% (baseline 506/2400 regions) to 31% (follow-up 734/2400 regions) (p<0.001). WSS angle, a measure of helicity near the aortic wall, decreased during follow-up. Magnitude WSS changes over time were associated with systolic blood pressure, peak aortic valve velocity, aortic valve regurgitation fraction, aortic stiffness indexes, and normalized flow displacement (all p<0.05).

CONCLUSIONS

An increase of regional WSS over time was observed in BAV patients, irrespective of aortic growth. The increasing WSSs comprising a larger area of the aorta warrants further research to investigate the possible predictive value for aortic dissection.

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.

Assessment of Complex Flow Patterns in Patients With Carotid Webs, Patients With Carotid Atherosclerosis, and Healthy Subjects Using 4D Flow MRI

BACKGROUND

Carotid webs (CaWs) are fibromuscular projections in the internal carotid artery (ICA) that cause mild luminal narrowing (<50%), but may be causative in up to one-third of seemingly cryptogenic strokes. Understanding hemodynamic alterations caused by CaWs is imperative to assessing stroke risk. Time-Average Wall Shear Stress (TAWSS) and Oscillatory Shear Index (OSI) are hemodynamic parameters linked to vascular dysfunction and thrombosis.

PURPOSE

To test the hypothesis: "CaWs are associated with lower TAWSS and higher OSI than mild atherosclerosis or healthy carotid bifurcation."

STUDY TYPE

Prospective study.

POPULATION

A total of 35 subjects (N = 14 bifurcations with CaW, 11F, age: 49 ± 10, 10 mild atherosclerosis 6F, age: 72 ± 9, 11 healthy 9F, age: 42 ± 13).

FIELD STRENGTH/SEQUENCE

4D flow/STAR-MATCH/3D TOF/3T MRI, CTA.

ASSESSMENT

4D Flow velocity data were analyzed in two ways: 1) 3D ROI in the ICA bulbar segment (complex flow patterns are expected) was used to quantify the regions with low TAWSS and high OSI. 2) 2D planes were placed perpendicular to the centerline of the carotid bifurcation for detailed analysis of TAWSS and OSI.

STATISTICAL TESTS

Independent-samples Kruskal-Wallis-H test with 0.05 used for statistical significance.

RESULTS

The percent surface area where low TAWSS was present in the ICA bulb was 12.3 ± 8.0% (95% CI: 7.6-16.9) in CaW subjects, 1.6 ± 1.9% (95% CI: 0.2-2.9) in atherosclerosis, and 8.5 ± 7.7% (95% CI: 3.6-13.4) in healthy subjects, all differences were statistically significant (ƞ2 = 0.3 [95% CI: 0.05-0.5], P-value CaW vs. healthy = 0.2). OSI had similar values in the CCA between groups (ƞ2 = 0.07 [95% CI: 0.0-0.2], P-value = 0.5), but OSI was significantly higher downstream of the bifurcation in CaW subjects compared to atherosclerosis and normal subjects. OSI returned to similar values between groups 1.5 diameters distal to the bifurcation (ƞ2 = 0.03 [95% CI: 0.0-0.2], P-value = 0.7).

CONCLUSION

Lower TAWSS and higher OSI are present in the ICA bulb in patients with CaW when compared to patients with atherosclerotic or healthy subjects.

EVIDENCE LEVEL

2 TECHNICAL EFFICACY: Stage 2.

Subjects with carotid webs demonstrate pro-thrombotic hemodynamics compared to subjects with carotid atherosclerosis

Carotid artery webs (CaW) are non-atherosclerotic projections into the vascular lumen and have been linked to up to one-third of cryptogenic strokes in younger patients. Determining how CaW affects local hemodynamics is essential for understanding clot formation and stroke risk. Computational fluid dynamics simulations were used to investigate patient-specific hemodynamics in carotid artery bifurcations with CaW, bifurcations with atherosclerotic lesions having a similar degree of lumen narrowing, and with healthy carotid bifurcations. Simulations were conducted using segmented computed tomography angiography geometries with inlet boundary conditions extracted from 2D phase contrast MRI scans. The study included carotid bifurcations with CaW (n = 13), mild atherosclerosis (n = 7), and healthy bifurcation geometries (n = 6). Hemodynamic parameters associated with vascular dysfunction and clot formation, including shear rate, oscillatory shear index (OSI), low velocity, and flow stasis were calculated and compared between the subject groups. Patients with CaW had significantly larger regions containing low shear rate, high OSI, low velocity, and flow stasis in comparison to subjects with mild atherosclerosis or normal bifurcations. These abnormal hemodynamic metrics in patients with CaW are associated with clot formation and vascular dysfunction and suggest that hemodynamic assessment may be a tool to assess stroke risk in these patients.

Abnormal flow pattern of low wall shear stress and high oscillatory shear index in spontaneous vertebral artery dissection with vertebral artery hypoplasia

Introduction

Spontaneous vertebral artery dissection (sVAD) might tend to develop in vertebral artery hypoplasia (VAH) with hemodynamic dysfunction and it is crucial to assess hemodynamics in sVAD with VAH to investigate this hypothesis. This retrospective study aimed to quantify hemodynamic parameters in patients with sVAD with VAH.

Methods

Patients who had suffered ischemic stroke due to an sVAD of VAH were enrolled in this retrospective study. The geometries of 14 patients (28 vessels) were reconstructed using Mimics and Geomagic Studio software from CT angiography (CTA). ANSYS ICEM and ANSYS FLUENT were utilized for mesh generation, set boundary conditions, solve governing equations, and perform numerical simulations. Slices were obtained at the upstream area, dissection or midstream area and downstream area of each VA. The blood flow patterns were visualized through instantaneous streamline and pressure at peak systole and late diastole. The hemodynamic parameters included pressure, velocity, time-averaged blood flow, time-averaged wall shear stress (TAWSS), oscillatory shear index (OSI), endothelial cell action potential (ECAP), relative residence time (RRT) and time-averaged nitric oxide production rate (TAR_NO).

Results

Significant focal increased velocity was present in the dissection area of steno-occlusive sVAD with VAH compared to other nondissected areas (0.910 m/s vs. 0.449 vs. 0.566, p < 0.001), while focal slow flow velocity was observed in the dissection area of aneurysmal dilatative sVAD with VAH according to velocity streamlines. Steno-occlusive sVAD with VAH arteries had a lower time-averaged blood flow (0.499 cm³/s vs. 2.268, p < 0.001), lower TAWSS (1.115 Pa vs. 2.437, p = 0.001), higher OSI (0.248 vs. 0.173, p = 0.006), higher ECAP (0.328 Pa^-1 vs. 0.094, p = 0.002), higher RRT (3.519 Pa^-1 vs. 1.044, p = 0.001) and deceased TAR_NO (104.014 nM/s vs. 158.195, p < 0.001) than the contralateral VAs.

Conclusion

Steno-occlusive sVAD with VAH patients had abnormal blood flow patterns of focal increased velocity, low time-averaged blood flow, low TAWSS, high OSI, high ECAP, high RRT and decreased TAR_NO. These results provide a good basis for further investigation of sVAD hemodynamics and support the applicability of the CFD method in testing the hemodynamic hypothesis of sVAD. More detailed hemodynamic conditions with different stages of sVAD are warranted in the future.

Advanced risk prediction for aortic dissection patients using imaging-based computational flow analysis

Patients with either a repaired or medically managed aortic dissection have varying degrees of risk of developing late complications. High-risk patients would benefit from earlier intervention to improve their long-term survival. Currently serial imaging is used for risk stratification, which is not always reliable. On the other hand, understanding aortic haemodynamics within a dissection is essential to fully evaluate the disease and predict how it may progress. In recent decades, computational fluid dynamics (CFD) has been extensively applied to simulate complex haemodynamics within aortic diseases, and more recently, four-dimensional (4D)-flow magnetic resonance imaging (MRI) techniques have been developed for in vivo haemodynamic measurement. This paper presents a comprehensive review on the application of image-based CFD simulations and 4D-flow MRI analysis for risk prediction in aortic dissection. The key steps involved in patient-specific CFD analyses are demonstrated. Finally, we propose a workflow incorporating computational modelling for personalised assessment to aid in risk stratification and treatment decision-making.

Flow-based method demonstrates improved accuracy for calculating wall shear stress in arterial flows from 4D flow MRI data

Four-dimensional flow magnetic resonance imaging (i.e., 4D flow MRI) has become a valuable tool for the in vivo assessment of blood flow within large vessels and cardiac chambers. As wall shear stress (WSS) has been correlated with the development and progression of cardiovascular disease, focus has been directed at developing techniques to quantify WSS directly from 4D flow MRI data. The goal of this study was to compare the accuracy of two such techniques - termed the velocity and flow-based methods - in the setting of simplified and complex flow scenarios. Synthetic MR data were created from exact solutions to the Navier-Stokes equations for the steady and pulsatile flow of an incompressible, Newtonian fluid through a rigid cylinder. In addition, synthetic MR data were created from the predicted velocity fields derived from a fluid-structure interaction (FSI) model of pulsatile flow through a thick-walled, multi-layered model of the carotid bifurcation. Compared to the analytical solutions for steady and pulsatile flow, the flow-based method demonstrated greater accuracy than the velocity-based method in calculating WSS across all changes in fluid velocity/flow rate, tube radius, and image signal-to-noise (p < 0.001). Furthermore, the velocity-based method was more sensitive to boundary segmentation than the flow-based method. When compared to results from the FSI model, the flow-based method demonstrated greater accuracy than the velocity-based method with average differences in time-averaged WSS of 0.31 ± 1.03 Pa and 0.45 ± 1.03 Pa, respectively (p <0.005). These results have implications on the utility, accuracy, and clinical translational of methods to determine WSS from 4D flow MRI.

Slice-based and time-specific hemodynamic measurements discriminate carotid artery vulnerable atherosclerotic plaques

BACKGROUND AND OBJECTIVE

Hemodynamic patterns play key roles in progression of carotid vulnerable plaques. However, most of previous studies utilized maximum or averaged value of hemodynamic measurements which is not an ideal representative of hemodynamic patterns. This study aimed to investigate the association of slice-based and time-specific hemodynamic measurements with carotid vulnerable plaque using magnetic resonance (MR) vessel wall imaging and histology.

METHODS

Thirty-two patients (mean age: 63.9±8.1 years; 25 males) with carotid atherosclerotic stenosis (≥50% stenosis) referred to carotid endarterectomy were recruited and underwent MR vessel wall imaging. Carotid plaque burden was evaluated on MR images and vulnerable plaque features including calcification, lipid-rich necrotic core, and intra-plaque hemorrhage (IPH) were identified by histology. The slice-based and time-specific hemodynamic measurements were extracted from computational fluid dynamics simulation of 3D carotid arterial model. Correlation coefficients between hemodynamic measurements and carotid plaque features were calculated and the logistic regressions with generalized estimating equation (GEE) were conducted. The value in discriminating carotid vulnerable plaque features was determined by receiver-operating-characteristic analysis.

RESULTS

Of 102 MR-histology matched slices from 32 patients, time-averaged wall shear stress (TAWSS) (r=0.263, p=0.008), oscillatory shear index (OSI) (r=-0.374, p<0.001), and peakWSS (r=0.232, p=0.019) were significantly associated with carotid IPH. The logistic regression with GEE revealed that peakWSS (OR, 1.206; 95% CI, 1.026-1.418; p, 0.023) and TAWSS (OR, 0.364, 95% CI, 0.138-0.959; p, 0.041) were significantly associated with presence of IPH after adjusting for age and BMI. In discriminating carotid IPH, the AUC of TAWSS, OSI, combined TAWSS with maximum wall thickness (MWT) and combined OSI with MWT was 0.656, 0.722, 0.761, and 0.764, respectively.

CONCLUSIONS

Slice-based and time-specific hemodynamic characteristics could effectively discriminate carotid IPH. Combination of hemodynamic measurements with carotid plaque burden might be a stronger indicator for carotid vulnerable plaque features than each measurement alone.

Hemodynamics in a three-dimensional printed aortic model: a comparison of four-dimensional phase-contrast magnetic resonance and image-based computational fluid dynamics

OBJECTIVE

This study aims to compare an electrocardiogram (ECG)-gated four-dimensional (4D) phase-contrast (PC) magnetic resonance imaging (MRI) technique and computational fluid dynamics (CFD) using variables controlled in a laboratory environment to minimize bias factors.

MATERIALS AND METHODS

Data from 4D PC-MRI were compared with computational fluid dynamics using steady and pulsatile flows at various inlet velocities. Anatomically realistic models for a normal aorta, a penetrating atherosclerotic ulcer, and an abdominal aortic aneurysm were constructed using a three-dimensional printer.

RESULTS

For the normal aorta model, the errors in the peak and the average velocities were within 5%. The peak velocities of the penetrating atherosclerotic ulcer and the abdominal aortic aneurysm models displayed a more extensive range of differences because of the high-speed and vortical fluid flows generated by the shape of the blood vessel. However, the average velocities revealed only relatively minor differences.

CONCLUSIONS

This study compared the characteristics of PC-MRI and CFD through a phantom study that only included controllable experimental parameters. Based on these results, 4D PC-MRI and CFD are powerful tools for analyzing blood flow patterns in vivo. However, there is room for future developments to improve velocity measurement accuracy.

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.

WSSNet: Aortic Wall Shear Stress Estimation Using Deep Learning on 4D Flow MRI

Wall shear stress (WSS) is an important contributor to vessel wall remodeling and atherosclerosis. However, image-based WSS estimation from 4D Flow MRI underestimates true WSS values, and the accuracy is dependent on spatial resolution, which is limited in 4D Flow MRI. To address this, we present a deep learning algorithm (WSSNet) to estimate WSS trained on aortic computational fluid dynamics (CFD) simulations. The 3D CFD velocity and coordinate point clouds were resampled into a 2D template of 48 × 93 points at two inward distances (randomly varied from 0.3 to 2.0 mm) from the vessel surface (“velocity sheets”). The algorithm was trained on 37 patient-specific geometries and velocity sheets. Results from 6 validation and test cases showed high accuracy against CFD WSS (mean absolute error 0.55 ± 0.60 Pa, relative error 4.34 ± 4.14%, 0.92 ± 0.05 Pearson correlation) and noisy synthetic 4D Flow MRI at 2.4 mm resolution (mean absolute error 0.99 ± 0.91 Pa, relative error 7.13 ± 6.27%, and 0.79 ± 0.10 Pearson correlation). Furthermore, the method was applied on in vivo 4D Flow MRI cases, effectively estimating WSS from standard clinical images. Compared with the existing parabolic fitting method, WSSNet estimates showed 2–3 × higher values, closer to CFD, and a Pearson correlation of 0.68 ± 0.12. This approach, considering both geometric and velocity information from the image, is capable of estimating spatiotemporal WSS with varying image resolution, and is more accurate than existing methods while still preserving the correct WSS pattern distribution.

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.

The impact of shape uncertainty on aortic-valve pressure-drop computations

Patient-specific image-based computational fluid dynamics (CFD) is widely adopted in the cardiovascular research community to study hemodynamics, and will become increasingly important for personalized medicine. However, segmentation of the flow domain is not exact and geometric uncertainty can be expected which propagates through the computational model, leading to uncertainty in model output. Seventy-four aortic-valves were segmented from computed tomography images at peak systole. Statistical shape modeling was used to obtain an approximate parameterization of the original segmentations. This parameterization was used to train a meta-model that related the first five shape mode coefficients and flowrate to the CFD-computed transvalvular pressure-drop. Consequently, shape uncertainty in the order of 0.5 and 1.0 mm was emulated by introducing uncertainty in the shape mode coefficients. A global variance-based sensitivity analysis was performed to quantify output uncertainty and to determine relative importance of the shape modes. The first shape mode captured the opening/closing behavior of the valve and uncertainty in this mode coefficient accounted for more than 90% of the output variance. However, sensitivity to shape uncertainty is patient-specific, and the relative importance of the fourth shape mode coefficient tended to increase with increases in valvular area. These results show that geometric uncertainty in the order of image voxel size may lead to substantial uncertainty in CFD-computed transvalvular pressure-drops. Moreover, this illustrates that it is essential to assess the impact of geometric uncertainty on model output, and that this should be thoroughly quantified for applications that wish to use image-based CFD models.

Experimental Validation of Enhanced Magnetic Resonance Imaging (EMRI) Using Particle Image Velocimetry (PIV)

Flow-sensitive four-dimensional Cardiovascular Magnetic Resonance Imaging (4D Flow CMR) has increasingly been utilised to characterise patients' blood flow, in association with patiens' state of health and disease, even though spatial and temporal resolutions still constitute a limit. Computational fluid dynamics (CFD) is a powerful tool that could expand these information and, if integrated with experimentally-obtained velocity fields, would enable to derive a large variety of the flow descriptors of interest. However, the accuracy of the flow parameters is highly influenced by the quality of the input data such as the anatomical model and boundary conditions typically derived from medical images including 4D Flow CMR. We previously proposed a novel approach in which 4D Flow CMR and CFD velocity fields are integrated to obtain an Enhanced 4D Flow CMR (EMRI), allowing to overcome the spatial-resolution limitation of 4D Flow CMR, and enable an accurate quantification of flow. In this paper, the proposed approach is validated in a U bend channel, an idealised model of the human aortic arch. The flow patterns were studied with 4D Flow CMR, CFD and EMRI, and compared with high resolution 2D PIV experiments obtained in pulsatile conditions. The main strengths and limitations of 4D Flow CMR and CFD were illustrated by exploiting the accuracy of PIV by comparing against PIV velocity fields. EMRI flow patterns showed a better qualitative and quantitative agreement with PIV results than the other techniques. EMRI enables to overcome the experimental limitations of MRI-based velocity measurements and the modelling simplifications of CFD, allowing an accurate prediction of complex flow patterns observed experimentally, while satisfying mass and momentum balance equations.

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.

Exploring the Relationships Between Hemodynamic Stresses in the Carotid Arteries

Background: Atherosclerosis manifests as a focal disease, often affecting areas with complex hemodynamics such as the carotid bifurcation. The magnitude and regularity of the hemodynamic shear stresses acting on the vessel wall are thought to generate risk patterns unique to each patient and play a role in the pathogenesis of atherosclerosis. The involvement of different expressions of shear stress in the pathogenesis of carotid atherosclerosis highlights the need to characterize and compare the differential impact of the various expressions of shear stress in the atherosclerotic carotid bifurcation. Therefore, the aim of this study is to characterize and compare hemodynamic wall shear stresses (WSS) in the carotid arteries of subjects with asymptomatic atherosclerotic plaques. Shear stresses were also compared against vessel diameter and bifurcation angle to examine the relationships with the geometry of the carotid bifurcation. Methods: 4D Flow MRI and contrast-enhanced MRA data were acquired for 245 subjects with atherosclerotic plaques of at least 2.7 mm in conjunction with the Swedish CArdioPulmonary bioImage Study (SCAPIS). Following automatic segmentation and geometric analysis, time-resolved WSS and near-wall turbulent kinetic energy (nwTKE) were derived from the 4D Flow data. Whole-cycle parameters including time-averaged WSS and nwTKE, and the oscillatory shear index (OSI) were calculated. Pairwise Spearman rank-correlation analyses were used to investigate relationships among the hemodynamic as well as geometric parameters. Results: One hundred and seventy nine subjects were successfully segmented using automated tools and subsequently geometric and hemodynamic analyses were performed. Temporally resolved WSS and nwTKE were strongly correlated, ρ = 0.64. Cycle-averaged WSS and nwTKE were moderately correlated, ρ = 0.57. Cycle-average nwTKE was weakly correlated to OSI (ρ = -0.273), revealing that nwTKE provides information about disturbed flow on the vessel wall that OSI does not. In this cohort, there was large inter-individual variation for both WSS and nwTKE. Both WSS and nwTKE varied most within the external carotid artery. WSS, nwTKE, and OSI were weakly correlated to vessel diameter and bifurcation angle. Conclusion: The turbulent and mean component of WSS were examined together in vivo for the first time, and a strong correlation was found between them. nwTKE presents the opportunity to quantify turbulent wall stresses in vivo and gain insight into the effects of disturbed flow on the vessel wall. Neither vessel diameter nor bifurcation angle were found to be strongly correlated to the turbulent or mean component of WSS in this cohort.

In vitro performance of echoPIV for assessment of laminar flow profiles in a carotid artery stent

Purpose: Detailed blood flow studies may contribute to improvements in carotid artery stenting. High-frame-rate contrast-enhanced ultrasound followed by particle image velocimetry (PIV), also called echoPIV, is a technique to study blood flow patterns in detail. The performance of echoPIV in presence of a stent has not yet been studied extensively. We compared the performance of echoPIV in stented and nonstented regions in an in vitro flow setup.

Approach: A carotid artery stent was deployed in a vessel-mimicking phantom. High-frame-rate contrast-enhanced ultrasound images were acquired with various settings. Signal intensities of the contrast agent, velocity values, and flow profiles were calculated.

Results: The results showed decreased signal intensities and correlation coefficients inside the stent, however, PIV analysis in the stent still resulted in plausible flow vectors.

Conclusions: Velocity values and laminar flow profiles can be measured in vitro in stented arteries using echoPIV.

Validation and Diagnostic Performance of a CFD-Based Non-invasive Method for the Diagnosis of Aortic Coarctation

Purpose: The clinical diagnosis of aorta coarctation (CoA) constitutes a challenge, which is usually tackled by applying the peak systolic pressure gradient (PSPG) method. Recent advances in computational fluid dynamics (CFD) have suggested that multi-detector computed tomography angiography (MDCTA)-based CFD can serve as a non-invasive PSPG measurement. The aim of this study was to validate a new CFD method that does not require any medical examination data other than MDCTA images for the diagnosis of CoA.

Materials and methods: Our study included 65 pediatric patients (38 with CoA, and 27 without CoA). All patients underwent cardiac catheterization to confirm if they were suffering from CoA or any other congenital heart disease (CHD). A series of boundary conditions were specified and the simulated results were combined to obtain a stenosis pressure-flow curve. Subsequently, we built a prediction model and evaluated its predictive performance by considering the AUC of the ROC by 5-fold cross-validation.

Results: The proposed MDCTA-based CFD method exhibited a good predictive performance in both the training and test sets (average AUC: 0.948 vs. 0.958; average accuracies: 0.881 vs. 0.877). It also had a higher predictive accuracy compared with the non-invasive criteria presented in the European Society of Cardiology (ESC) guidelines (average accuracies: 0.877 vs. 0.539).

Conclusion: The new non-invasive CFD-based method presented in this work is a promising approach for the accurate diagnosis of CoA, and will likely benefit clinical decision-making.

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.

A Dual-VENC Four-Dimensional Flow MRI Framework for Analysis of Subject-Specific Heterogeneous Nonlinear Vessel Deformation

Advancement of subject-specific in silico medicine requires new imaging protocols tailored to specific anatomical features, paired with new constitutive model development based on structure/function relationships. In this study, we develop a new dual-velocity encoding coefficient (VENC) 4D flow MRI protocol that provides unprecedented spatial and temporal resolution of in vivo aortic deformation. All previous dual-VENC 4D flow MRI studies in the literature focus on an isolated segment of the aorta, which fail to capture the full spectrum of aortic heterogeneity that exists along the vessel length. The imaging protocol developed provides high sensitivity to all blood flow velocities throughout the entire cardiac cycle, overcoming the challenge of accurately measuring the highly unsteady nonuniform flow field in the aorta. Cross-sectional area change, volumetric flow rate, and compliance are observed to decrease with distance from the heart, while pulse wave velocity (PWV) is observed to increase. A nonlinear aortic lumen pressure-area relationship is observed throughout the aorta such that a high vessel compliance occurs during diastole, and a low vessel compliance occurs during systole. This suggests that a single value of compliance may not accurately represent vessel behavior during a cardiac cycle in vivo. This high-resolution MRI data provide key information on the spatial variation in nonlinear aortic compliance, which can significantly advance the state-of-the-art of in-silico diagnostic techniques for the human aorta.

Highly accelerated aortic 4D flow MRI using compressed sensing: Performance at different acceleration factors in patients with aortic disease

PURPOSE

To systematically assess the feasibility and performance of a highly accelerated compressed sensing (CS) 4D flow MRI framework at three different acceleration factors (R) for the quantification of aortic flow dynamics and wall shear stress (WSS) in patients with aortic disease.

METHODS

Twenty patients with aortic disease (58 ± 15 y old; 19 M) underwent four 4D flow scans: one conventional (GRAPPA, R = 2) and three CS 4D flows with R = 5.7, 7.7, and 10.2. All scans were acquired with otherwise equivalent imaging parameters on a 1.5T scanner. Peak-systolic velocity (V_max ), peak flow (Q_max ), and net flow (Q_net ) were quantified at the ascending aorta (AAo), arch, and descending aorta (DAo). WSS was calculated at six regions within the AAo and arch.

RESULTS

Mean scan times for the conventional and CS 4D flows with R = 5.7, 7.7, and 10.2 were 9:58 ± 2:58 min, 3:40 ± 1:19 min, 2:50 ± 0:56 min, and 2:05 ± 0:42 min, respectively. V_max , Q_max , and Q_net were significantly underestimated by all CS protocols (underestimation ≤ -7%, -9%, and -10% by CS, R = 5.7, 7.7, and 10.2, respectively). WSS measurements showed the highest underestimation by all CS protocols (underestimation ≤ -9%, -12%, and -14% by CS, R = 5.7, 7.7, and 10.2).

CONCLUSIONS

Highly accelerated aortic CS 4D flow at R = 5.7, 7.7, and 10.2 showed moderate agreement with the conventional 4D flow, despite systematically underestimating various hemodynamic parameters. The shortened scan time may enable the clinical translation of CS 4D flow, although potential hemodynamic underestimation should be considered when interpreting the results.

Combining statistical shape modeling, CFD, and meta-modeling to approximate the patient-specific pressure-drop across the aortic valve in real-time

BACKGROUND

Advances in medical imaging, segmentation techniques, and high performance computing have stimulated the use of complex, patient-specific, three-dimensional Computational Fluid Dynamics (CFD) simulations. Patient-specific, CFD-compatible geometries of the aortic valve are readily obtained. CFD can then be used to obtain the patient-specific pressure-flow relationship of the aortic valve. However, such CFD simulations are computationally expensive, and real-time alternatives are desired.

AIM

The aim of this work is to evaluate the performance of a meta-model with respect to high-fidelity, three-dimensional CFD simulations of the aortic valve.

METHODS

Principal component analysis was used to build a statistical shape model (SSM) from a population of 74 iso-topological meshes of the aortic valve. Synthetic meshes were created with the SSM, and steady-state CFD simulations at flow-rates between 50 and 650 mL/s were performed to build a meta-model. The meta-model related the statistical shape variance, and flow-rate to the pressure-drop.

RESULTS

Even though the first three shape modes account for only 46% of shape variance, the features relevant for the pressure-drop seem to be captured. The three-mode shape-model approximates the pressure-drop with an average error of 8.8% to 10.6% for aortic valves with a geometric orifice area below 150 mm² . The proposed methodology was least accurate for aortic valve areas above 150 mm² . Further reduction to a meta-model introduces an additional 3% error.

CONCLUSIONS

Statistical shape modeling can be used to capture shape variation of the aortic valve. Meta-models trained by SSM-based CFD simulations can provide an estimate of the pressure-flow relationship in real-time.

Towards quantitative evaluation of wall shear stress from 4D flow imaging

Wall shear stress (WSS) is a relevant hemodynamic indicator of the local stress applied on the endothelium surface. More specifically, its spatiotemporal distribution reveals crucial in the evolution of many pathologies such as aneurysm, stenosis, and atherosclerosis. This paper introduces a new solution, called PaLMA, to quantify the WSS from 4D Flow MRI data. It relies on a two-step local parametric model, to accurately describe the vessel wall and the velocity-vector field in the neighborhood of a given point of interest. Extensive validations have been performed on synthetic 4D Flow MRI data, including four datasets generated from patient specific computational fluid dynamics simulations on carotids. The validation tests are focused on the impact of the noise component, of the resolution level, and of the segmentation accuracy concerning the vessel position in the context of complex flow patterns. In simulated cases aimed to reproduce clinical acquisition conditions, the WSS quantification performance reached by PaLMA is significantly higher (with a gain in RMSE of 12 to 27%) than the reference one obtained using the smoothing B-spline method proposed by Potters et al. (2015) method, while the computation time is equivalent for both WSS quantification methods.

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.

Bilateral Multiband 4D Flow MRI of the Carotid Arteries at 7T

PURPOSE

Simultaneous multislab (SMSb) 4D flow MRI was developed and implemented at 7T for accelerated acquisition of the 3D blood velocity vector field in both carotid bifurcations.

METHODS

SMSb was applied to 4D flow to acquire blood velocities in both carotid bifurcations in sagittal orientation using a local transmit/receive coil at 7T. B 1 + transmit efficiency was optimized by B 1 + shimming. SMSb 4D flow was obtained in 8 healthy subjects in single-band (SB) and multiband (MB) fashion. Additionally, MB data were retrospectively undersampled to simulate GRAPPA R = 2 (MB2_GRAPPA2), and both SB datasets were added to form an artificial MB dataset (SumSB). The band separation performance was quantified by signal leakage. Peak velocity and total flow values were calculated and compared to SB via intraclass correlation analysis (ICC).

RESULTS

Clean slab separation was achieved yielding a mean signal leakage of 13% above the mean SB noise level. Mean total flow for MB2, SumSB, and MB_GRAPPA2 deviated less than 9% from the SB values. Peak velocities averaged over all vessels and subjects were 0.48 ± 0.11 m/s for SB, 0.47 ± 0.12 m/s for SumSB, 0.50 ± 0.13 m/s for MB2, and 0.53 ± 0.13 m/s for MB2_GRAPPA2. ICC revealed excellent absolute agreement and consistency of total flow for all methods compared to SB2. Peak velocity showed good to excellent agreement and consistency for SumSB and MB2 and MB2_GRAPPA2 method showed poor to excellent agreement and good to excellent consistency.

CONCLUSION

Simultaneous multislab 4D Flow MRI allows accurate quantification of total flow and peak velocity while reducing scan times.

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.

An MRI-based method to register patient-specific wall shear stress data to histology

Wall shear stress (WSS), the frictional force exerted on endothelial cells by blood flow, is hypothesised to influence atherosclerotic plaque growth and composition. We developed a methodology for image registration of MR and histology images of advanced human carotid plaques and corresponding WSS data, obtained by MRI and computational fluid dynamics. The image registration method requires four types of input images, in vivo MRI, ex vivo MRI, photographs of transversally sectioned plaque tissue and histology images. These images are transformed to a shared 3D image domain by applying a combination of rigid and non-rigid registration algorithms. Transformation matrices obtained from registration of these images are used to transform subject-specific WSS data to the shared 3D image domain as well. WSS values originating from the 3D WSS map are visualised in 2D on the corresponding lumen locations in the histological sections and divided into eight radial segments. In each radial segment, the correlation between WSS values and plaque composition based on histological parameters can be assessed. The registration method was successfully applied to two carotid endarterectomy specimens. The resulting matched contours from the imaging modalities had Hausdorff distances between 0.57 and 0.70 mm, which is in the order of magnitude of the in vivo MRI resolution. We simulated the effect of a mismatch in the rigid registration of imaging modalities on WSS results by relocating the WSS data with respect to the stack of histology images. A 0.6 mm relocation altered the mean WSS values projected on radial bins on average by 0.59 Pa, compared to the output of original registration. This mismatch of one image slice did not change the correlation between WSS and plaque thickness. In conclusion, we created a method to investigate correlations between WSS and plaque composition.

Mechanotransmission of haemodynamic forces by the endothelial glycocalyx in a full-scale arterial model

The glycocalyx has been identified as a key mechano-sensor of the shear forces exerted by streaming blood onto the vascular endothelial lining. Although the biochemical reaction to the blood flow has been extensively studied, the mechanism of transmission of the haemodynamic shear forces to the endothelial transmembrane anchoring structures and, consequently, to the subcellular elements in the cytoskeleton, is still not fully understood. Here we apply a multiscale approach to elucidate how haemodynamic shear forces are transmitted to the transmembrane anchors of endothelial cells. Wall shear stress time histories, as obtained from image-based computational haemodynamics models of a carotid bifurcation, are used as a load and a continuum model is applied to obtain the mechanical response of the glycocalyx all along the cardiac cycle. The main findings of this in silico study are that: (1) the forces transmitted to the transmembrane anchors are in the range of 1-10 pN, which is in the order of magnitude reported for the different conformational states of transmembrane mechanotranductors; (2) locally, the forces transmitted to the anchors of the glycocalyx structure can be markedly different from the near-wall haemodynamic shear forces both in amplitude and frequency content. The findings of this in silico approach warrant future studies focusing on the actual forces transmitted to the transmembrane mechanotransductors, which might outperform haemodynamic descriptors of disturbed shear as localizing factors of vascular disease.

Modern Concepts in Regenerative Therapy for Ischemic Stroke: From Stem Cells for Promoting Angiogenesis to 3D-Bioprinted Scaffolds Customized via Carotid Shear Stress Analysis

Ischemic stroke is associated with a tremendous economic and societal burden, and only a few therapies are currently available for the treatment of this devastating disease. The main therapeutic approaches used nowadays for the treatment of ischemic brain injury aim to achieve reperfusion, neuroprotection and neurorecovery. Therapeutic angiogenesis also seems to represent a promising tool to improve the prognosis of cerebral ischemia. This review aims to present the modern concepts and the current status of regenerative therapy for ischemic stroke and discuss the main results of major clinical trials addressing the effectiveness of stem cell therapy for achieving neuroregeneration in ischemic stroke. At the same time, as a glimpse into the future, this article describes modern concepts for stroke prevention, such as the implantation of bioprinted scaffolds seeded with stem cells, whose 3D geometry is customized according to carotid shear stress.

Reconciling PC-MRI and CFD: An in-vitro study

Several well-resolved 4D Flow MRI acquisitions of an idealized rigid flow phantom featuring an aneurysm, a curved channel as well as a bifurcation were performed under pulsatile regime. The resulting hemodynamics were processed to remove MRI artifacts. Subsequently, they were compared with CFD predictions computed on the same flow domain, using an in-house high-order low dissipative flow solver. Results show that reaching a good agreement is not straightforward but requires proper treatments of both techniques. Several sources of discrepancies are highlighted and their impact on the final correlation evaluated. While a very poor correlation (r²  = 0.63) is found in the entire domain between raw MRI and CFD data, correlation as high as r²  = 0.97 is found when artifacts are removed by post-processing the MR data and down sampling the CFD results to match the MRI spatial and temporal resolutions. This work demonstrates that, in a well-controlled environment, both PC-MRI and CFD might bring reliable and correlated flow quantities when a proper methodology to reduce the errors is followed.

Spatial correlations between MRI-derived wall shear stress and vessel wall thickness in the carotid bifurcation

Background

To explore the possibility of creating three-dimensional (3D) estimation models for patient-specific wall thickness (WT) maps using patient-specific and cohort-averaged WT, wall shear stress (WSS), and vessel diameter maps in asymptomatic atherosclerotic carotid bifurcations.

Methods

Twenty subjects (aged 75 ± 6 years [mean ± standard deviation], eight women) underwent a 1.5-T MRI examination. Non-gated 3D phase-contrast gradient-echo images and proton density-weighted echo-planar images were retrospectively assessed for WSS, diameter estimation, and WT measurements. Spearman’s ρ and scatter plots were used to determine correlations between individual WT, WSS, and diameter maps. A bootstrapping technique was used to determine correlations between 3D cohort-averaged WT, WSS, and diameter maps. Linear regression between the cohort-averaged WT, WSS, and diameter maps was used to predict individual 3D WT.

Results

Spearman’s ρ averaged over the subjects was − 0.24 ± 0.18 (p < 0.001) and 0.07 ± 0.28 (p = 0.413) for WT versus WSS and for WT versus diameter relations, respectively. Cohort-averaged ρ, averaged over 1000 bootstraps, was − 0.56 (95% confidence interval [− 0.74,− 0.38]) for WT versus WSS and 0.23 (95% confidence interval [− 0.06, 0.52]) for WT versus diameter. Scatter plots did not reveal relationships between individual WT and WSS or between WT and diameter data. Linear relationships between these parameters became apparent after averaging over the cohort. Spearman’s ρ between the original and predicted WT maps was 0.21 ± 0.22 (p < 0.001).

Conclusions

With a combination of bootstrapping and cohort-averaging methods, 3D WT maps can be predicted from the individual 3D WSS and diameter maps. The methodology may help to elucidate pathological processes involving WSS in carotid atherosclerosis.

Electronic supplementary material

The online version of this article (10.1186/s41747-018-0058-1) contains supplementary material, which is available to authorized users.

Normal patterns of thoracic aortic wall shear stress measured using four-dimensional flow MRI in a large population

Wall shear stress (WSS) plays a governing role in vascular remodeling and a pathogenic role in vessel wall diseases. However, little is known of the normal WSS patterns in the aorta as there is currently no practical means to routinely measure WSS and no normal ranges derived from population data exist. WSS measurements were made on the aorta of 224 subjects with normal anatomy using four-dimensional flow MRI with multiple encoding velocities and an optimized postprocessing routine. The spatial and temporal variation in WSS and oscillatory shear index was analyzed using a flat map representation of the unfolded aorta. The influence of aortic shape and velocity on WSS was evaluated using regression analysis. WSS in the thoracic aorta is dominated by axial flow. Average peak systolic WSS was 1.79 ± 0.71 Pa in the aortic arch and was significantly higher at 2.23 ± 1.04 Pa in the descending aorta, with a strong negative correlation with advancing age. The spatial distribution of WSS is highly heterogeneous, with a localized region of elevated WSS along the length of the anterior wall seen across all individuals. Our data demonstrate that accurate four-dimensional flow-derived WSS measurement is feasible, and we further provide a standardized parametric approach for presentation and analysis. We present a normal range for WSS across the lifespan, demonstrating a decrease in WSS with advancing age as well as illustrating the high degree of spatial and temporal variation. NEW & NOTEWORTHY With the use of four-dimensional flow MRI and postprocessing, accurate direct measurement of wall shear stress (WSS) was performed in a population of normal thoracic aortas ( n = 224). WSS was higher in the descending aorta compared with the aortic arch and decreased with age. A heterogeneous pattern of elevated WSS along the length of the aorta anterior wall was consistent across the population. This work provides normal data across the adult age range, permitting comparison with pathology.

Hemodynamic analysis of carotid artery after endarterectomy: a preliminary and quantitative imaging study based on computational fluid dynamics and magnetic resonance angiography

Background

The carotid blood flow following carotid endarterectomy (CEA) is not fully understood. Computational fluid dynamics (CFD) is a promising method to study blood flow. This study is to investigate local hemodynamic characteristics after CEA via the use of unenhanced magnetic resonance angiography (MRA) and CFD.

Methods

Eight carotid arteries with atherosclerosis and sixteen normal carotid arteries were included in this study. Time-of-flight (TOF) and phase contrast (PC) MRA were applied for the measurement of three-dimensional artery geometries and velocity profile under CFD simulation. The hemodynamic parameters of the proximal internal carotid artery (ICA) including velocity, ICA/common carotid artery (CCA) velocity ratio, mean, maximum, minimum and gradient of wall shear stress (WSS_mean, WSS_max, WSS_min and WSSG) were calculated before and after CEA. Morphologic characteristics of the carotid including bifurcation angle, tortuosity and planarity were also analyzed.

Results

Compared with pre-CEA, there was a significant reduction in post-CEA velocity, WSS_max, WSS_mean, and WSSG, by 87.24%±13.38%, 86.86%±14.97%, 57.32%±56.71% and 69.74%±37.03% respectively, whereas WSS_min was almost unchanged. ICA/ CCA velocity ratios increased significantly after CEA. We also found that the post-CEA flow conditions were positively remodelled to approximate the conditions in normal arteries. The correlation between PC-MRA and CFD was excellent for the measurement of maximum velocity at the external carotid artery (r=0.846).

Conclusions

Our preliminary results indicated that major flow dynamics were restored shortly following CEA, and CFD based on MRA measurements could be useful for quantitative evaluation of hemodynamic outcomes after CEA.

Retrospective Study of Hemodynamic Changes Before and After Carotid Stenosis Formation by Vessel Surface Repairing

Prospective observation of hemodynamic changes before and after the formation of atherosclerotic stenosis in the carotid artery is difficult. Thus, a vessel surface repairing method was used for retrospective hemodynamic study before and after atherosclerotic stenosis formation in carotid artery. The three-dimensional geometry of sixteen sinus atherosclerotic stenosis carotid arteries were repaired and restored as normal arteries. Computational fluid dynamics analysis was performed to estimate wall shear stress (WSS), velocity and vortex in atherosclerosis-free areas and sinus in stenosis-repaired carotid artery. The analysis was also performed in the stenotic segment and upstream and downstream of stenosis in stenotic carotid artery. Compared to the atherosclerosis-free areas in stenosis-repaired carotid artery, sinus presented significantly lower WSS (P < 0.05), lower velocity (P < 0.05) and apparent vortex. Compared to the sinus, the WSS in the upstream of stenosis was lower (P < 0.05), while in the downstream area was similar (P = 0.87), both upstream and downstream of stenosis demonstrated similar velocity to sinus (P = 0.76 and P = 0.36, respectively) and apparent vortex. Atherosclerosis-prone areas including normal carotid sinus and upstream and downstream of stenosis in stenotic carotid artery were subjected to lower WSS and velocity as well as apparent vortex, thereby might be associated with the formation and progress of atherosclerosis.

Atherosclerosis at arterial bifurcations: evidence for the role of haemodynamics and geometry

Atherosclerotic plaques are found at distinct locations in the arterial system, despite the exposure to systemic risk factors of the entire vascular tree. From the study of arterial bifurcation regions, emerges ample evidence that haemodynamics are involved in the local onset and progression of the atherosclerotic disease. This observed co-localisation of disturbed flow regions and lesion prevalence at geometrically predisposed districts such as arterial bifurcations has led to the formulation of a ‘haemodynamic hypothesis’, that in this review is grounded to the most current research concerning localising factors of vascular disease. In particular, this review focuses on carotid and coronary bifurcations because of their primary relevance to stroke and heart attack. We highlight reported relationships between atherosclerotic plaque location, progression and composition, and fluid forces at vessel’s wall, in particular shear stress and its ‘easier-tomeasure’ surrogates, i.e. vascular geometric attributes (because geometry shapes the flow) and intravascular flow features (because they mediate disturbed shear stress), in order to give more insight in plaque initiation and destabilisation. Analogous to Virchow’s triad for thrombosis, atherosclerosis must be thought of as subject to a triad of, and especially interactions among, haemodynamic forces, systemic risk factors, and the biological response of the wall.

Influence of principal component analysis acceleration factor on velocity measurement in 2D and 4D PC-MRI

OBJECTIVE

The objective of the study was to determine how to optimize 2D and 4D phase-contrast magnetic resonance imaging (PC-MRI) acquisitions to acquire flow velocities in millimetric vessels. In particular, we search for the best compromise between acquisition time and accuracy and assess the influence of the principal component analysis (PCA).

MATERIALS AND METHODS

2D and 4D PC-MRI measurements are conducted within two in vitro vessel phantoms: a Y-bifurcation phantom, the branches of which range from 2 to 5 mm in diameter, and a physiological subject-specific phantom of the carotid bifurcation. The same sequences are applied in vivo in carotid vasculature.

RESULTS

For a vessel oriented in the axial direction, both 2D and axial 4D PC-MRI provided accuracy measurements regardless of the k-t PCA factor, while the acquisition time is reduced by a factor 6 for k-t PCA maximum value. The in vivo measurements show that the proposed sequences are adequate to acquire 2D and 4D velocity fields in millimetric vessels and with clinically realistic time durations.

CONCLUSION

The study shows the feasibility of conducting fast, high-resolution PC-MRI flow measurements in millimetric vessels and that it is worth maximizing the k-t PCA factor to reduce the acquisition time in the case of 2D and 4D axial acquisitions.

Phase-contrast MRI versus numerical simulation to quantify hemodynamical changes in cerebral aneurysms after flow diverter treatment

Cerebral aneurysms are a major risk factor for intracranial bleeding with devastating consequences for the patient. One recently established treatment is the implantation of flow-diverters (FD). Methods to predict their treatment success before or directly after implantation are not well investigated yet. The aim of this work was to quantitatively study hemodynamic parameters in patient-specific models of treated cerebral aneurysms and its correlation with the clinical outcome. Hemodynamics were evaluated using both computational fluid dynamics (CFD) and phase contrast (PC) MRI. CFD simulations and in vitro MRI measurements were done under similar flow conditions and results of both methods were comparatively analyzed. For preoperative and postoperative distribution of hemodynamic parameters, CFD simulations and PC-MRI velocity measurements showed similar results. In both cases where no occlusion of the aneurysm was observed after six months, a flow reduction of about 30-50% was found, while in the clinically successful case with complete occlusion of the aneurysm after 6 months, the flow reduction was about 80%. No vortex was observed in any of the three models after treatment. The results are in agreement with recent studies suggesting that CFD simulations can predict post-treatment aneurysm flow alteration already before implantation of a FD and PC-MRI could validate the predicted hemodynamic changes right after implantation of a FD.

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.

Vessel wall characterization using quantitative MRI: what's in a number?

The past decade has witnessed the rapid development of new MRI technology for vessel wall imaging. Today, with advances in MRI hardware and pulse sequences, quantitative MRI of the vessel wall represents a real alternative to conventional qualitative imaging, which is hindered by significant intra- and inter-observer variability. Quantitative MRI can measure several important morphological and functional characteristics of the vessel wall. This review provides a detailed introduction to novel quantitative MRI methods for measuring vessel wall dimensions, plaque composition and permeability, endothelial shear stress and wall stiffness. Together, these methods show the versatility of non-invasive quantitative MRI for probing vascular disease at several stages. These quantitative MRI biomarkers can play an important role in the context of both treatment response monitoring and risk prediction. Given the rapid developments in scan acceleration techniques and novel image reconstruction, we foresee the possibility of integrating the acquisition of multiple quantitative vessel wall parameters within a single scan session.

The Effect of Spatial and Temporal Resolution of Cine Phase Contrast MRI on Wall Shear Stress and Oscillatory Shear Index Assessment

Introduction

Wall shear stress (WSS) and oscillatory shear index (OSI) are associated with atherosclerotic disease. Both parameters are derived from blood velocities, which can be measured with phase-contrast MRI (PC-MRI). Limitations in spatiotemporal resolution of PC-MRI are known to affect these measurements. Our aim was to investigate the effect of spatiotemporal resolution using a carotid artery phantom.

Methods

A carotid artery phantom was connected to a flow set-up supplying pulsatile flow. MRI measurement planes were placed at the common carotid artery (CCA) and internal carotid artery (ICA). Two-dimensional PC-MRI measurements were performed with thirty different spatiotemporal resolution settings. The MRI flow measurement was validated with ultrasound probe measurements. Mean flow, peak flow, flow waveform, WSS and OSI were compared for these spatiotemporal resolutions using regression analysis. The slopes of the regression lines were reported in %/mm and %/100ms. The distribution of low and high WSS and OSI was compared between different spatiotemporal resolutions.

Results

The mean PC-MRI CCA flow (2.5±0.2mL/s) agreed with the ultrasound probe measurements (2.7±0.02mL/s). Mean flow (mL/s) depended only on spatial resolution (CCA:-13%/mm, ICA:-49%/mm). Peak flow (mL/s) depended on both spatial (CCA:-13%/mm, ICA:-17%/mm) and temporal resolution (CCA:-19%/100ms, ICA:-24%/100ms). Mean WSS (Pa) was in inverse relationship only with spatial resolution (CCA:-19%/mm, ICA:-33%/mm). OSI was dependent on spatial resolution for CCA (-26%/mm) and temporal resolution for ICA (-16%/100ms). The regions of low and high WSS and OSI matched for most of the spatiotemporal resolutions (CCA:30/30, ICA:28/30 cases for WSS; CCA:23/30, ICA:29/30 cases for OSI).

Conclusion

We show that both mean flow and mean WSS are independent of temporal resolution. Peak flow and OSI are dependent on both spatial and temporal resolution. However, the magnitude of mean and peak flow, WSS and OSI, and the spatial distribution of OSI and WSS did not exhibit a strong dependency on spatiotemporal resolution.

High shear stress relates to intraplaque haemorrhage in asymptomatic carotid plaques

BACKGROUND AND AIMS

Carotid artery plaques with vulnerable plaque components are related to a higher risk of cerebrovascular accidents. It is unknown which factors drive vulnerable plaque development. Shear stress, the frictional force of blood at the vessel wall, is known to influence plaque formation. We evaluated the association between shear stress and plaque components (intraplaque haemorrhage (IPH), lipid rich necrotic core (LRNC) and/or calcifications) in relatively small carotid artery plaques in asymptomatic persons.

METHODS

Participants (n = 74) from the population-based Rotterdam Study, all with carotid atherosclerosis assessed on ultrasound, underwent carotid MRI. Multiple MRI sequences were used to evaluate the presence of IPH, LRNC and/or calcifications in plaques in the carotid arteries. Images were automatically segmented for lumen and outer wall to obtain a 3D reconstruction of the carotid bifurcation. These reconstructions were used to calculate minimum, mean and maximum shear stresses by applying computational fluid dynamics with subject-specific inflow conditions. Associations between shear stress measures and plaque composition were studied using generalized estimating equations analysis, adjusting for age, sex and carotid wall thickness.

RESULTS

The study group consisted of 93 atherosclerotic carotid arteries of 74 participants. In plaques with higher maximum shear stresses, IPH was more often present (OR per unit increase in maximum shear stress (log transformed) = 12.14; p = 0.001). Higher maximum shear stress was also significantly associated with the presence of calcifications (OR = 4.28; p = 0.015).

CONCLUSIONS

Higher maximum shear stress is associated with intraplaque haemorrhage and calcifications.

Age-related changes in aortic 3D blood flow velocities and wall shear stress: Implications for the identification of altered hemodynamics in patients with aortic valve disease

PURPOSE

To investigate age-related changes in peak systolic aortic 3D velocity and wall shear stress (WSS) in healthy controls and to investigate the importance of age-matching for 3D mapping of abnormal aortic hemodynamics in bicuspid aortic valve disease (BAV).

MATERIALS AND METHODS

4D flow MRI (fields strengths = 1.5-3T; resolution = 2.2-3.9 × 1.7-2.6 × 2.2-4.0 mm(3) ; Venc = 150-250 cm/s; TE/TR/FA = 2.3-2.8/4.7-5.4msec/7-15°) was performed in 56 controls (age range: 19-78 years) and in two BAV patient groups each consisting of 10 subjects (group 1: 20-29 years, group 2: 52-57 years). Heat maps showing abnormal 3D velocity and WSS were created for the BAV patients by comparison with an age-matched and with an unmatched control group. The fraction of the aorta exposed to abnormal velocity/WSS was calculated relative to the total aortic volume/surface.

RESULTS

Significant inverse relationships between age and healthy velocity/WSS were found (R(2)  = 0.32/0.39, P < 0.001). For BAV group 1, abnormally elevated velocity/WSS was overestimated when compared with older controls (51-60 years) than when correctly age-matched (∼25 ± 14% vs. ∼8 ± 5%). For BAV group 2, abnormally decreased velocity/WSS was overestimated when compared with younger controls (21-30 years) than when correctly age-matched (∼9 ± 7% vs. 1 ± 1%).

CONCLUSION

Significant correlations exist between age and peak systolic velocity and WSS. Therefore, robust age-matching is important when creating abnormal 3D aortic velocity and WSS maps for patients with BAV.

Computational fluid dynamic simulation of human carotid artery bifurcation based on anatomy and volumetric blood flow rate measured with magnetic resonance imaging

Blood flow patterns and local hemodynamic parameters have been widely associated with the onset and progression of atherosclerosis in the carotid artery. Assessment of these parameters can be performed noninvasively using cine phase-contrast (PC) magnetic resonance imaging (MRI). In addition, in the last two decades, computational fluid dynamics (CFD) simulation in three dimensional models derived from anatomic medical images has been employed to investigate the blood flow in the carotid artery. This study developed a workflow of a subject-specific CFD analysis using MRI to enhance estimating hemodynamics of the carotid artery. Time-of-flight (TOF) MRI scans were used to construct three-dimensional computational models. PC-MRI measurements were utilized to impose the boundary condition at the inlet and a 0-dimensional lumped parameter model was employed for the outflow boundary condition. The choice of different viscosity models of blood flow as a source of uncertainty was studied, by means of the axial velocity, wall shear stress, and oscillatory shear index. The sequence of workflow in CFD analysis was optimized for a healthy subject using PC-MRI. Then, a patient with carotid artery stenosis and its hemodynamic parameters were examined. The simulations indicated that the lumped parameter model used at the outlet gives physiologically reasonable values of hemodynamic parameters. Moreover, the dependence of hemodynamics parameters on the viscosity models was observed to vary for different geometries. Other factors, however, may be required for a more accurate CFD analysis, such as the segmentation and smoothness of the geometrical model, mechanical properties of the artery's wall, and the prescribed velocity profile at the inlet.

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.

MRI-based biomechanical parameters for carotid artery plaque vulnerability assessment

Carotid atherosclerotic plaques are a major cause of ischaemic stroke. The biomechanical environment to which the arterial wall and plaque is subjected to plays an important role in the initiation, progression and rupture of carotid plaques. MRI is frequently used to characterize the morphology of a carotid plaque, but new developments in MRI enable more functional assessment of carotid plaques. In this review, MRI based biomechanical parameters are evaluated on their current status, clinical applicability, and future developments. Blood flow related biomechanical parameters, including endothelial wall shear stress and oscillatory shear index, have been shown to be related to plaque formation. Deriving these parameters directly from MRI flow measurements is feasible and has great potential for future carotid plaque development prediction. Blood pressure induced stresses in a plaque may exceed the tissue strength, potentially leading to plaque rupture. Multi-contrast MRI based stress calculations in combination with tissue strength assessment based on MRI inflammation imaging may provide a plaque stress-strength balance that can be used to assess the plaque rupture risk potential. Direct plaque strain analysis based on dynamic MRI is already able to identify local plaque displacement during the cardiac cycle. However, clinical evidence linking MRI strain to plaque vulnerability is still lacking. MRI based biomechanical parameters may lead to improved assessment of carotid plaque development and rupture risk. However, better MRI systems and faster sequences are required to improve the spatial and temporal resolution, as well as increase the image contrast and signal-to-noise ratio.

Assessment of Hemodynamic Conditions in the Aorta Following Root Replacement with Composite Valve-Conduit Graft

This paper presents the analysis of detailed hemodynamics in the aortas of four patients following replacement with a composite bio-prosthetic valve-conduit. Magnetic resonance image-based computational models were set up for each patient with boundary conditions comprising subject-specific three-dimensional inflow velocity profiles at the aortic root and central pressure waveform at the model outlet. Two normal subjects were also included for comparison. The purpose of the study was to investigate the effects of the valve-conduit on flow in the proximal and distal aorta. The results suggested that following the composite valve-conduit implantation, the vortical flow structure and hemodynamic parameters in the aorta were altered, with slightly reduced helical flow index, elevated wall shear stress and higher non-uniformity in wall shear compared to normal aortas. Inter-individual analysis revealed different hemodynamic conditions among the patients depending on the conduit configuration in the ascending aorta, which is a key factor in determining post-operative aortic flow. Introducing a natural curvature in the conduit to create a smooth transition between the conduit and native aorta may help prevent the occurrence of retrograde and recirculating flow in the aortic arch, which is particularly important when a large portion or the entire ascending aorta needs to be replaced.

Advanced MRI for carotid plaque imaging

Atherosclerosis is the ubiquitous underling pathological process that manifests in heart attack and stroke, cumulating in the death of one in three North American adults. High-resolution magnetic resonance imaging (MRI) is able to delineate atherosclerotic plaque components and total plaque burden within the carotid arteries. Using dedicated hardware, high resolution images can be obtained. Combining pre- and post-contrast T1, T2, proton-density, and magnetization-prepared rapid acquisition gradient echo weighted fat-saturation imaging, plaque components can be defined. Post-processing software allows for semi- and fully automated quantitative analysis. Imaging correlation with surgical specimens suggests that this technique accurately differentiates plaque features. Total plaque burden and specific plaque components such as a thin fibrous cap, large fatty or necrotic core and intraplaque hemorrhage are accepted markers of neuroischemic events. Given the systemic nature of atherosclerosis, emerging science suggests that the presence of carotid plaque is also an indicator of coronary artery plaque burden, although the preliminary data primarily involves patients with stable coronary disease. While the availability and cost-effectiveness of MRI will ultimately be important determinants of whether carotid MRI is adopted clinically in cardiovascular risk assessment, the high accuracy and reliability of this technique suggests that it has potential as an imaging biomarker of future risk.

The effect of resolution on viscous dissipation measured with 4D flow MRI in patients with Fontan circulation: Evaluation using computational fluid dynamics

Viscous dissipation inside Fontan circulation, a parameter associated with the exercise intolerance of Fontan patients, can be derived from computational fluid dynamics (CFD) or 4D flow MRI velocities. However, the impact of spatial resolution and measurement noise on the estimation of viscous dissipation is unclear. Our aim was to evaluate the influence of these parameters on viscous dissipation calculation. Six Fontan patients underwent whole heart 4D flow MRI. Subject-specific CFD simulations were performed. The CFD velocities were down-sampled to isotropic spatial resolutions of 0.5mm, 1mm, 2mm and to MRI resolution. Viscous dissipation was compared between (1) high resolution CFD velocities, (2) CFD velocities down-sampled to MRI resolution, (3) down-sampled CFD velocities with MRI mimicked noise levels, and (4) in-vivo 4D flow MRI velocities. Relative viscous dissipation between subjects was also calculated. 4D flow MRI velocities (15.6 ± 3.8 cm/s) were higher, although not significantly different than CFD velocities (13.8 ± 4.7 cm/s, p=0.16), down-sampled CFD velocities (12.3 ± 4.4 cm/s, p=0.06) and the down-sampled CFD velocities with noise (13.2 ± 4.2 cm/s, p=0.06). CFD-based viscous dissipation (0.81 ± 0.55 mW) was significantly higher than those based on down-sampled CFD (0.25 ± 0.19 mW, p=0.03), down-sampled CFD with noise (0.49 ± 0.26 mW, p=0.03) and 4D flow MRI (0.56 ± 0.28 mW, p=0.06). Nevertheless, relative viscous dissipation between different subjects was maintained irrespective of resolution and noise, suggesting that comparison of viscous dissipation between patients is still possible.

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.