Aortic valve disease is often treated with bioprosthetic valves. An alternative treatment is aortic valve neocuspidization which is a relatively new reparative procedure whereby the three aortic cusps are replaced with patient pericardium or bovine tissues. Recent research indicates that aortic blood flow is disturbed, and turbulence effects have yet to be evaluated in either bioprosthetic or aortic valve neocuspidization valve types in patient-specific settings. The aim of this study is to better understand turbulence production in the aorta and evaluate its effects on laminar and turbulent wall shear stress. Four patients with aortic valve disease were treated with either bioprosthetic valves (n=2) or aortic valve neocuspidization valvular repair (n=2). Aortic geometries were segmented from magnetic resonance images (MRI), and 4D flow MRI was used to derive physiological inlet and outlet boundary conditions. Pulsatile large-eddy simulations were performed to capture the full range of laminar, transitional and turbulence characteristics in the aorta. Turbulence was produced in all aortas with highest levels occurring during systolic deceleration. In the ascending aorta, turbulence production is attributed to a combination of valvular skew, valvular eccentricity, and ascending aortic dilation. In the proximal descending thoracic aorta, turbulence production is dependent on the type of arch-descending aorta connection (e.g., a narrowing or sharp bend) which induces flow separation. Laminar and turbulent wall shear stresses are of similar magnitude throughout late systolic deceleration and diastole, although turbulent wall shear stress magnitudes exceed laminar wall shear stresses between 27.3% and 61.1% of the cardiac cycle. This emphasises the significance of including turbulent wall shear stress to improve our comprehension of progressive arterial wall diseases. The findings of this study recommend that aortic valve treatments should prioritise minimising valvular eccentricity and skew in order to mitigate turbulence generation.

The opening and closing dynamics of the aortic valve (AV) has a strong influence on haemodynamics in the aortic root, and both play a pivotal role in maintaining normal physiological functions of the valve. The aim of this study was to establish a subject-specific fluid–structure interaction (FSI) workflow capable of simulating the motion of a tricuspid healthy valve and the surrounding haemodynamics under physiologically realistic conditions. A subject-specific aortic root was reconstructed from magnetic resonance (MR) images acquired from a healthy volunteer, whilst the valve leaflets were built using a parametric model fitted to the subject-specific aortic root geometry. The material behaviour of the leaflets was described using the isotropic hyperelastic Ogden model, and subject-specific boundary conditions were derived from 4D-flow MR imaging (4D-MRI). Strongly coupled FSI simulations were performed using a finite volume-based boundary conforming method implemented in FlowVision. Our FSI model was able to simulate the opening and closing of the AV throughout the entire cardiac cycle. Comparisons of simulation results with 4D-MRI showed a good agreement in key haemodynamic parameters, with stroke volume differing by 7.5% and the maximum jet velocity differing by less than 1%. Detailed analysis of wall shear stress (WSS) on the leaflets revealed much higher WSS on the ventricular side than the aortic side and different spatial patterns amongst the three leaflets.

Purpose: Finite element analysis (FEA) has been used to predict wall stress in ascending thoracic aortic aneurysm (ATAA) in order to evaluate risk of dissection or rupture. Patient-specific FEA requires detailed information on ATAA geometry, loading conditions, material properties, and wall thickness. Unfortunately, measuring aortic wall thickness and mechanical properties non-invasively poses a significant challenge, necessitating the use of non-patient-specific data in most FE simulations. This study aimed to assess the impact of employing non-patient-specific material properties and wall thickness on ATAA wall stress predictions. Methods: FE simulations were performed on 13 ATAA geometries reconstructed from computed tomography angiography (CTA) images. Patient-specific material properties and wall thicknesses were made available from a previous study where uniaxial tensile testing was performed on tissue samples obtained from the same patients. The ATAA wall models were discretised with hexahedral elements and prestressed. For each ATAA model, FE simulations were conducted using patient-specific material properties and wall thicknesses, and group-mean values derived from all tissue samples included in the same experimental study. Literature-based material property and wall thickness were also obtained from the literature and applied to 4 representative cases. Additional FE simulations were performed on these 4 cases by employing group-mean and literature-based wall thicknesses. Results: FE simulations using the group-mean material property produced peak wall stresses comparable to those obtained using patient-specific material properties, with a mean deviation of 7.8%. Peak wall stresses differed by 20.8% and 18.7% in patients with exceptionally stiff or compliant walls, respectively. Comparison to results using literature-based material properties revealed larger discrepancies, ranging from 5.4% to 28.0% (mean 20.1%). Bland-Altman analysis showed significant discrepancies in areas of high wall stress, where wall stress obtained using patient-specific and literature-based properties differed by up to 674 kPa, compared to 227 kPa between patient-specific and group-mean properties. Regarding wall thickness, using the literature-based value resulted in even larger discrepancies in predicted peak stress, ranging from 24.2% to 30.0% (mean 27.3%). Again, using the group-mean wall thickness offered better predictions with a difference less than 5% in three out of four cases. While peak wall stresses were most affected by the choice of mechanical properties or wall thickness, the overall distribution of wall stress hardly changed. Conclusions: Our study demonstrated the importance of incorporating patient-specific material properties and wall thickness in FEA for risk prediction of aortic dissection or rupture. Our future efforts will focus on developing inverse methods for non-invasive determination of patient-specific wall material parameters and wall thickness.

Background and Objective: Numerical simulations of blood flow are a valuable tool to investigate the pathophysiology of ascending thoratic aortic aneurysms (ATAA). To accurately reproduce in vivo hemodynamics, computational fluid dynamics (CFD) models must employ realistic inflow boundary conditions (BCs). However, the limited availability of in vivo velocity measurements, still makes researchers resort to idealized BCs. The aim of this study was to generate and thoroughly characterize a large dataset of synthetic 4D aortic velocity profiles sampled on a 2D cross-section along the ascending aorta with features similar to clinical cohorts of patients with ATAA. Methods: Time-resolved 3D phase contrast magnetic resonance (4D flow MRI) scans of 30 subjects with ATAA were processed through in-house code to extract anatomically consistent cross-sectional planes along the ascending aorta, ensuring spatial alignment among all planes and interpolating all velocity fields to a reference configuration. Velocity profiles of the clinical cohort were extensively characterized by computing flow morphology descriptors of both spatial and temporal features. By exploiting principal component analysis (PCA), a statistical shape model (SSM) of 4D aortic velocity profiles was built and a dataset of 437 synthetic cases with realistic properties was generated. Results: Comparison between clinical and synthetic datasets showed that the synthetic data presented similar characteristics as the clinical population in terms of key morphological parameters. The average velocity profile qualitatively resembled a parabolic-shaped profile, but was quantitatively characterized by more complex flow patterns which an idealized profile would not replicate. Statistically significant correlations were found between PCA principal modes of variation and flow descriptors. Conclusions: We built a data-driven generative model of 4D aortic inlet velocity profiles, suitable to be used in computational studies of blood flow. The proposed software system also allows to map any of the generated velocity profiles to the inlet plane of any virtual subject given its coordinate set.

Introduction: Thoracic endovascular aortic repair (TEVAR) of the arch is challenging given its complex geometry and the involvement of supra-aortic arteries. Different branched endografts have been designed for use in this region, but their haemodynamic performance and the risk for post-intervention complications are not yet clear. This study aims to examine aortic haemodynamics and biomechanical conditions following TVAR treatment of an aortic arch aneurysm with a two-component single-branched endograft. Methods: Computational fluid dynamics and finite element analysis were applied to a patient-specific case at different stages: pre-intervention, post-intervention and follow-up. Physiologically accurate boundary conditions were used based on available clinical information. Results: Computational results obtained from the post-intervention model confirmed technical success of the procedure in restoring normal flow to the arch. Simulations of the follow-up model, where boundary conditions were modified to reflect change in supra-aortic vessel perfusion observed on the follow-up scan, predicted normal flow patterns but high levels of wall stress (up to 1.3M MPa) and increased displacement forces in regions at risk of compromising device stability. This might have contributed to the suspected endoleaks or device migration identified at the final follow up. Discussion: Our study demonstrated that detailed haemodynamic and biomechanical analysis can help identify possible causes for post-TEVAR complications in a patient-specific setting. Further refinement and validation of the computational workflow will allow personalised assessment to aid in surgical planning and clinical decision making.