No abstract

Bicuspid aortic valve (BAV) is the most common congenital cardiac disease and is a foremost risk factor for aortopathies. Despite the genetic basis of BAV and of the associated aortopathies, BAV-related alterations in aortic fluid-dynamics, and particularly in wall shear stresses (WSSs), likely play a role in the progression of aortopathy, and may contribute to its pathogenesis. To test whether WSS may trigger aortopathy, in this study we used 4D Flow sequences of phase-contrast cardiac magnetic resonance imaging (CMR) to quantitatively compare the in vivo fluid dynamics in the thoracic aorta of two groups of subjects: (i) five prospectively enrolled young patients with normo-functional BAV and with no aortic dilation and (ii) ten age-matched healthy volunteers. Through the semi-automated processing of 4D Flow data, the aortic bulk flow at peak systole was quantified, and WSSs acting on the endothelium of the ascending aorta were characterized throughout the systolic phase in terms of magnitude and time-dependency through a method recently developed by our group. Variables computed for each BAV patient were compared vs. the corresponding distribution of values obtained for healthy controls. In BAV patients, ascending aorta diameter was measured on cine-CMR images at baseline and at 3-year follow-up. As compared to controls, normo-functional BAV patients were characterized by minor bulk flow disturbances at peak systole. However, they were characterized by evident alterations of WSS distribution and peak values in the ascending aorta. In particular, in four BAV patients, who were characterized by right-left leaflet fusion, WSS peak values exceeded by 27–46% the 90th percentile of the distribution obtained for healthy volunteers. Only in the BAV patient with right-non-coronary leaflet fusion the same threshold was exceeded by 132%. Also, evident alterations in the time-dependency of WSS magnitude and direction were observed. Despite, these fluid-dynamic alterations, no clinically relevant anatomical remodeling was observed in the BAV patients at 3-year follow-up. In light of previous evidence from the literature, our results suggest that WSS alterations may precede the onset of aortopathy and may contribute to its triggering, but WSS-driven anatomical remodeling, if any, is a very slow process.

Blood flow in the aorta is helical, but most computational studies ignore the presence of secondary flow components at the ascending aorta (AAo) inlet. The aim of this study is to ascertain the importance of inlet boundary conditions (BCs) in computational analysis of flow patterns in the thoracic aorta based on patient-specific images, with a particular focus on patients with an abnormal aortic valve. Two cases were studied: one presenting a severe aortic valve stenosis and the other with a mechanical valve. For both aorta models, three inlet BCs were compared; these included the flat profile and 1D through-plane velocity and 3D phase-contrast magnetic resonance imaging derived velocity profiles, with the latter being used for benchmarking. Our results showed that peak and mean velocities at the proximal end of the ascending aorta were underestimated by up to 41% when the secondary flow components were neglected. The results for helical flow descriptors highlighted the strong influence of secondary velocities on the helical flow structure in the AAo. Differences in all wall shear stress (WSS)-derived indices were much more pronounced in the AAo and aortic arch (AA) than in the descending aorta (DAo). Overall, this study demonstrates that using 3D velocity profiles as inlet BC is essential for patient-specific analysis of hemodynamics and WSS in the AAo and AA in the presence of an abnormal aortic valve. However, predicted flow in the DAo is less sensitive to the secondary velocities imposed at the inlet; hence, the 1D through-plane profile could be a sufficient inlet BC for studies focusing on distal regions of the thoracic aorta.

Computational hemodynamics studies of aortic dissections usually combine patient-specific geometries with idealized or generic boundary conditions. In this study we present a comprehensive methodology for simulations of hemodynamics in type B aortic dissection (TBAD) based on fully patient-specific boundary conditions. Methods: Pre-operative 4D flow magnetic resonance imaging (MRI) and Doppler-wire pressure measurements (pre-and post-operative) were acquired from a TBAD patient. These data were used to derive boundary conditions for computational modelling of flow before and after thoracic endovascular repair (TEVAR). Validations of the computational results were performed by comparing predicted flow patterns with pre-TEVAR 4D flow MRI, as well as pressures with in vivo measurements. Results and Conclusion: Comparison of instantaneous velocity streamlines showed a good qualitative agreement with 4D flow MRI. Quantitative comparison of predicted pressures with pressure measurements revealed a maximum difference of 11 mmHg (-9.7%). Furthermore, our model correctly predicted the reduction of true lumen pressure from 74/115 mmHg pre-TEVAR to 64/107 mmHg post-TEVAR (diastolic/systolic pressures at entry tear level), compared to the corresponding measurements of 72/118 mmHg and 64/114 mmHg. This demonstrates that pre-TEVAR 4D flow MRI can be used to tune boundary conditions for post-TEVAR hemodynamic analyses.

Bicuspid aortic valve (BAV), i.e. the fusion of two aortic valve cusps, is the most frequent congenital cardiac malformation. Its progression is often characterized by accelerated leaflet calcification and aortic wall dilation. These processes are likely enhanced by altered biomechanical stimuli, including fluid-dynamic wall shear stresses (WSS) acting on both the aortic wall and the aortic valve. Several studies have proposed the exploitation of 4D-flow magnetic resonance imaging sequences to characterize abnormal in vivo WSS in BAV-affected patients, to support prognosis and timing of intervention. However, current methods fail to quantify WSS peak values. On this basis, we developed two new methods for the improved quantification of in vivo WSS acting on the aortic wall based on 4D-flow data. We tested both methods separately and in combination on synthetic datasets obtained by two computational fluid-dynamics (CFD) models of the aorta with healthy and bicuspid aortic valve. Tests highlighted the need for data spatial resolution at least comparable to current clinical guidelines, the low sensitivity of the methods to data noise, and their capability, when used jointly, to compute more realistic peak WSS values as compared to state-of-the-art methods. The integrated application of the two methods on the real 4D-flow data from a preliminary cohort of three healthy volunteers and three BAV-affected patients confirmed these indications. In particular, quantified WSS peak values were one order of magnitude higher than those reported in previous 4D-flow studies, and much closer to those computed by highly time- and space-resolved CFD simulations.

In order for computational fluid dynamics to provide quantitative parameters to aid in the clinical assessment of type B aortic dissection, the results must accurately mimic the hemodynamic environment within the aorta. The choice of inlet velocity profile (IVP) therefore is crucial; however, idealised profiles are often adopted, and the effect of IVP on hemodynamics in a dissected aorta is unclear. This study examined two scenarios with respect to the influence of IVP—using (a) patient-specific data in the form of a three-directional (3D), through-plane (TP) or flat IVP; and (b) non-patient-specific flow waveform. The results obtained from nine simulations using patient-specific data showed that all forms of IVP were able to reproduce global flow patterns as observed with 4D flow magnetic resonance imaging. Differences in maximum velocity and time-averaged wall shear stress near the primary entry tear were up to 3% and 6%, respectively, while pressure differences across the true and false lumen differed by up to 6%. More notable variations were found in regions of low wall shear stress when the primary entry tear was close to the left subclavian artery. The results obtained with non-patient-specific waveforms were markedly different. Throughout the aorta, a 25% reduction in stroke volume resulted in up to 28% and 35% reduction in velocity and wall shear stress, respectively, while the shape of flow waveform had a profound influence on the predicted pressure. The results of this study suggest that 3D, TP and flat IVPs all yield reasonably similar velocity and time-averaged wall shear stress results, but TP IVPs should be used where possible for better prediction of pressure. In the absence of patient-specific velocity data, effort should be made to acquire patient’s stroke volume and adjust the applied IVP accordingly.

Blood flow in the aorta is often assumed laminar, however aortic valve pathologies may induce transition to turbulence and our understanding of turbulence effects is incomplete. The aim of the study was to provide a detailed analysis of turbulence effects in aortic valve stenosis (AVS). Methods Large-eddy simulation (LES) of flow through a patient-specific aorta with AVS was conducted. Magnetic resonance imaging (MRI) was performed and used for geometric reconstruction and patient-specific boundary conditions. Computed velocity field was compared with 4D flow MRI to check qualitative and quantitative consistency. The effect of turbulence was evaluated in terms of fluctuating kinetic energy, turbulence-related wall shear stress (WSS) and energy loss. Results Our analysis suggested that turbulence was induced by a combination of a high velocity jet impinging on the arterial wall and a dilated ascending aorta which provided sufficient space for turbulence to develop. Turbulent WSS contributed to 40% of the total WSS in the ascending aorta and 38% in the entire aorta. Viscous and turbulent irreversible energy losses accounted for 3.9 and 2.7% of the total stroke work, respectively. Conclusions This study demonstrates the importance of turbulence in assessing aortic haemodynamics in a patient with AVS. Neglecting the turbulent contribution to WSS could potentially result in a significant underestimation of the total WSS. Further work is warranted to extend the analysis to more AVS cases and patients with other aortic valve diseases.

Severity of aortic coarctation (CoA) is currently assessed by estimating trans-coarctation pressure drops through cardiac catheterization or echocardiography. In principle, more detailed information could be obtained non-invasively based on space-and time-resolved magnetic resonance imaging (4D flow) data. Yet the limitations of this imaging technique require testing the accuracy of 4D flow-derived hemodynamic quantities against other methodologies. With the objective of assessing the feasibility and accuracy of this non-invasive method to support the clinical diagnosis of CoA, we developed an algorithm (4DF-FEPPE) to obtain relative pressure distributions from 4D flow data by solving the Poisson pressure equation. 4DF-FEPPE was tested against results from a patient-specific fluid-structure interaction (FSI) simulation, whose patient-specific boundary conditions were prescribed based on 4D flow data. Since numerical simulations provide noise-free pressure fields on fine spatial and temporal scales, our analysis allowed to assess the uncertainties related to 4D flow noise and limited resolution. 4DF-FEPPE and FSI results were compared on a series of cross-sections along the aorta. Bland-Altman analysis revealed very good agreement between the two methodologies in terms of instantaneous data at peak systole, end-diastole and time-averaged values: biases (means of differences) were +0.4 mmHg,-1.1 mmHg and +0.6 mmHg, respectively. Limits of agreement (2 SD) were ±0.978 mmHg, ±1.06 mmHg and ±1.97 mmHg, respectively. Peak-to-peak and maximum trans-coarctation pressure drops obtained with 4DF-FEPPE differed from FSI results by 0.75 mmHg and-1.34 mmHg respectively. The present study considers important validation aspects of non-invasive pressure difference estimation based on 4D flow MRI, showing the potential of this technology to be more broadly applied to the clinical practice.