In vivo hemodynamic analysis of intracranial aneurysms obtained by magnetic resonance fluid dynamics (MRFD) based on time-resolved three-dimensional phase-contrast MRI

Rupture-Risk Stratifying Patients with Cerebral Arteriovenous Malformations Using Quantitative Hemodynamic Flow Measurements

Arteriovenous malformations (AVMs) are high-pressure, low-resistance arterial-venous shunts without intervening capillaries. Up to 60% of AVMs present with an intracranial hemorrhage; however, noninvasive neuroimaging has increasingly diagnosed incidental AVMs. AVM management depends on weighing the lifetime rupture risk against the risks of intervention. Although AVM rupture risk relies primarily on angioarchitectural features, measuring hemodynamic flow is gaining traction. Accurate understanding of AVM hemodynamic flow parameters will help endovascular neurosurgeons and interventional neuroradiologists stratify patients by rupture risk and select treatment plans. This review examines various neuroimaging modalities and their capabilities to quantify AVM flow, as well as the relationship between AVM flow and rupture risk. Quantitative hemodynamic studies on the relationship between AVM flow and rupture risk have not reached a clear consensus; however, the preponderance of data suggests that higher arterial inflow and lower venous outflow in the AVM nidus contribute to increased hemorrhagic risk. Future studies should consider using larger sample sizes and standardized definitions of hemodynamic parameters to reach a consensus. In the meantime, classic angioarchitectural features may be more strongly correlated with AVM rupture than the amount of blood flow.

Automatic 4D Flow MRI Segmentation Using the Standardized Difference of Means Velocity

We present a method to automatically segment 4D flow magnetic resonance imaging (MRI) by identifying net flow effects using the standardized difference of means (SDM) velocity. The SDM velocity quantifies the ratio between the net flow and observed flow pulsatility in each voxel. Vessel segmentation is performed using an F-test, identifying voxels with significantly higher SDM velocity values than background voxels. We compare the SDM segmentation algorithm against pseudo-complex difference (PCD) intensity segmentation of 4D flow measurements in in vitro cerebral aneurysm models and 10 in vitro Circle of Willis (CoW) datasets. We also compared the SDM algorithm to convolutional neural network (CNN) segmentation in 5 thoracic vasculature datasets. The in vitro flow phantom geometry is known, while the ground truth geometries for the CoW and thoracic aortas are derived from high-resolution time-of-flight (TOF) magnetic resonance angiography and manual segmentation, respectively. The SDM algorithm demonstrates greater robustness than PCD and CNN approaches and can be applied to 4D flow data from other vascular territories. The SDM to PCD comparison demonstrated an approximate 48% increase in sensitivity in vitro and 70% increase in the CoW, respectively; the SDM and CNN sensitivities were similar. The vessel surface derived from the SDM method was 46% closer to the in vitro surfaces and 72% closer to the in vitro TOF surfaces than the PCD approach. The SDM and CNN approaches both accurately identify vessel surfaces. The SDM algorithm is a repeatable segmentation method, enabling reliable computation of hemodynamic metrics associated with cardiovascular disease.

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.

Why 4D Flow MRI? Real Advantages

This special issue of Magnetic Resonance in Medical Sciences features the most recent reviews on 4D Flow MRI. These reviews deal with the current status of the emerging technique of 4D Flow MRI facilitated in various areas that are difficult to obtain with conventional flowmetry. MR signals inherently contain flow velocity information. In previous decades, in vivo blood flow measurement was traditionally performed by 2D methods, such as Doppler ultrasonography and 2D phase-contrast MRI, which have long been regarded as mature techniques in hemodynamic flowmetry. Although 2D velocimetries have many advantages over 4D Flow MRI in terms of cost and accessibility, and provide excellent temporal and in-plane spatial resolutions, they also have some disadvantages. The emerging technology of 4D Flow MRI can overcome the shortcomings of conventional 2D imaging. In recent years, hemodynamic analysis has witnessed significant progress that is primarily attributable to advances in 4D Flow MRI.

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.

Dynamic modes of inflow jet in brain aneurysms

The transition of the inflow jet to turbulence is crucial in understanding the pathology of brain aneurysms. Previous works Le et al. (2010, 2013) have shown evidence for a highly dynamic inflow jet in the ostium of brain aneurysms. While it is highly desired to investigate this inflow jet dynamics in clinical practice, the constraints on spatial and temporal resolutions of in vivo data do not allow a detailed analysis of this transition. In this work, Dynamic Mode Decomposition (DMD) is used to identify the most energetic modes of the inflow jet in patient-specific models of internal carotid aneurysms via the utilization of high-resolution simulation data. It is hypothesized that dynamic modes are not solely controlled by the blood flow waveform at the parent artery. They are also dependent on jet-wall interaction phenomena. DMD analysis shows that the spatial extent of low- frequency modes corresponds well to the most energetic areas of the inflow jet. The high-frequency modes are short-lived and correspond to the flow separation at the proximal neck and the jet's impingement onto the aneurysmal wall. Low-frequency modes can be reconstructed at relatively low spatial and temporal resolutions comparable to ones of in vivo data. The current results suggest that DMD can be practically useful in analyzing blood flow patterns of brain aneurysms with in vivo data.

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.

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.

What does computational fluid dynamics tell us about intracranial aneurysms? A meta-analysis and critical review

Despite the plethora of published studies on intracranial aneurysms (IAs) hemodynamic using computational fluid dynamics (CFD), limited progress has been made towards understanding the complex physics and biology underlying IA pathophysiology. Guided by 1733 published papers, we review and discuss the contemporary IA hemodynamics paradigm established through two decades of IA CFD simulations. We have traced the historical origins of simplified CFD models which impede the progress of comprehending IA pathology. We also delve into the debate concerning the Newtonian fluid assumption used to represent blood flow computationally. We evidently demonstrate that the Newtonian assumption, used in almost 90% of studies, might be insufficient to describe IA hemodynamics. In addition, some fundamental properties of the Navier-Stokes equation are revisited in supplementary material to highlight some widely spread misconceptions regarding wall shear stress (WSS) and its derivatives. Conclusively, our study draws a roadmap for next-generation IA CFD models to help researchers investigate the pathophysiology of IAs.

Factors influencing blood flow resistance from a large internal carotid artery aneurysm revealed by a computational fluid dynamics model

Hyperperfusion syndrome occurs after treatment of a large or giant cerebral aneurysm. Recently, flow-diverter stent placement has emerged as an effective treatment method for a large cerebral aneurysm, but postoperative ipsilateral delayed intraparenchymal hemorrhage occurs in a minority of cases. The mechanism underlying delayed intraparenchymal hemorrhage is not established, but one possibility is hyperperfusion syndrome. The incidence of delayed intraparenchymal hemorrhage appears to be higher for giant aneurysms; hence, we speculated that large/giant aneurysms may create flow resistance, and mitigation by flow-diverter stent deployment leads to hyperperfusion syndrome and delayed intraparenchymal hemorrhage. The purpose of this study was to identify aneurysm characteristics promoting flow resistance by the analysis of pressure loss in an internal carotid artery paraclinoid aneurysm model using computational fluid dynamics. A virtual U-shaped model of the internal carotid artery siphon portion was created with a spherical aneurysm of various angles, body diameters, and neck diameters. Visualization of streamlines, were calculated of pressure loss between proximal and distal sides of the aneurysm, and vorticity within the aneurysm were calculated. The pressure loss and vorticity demonstrated similar changes according to angle, peaking at 60°. In contrast, aneurysm diameter had little influence on pressure loss. Larger neck width, however, increases pressure loss. Our model predicts that aneurysm location and neck diameter can increase the flow resistance from a large internal carotid artery aneurysm. Patients with large aneurysm angles and neck diameters may be at increased risk of hyperperfusion syndrome and ensuing delayed intraparenchymal hemorrhage following flow-diverter stent treatment.

Hemodynamic Characteristics Associated With Paraclinoid Aneurysm Recurrence in Patients After Embolization

Objective: To investigate the hemodynamic features before and after embolization of paraclinoidal aneurysms using hemodynamic numerical simulation and the influence of embolization on recurrence after embolization.

Methods: From January 2016 to December 2017, we enrolled a total of 113 paraclinoidal aneurysms treated with embolization. They were divided into recurrent group and stable group depending on follow-up results. An aneurysm model was generated based on 3D-DSA before and after embolization. The hemodynamic characteristics were analyzed between two groups using Computational fluid dynamic (CFD).

Results: In the recurrent group, the peak systolic WSS, OSI and velocity around the aneurysm neck areas prior to embolization were 20.47 ± 3.04 Pa, 0.06 ± 0.02 and 0.07 ± 0.03 m/s, respectively. These values were 23.50 ± 4.11 Pa, 0.06 ± 0.01 and 0.11 ± 0.02 m/s, respectively in the stable group (P > 0.05). The WSS, OSI, velocity around the same areas in the recurrent group after embolization were 35.59 ± 8.75 Pa, 0.07 ± 0.02 and 0.12 ± 0.03 m/s, respectively (P < 0.01). In the stable group, the WSS, OSI and velocity were 13.08 ± 2.89 Pa, 0.04 ± 0.01 and 0.07 ± 0.02 m/s, respectively (P < 0.01). After embolization, the WSS, OSI and velocity around the aneurysm neck areas in the recurrent group were significantly higher than those in the stable group.

Conclusions: High peak systolic WSS, OSI and velocity around aneurysm neck areas after embolization of paraclinoidal aneurysms may be important factors leading to recurrence.

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.

Novel and facile criterion to assess the accuracy of WSS estimation by 4D flow MRI

Four-dimensional flow magnetic resonance imaging (4D flow MRI) is a versatile tool to obtain hemodynamic information and anatomic information simultaneously. The wall shear stress (WSS), a force exerted on a vessel wall in parallel, is one of the hemodynamic parameters available with 4D flow MRI and is thought to play an important role in clinical applications such as assessing the development of atherosclerosis. Nevertheless, the accuracy of WSS obtained with 4D flow MRI is rarely evaluated or reported in literature, especially in the in vivo studies. We propose a novel and facile criterion called Reynolds resolution to assess the accuracy of WSS estimation in 4D flow MRI studies. Reynolds resolution consists of a spatial resolution, encoding velocity, kinematic viscosity of a working fluid, and signal-to-noise ratio, which are readily accessible information in 4D flow MRI measurements. We explored the relationship between Reynolds resolution and the WSS error. To include diverse and extensive cases, we measured three circular tubing flows with a diameter of 40, 8, and 2 mm. The 40 mm tubing flow was measured by 3 Tesla (T) human MR scanner with a knee coil and spatial resolution of 0.5 mm. The 8 and 2 mm tubing flows were both measured by 4.7 T MR scanner, but the scans were performed with a conventional birdcage coil (8 mm tubing) and a custom-made solenoid coil (2 mm tubing), respectively. The spatial resolution was varied from 0.2, 0.4 or 0.8 mm for the 8 mm tubing flow, but was fixed at 0.090 mm for 2 mm tubing flow. In addition, the near-wall velocity gradient, required to be determined prior to the WSS, was calculated using two methods; these included assuming a linear velocity profile or quadratic velocity profile near wall. The accuracy of WSS obtained using each method and tubing flow was evaluated against the theoretical WSS value. As a result, we found that Reynolds resolution is in logarithmic relation to the WSS error.

Highly accelerated intracranial 4D flow MRI: evaluation of healthy volunteers and patients with intracranial aneurysms

OBJECTIVES

To evaluate an accelerated 4D flow MRI method that provides high temporal resolution in a clinically feasible acquisition time for intracranial velocity imaging.

MATERIALS AND METHODS

Accelerated 4D flow MRI was developed by using a pseudo-random variable-density Cartesian undersampling strategy (CIRCUS) with the combination of k-t, parallel imaging and compressed sensing image reconstruction techniques (k-t SPARSE-SENSE). Four-dimensional flow data were acquired on five healthy volunteers and eight patients with intracranial aneurysms using CIRCUS (acceleration factor of R = 4, termed CIRCUS4) and GRAPPA (R = 2, termed GRAPPA2) as the reference method. Images with three times higher temporal resolution (R = 12, CIRCUS12) were also reconstructed from the same acquisition as CIRCUS4. Qualitative and quantitative image assessment was performed on the images acquired with different methods, and complex flow patterns in the aneurysms were identified and compared.

RESULTS

Four-dimensional flow MRI with CIRCUS was achieved in 5 min and allowed further improved temporal resolution of <30 ms. Volunteer studies showed similar qualitative and quantitative evaluation obtained with the proposed approach compared to the reference (overall image scores: GRAPPA2 3.2 ± 0.6; CIRCUS4 3.1 ± 0.7; CIRCUS12 3.3 ± 0.4; difference of the peak velocities: -3.83 ± 7.72 cm/s between CIRCUS4 and GRAPPA2, -1.72 ± 8.41 cm/s between CIRCUS12 and GRAPPA2). In patients with intracranial aneurysms, the higher temporal resolution improved capturing of the flow features in intracranial aneurysms (pathline visualization scores: GRAPPA2 2.2 ± 0.2; CIRCUS4 2.5 ± 0.5; CIRCUS12 2.7 ± 0.6).

CONCLUSION

The proposed rapid 4D flow MRI with a high temporal resolution is a promising tool for evaluating intracranial aneurysms in a clinically feasible acquisition time.

Hemodynamic vascular biomarkers for initiation of paraclinoid internal carotid artery aneurysms using patient-specific computational fluid dynamic simulation based on magnetic resonance imaging

PURPOSE

We performed computational fluid dynamics (CFD) for patients with and without paraclinoid internal carotid artery (ICA) aneurysms to evaluate the distribution of vascular biomarkers at the aneurysm initiation sites of the paraclinoid ICA.

METHODS

This study included 35 patients who were followed up for aneurysms using 3D time of flight (TOF) magnetic resonance angiography (MRA) and 3D cine phase-contrast MR imaging. Fifteen affected ICAs were included in group A with the 15 unaffected contralateral ICAs in group B. Thirty-three out of 40 paraclinoid ICAs free of aneurysms and arteriosclerotic lesions were included in group C. We deleted the aneurysms in group A based on the 3D TOF MRA dataset. We performed CFD based on MR data set and obtained wall shear stress (WSS), its derivatives, and streamlines. We qualitatively evaluated their distributions at and near the intracranial aneurysm initiation site among three groups. We also calculated and compared the normalized highest (nh-) WSS and nh-spatial WSS gradient (SWSSG) around the paraclinoid ICA among three groups.

RESULTS

High WSS and SWSSG distribution were observed at and near the aneurysm initiation site in group A. High WSS and SWSSG were also observed at similar locations in group B and group C. However, nh-WSS and nh-SWSSG were significantly higher in group A than in group C, and nh-SWSSG was significantly higher in group A than in group B.

CONCLUSION

Our findings indicated that nh-WSS and nh-SWSSG were good biomarkers for aneurysm initiation in the paraclinoid ICA.

Accelerated dual-venc 4D flow MRI for neurovascular applications

PURPOSE

To improve velocity-to-noise ratio (VNR) and dynamic velocity range of 4D flow magnetic resonance imaging (MRI) by using dual-velocity encoding (dual-venc) with k-t generalized autocalibrating partially parallel acquisition (GRAPPA) acceleration.

MATERIALS AND METHODS

A dual-venc 4D flow MRI sequence with k-t GRAPPA acceleration was developed using a shared reference scan followed by three-directional low- and high-venc scans (repetition time / echo time / flip angle = 6.1 msec / 3.4 msec / 15°, temporal/spatial resolution = 43.0 msec/1.2 × 1.2 × 1.2 mm³ ). The high-venc data were used to correct for aliasing in the low-venc data, resulting in a single dataset with the favorable VNR of the low-venc but without velocity aliasing. The sequence was validated with a 3T MRI scanner in phantom experiments and applied in 16 volunteers to investigate its feasibility for assessing intracranial hemodynamics (net flow and peak velocity) at the major intracranial vessels. In addition, image quality and image noise were assessed in the in vivo acquisitions.

RESULTS

All 4D flow MRI scans were acquired successfully with an acquisition time of 20 ± 4 minutes. The shared reference scan reduced the total acquisition time by 12.5% compared to two separate scans. Phantom experiments showed 51.4% reduced noise for dual-venc compared to high-venc and an excellent agreement of velocities (ρ = 0.8, P < 0.001). The volunteer data showed decreased noise in dual-venc data (54.6% lower) compared to high-venc, and improved image quality, as graded by two observers: fewer artifacts (P < 0.0001), improved vessel conspicuity (P < 0.0001), and reduced noise (P < 0.0001).

CONCLUSION

Dual-venc 4D flow MRI exhibits the superior VNR of the low-venc acquisition and reliably incorporates low- and high-velocity fields simultaneously. In vitro and in vivo data demonstrate improved flow visualization, image quality, and image noise.

LEVEL OF EVIDENCE

2 Technical Efficacy: Stage 1 J. MAGN. RESON. IMAGING 2017;46:102-114.

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.

Spatial phenotyping of the endocardial endothelium as a function of intracardiac hemodynamic shear stress

Despite substantial evidence for the central role of hemodynamic shear stress in the functional integrity of vascular endothelial cells, hemodynamic and molecular regulation of the endocardial endothelium lining the heart chambers remains understudied. We propose that regional differences in intracardiac hemodynamics influence differential endocardial gene expression leading to phenotypic heterogeneity of this cell layer. Measurement of intracardiac hemodynamics was performed using 4-dimensional flow MRI in healthy humans (n=8) and pigs (n=5). Local wall shear stress (WSS) and oscillatory shear indices (OSI) were calculated in three distinct regions of the LV - base, mid-ventricle (midV), and apex. In both the humans and pigs, WSS values were significantly lower in the apex and midV relative to the base. Additionally, both the apex and midV had greater oscillatory shear indices (OSI) than the base. To investigate regional phenotype, endocardial endothelial cells (EEC) were isolated from an additional 8 pigs and RNA sequencing was performed. A false discovery rate of 0.10 identified 1051 differentially expressed genes between the base and apex, and 321 between base and midV. Pathway analyses revealed apical upregulation of genes associated with translation initiation. Furthermore, tissue factor pathway inhibitor (TFPI; mean 50-fold) and prostacyclin synthase (PTGIS; 5-fold), genes prominently associated with antithrombotic protection, were consistently upregulated in LV apex. These spatio-temporal WSS values in defined regions of the left ventricle link local hemodynamics to regional heterogeneity in endocardial gene expression.

Four-dimensional MRI flow examinations in cerebral and extracerebral vessels - ready for clinical routine?

PURPOSE OF REVIEW

To evaluate the feasibility of 4-dimensional (4D) flow MRI for the clinical assessment of cerebral and extracerebral vascular hemodynamics in patients with neurovascular disease.

RECENT FINDINGS

4D flow MRI has been applied in multiple studies to qualitatively and quantitatively study intracranial aneurysm blood flow for potential risk stratification and to assess treatment efficacy of various neurovascular lesions, including intraaneurysmal and parent artery blood flow after flow diverter stent placement and staged embolizations of arteriovenous malformations and vein of Galen aneurysmal malformations. Recently, the technique has been utilized to characterize age-related changes of normal cerebral hemodynamics in healthy individuals over a broad age range.

SUMMARY

4D flow MRI is a useful tool for the noninvasive, volumetric and quantitative hemodynamic assessment of neurovascular disease without the need for gadolinium contrast agents. Further improvements are warranted to overcome technical limitations before broader clinical implementation. Current developments, such as advanced acceleration techniques (parallel imaging and compressed sensing) for faster data acquisition, dual or multiple velocity encoding strategies for more accurate arterial and venous flow quantification, ultrahigh-field strengths to achieve higher spatial resolution and streamlined postprocessing workflow for more efficient and standardized flow analysis, are promising advancements in 4D flow MRI.

Hemodynamic Measurement Using Four-Dimensional Phase-Contrast MRI: Quantification of Hemodynamic Parameters and Clinical Applications

Recent improvements have been made to the use of time-resolved, three-dimensional phase-contrast (PC) magnetic resonance imaging (MRI), which is also named four-dimensional (4D) PC-MRI or 4D flow MRI, in the investigation of spatial and temporal variations in hemodynamic features in cardiovascular blood flow. The present article reviews the principle and analytical procedures of 4D PC-MRI. Various fluid dynamic biomarkers for possible clinical usage are also described, including wall shear stress, turbulent kinetic energy, and relative pressure. Lastly, this article provides an overview of the clinical applications of 4D PC-MRI in various cardiovascular regions.

Predictive modeling and in vivo assessment of cerebral blood flow in the management of complex cerebral aneurysms

Cerebral aneurysms are weakened blood vessel dilatations that can result in spontaneous, devastating hemorrhage events. Aneurysm treatment aims to reduce hemorrhage events, and strategies for complex aneurysms often require surgical bypass or endovascular stenting for blood flow diversion. Interventions that divert blood flow from their normal circulation patterns have the potential to result in unintentional ischemia. Recent developments in computational modeling and in vivo assessment of hemodynamics for cerebral aneurysm treatment have entered into clinical practice. Herein, we review how these techniques are currently utilized to improve risk stratification and treatment planning.

Intra-left ventricular flow dynamics in patients with preserved and impaired left ventricular function: Analysis with 3D cine phase contrast MRI (4D-Flow)

PURPOSE

To examine how left ventricular (LV) volume and function affect flow dynamics by analyzing 3D intra-LV vortex features using 4D-Flow.

MATERIALS AND METHODS

Twenty-one patients with preserved (LVEF > 60%) and 14 with impaired LV function (LVEF < 40%) underwent 4D-Flow (at 3T).

RESULTS

In patients with preserved LV function, the intra-LV vortices developed in both the early and late diastolic phases. The shift of inflow vectors at the basal LV toward the posterior-lateral side of the LV and the mid-ventricular turn of inflow vectors toward the LV outflow could explain clearer vortex formation in the late diastolic phase. In patients with impaired LV function, the intra-LV vortices during the diastolic phase located at the more apical LV were larger and more spherically shaped. Both the distance to the vortex core and the vortex area correlated significantly with LV end-diastolic volume (r = 0.66 and 0.73), LVEF (r = -0.74 and -0.68), LV sphericity index (r = -0.60 and -0.65), and peak filling rate (r = -0.61 and -0.64), respectively (P < 0.01). The intra-LV vortices developed during the systolic phase in 10 cases. In those, some of the particles at the apical LV rotated within the LV, whereas in patients with preserved LV function, all of the particles were directed straight to the ascending aorta with accelerated flow velocity (256.8 ± 120.2 cm/s vs. 414.3 ± 88.2 cm/s, P < 0.01).

CONCLUSION

Vortex formation during the diastolic phase may be critical for both LV filling and ejection. 4D-Flow showed the 3D alterations of intra-LV flow dynamics by LV dilatation and dysfunction in a noninvasive and comprehensive manner. J. Magn. Reson. Imaging 2016;44:1493-1503.

Integrated Regional Cardiac Hemodynamic Imaging and RNA Sequencing Reveal Corresponding Heterogeneity of Ventricular Wall Shear Stress and Endocardial Transcriptome

Background

Unlike arteries, in which regionally distinct hemodynamics are associated with phenotypic heterogeneity, the relationships between endocardial endothelial cell phenotype and intraventricular flow remain largely unexplored. We investigated regional differences in left ventricular wall shear stress and their association with endocardial endothelial cell gene expression.

Methods and Results

Local wall shear stress was calculated from 4‐dimensional flow magnetic resonance imaging in 3 distinct regions of human (n=8) and pig (n=5) left ventricle: base, adjacent to the outflow tract; midventricle; and apex. In both species, wall shear stress values were significantly lower in the apex and midventricle relative to the base; oscillatory shear index was elevated in the apex. RNA sequencing of the endocardial endothelial cell transcriptome in pig left ventricle (n=8) at a false discovery rate ≤10% identified 1051 genes differentially expressed between the base and the apex and 327 between the base and the midventricle; no differentially expressed genes were detected at this false discovery rate between the apex and the midventricle. Enrichment analyses identified apical upregulation of genes associated with translation initiation including mammalian target of rapamycin, and eukaryotic initiation factor 2 signaling. Genes of mitochondrial dysfunction and oxidative phosphorylation were also consistently upregulated in the left ventricular apex, as were tissue factor pathway inhibitor (mean 50‐fold) and prostacyclin synthase (5‐fold)—genes prominently associated with antithrombotic protection.

Conclusions

We report the first spatiotemporal measurements of wall shear stress within the left ventricle and linked regional hemodynamics to heterogeneity in ventricular endothelial gene expression, most notably to translation initiation and anticoagulation properties in the left ventricular apex, in which oscillatory shear index is increased and wall shear stress is decreased.

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.

4D flow cardiovascular magnetic resonance consensus statement

Pulsatile blood flow through the cavities of the heart and great vessels is time-varying and multidirectional. Access to all regions, phases and directions of cardiovascular flows has formerly been limited. Four-dimensional (4D) flow cardiovascular magnetic resonance (CMR) has enabled more comprehensive access to such flows, with typical spatial resolution of 1.5×1.5×1.5 - 3×3×3 mm(3), typical temporal resolution of 30-40 ms, and acquisition times in the order of 5 to 25 min. This consensus paper is the work of physicists, physicians and biomedical engineers, active in the development and implementation of 4D Flow CMR, who have repeatedly met to share experience and ideas. The paper aims to assist understanding of acquisition and analysis methods, and their potential clinical applications with a focus on the heart and greater vessels. We describe that 4D Flow CMR can be clinically advantageous because placement of a single acquisition volume is straightforward and enables flow through any plane across it to be calculated retrospectively and with good accuracy. We also specify research and development goals that have yet to be satisfactorily achieved. Derived flow parameters, generally needing further development or validation for clinical use, include measurements of wall shear stress, pressure difference, turbulent kinetic energy, and intracardiac flow components. The dependence of measurement accuracy on acquisition parameters is considered, as are the uses of different visualization strategies for appropriate representation of time-varying multidirectional flow fields. Finally, we offer suggestions for more consistent, user-friendly implementation of 4D Flow CMR acquisition and data handling with a view to multicenter studies and more widespread adoption of the approach in routine clinical investigations.

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.

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.

Characteristics of intra-left atrial flow dynamics and factors affecting formation of the vortex flow – analysis with phase-resolved 3-dimensional cine phase contrast magnetic resonance imaging

BACKGROUND

The intra-left atrial (LA) blood flow from pulmonary veins (PVs) to the left ventricle (LV) changes under various conditions and might affect global cardiac function. By using phase-resolved 3-dimensional cine phase contrast magnetic resonance imaging (4D-Flow), the intra-LA vortex formation was visualized and the factors affecting the intra-LA flow dynamics were examined.

METHODS AND RESULTS

Thirty-two patients with or without organic heart diseases underwent 4D-Flow and transthoracic echocardiography. The intra-LA velocity vectors from each PV were post-processed to delineate streamline and pathline images. The vector images revealed intra-LA vortex formation in 20 of 32 patients. All the vortices developed during the late systolic and early diastolic phases and were directed counter-clockwise when viewed from the subjects' cranial side. The flow vectors from the right PVs lengthened predominantly toward the mitral valves and partly toward the LA appendage, whereas those from the left PVs directed rightward along the posterior wall and joined the vortex. Patients with vortex had less organic heart diseases, smaller LV and LA volume, and greater peak flow velocity and volume mainly in the left PVs, although the flow directions from each PV or PV areas did not differ.

CONCLUSIONS

4D-Flow can clearly visualize the intra-LA vortex formation and analyze its characteristic features. The vortex formation might depend on LV and LA volume and on flow velocity and volume from PVs.

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.

Identification of the inflow zone of unruptured cerebral aneurysms: comparison of 4D flow MRI and 3D TOF MRA data

BACKGROUND AND PURPOSE

The hemodynamics of the inflow zone of cerebral aneurysms may be a key factor in coil compaction and recanalization after endovascular coil embolization. We performed 4D flow MR imaging in conjunction with 3D TOF MRA and compared their ability to identify the inflow zone of unruptured cerebral aneurysms.

MATERIALS AND METHODS

This series comprised 50 unruptured saccular cerebral aneurysms in 44 patients. Transluminal color-coded 3D MRA images were created by selecting the signal-intensity ranges on 3D TOF MRA images that corresponded with both the luminal margin and the putative inflow.

RESULTS

4D flow MR imaging demonstrated the inflow zone and yielded inflow velocity profiles for all 50 aneurysms. In 18 of 24 lateral-projection aneurysms (75%), the inflow zone was located distally on the aneurysmal neck. The maximum inflow velocity ranged from 285 to 922 mm/s. On 4D flow MR imaging and transluminal color-coded 3D MRA studies, the inflow zone of 32 aneurysms (64%) was at a similar location. In 91% of aneurysms whose neck section plane angle was <30° with respect to the imaging section direction on 3D TOF MRA, depiction of the inflow zone was similar on transluminal color-coded 3D MRA and 4D flow MR images.

CONCLUSIONS

4D flow MR imaging can demonstrate the inflow zone and provide inflow velocity profiles. In aneurysms whose angle of the neck-section plane is obtuse vis-a-vis the imaging section on 3D TOF MRA scans, transluminal color-coded 3D MRA may depict the inflow zone reliably.

Three-dimensional hemodynamics in intracranial aneurysms: influence of size and morphology

PURPOSE

To use four-dimensional (4D)-flow MRI for the comprehensive in vivo analysis of hemodynamics and its relationship to size and morphology of different intracranial aneurysms (IA). We hypothesize that different IA groups, defined by size and morphology, exhibit different velocity fields, wall shear stress, and vorticity.

MATERIALS AND METHODS

The 4D-flow MRI (spatial resolution = 0.99-1.8 × 0.78-1.46 × 1.2-1.4 mm(3) , temporal resolution = 44-48 ms) was performed in 19 IAs (18 patients, age = 55.4 ± 13.8 years) with saccular (n = 16) and fusiform (n = 3) morphology and different sizes ranging from small (n = 8; largest dimension = 6.2 ± 0.4 mm) to large and giant (n = 11; 25 ± 7 mm). Analysis included quantification of volumetric spatial-temporal velocity distribution, vorticity, and wall shear stress (WSS) along the aneurysm's 3D surface.

RESULTS

The 4D-flow MRI revealed distinct hemodynamic patterns for large/giant saccular aneurysms (Group 1), small saccular aneurysms (Group 2), and large/giant fusiform aneurysms (Group 3). Saccular IA (Groups 1, 2) demonstrated significantly higher peak velocities (P < 0.002) and WSS (P < 0.001) compared with fusiform aneurysms. Although intra-aneurysmal 3D velocity distributions were similar for Group 1 and 2, vorticity and WSS was significantly (P < 0.001) different (increased in Group 1 by 54%) indicating a relationship between IA size and hemodynamics. Group 3 showed reduced velocities (P < 0.001) and WSS (P < 0.001).

CONCLUSION

The 4D-flow MRI demonstrated the influence of lesion size and morphology on aneurysm hemodynamics suggesting the potential of 4D-flow MRI to assist in the classification of individual aneurysms.

Accurate determination of patient-specific boundary conditions in computational vascular hemodynamics using 3D cine phase-contrast MRI

In the patient-specific vascular CFD, determination of the inlet and outlet boundary conditions (BCs) is an important issue for a valid diagnosis. The 3D cine phase-contrast MRI (4D Flow) velocimetry is promising for this issue; yet, its measured velocities contain relatively large error and are not admissible as the BCs without any correction. This paper proposes a novel correction method for determining the BCs accurately using the 4D Flow velocimetry. First, we reveal that the error of the velocity measured by the 4D Flow at each measurement voxel is large but is distributed symmetrically. Secondly, our method pays attention to the incompressibility of the blood and the fact that the volume flow rate (VFR) in each vessel is constant on any cross sections. We reveal that the average of the cross-sectional VFRs integrated from many measurement voxel in each vessel is accurate despite the large error. Finally, we propose the novel correction method, which applies a smoothing to the measured velocities on each inlet or outlet boundary with a low-pass filter and then corrects them with the VFR. The results of the several phantom studies are presented to validate the accuracy of our method. A demonstrative analysis for an actual aneurysm is also presented to show the feasibility and effectiveness of our method.

Image-Based Evaluation of Vascular Function and Hemodynamics

The noticeable characteristics of the blood vascular structure are the inconsistent viscosity of blood and the stiffness of the vascular wall. If we can control these two factors, we can solve more problems related to hemodynamics and vascular wall function. Understanding the properties of hemodynamics and vascular wall function may provide more information applicable to clinical practice for cardiovascular disease. The bedside techniques evaluating vascular function usually measure indirect parameters. In contrast, some medical imaging techniques provide clear and direct depictions of functional cardiovascular characteristics. In this review, image-based evaluation of hemodynamic and vascular wall functions is discussed from the perspective of blood flow velocity, flow volume, flow pattern, peripheral vascular resistance, intraluminal pressure, vascular wall stress, and wall stiffness.

Intracranial 4D flow MRI: toward individualized assessment of arteriovenous malformation hemodynamics and treatment-induced changes

BACKGROUND AND PURPOSE

Arteriovenous malformations are an important etiology of hemorrhagic stroke. However, current imaging modalities and risk do not provide insights into individual AVM hemodynamics and its role in pathophysiology. The aims of this study are to determine whether intracranial 4D flow MR imaging can provide insights into arteriovenous malformation hemodynamics independent of the Spetzler-Martin grade and to report the changes in flow observed during staged embolization.

MATERIALS AND METHODS

Intracranial 3D blood flow was assessed in 20 patients with AVM (age = 39 ± 15 years, Spetzler-Martin grade ranging from 1-4) with the use of 4D flow MR imaging (temporal resolution = 45 ms, spatial resolution = [1.2-1.6mm](3)). AVM hemodynamics were visualized by means of time-integrated 3D pathlines depicting the AVM arterial feeding and venous draining patterns over the cardiac cycle. Analysis included the grading of feeding and draining velocities on a 3-point scale (0 = low <25 cm/s, 1 = medium <50 cm/s, 2 = high >50 cm/s). For 4 of 20 patients undergoing 4D flow MR imaging follow-up after staged embolization, peak velocities were quantified in arterial feeders, draining veins, the sagittal sinus, and contralateral arteries.

RESULTS

In 50% of the cases with Spetzler-Martin grade >2, heterogeneous flow (velocity grade differences >1) was found across arteries and veins. Velocities in draining veins increased from Spetzler-Martin grade = 1 (grading = 0.5 ± 0.6) to Spetzler-Martin grade ≥3 (1.1 ± 0.6), whereas arterial velocities were similar (1.7 ± 0.6 versus 1.5 ± 0.6). In the postembolization subgroup of 4 patients, 4D flow MR imaging demonstrated successively more compact AVM and redistribution of velocities. Changes in arterial and venous velocities during treatment were highly different among individuals.

CONCLUSIONS

Spetzler-Martin grade does not reflect differences in 3D AVM arterial and venous hemodynamics, and an individual assessment of AVM hemodynamics may be needed for improved lesion characterization. Four-dimensional flow MR imaging may have the potential to monitor and guide embolization treatment planning.

Physical factors effecting cerebral aneurysm pathophysiology

Many factors that are either blood-, wall-, or hemodynamics-borne have been associated with the initiation, growth, and rupture of intracranial aneurysms. The distribution of cerebral aneurysms around the bifurcations of the circle of Willis has provided the impetus for numerous studies trying to link hemodynamic factors (flow impingement, pressure, and/or wall shear stress) to aneurysm pathophysiology. The focus of this review is to provide a broad overview of such hemodynamic associations as well as the subsumed aspects of vascular anatomy and wall structure. Hemodynamic factors seem to be correlated to the distribution of aneurysms on the intracranial arterial tree and complex, slow flow patterns seem to be associated with aneurysm growth and rupture. However, both the prevalence of aneurysms in the general population and the incidence of ruptures in the aneurysm population are extremely low. This suggests that hemodynamic factors and purely mechanical explanations by themselves may serve as necessary, but never as necessary and sufficient conditions of this disease's causation. The ultimate cause is not yet known, but it is likely an additive or multiplicative effect of a handful of biochemical and biomechanical factors.

4D flow MRI

Traditionally, magnetic resonance imaging (MRI) of flow using phase contrast (PC) methods is accomplished using methods that resolve single-directional flow in two spatial dimensions (2D) of an individual slice. More recently, three-dimensional (3D) spatial encoding combined with three-directional velocity-encoded phase contrast MRI (here termed 4D flow MRI) has drawn increased attention. 4D flow MRI offers the ability to measure and to visualize the temporal evolution of complex blood flow patterns within an acquired 3D volume. Various methodological improvements permit the acquisition of 4D flow MRI data encompassing individual vascular structures and entire vascular territories such as the heart, the adjacent aorta, the carotid arteries, abdominal, or peripheral vessels within reasonable scan times. To subsequently analyze the flow data by quantitative means and visualization of complex, three-directional blood flow patterns, various tools have been proposed. This review intends to introduce currently used 4D flow MRI methods, including Cartesian and radial data acquisition, approaches for accelerated data acquisition, cardiac gating, and respiration control. Based on these developments, an overview is provided over the potential this new imaging technique has in different parts of the body from the head to the peripheral arteries.

Experimental insights into flow impingement in cerebral aneurysm by stereoscopic particle image velocimetry: transition from a laminar regime

This study experimentally investigated the instability of flow impingement in a cerebral aneurysm, which was speculated to promote the degradation of aneurysmal wall. A patient-specific, full-scale and elastic-wall replica of cerebral artery was fabricated from transparent silicone rubber. The geometry of the aneurysm corresponded to that found at 9 days before rupture. The flow in a replica was analysed by quantitative flow visualization (stereoscopic particle image velocimetry) in a three-dimensional, high-resolution and time-resolved manner. The mid-systolic and late-diastolic flows with a Reynolds number of 450 and 230 were compared. The temporal and spatial variations of near-wall velocity at flow impingement delineated its inherent instability at a low Reynolds number. Wall shear stress (WSS) at that site exhibited a combination of temporal fluctuation and spatial divergence. The frequency range of fluctuation was found to exceed significantly that of the heart rate. The high-frequency-fluctuating WSS appeared only during mid-systole and disappeared during late diastole. These results suggested that the flow impingement induced a transition from a laminar regime. This study demonstrated that the hydrodynamic instability of shear layer could not be neglected even at a low Reynolds number. No assumption was found to justify treating the aneurysmal haemodynamics as a fully viscous laminar flow.

Wall shear stress estimated with phase contrast MRI in an in vitro and in vivo intracranial aneurysm

PURPOSE

To evaluate wall shear stress (WSS) estimations in an in vitro and in vivo intracranial aneurysm, WSS was estimated from phase contrast magnetic resonance imaging (PC-MRI) and compared with computational fluid dynamics (CFD).

MATERIALS AND METHODS

First, WSS was estimated using a high-resolution in vitro PC-MRI measurement under steady and pulsatile flow conditions and compared with CFD simulations. Second, WSS was estimated in steady PC-MRI data acquired at different spatial resolutions. Third, WSS estimations in pulsatile in vivo data were compared with CFD. The direction and magnitude of WSS vectors were computed and compared.

RESULTS

Quantitative agreement between PC-MRI and CFD-based WSS estimations was moderate for the phantom (Spearman ρ = 0.69). The WSS magnitude derived from PC-MRI data was lower than CFD for both the in vitro and in vivo case. However, there was qualitative agreement between PC-MRI and CFD, i.e. WSS vector direction was similar for both modalities. Circular WSS patterns were found both in vitro and in vivo for PC-MRI and CFD. Increasing PC-MRI resolution increased mean WSS magnitude and uncovered complex WSS patterns.

CONCLUSION

WSS patterns can be estimated based on PC-MRI data in in vitro and in vivo aneurysm geometries. Similar WSS directions as CFD can be discerned.

Left atrial flow velocity distribution and flow coherence using four-dimensional FLOW MRI: a pilot study investigating the impact of age and Pre- and Postintervention atrial fibrillation on atrial hemodynamics

PURPOSE

To use four-dimensional (4D) flow MRI to characterize and quantify 3D blood flow in the left atria (LA) of patients with a history of atrial fibrillation (AF).

MATERIALS AND METHODS

The 4D flow MRI was acquired in 19 volunteers (n = 9<30 years, n = 10>50 years) and 10 patients with AF (62 ± 9.6 years; n = 4 in persistent AF, n = 6 postintervention). The LA in each dataset was segmented, and intra-atrial blood flow velocity was quantified. Flow coherence was measured as the consistency of the net blood flow vector.

RESULTS

Quantification of atrial flow revealed significant differences in atrial hemodynamics between age groups. Postintervention AF patients had a mean blood flow of 0.22 ± 0.04 m/s, which was not significantly different than age-matched volunteers (0.21 ± 0.03 m/s). Patients with persistent AF had a mean blood flow of 0.13 ± 0.01 m/s, lower than AF patients in sinus rhythm (0.22 ± 0.04 m/s, P = 0.005), or age-matched volunteers (0.21 ± 0.03 m/s, P < 0.001). Flow coherence was significantly impaired in patients in AF.

CONCLUSION

Flow-sensitive MRI shows that patients with a history of AF had global hemodynamics in the LA similar to those of age-matched volunteers. Additional studies with larger cohorts of AF patients and correlation with outcome are needed to further investigate the potential of atrial 4D flow MRI to flow patterns indicative of stroke risk in AF.

Experimental validation of numerical simulations on a cerebral aneurysm phantom model

The treatment of cerebral aneurysms, found in roughly 5% of the population and associated in case of rupture to a high mortality rate, is a major challenge for neurosurgery and neuroradiology due to the complexity of the intervention and to the resulting, high hazard ratio. Improvements are possible but require a better understanding of the associated, unsteady blood flow patterns in complex 3D geometries. It would be very useful to carry out such studies using suitable numerical models, if it is proven that they reproduce accurately enough the real conditions. This validation step is classically based on comparisons with measured data. Since in vivo measurements are extremely difficult and therefore of limited accuracy, complementary model-based investigations considering realistic configurations are essential. In the present study, simulations based on computational fluid dynamics (CFD) have been compared with in situ, laser-Doppler velocimetry (LDV) measurements in the phantom model of a cerebral aneurysm. The employed 1:1 model is made from transparent silicone. A liquid mixture composed of water, glycerin, xanthan gum and sodium chloride has been specifically adapted for the present investigation. It shows physical flow properties similar to real blood and leads to a refraction index perfectly matched to that of the silicone model, allowing accurate optical measurements of the flow velocity. For both experiments and simulations, complex pulsatile flow waveforms and flow rates were accounted for. This finally allows a direct, quantitative comparison between measurements and simulations. In this manner, the accuracy of the employed computational model can be checked.