3D Quantification of Wall Shear Stress and Oscillatory Shear Index Using a Finite-Element Method in 3D CINE PC-MRI Data of the Thoracic Aorta

Simulating Subject-Specific Aortic Hemodynamic Effects of Valvular Lesions in Rheumatic Heart Disease

Rheumatic heart disease (RHD) is a neglected tropical disease despite the substantial global health burden. In this study, we aimed to develop a lower cost method of modeling aortic blood flow using subject-specific velocity profiles, aiding our understanding of RHD's consequences on the structure and function of the ascending aorta. Echocardiography and cardiovascular magnetic resonance (CMR) are often used for diagnosis, including valve dysfunction assessments. However, there is a need to further characterize aortic valve lesions to improve treatment options and timing for patients, while using accessible and affordable imaging strategies. Here, we simulated effects of RHD aortic valve lesions on the aorta using computational fluid dynamics (CFD). We hypothesized that inlet velocity distribution and wall shear stress (WSS) will differ between RHD and non-RHD individuals, as well as between subject-specific and standard Womersley velocity profiles. Phase-contrast CMR data from South Africa of six RHD subjects with aortic stenosis and/or regurgitation and six matched controls were used to estimate subject-specific velocity inlet profiles and the mean velocity for Womersley profiles. Our findings were twofold. First, we found WSS in subject-specific RHD was significantly higher (p < 0.05) than control subject simulations, while Womersley simulation groups did not differ. Second, evaluating spatial velocity differences (ΔSV) between simulation types revealed that simulations of RHD had significantly higher ΔSV than non-RHD (p < 0.05), these results highlight the need for implementing subject-specific input into RHD CFD, which we demonstrate how to accomplish through accessible methods.

Influence of Non-Newtonian Viscosity on Flow Structures and Wall Deformation in Compliant Serpentine Microchannels: A Numerical Study

The viscosity of fluid plays a major role in the flow dynamics of microchannels. Viscous drag and shear forces are the primary tractions for microfluidic fluid flow. Capillary blood vessels with a few microns diameter are impacted by the rheology of blood flowing through their conduits. Hence, regenerated capillaries should be able to withstand such impacts. Consequently, there is a need to understand the flow physics of culture media through the lumen of the substrate as it is one of the vital promoting factors for vasculogenesis under optimal shear conditions at the endothelial lining of the regenerated vessel. Simultaneously, considering the diffusive role of capillaries for ion exchange with the surrounding tissue, capillaries have been found to reorient themselves in serpentine form for modulating the flow conditions while developing sustainable shear stress. In the current study, S-shaped (S1) and delta-shaped (S2) serpentine models of capillaries were considered to evaluate the shear stress distribution and the oscillatory shear index (OSI) and relative residual time (RRT) of the derivatives throughout the channel (due to the phenomena of near-wall stress fluctuation), along with the influence of culture media rheology on wall stress parameters. The non-Newtonian power-law formulation was implemented for defining rheological viscosity of the culture media. The flow actuation of the media was considered to be sinusoidal and physiological, realizing the pulsatile blood flow behavior in the circulatory network. A distinct difference in shear stress distributions was observed in both the serpentine models. The S1 model showed higher change in shear stress in comparison to the S2 model. Furthermore, the non-Newtonian viscosity formulation was found to produce more sustainable shear stress near the serpentine walls compared to the Newtonian formulation fluid, emphasizing the influence of rheology on stress generation. Further, cell viability improved in the bending regions of serpentine channels compared to the long run section of the same channel.

Physics-informed compressed sensing for PC-MRI: an inverse Navier-Stokes problem

We formulate a physics-informed compressed sensing (PICS) method for the reconstruction of velocity fields from noisy and sparse phase-contrast magnetic resonance signals. The method solves an inverse Navier-Stokes boundary value problem, which permits us to jointly reconstruct and segment the velocity field, and at the same time infer hidden quantities such as the hydrodynamic pressure and the wall shear stress. Using a Bayesian framework, we regularize the problem by introducing a priori information about the unknown parameters in the form of Gaussian random fields. This prior information is updated using the Navier-Stokes problem, an energy-based segmentation functional, and by requiring that the reconstruction is consistent with the k-space signals. We create an algorithm that solves this inverse problem, and test it for noisy and sparse k-space signals of the flow through a converging nozzle. We find that the method is capable of reconstructing and segmenting the velocity fields from sparsely-sampled (15% k-space coverage), low (~10) signal-to-noise ratio (SNR) signals, and that the reconstructed velocity field compares well with that derived from fully-sampled (100% k-space coverage) high (>40) SNR signals of the same flow.

The role of 4-dimensional flow in the assessment of bicuspid aortic valve and its valvulo-aortopathies

Bicuspid aortic valve is the most common congenital cardiac malformation and the leading cause of aortopathy and aortic stenosis in younger patients. Aortic wall remodelling secondary to altered haemodynamic flow patterns, changes in peak velocity, and wall shear stress may be implicated in the development of aortopathy in the presence of bicuspid aortic valve and dysfunction. Assessment of these parameters as potential predictors of disease severity and progression is thus desirable. The anatomic and functional information acquired from 4D flow MRI can allow simultaneous visualisation and quantification of the pathological geometric and haemodynamic changes of the aorta. We review the current clinical utility of haemodynamic quantities including velocity, wall sheer stress and energy losses, as well as visual descriptors such as vorticity and helicity, and flow direction in assessing the aortic valve and associated aortopathies.

The impact of 4D-Flow MRI spatial resolution on patient-specific CFD simulations of the thoracic aorta

Magnetic Resonance Imaging (MRI) is considered the gold standard of medical imaging technologies as it allows for accurate imaging of blood vessels. 4-Dimensional Flow Magnetic Resonance Imaging (4D-Flow MRI) is built on conventional MRI, and provides flow data in the three vector directions and a time resolved magnitude data set. As such it can be used to retrospectively calculate haemodynamic parameters of interest, such as Wall Shear Stress (WSS). However, multiple studies have indicated that a significant limitation of the imaging technique is the spatiotemporal resolution that is currently available. Recent advances have proposed and successfully integrated 4D-Flow MRI imaging techniques with Computational Fluid Dynamics (CFD) to produce patient-specific simulations that have the potential to aid in treatments,surgical decision making, and risk stratification. However, the consequences of using insufficient 4D-Flow MRI spatial resolutions on any patient-specific CFD simulations is currently unclear, despite being a recognised limitation. The research presented in this study aims to quantify the inaccuracies in patient-specific 4D-Flow MRI based CFD simulations that can be attributed to insufficient spatial resolutions when acquiring 4D-Flow MRI data. For this research, a patient has undergone four 4D-Flow MRI scans acquired at various isotropic spatial resolutions and patient-specific CFD simulations have subsequently been run using geometry and velocity data produced from each scan. It was found that compared to CFD simulations based on a \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1.5\,{\text {mm}} \times 1.5\,{\text {mm}} \times 1.5\,{\text {mm}}$$\end{document}1.5mm×1.5mm×1.5mm, using a spatial resolution of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$4\,{\text {mm}} \times 4\,{\text {mm}} \times 4\,{\text {mm}}$$\end{document}4mm×4mm×4mm substantially underestimated the maximum velocity magnitude at peak systole by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$110.55\%$$\end{document}110.55%. The impacts of 4D-Flow MRI spatial resolution on WSS calculated from CFD simulations have been investigated and it has been shown that WSS is underestimated in CFD simulations that are based on a coarse 4D-Flow MRI spatial resolution. The authors have concluded that a minimum 4D-Flow MRI spatial resolution of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1.5\,{\text {mm}} \times 1.5\,{\text {mm}} \times 1.5\,{\text {mm}}$$\end{document}1.5mm×1.5mm×1.5mm must be used when acquiring 4D-Flow MRI data to perform patient-specific CFD simulations. A coarser spatial resolution will produce substantial differences within the flow field and geometry.

Fully Three-Dimensional Hemodynamic Characterization of Altered Blood Flow in Bicuspid Aortic Valve Patients With Respect to Aortic Dilatation: A Finite Element Approach

Background and Purpose

Prognostic models based on cardiovascular hemodynamic parameters may bring new information for an early assessment of patients with bicuspid aortic valve (BAV), playing a key role in reducing the long-term risk of cardiovascular events. This work quantifies several three-dimensional hemodynamic parameters in different patients with BAV and ranks their relationships with aortic diameter.

Materials and Methods

Using 4D-flow CMR data of 74 patients with BAV (49 right-left and 25 right-non-coronary) and 48 healthy volunteers, aortic 3D maps of seventeen 17 different hemodynamic parameters were quantified along the thoracic aorta. Patients with BAV were divided into two morphotype categories, BAV-Non-AAoD (where we include 18 non-dilated patients and 7 root-dilated patients) and BAV-AAoD (where we include the 49 patients with dilatation of the ascending aorta). Differences between volunteers and patients were evaluated using MANOVA with Pillai's trace statistic, Mann–Whitney U test, ROC curves, and minimum redundancy maximum relevance algorithm. Spearman's correlation was used to correlate the dilation with each hemodynamic parameter.

Results

The flow eccentricity, backward velocity, velocity angle, regurgitation fraction, circumferential wall shear stress, axial vorticity, and axial circulation allowed to discriminate between volunteers and patients with BAV, even in the absence of dilation. In patients with BAV, the diameter presented a strong correlation (&gt; |+/−0.7|) with the forward velocity and velocity angle, and a good correlation (&gt; |+/−0.5|) with regurgitation fraction, wall shear stress, wall shear stress axial, and vorticity, also for morphotypes and phenotypes, some of them are correlated with the diameter. The velocity angle proved to be an excellent biomarker in the differentiation between volunteers and patients with BAV, BAV morphotypes, and BAV phenotypes, with an area under the curve bigger than 0.90, and higher predictor important scores.

Conclusions

Through the application of a novel 3D quantification method, hemodynamic parameters related to flow direction, such as flow eccentricity, velocity angle, and regurgitation fraction, presented the best relationships with a local diameter and effectively differentiated patients with BAV from healthy volunteers.

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.

Numerical simulation in the abdominal aorta and the visceral arteries with or without stenosis based on 2D PCMRI

It is important to obtain accurate boundary conditions (BCs) in hemodynamic simulations. This article aimed to improve the accuracy of BCs in computational fluid dynamics (CFD) simulation and analyze the differences in hemodynamics between healthy volunteers and patients with visceral arterial stenosis (VAS). The geometric models of seven cases were reconstructed using the magnetic resonance angiogram (MRA) or computed tomography angiogram (CTA) imaging data. The physiological flow waveforms obtained from 2D Phase Contrast Magnetic Resonance Imaging (PCMRI) were imposed on the aortic inlet and the visceral arteries' outlets. The individualized RCR values of the three-element Windkessel model were imposed on the aortic outlet. CFD simulations were run in the open-source software: svSolver. Two specific time points were selected to compare the hemodynamics of healthy volunteers and patients with VAS. The results suggested that blood in the stenotic visceral arteries flowed at high speed throughout the cardiac cycle. The low pressure is distributed at stenotic lesions. The wall shear stress (WSS) reached 4 Pa near stenotic locations. The low time average wall shear stress (TAWSS), high oscillatory shear index (OSI), and high relative residence time (RRT) concentrated in the abdominal aorta. Besides, the ratios of the areas with low TAWSS, high OSI, and high RRT to the computational domain were higher in patients with VAS than which in the healthy volunteers. The individualized BCs were used for hemodynamic simulations and results suggest that patients with stenosis have a higher risk of blood retention and atherosclerosis formation in the abdominal aorta.

Tissue engineered in-vitro vascular patch fabrication using hybrid 3D printing and electrospinning

Three-dimensional (3D) engineered cardiovascular tissues have shown great promise to replace damaged structures. Specifically, tissue engineering vascular grafts (TEVG) have the potential to replace biological and synthetic grafts. We aimed to design an in-vitro patient-specific patch based on a hybrid 3D print combined with vascular smooth muscle cells (VSMC) differentiation. Based on the medical images of a 2 months-old girl with aortic arch hypoplasia and using computational modelling, we evaluated the most hemodynamically efficient aortic patch surgical repair. Using the designed 3D patch geometry, the scaffold was printed using a hybrid fused deposition modelling (FDM) and electrospinning techniques. The scaffold was seeded with multipotent mesenchymal stem cells (MSC) for later maturation to derived VSMC (dVSMC). The graft showed adequate resistance to physiological aortic pressure (burst pressure 101 ​± ​15 ​mmHg) and a porosity gradient ranging from 80 to 10 ​μm allowing cells to infiltrate through the entire thickness of the patch. The bio-scaffolds showed good cell viability at days 4 and 12 and adequate functional vasoactive response to endothelin-1. In summary, we have shown that our method of generating patient-specific patch shows adequate hemodynamic profile, mechanical properties, dVSMC infiltration, viability and functionality. This innovative 3D biotechnology has the potential for broad application in regenerative medicine and potentially in heart disease prevention.

Influence of Spatial Resolution and Compressed SENSE Acceleration Factor on Flow Quantification with 4D Flow MRI at 3 Tesla

Four-dimensional (4D) flow MRI allows quantifying flow in blood vessels–non invasively and in vivo. The clinical use of 4D flow MRI in small vessels, however, is hampered by long examination times and limited spatial resolution. Compressed SENSE (CS-SENSE) is a technique that can accelerate 4D flow dramatically. Here, we investigated the effect of spatial resolution and CS acceleration on flow measurements by using 4D flow MRI in small vessels in vitro at 3 T. We compared the flow in silicon tubes (inner diameters of 2, 3, 4, and 5 mm) measured with 4D flow MRI, accelerated with four CS factors (CS = 2.5, 4.5, 6.5, and 13) and three voxel sizes (0.5, 1, and 1.5 mm³) to 2D flow MRI and a flow sensor. Additionally, the velocity field in an aneurysm model acquired with 4D flow MRI was compared to the one simulated with computational fluid dynamics (CFD). A strong correlation was observed between flow sensor, 2D flow MRI, and 4D flow MRI (rho > 0.94). The use of fewer than seven voxels per vessel diameter (nROI) resulted in an overestimation of flow in more than 5% of flow measured with 2D flow MRI. A negative correlation (rho = −0.81) between flow error and nROI were found for CS = 2.5 and 4.5. No statistically significant impact of CS factor on differences in flow rates was observed. However, a trend of increased flow error with increased CS factor was observed. In an aneurysm model, the peak velocity and stagnation zone were detected by CFD and all 4D flow MRI variants. The velocity difference error in the aneurysm sac did not exceed 11% for CS = 4.5 in comparison to CS = 2.5 for all spatial resolutions. Therefore, CS factors from 2.5–4.5 can appear suitable to improve spatial or temporal resolution for accurate quantification of flow rate and velocity. We encourage reporting the number of voxels per vessel diameter to standardize 4D flow MRI protocols.

Identification of hemodynamic biomarkers for bicuspid aortic valve induced aortic dilation using machine learning

Recent advances in medical imaging have confirmed the presence of altered hemodynamics in bicuspid aortic valve (BAV) patients. Therefore, there is a need for new hemodynamic biomarkers to refine disease monitoring and improve patient risk stratification. This research aims to analyze and extract multiple correlation patterns of hemodynamic parameters from 4D Flow MRI data and find which parameters allow an accurate classification between healthy volunteers (HV) and BAV patients with dilated and non-dilated ascending aorta using machine learning. Sixteen hemodynamic parameters were calculated in the ascending aorta (AAo) and aortic arch (AArch) at peak systole from 4D Flow MRI. We used sequential forward selection (SFS) and principal component analysis (PCA) as feature selection algorithms. Then, eleven machine-learning classifiers were implemented to separate HV and BAV patients (non- and dilated ascending aorta). Multiple correlation patterns from hemodynamic parameters were extracted using hierarchical clustering. The linear discriminant analysis and random forest are the best performing classifiers, using five hemodynamic parameters selected with SFS (velocity angle, forward velocity, vorticity, and backward velocity in AAo; and helicity density in AArch) a 96.31 ± 1.76% and 96.00 ± 0.83% accuracy, respectively. Hierarchical clustering revealed three groups of correlated features. According to this analysis, we observed that features selected by SFS have a better performance than those selected by PCA because the five selected parameters were distributed according to 3 different clusters. Based on the proposed method, we concluded that the feature selection method found five potentially hemodynamic biomarkers related to this disease.

2D Projection Maps of WSS and OSI Reveal Distinct Spatiotemporal Changes in Hemodynamics in the Murine Aorta during Ageing and Atherosclerosis

Growth, ageing and atherosclerotic plaque development alter the biomechanical forces acting on the vessel wall. However, monitoring the detailed local changes in wall shear stress (WSS) at distinct sites of the murine aortic arch over time has been challenging. Here, we studied the temporal and spatial changes in flow, WSS, oscillatory shear index (OSI) and elastic properties of healthy wildtype (WT, n = 5) and atherosclerotic apolipoprotein E-deficient (Apoe^-/-, n = 6) mice during ageing and atherosclerosis using high-resolution 4D flow magnetic resonance imaging (MRI). Spatially resolved 2D projection maps of WSS and OSI of the complete aortic arch were generated, allowing the pixel-wise statistical analysis of inter- and intragroup hemodynamic changes over time and local correlations between WSS, pulse wave velocity (PWV), plaque and vessel wall characteristics. The study revealed converse differences of local hemodynamic profiles in healthy WT and atherosclerotic Apoe^-/- mice, and we identified the circumferential WSS as potential marker of plaque size and composition in advanced atherosclerosis and the radial strain as a potential marker for vascular elasticity. Two-dimensional (2D) projection maps of WSS and OSI, including statistical analysis provide a powerful tool to monitor local aortic hemodynamics during ageing and atherosclerosis. The correlation of spatially resolved hemodynamics and plaque characteristics could significantly improve our understanding of the impact of hemodynamics on atherosclerosis, which may be key to understand plaque progression towards vulnerability.

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.

Hemodynamic Comparison of Stent-Grafts for the Treatment of Aortoiliac Occlusive Disease

Purpose:

To compare the flow patterns and hemodynamics of the AFX stent-graft and the covered endovascular reconstruction of aortic bifurcation (CERAB) configuration using laser particle image velocimetry (PIV) experiments.

Materials and Methods:

Two anatomically realistic aortoiliac phantoms were constructed using polydimethylsiloxane polymer. An AFX stent-graft with a transparent cover made with a new method was inserted into one phantom. A CERAB configuration using Atrium’s Avanta V12 with transparent covers made with a previously established method was inserted into the other phantom, both modified stent-grafts were suitable for laser PIV, enabling visualization of the flow fields and quantification of time average wall shear stress (TAWSS), oscillatory shear index (OSI), and relative residence time (RRT).

Results:

Disturbed flow was observed at the bifurcation region of the AFX, especially at the end systolic velocity (ESV) time-point where recirculation was noticeable due to vortical flow. In contrast, predominantly unidirectional flow was observed at the CERAB bifurcation. These observations were confirmed by the quantified hemodynamic results from PIV analysis where mean TAWSS of 0.078 Pa (range: 0.009–0.242 Pa) was significantly lower in AFX as compared with 0.229 Pa (range: 0.013–0.906 Pa) for CERAB (p<0.001). Mean OSI of 0.318 (range: 0.123–0.496) in AFX was significantly higher than 0.252 (range: 0.055–0.472) in CERAB (p<0.001). Likewise, mean RRT of 180 Pa−1 (range: 9–3603 Pa−1) in AFX was also significantly higher than 88 Pa−1 (range: 2–840 Pa−1) in CERAB (p=0.0086).

Conclusion:

In this in vitro study, the flow pattern of a modified AFX stent-graft was found to be more disturbed especially at the end systolic phase, its hemodynamic outcomes less desirable than CERAB configuration.

Clinical Relevance:

While the AFX stent-graft has an advantage over the CERAB configuration in eliminating radial mismatch, and maintaining the anatomical bifurcation for future endovascular intervention, this in vitro study revealed that the associated lower TAWSS, higher OSI and RRT may predispose to thrombosis and are, thus, less desirable as compared to a CERAB configuration. Further investigation is warranted to confirm whether these findings translate into the clinical setting.

Hemodynamic Comparison of Stent-Grafts for the Treatment of Aortoiliac Occlusive Disease

PURPOSE

To compare the flow patterns and hemodynamics of the AFX stent-graft and the covered endovascular reconstruction of aortic bifurcation (CERAB) configuration using laser particle image velocimetry (PIV) experiments.

MATERIALS AND METHODS

Two anatomically realistic aortoiliac phantoms were constructed using polydimethylsiloxane polymer. An AFX stent-graft with a transparent cover made with a new method was inserted into one phantom. A CERAB configuration using Atrium's Avanta V12 with transparent covers made with a previously established method was inserted into the other phantom, both modified stent-grafts were suitable for laser PIV, enabling visualization of the flow fields and quantification of time average wall shear stress (TAWSS), oscillatory shear index (OSI), and relative residence time (RRT).

RESULTS

Disturbed flow was observed at the bifurcation region of the AFX, especially at the end systolic velocity (ESV) time-point where recirculation was noticeable due to vortical flow. In contrast, predominantly unidirectional flow was observed at the CERAB bifurcation. These observations were confirmed by the quantified hemodynamic results from PIV analysis where mean TAWSS of 0.078 Pa (range: 0.009-0.242 Pa) was significantly lower in AFX as compared with 0.229 Pa (range: 0.013-0.906 Pa) for CERAB (p<0.001). Mean OSI of 0.318 (range: 0.123-0.496) in AFX was significantly higher than 0.252 (range: 0.055-0.472) in CERAB (p<0.001). Likewise, mean RRT of 180 Pa^-1 (range: 9-3603 Pa^-1) in AFX was also significantly higher than 88 Pa^-1 (range: 2-840 Pa^-1) in CERAB (p=0.0086).

CONCLUSION

In this in vitro study, the flow pattern of a modified AFX stent-graft was found to be more disturbed especially at the end systolic phase, its hemodynamic outcomes less desirable than CERAB configuration.

CLINICAL RELEVANCE

While the AFX stent-graft has an advantage over the CERAB configuration in eliminating radial mismatch, and maintaining the anatomical bifurcation for future endovascular intervention, this in vitro study revealed that the associated lower TAWSS, higher OSI and RRT may predispose to thrombosis and are, thus, less desirable as compared to a CERAB configuration. Further investigation is warranted to confirm whether these findings translate into the clinical setting.

False lumen pressure estimation in type B aortic dissection using 4D flow cardiovascular magnetic resonance: comparisons with aortic growth

BACKGROUND

Chronic type B aortic dissection (TBAD) is associated with poor long-term outcome, and accurate risk stratification tools remain lacking. Pressurization of the false lumen (FL) has been recognized as central in promoting aortic growth. Several surrogate imaging-based metrics have been proposed to assess FL hemodynamics; however, their relationship to enlarging aortic dimensions remains unclear. We investigated the association between aortic growth and three cardiovascular magnetic resonance (CMR)-derived metrics of FL pressurization: false lumen ejection fraction (FLEF), maximum systolic deceleration rate (MSDR), and FL relative pressure (FL ΔPmax).

METHODS

CMR/CMR angiography was performed in 12 patients with chronic dissection of the descending thoracoabdominal aorta, including contrast-enhanced CMR angiography and time-resolved three-dimensional phase-contrast CMR (4D Flow CMR). Aortic growth rate was calculated as the change in maximal aortic diameter between baseline and follow-up imaging studies over the time interval, with patients categorized as having either 'stable' (< 3 mm/year) or 'enlarging' (≥ 3 mm/year) growth. Three metrics relating to FL pressurization were defined as: (1) FLEF: the ratio between retrograde and antegrade flow at the TBAD entry tear, (2) MSDR: the absolute difference between maximum and minimum systolic acceleration in the proximal FL, and (3) FL ΔPmax: the difference in absolute pressure between aortic root and distal FL.

RESULTS

FLEF was higher in enlarging TBAD (49.0 ± 17.9% vs. 10.0 ± 11.9%, p = 0.002), whereas FL ΔPmax was lower (32.2 ± 10.8 vs. 57.2 ± 12.5 mmHg/m, p = 0.017). MSDR and conventional anatomic variables did not differ significantly between groups. FLEF showed positive (r = 0.78, p = 0.003) correlation with aortic growth rate whereas FL ΔPmax showed negative correlation (r = - 0.64, p = 0.026). FLEF and FL ΔPmax remained as independent predictors of aortic growth rate after adjusting for baseline aortic diameter.

CONCLUSION

Comparative analysis of three 4D flow CMR metrics of TBAD FL pressurization demonstrated that those that focusing on retrograde flow (FLEF) and relative pressure (FL ΔPmax) independently correlated with growth and differentiated patients with enlarging and stable descending aortic dissections. These results emphasize the highly variable nature of aortic hemodynamics in TBAD patients, and suggest that 4D Flow CMR derived metrics of FL pressurization may be useful to separate patients at highest and lowest risk for progressive aortic growth and complications.

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.

An LDV based method to quantify the error of PC-MRI derived Wall Shear Stress measurement

Wall Shear Stress (WSS) has been demonstrated to be a biomarker of the development of atherosclerosis. In vivo assessment of WSS is still challenging, but 4D Flow MRI represents a promising tool to provide 3D velocity data from which WSS can be calculated. In this study, a system based on Laser Doppler Velocimetry (LDV) was developed to validate new improvements of 4D Flow MRI acquisitions and derived WSS computing. A hydraulic circuit was manufactured to allow both 4D Flow MRI and LDV velocity measurements. WSS profiles were calculated with one 2D and one 3D method. Results indicated an excellent agreement between MRI and LDV velocity data, and thus the set-up enabled the evaluation of the improved performances of 3D with respect to the 2D-WSS computation method. To provide a concrete example of the efficacy of this method, the influence of the spatial resolution of MRI data on derived 3D-WSS profiles was investigated. This investigation showed that, with acquisition times compatible with standard clinical conditions, a refined MRI resolution does not improve WSS assessment, if the impact of noise is unreduced. This study represents a reliable basis to validate with LDV WSS calculation methods based on 4D Flow MRI.

Impact of aortic arch curvature in flow haemodynamics in patients with transposition of the great arteries after arterial switch operation

AIMS 

In this study, we will describe a comprehensive haemodynamic analysis and its relationship to the dilation of the aorta in transposition of the great artery (TGA) patients post-arterial switch operation (ASO) and controls using 4D-flow magnetic resonance imaging (MRI) data.

METHODS AND RESULTS 

Using 4D-flow MRI data of 14 TGA young patients and 8 age-matched normal controls obtained with 1.5 T GE-MR scanner, we evaluate 3D maps of 15 different haemodynamics parameters in six regions; three of them in the aortic root and three of them in the ascending aorta (anterior-left, -right, and posterior for both cases) to find its relationship with the aortic arch curvature and root dilation. Differences between controls and patients were evaluated using Mann-Whitney U test, and the relationship with the curvature was accessed by unpaired t-test. For statistical significance, we consider a P-value of 0.05. The aortic arch curvature was significantly different between patients 46.238 ± 5.581 m-1 and controls 41.066 ± 5.323 m-1. Haemodynamic parameters as wall shear stress circumferential (WSS-C), and eccentricity (ECC), were significantly different between TGA patients and controls in both the root and ascending aorta regions. The distribution of forces along the ascending aorta is highly inhomogeneous in TGA patients. We found that the backward velocity (B-VEL), WSS-C, velocity angle (VEL-A), regurgitation fraction (RF), and ECC are highly correlated with the aortic arch curvature and root dilatation.

CONCLUSION 

We have identified six potential biomarkers (B-VEL, WSS-C, VEL-A, RF, and ECC), which may be helpful for follow-up evaluation and early prediction of aortic root dilatation in this patient population.

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

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

Low and Oscillatory Wall Shear Stress Is Not Related to Aortic Dilation in Patients With Bicuspid Aortic Valve: A Time-Resolved 3-Dimensional Phase-Contrast Magnetic Resonance Imaging Study

Supplemental Digital Content is available in the text.

Objective:

To assess the relationship between regional wall shear stress (WSS) and oscillatory shear index (OSI) and aortic dilation in patients with bicuspid aortic valve (BAV).

Approach and Results:

Forty-six consecutive patients with BAV (63% with right-left-coronary-cusp fusion, aortic diameter ≤ 45 mm and no severe valvular disease) and 44 healthy volunteers were studied by time-resolved 3-dimensional phase-contrast magnetic resonance imaging. WSS and OSI were quantified at different levels of the ascending aorta and the aortic arch, and regional WSS and OSI maps were obtained. Seventy percent of BAV had ascending aorta dilation. Compared with healthy volunteers, patients with BAV had increased WSS and decreased OSI in most of the ascending aorta and the aortic arch. In both BAV and healthy volunteers, regions of high WSS matched regions of low OSI and vice versa. No regions of both low WSS and high OSI were identified in BAV compared with healthy volunteers. Patients with BAV with dilated compared with nondilated aorta presented low and oscillatory WSS in the aortic arch, but not in the ascending aorta where dilation is more prevalent. Furthermore, no regions of concomitant low WSS and high OSI were identified when BAV were compared according to leaflet fusion pattern, despite the well-known differences in regional dilation prevalence.

Conclusions:

Regions with low WSS and high OSI do not match those with the highest prevalence of dilation in patients with BAV, thus providing no evidence to support the low and oscillatory shear stress theory in the pathogenesis of proximal aorta dilation in the presence of BAV.

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.

Aortic flow patterns and wall shear stress maps by 4D-flow cardiovascular magnetic resonance in the assessment of aortic dilatation in bicuspid aortic valve disease

BACKGROUND

In patients with bicuspid valve (BAV), ascending aorta (AAo) dilatation may be caused by altered flow patterns and wall shear stress (WSS). These differences may explain different aortic dilatation morphotypes. Using 4D-flow cardiovascular magnetic resonance (CMR), we aimed to analyze differences in flow patterns and regional axial and circumferential WSS maps between BAV phenotypes and their correlation with ascending aorta dilatation morphotype.

METHODS

One hundred and one BAV patients (aortic diameter ≤ 45 mm, no severe valvular disease) and 20 healthy subjects were studied by 4D-flow CMR. Peak velocity, flow jet angle, flow displacement, in-plane rotational flow (IRF) and systolic flow reversal ratio (SFRR) were assessed at different levels of the AAo. Peak-systolic axial and circumferential regional WSS maps were also estimated. Unadjusted and multivariable adjusted linear regression analyses were used to identify independent correlates of aortic root or ascending dilatation. Age, sex, valve morphotype, body surface area, flow derived variables and WSS components were included in the multivariable models.

RESULTS

The AAo was non-dilated in 24 BAV patients and dilated in 77 (root morphotype in 11 and ascending in 66). BAV phenotype was right-left (RL-) in 78 patients and right-non-coronary (RN-) in 23. Both BAV phenotypes presented different outflow jet direction and velocity profiles that matched the location of maximum systolic axial WSS. RL-BAV velocity profiles and maximum axial WSS were homogeneously distributed right-anteriorly, however, RN-BAV showed higher variable profiles with a main proximal-posterior distribution shifting anteriorly at mid-distal AAo. Compared to controls, BAV patients presented similar WSS magnitude at proximal, mid and distal AAo (p = 0.764, 0.516 and 0.053, respectively) but lower axial and higher circumferential WSS components (p < 0.001 for both, at all aortic levels). Among BAV patients, RN-BAV presented higher IRF at all levels (p = 0.024 proximal, 0.046 mid and 0.002 distal AAo) and higher circumferential WSS at mid and distal AAo (p = 0.038 and 0.046, respectively) than RL-BAV. However, axial WSS was higher in RL-BAV compared to RN-BAV at proximal and mid AAo (p = 0.046, 0.019, respectively). Displacement and axial WSS were independently associated with the root-morphotype, and circumferential WSS and SFRR with the ascending-morphotype.

CONCLUSIONS

Different BAV-phenotypes present different flow patterns with an anterior distribution in RL-BAV, whereas, RN-BAV patients present a predominant posterior outflow jet at the sinotubular junction that shifts to anterior or right anterior in mid and distal AAo. Thus, RL-BAV patients present a higher axial WSS at the aortic root while RN-BAV present a higher circumferential WSS in mid and distal AAo. These results may explain different AAo dilatation morphotypes in the BAV population.

Wall shear stress estimation in the aorta: Impact of wall motion, spatiotemporal resolution, and phase noise

BACKGROUND

Wall shear stress (WSS) presents an important parameter for assessing blood flow characteristics and evaluating flow-mediated lesions in the aorta.

PURPOSE

To investigate the robustness of WSS and oscillatory shear index (OSI) estimation based on 4D flow MRI against vessel wall motion, spatiotemporal resolution, and velocity encoding (VENC).

STUDY TYPE

Simulated and prospective.

POPULATION

Synthetic 4D flow MRI data of the aorta, simulated using the Lattice-Boltzmann method; in vivo 4D flow MRI data of the aorta from healthy volunteers (n = 11) and patients with congenital heart defects (n = 17).

FIELD STRENGTH/SEQUENCE

1.5T; 4D flow MRI with PEAK-GRAPPA acceleration and prospective electrocardiogram triggering.

ASSESSMENT

Predicated upon 3D cubic B-splines interpolation of the image velocity field, WSS was estimated in mid-systole, early-diastole, and late-diastole and OSI was derived. We assessed the impact of spatiotemporal resolution and phase noise, and compared results based on tracked-using deformable registration-and static vessel wall location.

STATISTICAL TESTS

Bland-Altman analysis to assess WSS/OSI differences; Hausdorff distance (HD) to assess wall motion; and Pearson's correlation coefficient (PCC) to assess correlation of HD with WSS.

RESULTS

Synthetic data results show systematic over-/underestimation of WSS when different spatial resolution (mean ± 1.96 SD up to -0.24 ± 0.40 N/m2 and 0.5 ± 1.38 N/m2 for 8-fold and 27-fold voxel size, respectively) and VENC-depending phase noise (mean ± 1.96 SD up to 0.31 ± 0.12 N/m2 and 0.94 ± 0.28 N/m2 for 2-fold and 4-fold VENC increase, respectively) are given. Neglecting wall motion when defining the vessel wall perturbs WSS estimates to a considerable extent (1.96 SD up to 1.21 N/m2 ) without systematic over-/underestimation (Bland-Altman mean range -0.06 to 0.05).

DATA CONCLUSION

In addition to sufficient spatial resolution and velocity to noise ratio, accurate tracking of the vessel wall is essential for reliable image-based WSS estimation and should not be neglected if wall motion is present.

LEVEL OF EVIDENCE

2 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2018.

A FSI computational framework for vascular physiopathology: A novel flow-tissue multiscale strategy

A novel fluid-structure computational framework for vascular applications is herein presented. It is developed by combining the double multi-scale nature of vascular physiopathology in terms of both tissue properties and blood flow. Addressing arterial tissues, they are modelled via a nonlinear multiscale constitutive rationale, based only on parameters having a clear histological and biochemical meaning. Moreover, blood flow is described by coupling a three-dimensional fluid domain (undergoing physiological inflow conditions) with a zero-dimensional model, which allows to reproduce the influence of the downstream vasculature, furnishing a realistic description of the outflow proximal pressure. The fluid-structure interaction is managed through an explicit time-marching approach, able to accurately describe tissue nonlinearities within each computational step for the fluid problem. A case study associated to a patient-specific aortic abdominal aneurysmatic geometry is numerically investigated, highlighting advantages gained from the proposed multiscale strategy, as well as showing soundness and effectiveness of the established framework for assessing useful clinical quantities and risk indexes.

Enhancing the Velocity Data From 4D Flow MR Images by Reducing its Divergence

Velocity measurements from 4D flow MRI are prone to be affected by several imperfections of the MR system. Assuming that blood is incompressible, we propose a novel method for enhancing the velocity field by reducing its divergence. To enhance the velocity data, we added a corrector velocity to each voxel such that the divergence is minimized. The method was validated using an analytical Womersley flow model for different settings of resolution and noise levels. The performance of the proposed method was also assessed in volunteers and patients. Results demonstrated a significant reduction of the divergence depending on the size of the regularization term, obtaining a reduction close to 50% of the mean divergence with negligible modification of flow parameters. Remarkably, we found that the reduction of the divergence, in percentage, was independent of volunteers, resolution or noise.