Variability of 4D flow parameters when subjected to changes in MRI acquisition parameters using a realistic thoracic aortic phantom

Optimization of 4D Flow MRI Spatial and Temporal Resolution for Examining Complex Hemodynamics in the Carotid Artery Bifurcation

BACKGROUND

Three-dimensional, ECG-gated, time-resolved, three-directional, velocity-encoded phase-contrast MRI (4D flow MRI) has been applied extensively to measure blood velocity in great vessels but has been much less used in diseased carotid arteries. Carotid artery webs (CaW) are non-inflammatory intraluminal shelf-like projections into the internal carotid artery (ICA) bulb that are associated with complex flow and cryptogenic stroke.

PURPOSE

Optimize 4D flow MRI for measuring the velocity field of complex flow in the carotid artery bifurcation model that contains a CaW.

METHODS

A 3D printed phantom model created from computed tomography angiography (CTA) of a subject with CaW was placed in a pulsatile flow loop within the MRI scanner. 4D Flow MRI images of the phantom were acquired with five different spatial resolutions (0.50-2.00  mm3) and four different temporal resolutions (23-96 ms) and compared to a computational fluid dynamics (CFD) solution of the flow field as a reference. We examined four planes perpendicular to the vessel centerline, one in the common carotid artery (CCA) and three in the internal carotid artery (ICA) where complex flow was expected. At these four planes pixel-by-pixel velocity values, flow, and time average wall shear stress (TAWSS) were compared between 4D flow MRI and CFD.

HYPOTHESIS

An optimized 4D flow MRI protocol will provide a good correlation with CFD velocity and TAWSS values in areas of complex flow within a clinically feasible scan time (~ 10 min).

RESULTS

Spatial resolution affected the velocity values, time average flow, and TAWSS measurements. Qualitatively, a spatial resolution of 0.50  mm3 resulted in higher noise, while a lower spatial resolution of 1.50-2.00  mm3 did not adequately resolve the velocity profile. Isotropic spatial resolutions of 0.50-1.00  mm3 showed no significant difference in total flow compared to CFD. Pixel-by-pixel velocity correlation coefficients between 4D flow MRI and CFD were > 0.75 for 0.50-1.00  mm3 but were < 0.5 for 1.50 and 2.00  mm3. Regional TAWSS values determined from 4D flow MRI were generally lower than CFD and decreased at lower spatial resolutions (larger pixel sizes). TAWSS differences between 4D flow and CFD were not statistically significant at spatial resolutions of 0.50-1.00  mm3 but were different at 1.50 and 2.00 mm3. Differences in temporal resolution only affected the flow values when temporal resolution was > 48.4 ms; temporal resolution did not affect TAWSS values.

CONCLUSION

A spatial resolution of 0.74-1.00  mm3 and a temporal resolution of 23-48 ms (1-2 k-space segments) provides a 4D flow MRI protocol capable of imaging velocity and TAWSS in regions of complex flow within the carotid bifurcation at a clinically acceptable scan time.

Repeatability and reproducibility of various 4D Flow MRI postprocessing software programs in a multi-software and multi-vendor cross-over comparison study

BACKGROUND

Different software programs are available for the evaluation of 4D Flow cardiovascular magnetic resonance (CMR). A good agreement of the results between programs is a prerequisite for the acceptance of the method. Therefore, the goal was to compare quantitative results from a cross-over comparison in individuals examined on two scanners of different vendors analyzed with four postprocessing software packages.

METHODS

Eight healthy subjects (27 ± 3 years, 3 women) were each examined on two 3T CMR systems (Ingenia, Philips Healthcare; MAGNETOM Skyra, Siemens Healthineers) with a standardized 4D Flow CMR sequence. Six manually placed aortic contours were evaluated with Caas (Pie Medical Imaging, SW-A), cvi42 (Circle Cardiovascular Imaging, SW-B), GTFlow (GyroTools, SW-C), and MevisFlow (Fraunhofer Institute MEVIS, SW-D) to analyze seven clinically used parameters including stroke volume, peak flow, peak velocity, and area as well as typically scientifically used wall shear stress values. Statistical analysis of inter- and intrareader variability, inter-software and inter-scanner comparison included calculation of absolute and relative error (ER), intraclass correlation coefficient (ICC), Bland-Altman analysis, and equivalence testing based on the assumption that inter-software differences needed to be within 80% of the range of intrareader differences.

RESULTS

SW-A and SW-C were the only software programs showing agreement for stroke volume (ICC = 0.96; ER = 3 ± 8%), peak flow (ICC: 0.97; ER = -1 ± 7%), and area (ICC = 0.81; ER = 2 ± 22%). Results from SW-A/D and SW-C/D were equivalent only for area and peak flow. Other software pairs did not yield equivalent results for routinely used clinical parameters. Especially peak maximum velocity yielded poor agreement (ICC ≤ 0.4) between all software packages except SW-A/D that showed good agreement (ICC = 0.80). Inter- and intrareader consistency for clinically used parameters was best for SW-A and SW-D (ICC = 0.56-97) and worst for SW-B (ICC = -0.01-0.71). Of note, inter-scanner differences per individual tended to be smaller than inter-software differences.

CONCLUSIONS

Of all tested software programs, only SW-A and SW-C can be used equivalently for determination of stroke volume, peak flow, and vessel area. Irrespective of the applied software and scanner, high intra- and interreader variability for all parameters have to be taken into account before introducing 4D Flow CMR in clinical routine. Especially in multicenter clinical trials a single image evaluation software should be applied.

Comparison of Improved Unidirectional Dual Velocity-Encoding MRI Methods

BACKGROUND

In phase-contrast (PC) MRI, several dual velocity encoding methods have been proposed to robustly increase velocity-to-noise ratio (VNR), including a standard dual-VENC (SDV), an optimal dual-VENC (ODV), and bi- and triconditional methods.

PURPOSE

To develop a correction method for the ODV approach and to perform a comparison between methods.

STUDY TYPE

Case-control study.

POPULATION

Twenty-six volunteers.

FIELD STRENGTH/SEQUENCE

1.5 T phase-contrast MRI with VENCs of 50, 75, and 150 cm/second.

ASSESSMENT

Since we acquired single-VENC protocols, we used the background phase from high-VENC (VENC_H ) to reconstruct the low-VENC (VENC_L ) phase. We implemented and compared the unwrapping methods for different noise levels and also developed a correction of the ODV method.

STATISTICAL TESTS

Shapiro-Wilk's normality test, two-way analysis of variance with homogeneity of variances was performed using Levene's test, and the significance level was adjusted by Tukey's multiple post hoc analysis with Bonferroni (P < 0.05).

RESULTS

Statistical analysis revealed no extreme outliers, normally distributed residuals, and homogeneous variances. We found statistically significant interaction between noise levels and the unwrapping methods. This implies that the number of non-unwrapped pixels increased with the noise level. We found that for β = VENC_L /VENC_H  = 1/2, unwrapping methods were more robust to noise. The post hoc test showed a significant difference between the ODV corrected and the other methods, offering the best results regarding the number of unwrapped pixels.

DATA CONCLUSIONS

All methods performed similarly without noise, but the ODV corrected method was more robust to noise at the price of a higher computational time.

LEVEL OF EVIDENCE

4 TECHNICAL EFFICACY STAGE: 1.

Assessment of 4D flow MRI's quality by verifying its Navier-Stokes compatibility

4D Flow Magnetic Resonance Imaging (MRI) is the state-of-the-art technique to comprehensively measure the complex spatio-temporal and multidirectional patterns of blood flow. However, it is subject to artifacts such as noise and aliasing, which due to the 3D and dynamic structure is difficult to detect in clinical practice. In this work, a new mathematical and computational model to determine the quality of 4D Flow MRI is presented. The model is derived by assuming the true velocity satisfies the incompressible Navier-Stokes equations and that can be decomposed by the measurements u→meas plus an extra field w→ . Therefore, a non-linear problem with w→ as unknown arises, which serves as a measure of data quality. A stabilized finite element formulation tailored to this problem is proposed and analyzed. Then, extensive numerical examples-using synthetic 4D Flow MRI data as well as real measurements on experimental phantom and subjects-illustrate the ability to use w→ for assessing the quality of 4D Flow MRI measurements over space and time.

Impact of Respiratory Gating on Hemodynamic Parameters from 4D Flow MRI

The hemodynamic parameters from 4D flow datasets have shown promising diagnostic value in different cardiovascular pathologies. However, the behavior of these parameters can be affected when the 4D flow data are corrupted by respiratory motion. The purpose of this work was to perform a quantitative comparison between hemodynamic parameters computed from 4D flow cardiac MRI both with and without respiratory self-gating. We considered 4D flow MRI data from 15 healthy volunteers (10 men and 5 women, 30.40 ± 6.23 years of age) that were acquired at 3T. Using a semiautomatic segmentation process of the aorta, we obtained the hemodynamic parameters from the 4D flow MRI, with and without respiratory self-gating. A statistical analysis, using the Wilcoxon signed-rank test and Bland–Altman, was performed to compare the hemodynamic parameters from both acquisitions. We found that the calculations of the hemodynamic parameters from 4D flow data that were acquired without respiratory self-gating showed underestimated values in the aortic arch, and the descending and diaphragmatic aorta. We also found a significant variability of the hemodynamic parameters in the ascending aorta of healthy volunteers when comparing both methods. The 4D flow MRI requires respiratory compensation to provide reliable calculations of hemodynamic parameters.

Valvular regurgitation flow jet assessment using in vitro 4D flow MRI: Implication for mitral regurgitation

PURPOSE

The purpose of this study was to evaluate the accuracy of four-dimensional (4D) flow MRI for direct assessment of peak velocity, flow volume, and momentum of a mitral regurgitation (MR) flow jets using an in vitro pulsatile jet flow phantom. We systematically investigated the impact of spatial resolution and quantification location along the jet on flow quantities with Doppler ultrasound as a reference for peak velocity.

METHODS

Four-dimensional flow MRI data of a pulsatile jet through a circular, elliptical, and 3D-printed patient-specific MR orifice model was acquired with varying spatial resolution (1.5-5 mm isotropic voxel). Flow rate and momentum of the jet were quantified at various axial distances (x = 0-50 mm) and integrated over time to calculate Vol_jet and MTI_jet . In vivo assessment of Vol_jet and MTI_jet was performed on 3 MR patients.

RESULTS

Peak velocities were comparable to Doppler ultrasound (3% error, 1.5 mm voxel), but underestimated with decreasing spatial resolution (-40% error, 5 mm voxel). Vol_jet was similar to regurgitant volume (RVol) within 5 mm, and then increased linearly with the axial distance (19%/cm) because of flow entrainment. MTI_jet remained steady throughout the jet (2%/cm) as theoretically predicted. Four and 9 voxels across the jet were required to measure flow volume and momentum-time-integral within 10% error, respectively.

CONCLUSION

Four-dimensional flow MRI detected accurate peak velocity, flow rate, and momentum for in vitro MR-mimicking flow jets. Spatial resolution significantly impacted flow quantitation, which otherwise followed predictions of flow entrainment and momentum conservation. This study provides important preliminary information for accurate in vivo MR assessment using 4D flow MRI.

A novel MRI-based data fusion methodology for efficient, personalised, compliant simulations of aortic haemodynamics

We present a novel, cost-efficient methodology to simulate aortic haemodynamics in a patient-specific, compliant aorta using an MRI data fusion process. Based on a previously-developed Moving Boundary Method, this technique circumvents the high computational cost and numerous structural modelling assumptions required by traditional Fluid-Structure Interaction techniques. Without the need for Computed Tomography (CT) data, the MRI images required to construct the simulation can be obtained during a single imaging session. Black Blood MR Angiography and 2D Cine-MRI data were used to reconstruct the luminal geometry and calibrate wall movement specifically to each region of the aorta. 4D-Flow MRI and non-invasive pressure measurements informed patient-specific inlet and outlet boundary conditions. Luminal area closely matched 2D Cine-MRI measurements with a mean error of less than 4.6% across the cardiac cycle, while physiological pressure and flow distributions were simulated to within 3.3% of patient-specific targets. Moderate agreement with 4D-Flow MRI velocity data was observed. Despite lower peak velocity, an equivalent rigid-wall simulation predicted a mean Time-Averaged Wall Shear Stress (TAWSS) 13% higher than the compliant simulation. The agreement observed between compliant simulation results and MRI data is testament to the accuracy and efficiency of this MRI-based simulation technique.

Validation of 4D Flow based relative pressure maps in aortic flows

While the clinical gold standard for pressure difference measurements is invasive catheterization, 4D Flow MRI is a promising tool for enabling a non-invasive quantification, by linking highly spatially resolved velocity measurements with pressure differences via the incompressible Navier-Stokes equations. In this work we provide a validation and comparison with phantom and clinical patient data of pressure difference maps estimators. We compare the classical Pressure Poisson Estimator (PPE) and the new Stokes Estimator (STE) against catheter pressure measurements under a variety of stenosis severities and flow intensities. Specifically, we use several 4D Flow data sets of realistic aortic phantoms with different anatomic and hemodynamic severities and two patients with aortic coarctation. The phantom data sets are enriched by subsampling to lower resolutions, modification of the segmentation and addition of synthetic noise, in order to study the sensitivity of the pressure difference estimators to these factors. Overall, the STE method yields more accurate results than the PPE method compared to catheterization data. The superiority of the STE becomes more evident at increasing Reynolds numbers with a better capacity of capturing pressure gradients in strongly convective flow regimes. The results indicate an improved robustness of the STE method with respect to variation in lumen segmentation. However, with heuristic removal of the wall-voxels, the PPE can reach a comparable accuracy for lower Reynolds' numbers.

Traveling Volunteers: A Multi-Vendor, Multi-Center Study on Reproducibility and Comparability of 4D Flow Derived Aortic Hemodynamics in Cardiovascular Magnetic Resonance

BACKGROUND

Implementation of four-dimensional flow magnetic resonance (4D Flow MR) in clinical routine requires awareness of confounders.

PURPOSE

To investigate inter-vendor comparability of 4D Flow MR derived aortic hemodynamic parameters, assess scan-rescan repeatability, and intra- and interobserver reproducibility.

STUDY TYPE

Prospective multicenter study.

POPULATION

Fifteen healthy volunteers (age 24.5 ± 5.3 years, 8 females).

FIELD STRENGTH/SEQUENCE

3 T, vendor-provided and clinically used 4D Flow MR sequences of each site.

ASSESSMENT

Forward flow volume, peak velocity, average, and maximum wall shear stress (WSS) were assessed via nine planes (P1-P9) throughout the thoracic aorta by a single observer (AD, 2 years of experience). Inter-vendor comparability as well as scan-rescan, intra- and interobserver reproducibility were examined.

STATISTICAL TESTS

Equivalence was tested setting the 95% confidence interval of intraobserver and scan-rescan difference as the limit of clinical acceptable disagreement. Intraclass correlation coefficient (ICC) and Bland-Altman plots were used for scan-rescan reproducibility and intra- and interobserver agreement. A P-value <0.05 was considered statistically significant. ICCs ≥ 0.75 indicated strong correlation (>0.9: excellent, 0.75-0.9: good).

RESULTS

Ten volunteers finished the complete study successfully. 4D flow derived hemodynamic parameters between scanners of three different vendors are not equivalent exceeding the equivalence range. P3-P9 differed significantly between all three scanners for forward flow (59.1 ± 13.1 mL vs. 68.1 ± 12.0 mL vs. 55.4 ± 13.1 mL), maximum WSS (1842.0 ± 190.5 mPa vs. 1969.5 ± 398.7 mPa vs. 1500.6 ± 247.2 mPa), average WSS (1400.0 ± 149.3 mPa vs. 1322.6 ± 211.8 mPa vs. 1142.0 ± 198.5 mPa), and peak velocity between scanners I vs. III (114.7 ± 12.6 cm/s vs. 101.3 ± 15.6 cm/s). Overall, the plane location at the sinotubular junction (P1) presented most inter-vendor stability (forward: 78.5 ± 15.1 mL vs. 80.3 ± 15.4 mL vs. 79.5 ± 19.9 mL [P = 0.368]; peak: 126.4 ± 16.7 cm/s vs. 119.7 ± 13.6 cm/s vs. 111.2 ± 22.6 cm/s [P = 0.097]). Scan-rescan reproducibility and intra- and interobserver variability were good to excellent (ICC ≥ 0.8) with best agreement for forward flow (ICC ≥ 0.98).

DATA CONCLUSION

The clinical protocol used at three different sites led to differences in hemodynamic parameters assessed by 4D flow.

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY STAGE: 2.

Parametric Sequential Method for MRI-Based Wall Shear Stress Quantification

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

Bicuspid aortic valve morphology and aortic valvular outflow jets: an experimental analysis using an MRI-compatible pulsatile flow circulation system

The characteristics of aortic valvular outflow jet affect aortopathy in the bicuspid aortic valve (BAV). This study aimed to elucidate the effects of BAV morphology on the aortic valvular outflow jets. Morphotype-specific valve-devising apparatuses were developed to create aortic valve models. A magnetic resonance imaging-compatible pulsatile flow circulation system was developed to quantify the outflow jet. The eccentricity and circulation values of the peak systolic jet were compared among tricuspid aortic valve (TAV), three asymmetric BAVs, and two symmetric BAVs. The results showed mean aortic flow and leakage did not differ among the five BAVs (six samples, each). Asymmetric BAVs demonstrated the eccentric outflow jets directed to the aortic wall facing the smaller leaflets. In the asymmetric BAV with the smaller leaflet facing the right-anterior, left-posterior, and left-anterior quadrants of the aorta, the outflow jets exclusively impinged on the outer curvature of the ascending aorta, proximal arch, and the supra-valvular aortic wall, respectively. Symmetric BAVs demonstrated mildly eccentric outflow jets that did not impinge on the aortic wall. The circulation values at peak systole increased in asymmetric BAVs. The bicuspid symmetry and the position of smaller leaflet were determinant factors of the characteristics of aortic valvular outflow jet.

Impact of sequence type and field strength (1.5, 3, and 7T) on 4D flow MRI hemodynamic aortic parameters in healthy volunteers

PURPOSE

4D flow magnetic resonance imaging (4D-MRI) allows time-resolved visualization of blood flow patterns, quantification of volumes, velocities, and advanced parameters, such as wall shear stress (WSS). As 4D-MRI enters the clinical arena, standardization and awareness of confounders are important. Our aim was to evaluate the equivalence of 4D flow-derived aortic hemodynamics in healthy volunteers using different sequences and field strengths.

METHODS

4D-MRI was acquired in 10 healthy volunteers at 1.5T using three different prototype sequences, at 3T and at 7T (Siemens Healthineers). After evaluation of diagnostic quality in three segments (ascending-, descending aorta, aortic arch), peak velocity, flow volumes, and WSS were investigated. Equivalence limits for comparison of field strengths/sequences were based on the limits of Bland-Altman analyses of the intraobserver variability.

RESULTS

Non-diagnostic quality was found in 10/144 segments, 9/10 were obtained at 7T. Apart for the comparison of forward flow between sequence 1 and 3, the differences in measurements between field strengths/sequences exceeded the range of agreement. Significant differences were found between field strengths/sequences for forward flow (1.5T vs. 3T, 3T vs. 7T, sequence 1 vs. 3, 2 vs. 3 [P < .001]), WSS (1.5T vs. 3T [P < .05], sequence 1 vs. 2, 1 vs. 3, 2 vs. 3 [P < .001]), and peak velocity (1.5T vs. 7T, sequence 1 vs. 3 [P > .001]). All parameters at all field strengths/with all sequences correlated moderately to strongly (r ≥ 0.5).

CONCLUSION

Data from all sequences could be acquired and resulting images showed sufficient quality for further analysis. However, the variability of the measurements of peak velocity, flow volumes, and WSS was higher when comparing field strengths/sequences as the equivalence limits defined by the intraobserver assessments.

4D flow imaging with UNFOLD in a reduced FOV

PURPOSE

Two-dimensional selective excitation (2DRF) allows shortening 4D flow scan times by reducing the FOV, but the longer 2DRF pulse duration decreases the temporal resolution, yielding underestimated peak flow values. Multiple k-space lines per cardiac phase, n_l ≥ 2, are commonly applied in 4D flow MRI to shorten the inherent long scan times. We demonstrate that 2DRF 4D flow with n_l ≥ 2 can be easily combined with UNFOLD (UNaliasing by Fourier-encoding the Overlaps using the temporaL Dimension), a technique that allows regaining nominally the temporal resolution of the respective acquisition with n_l = 1, to assure peak flow quantification.

METHODS

Two different 2DRF pulses with spiral k-space trajectories were designed and integrated into a 4D flow sequence. Flow phantom experiments and 7 healthy control 4D flow in vivo measurements, with and without UNFOLD reconstructions, were compared with conventional reconstruction and 1D slab-selective excitation (1DRF) by evaluating time-resolved flow curves, peak flow, peak velocity, blood flow volume per cardiac cycle, and spatial aliasing.

RESULTS

Applying UNFOLD to 4D flow imaging with 2DRF and reduced FOV increased the quantified in vivo peak flow values significantly by 3.7% ± 2.3% to 5.2% ± 2.4% (P < .05). Accordingly, the peak flow underestimation of 2DRF scans compared with conventional 1DRF scans decreased with UNFOLD. Finally, 2DRF combined with UNFOLD accelerated the 4D flow acquisition 3.5 ± 1.4 fold by reducing the FOV and increasing the effective temporal resolution by 6.7% compared with conventional 1D selective excitation, with 2 k-space lines per cardiac phase.

CONCLUSION

Two-dimensional selective excitation combined with UNFOLD allows limiting the FOV to shorten 4D flow scan times and compensates for the loss in temporal resolution with 2DRF (Δt = 64.8 ms) compared with 1DRF (Δt = 43.2 ms), yielding an effective resolution of Δt_eff = 40.5 ms to enhance peak flow quantification.

MR-based wall shear stress measurements in fully developed turbulent flow using the Clauser plot method

In arterial blood flow wall shear stress (WSS) quantifies the frictional force that flowing blood exerts on a vessel wall. WSS can be directly estimated from phase-contrast (PC) MR velocity measurements and has been suggested as a biomarker in cardio-vascular diseases. We present and investigate the application of the Clauser plot method for estimating WSS in fully developed turbulent stationary flow using PC velocity measurements. The Clauser plot method estimates WSS from the logarithmic region of boundary layer in fully developed turbulent stationary flow. The Clauser plot method was evaluated using 2D PC-MR phantom measurements at 3 T for different in-plane resolutions at various Reynolds numbers. WSS values derived from the Clauser plot were compared to results from Laser Doppler Velocimetry (LDV) measurements and theoretical results calculated using the friction factor formula for smooth pipe flow. For all Reynolds numbers, WSS values derived from the Clauser plot were in good agreement with results from LDV measurements and values using the friction factor formula (relative deviations ∼5%). Furthermore, Clauser plot derived results were almost independent of spatial resolution, in contrast to WSS results obtained with our in-house software tool for MR-based WSS quantification showing relative deviations of more than 100%. In fully developed turbulent flow, the Clauser plot method provides highly consistent WSS independent of the underlying spatial resolution. Therefore, it renders a valuable approach for MR-based WSS estimates in controllable flow settings. Although its direct in vivo applicability is severely limited because of the different flow character, it may serve as helpful approach for validation of MR-based WSS quantification algorithms prior to their clinical application.

Two-Minute k-Space and Time-accelerated Aortic Four-dimensional Flow MRI: Dual-Center Study of Feasibility and Impact on Velocity and Wall Shear Stress Quantification

Purpose

To investigate the two-center feasibility of highly k-space and time (k-t)-accelerated 2-minute aortic four-dimensional (4D) flow MRI and to evaluate its performance for the quantification of velocities and wall shear stress (WSS).

Materials and Methods

This cross-sectional study prospectively included 68 participants (center 1, 11 healthy volunteers [mean age ± standard deviation, 61 years ± 15] and 16 patients with aortic disease [mean age, 60 years ± 10]; center 2, 14 healthy volunteers [mean age, 38 years ± 13] and 27 patients with aortic or cardiac disease [mean age, 78 years ± 18]). Each participant underwent highly accelerated 4D flow MRI (k-t acceleration, acceleration factor of 5) of the thoracic aorta. For comparison, conventional 4D flow MRI (acceleration factor of 2) was acquired in the participants at center 1 (n = 27). Regional aortic peak systolic velocities and three-dimensional WSS were quantified.

Results

k-t-accelerated scan times (center 1, 2:03 minutes ± 0:29; center 2, 2:06 minutes ± 0:20) were significantly reduced compared with conventional 4D flow MRI (center 1, 12:38 minutes ± 2:25; P < .0001). Overall good agreement was found between the two techniques (absolute differences ≤15%), but proximal aortic WSS was significantly underestimated in patients by using k-t-accelerated 4D flow when compared with conventional 4D flow (P ≤ .03). k-t-accelerated 4D flow MRI was reproducible (intra- and interobserver intraclass correlation coefficient ≥0.98) and identified significantly increased peak velocities and WSS in patients with stenotic (P ≤ .003) or bicuspid (P ≤ .04) aortic valves compared with healthy volunteers. In addition, k-t-accelerated 4D flow MRI-derived velocities and WSS were inversely related to age (r ≥-0.53; P ≤ .03) over all healthy volunteers.

Conclusion

k-t-accelerated aortic 4D flow MRI providing 2-minute scan times was feasible and reproducible at two centers. Although consistent healthy aging- and disease-related changes in aortic hemodynamics were observed, care should be taken when considering WSS, which can be underestimated in patients.© RSNA, 2019See also the commentary by François in this issue.

Optimal Dual-VENC Unwrapping in Phase-Contrast MRI

Dual-VENC strategies have been proposed to improve the velocity-to-noise ratio in phase-contrast MRI. However, they are based on aliasing-free high-VENC data. The aim of this paper is to propose a dual-VENC velocity estimation method allowing high-VENC aliased data. For this purpose, we reformulate the phase-contrast velocity as a least squares estimator, providing a natural framework for including multiple encoding gradient measurements. By analyzing the mathematical properties of both single- and dual-VENC problems, we can justify theoretically high/low-VENC ratios such that the aliasing velocity can be minimized. The resulting reconstruction algorithm was assessed using three types of data: numerical, experimental, and volunteers. In clinical practice, this method would allow shorter examination times by avoiding tedious adaptation of VENC values by repeated scans.

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.