Computational Analysis of Balloon Catheter Behaviour at Variable Inflation Levels

Aortic valvuloplasty is a minimally invasive procedure for the dilatation of stenotic aortic valves. Rapid ventricular pacing is an established technique for balloon stabilization during this procedure. However, low cardiac output due to the pacing is one of the inherent risks, which is also associated with several potential complications. This paper proposes a numerical modelling approach to understand the effect of different inflation levels of a valvuloplasty balloon catheter on the positional instability caused by a pulsating blood flow. An unstretched balloon catheter model was crimped into a tri-folded configuration and inflated to several levels. Ten different inflation levels were then tested, and a Fluid-Structure Interaction model was built to solve interactions between the balloon and the blood flow modelled in an idealised aortic arch. Our computational results show that the maximum displacement of the balloon catheter increases with the inflation level, with a small step at around 50% inflation and a sharp increase after reaching 85% inflation. This work represents a substantial progress towards the use of simulations to solve the interactions between a balloon catheter and pulsating blood flow.

Effect of Valve Height on the Opening and Closing Performance of the Aortic Valve Under Aortic Root Dilatation

Patients with aortic valve disease can suffer from valve insufficiency after valve repair surgery due to aortic root dilatation. The paper investigates the effect of valve height (Hv) on the aortic valve opening and closing in order to select the appropriate range of Hv for smoother blood flow through the aortic valve and valve closure completely in the case of continuous aortic root dilatation. A total of 20 parameterized three-dimensional models of the aortic root were constructed following clinical surgical guidance. Aortic annulus diameter (DAA) was separately set to 26, 27, 28, 29, and 30 mm to simulate aortic root dilatation. H_V value was separately set to 13.5, 14, 14.5, and 15 mm to simulate aortic valve alterations in surgery. Time-varying pressure loads were applied to the valve, vessel wall of the ascending aorta, and left ventricle. Then, finite element analysis software was employed to simulate the movement and mechanics of the aortic root. The feasible design range of the valve size was evaluated using maximum stress, geometric orifice area (GOA), and leaflet contact force. The results show that the valve was incompletely closed when H_V was 13.5 mm and D_AA was 29 or 30 mm. The GOA of the valve was small when H_V was 15 mm and D_AA was 26 or 27 mm. The corresponding values of the other models were within the normal range. Compared with the model with an H_V of 14 mm, the model with an H_V of 14.5 mm could effectively reduce maximum stress and had relatively larger GOA and less change in contact force. As a result, valve height affects the performance of aortic valve opening and closing. Smaller H_V is adapted to smaller D_AA and vice versa. When H_V is 14.5 mm, the valve is well adapted to the dilatation of the aortic root to enhance repair durability. Therefore, more attention should be paid to H_V in surgical planning.

A Fluid-Structure Interaction Study of Different Bicuspid Aortic Valve Phenotypes Throughout the Cardiac Cycle

The bicuspid aortic valve (BAV) is a congenital malformation of the aortic valve with a variety of structural features. The current research on BAV mainly focuses on the systolic phase, while ignoring the diastolic hemodynamic characteristics and valve mechanics. The purpose of this study is to compare the differences in hemodynamics and mechanical properties of BAV with different phenotypes throughout the cardiac cycle by means of numerical simulation. Based on physiological anatomy, we established an idealized tricuspid aortic valve (TAV) model and six phenotypes of BAV models (including Type 0 a-p, Type 0 lat, Type 1 L-R, Type 1 N-L, Type 1 R-N, and Type 2), and simulated the dynamic changes of the aortic valve during the cardiac cycle using the fluid-structure interaction method. The morphology of the leaflets, hemodynamic parameters, flow patterns, and strain were analyzed. Compared with TAV, the cardiac output and effective orifice area of different BAV phenotypes decreased certain degree, along with the peak velocity and mean pressure difference increased both. Among all BAV models, Type 2 exhibited the worst hemodynamic performance. During the systole, obvious asymmetric flow field was observed in BAV aorta, which was related to the orientation of BAV. Higher strain was generated in diastole for BAV models. The findings of this study suggests specific differences in the hemodynamic characteristics and valve mechanics of different BAV phenotypes, including different severity of stenosis, flow patterns, and leaflet strain, which may be critical for prediction of other subsequent aortic diseases and differential treatment strategy for certain BAV phenotype.

Biomechanical Assessment of Bicuspid Aortic Valve Phenotypes: A Fluid-Structure Interaction Modelling Approach

PURPOSE

Bicuspid aortic valve (BAV) is a congenital heart malformation with phenotypic heterogeneity. There is no prior computational study that assesses the haemodynamic and valve mechanics associated with BAV type 2 against a healthy tricuspid aortic valve (TAV) and other BAV categories.

METHODS

A proof-of-concept study incorporating three-dimensional fluid-structure interaction (FSI) models with idealised geometries (one TAV and six BAVs, namely type 0 with lateral and anterior-posterior orientations, type 1 with R-L, N-R and N-L leaflet fusion and type 2) has been developed. Transient physiological boundary conditions have been applied and simulations were run using an Arbitrary Lagrangian-Eulerian formulation.

RESULTS

Our results showed the presence of abnormal haemodynamics in the aorta and abnormal valve mechanics: type 0 BAVs yielded the best haemodynamical and mechanical outcomes, but cusp stress distribution varied with valve orifice orientation, which can be linked to different cusp calcification location onset; type 1 BAVs gave rise to similar haemodynamics and valve mechanics, regardless of raphe position, but this position altered the location of abnormal haemodynamic features; finally, type 2 BAV constricted the majority of blood flow, exhibiting the most damaging haemodynamic and mechanical repercussions when compared to other BAV phenotypes.

CONCLUSION

The findings of this proof-of-concept work suggest that there are specific differences across haemodynamics and valve mechanics associated with BAV phenotypes, which may be critical to subsequent processes associated with their pathophysiology processes.

NUMERICAL SIMULATION OF DILATION PATTERNS OF THE ASCENDING AORTA IN AORTOPATHIES

Aortic dilation is associated with congenital bicuspid aortic valve (BAV) disease, and its etiology is still not completely understood. The aim of this study is to provide further insight into aortic hemodynamics in a BAV population with different degrees of aortic dilation and regurgitation in comparison with a patient without pathology. A fluid–structure interaction (FSI) numerical approach is implemented regarding patient-specific geometries, where the aortic valves are defined by analytical orifices. Results show that, while the patient without pathology displays a typical hemodynamic behavior of flows in bends, BAV-related aortas present an accelerated flow along the outer aortic wall. Wall shear stress (WSS) overload in the outer curvature is observed, more marked in more dilated aortas. Moreover, helices in the ascending aorta are present in these patients, enhanced with greater dilation. These findings support the fact that hemodynamic factors play an important role in aortic dilation onset and development in BAV patients, caused by a prolonged exposure of the outer ascending aortic curvature to altered WSS. Besides, our results suggest that greater aortic regurgitation may be associated with abnormal WSS distributions in the ascending aorta during diastole, which can facilitate aortic root dilation.

Design and Computational Validation of a Novel Bioreactor for Conditioning Vascular Tissue to Time-Varying Multidirectional Fluid Shear Stress

PURPOSE

The cardiovascular endothelium experiences pulsatile and multidirectional fluid wall shear stress (WSS). While the effects of non-physiologic WSS magnitude and pulsatility on cardiovascular function have been studied extensively, the impact of directional abnormalities remains unknown due to the challenge to replicate this characteristic in vitro. To address this gap, this study aimed at designing a bioreactor capable of subjecting cardiovascular tissue to time-varying WSS magnitude and directionality.

METHODS

The device consisted of a modified cone-and-plate bioreactor. The cone rotation generates a fluid flow subjecting tissue to desired WSS magnitude, while WSS directionality is achieved by altering the alignment of the tissue relative to the flow at each instant of time. Computational fluid dynamics was used to verify the device ability to replicate the native WSS of the proximal aorta. Cone and tissue mount velocities were determined using an iterative optimization procedure.

RESULTS

Using conditions derived from cone-and-plate theory, the initial simulations yielded root-mean-square errors of 22.8 and 8.4% in WSS magnitude and angle, respectively, between the predicted and the target signals over one cycle, relative to the time-averaged target values. The conditions obtained after two optimization iterations reduced those errors to 3.5 and 0.5%, respectively, and generated 0.2% and 0.01% difference in time-averaged WSS magnitude and angle, respectively, relative to the target waveforms.

CONCLUSIONS

A bioreactor capable of generating simultaneously desired time-varying WSS magnitude and directionality was designed and validated computationally. The ability to subject tissue to in vivo-like WSS will provide new insights into cardiovascular mechanobiology and disease.

Fluid-Structure Interaction Models of Bicuspid Aortic Valves: The Effects of Nonfused Cusp Angles

Bicuspid aortic valve (BAV) is the most common type of congenital heart disease, occurring in 0.5-2% of the population, where the valve has only two rather than the three normal cusps. Valvular pathologies, such as aortic regurgitation and aortic stenosis, are associated with BAVs, thereby increasing the need for a better understanding of BAV kinematics and geometrical characteristics. The aim of this study is to investigate the influence of the nonfused cusp (NFC) angle in BAV type-1 configuration on the valve's structural and hemodynamic performance. Toward that goal, a parametric fluid-structure interaction (FSI) modeling approach of BAVs is presented. Four FSI models were generated with varying NFC angles between 120 deg and 180 deg. The FSI simulations were based on fully coupled structural and fluid dynamic solvers and corresponded to physiologic values, including the anisotropic hyper-elastic behavior of the tissue. The simulated angles led to different mechanical behavior, such as eccentric jet flow direction with a wider opening shape that was found for the smaller NFC angles, while a narrower opening orifice followed by increased jet flow velocity was observed for the larger NFC angles. Smaller NFC angles led to higher concentrated flow shear stress (FSS) on the NFC during peak systole, while higher maximal principal stresses were found in the raphe region during diastole. The proposed biomechanical models could explain the early failure of BAVs with decreased NFC angles, and suggests that a larger NFC angle is preferable in suture annuloplasty BAV repair surgery.

Impact of the bicuspid aortic valve on aortic root haemodynamics: three-dimensional computed fluid dynamics simulation

OBJECTIVES

The aim was to evaluate the impact of a bicuspid aortic valve (BAV) on local shear stress and on the pressure profile on the elements of the aortic root (AoR).

METHODS

The experiment setup included a BAV with aortic valve stenosis (n = 5 pigs, 67 ± 3.5 kg) and insufficiency (n = 5 pigs, 66.7 ± 4.4 kg). By implanting 6 high-fidelity microsonometric crystals in each AoR, we determined the 3-dimensional (3D) geometry of the AoR. Experimental and geometric data were used to create a 3D time- and pressure-dependent computed fluid dynamic model of the AoR with the BAV.

RESULTS

3D AoR geometry was determined by AoR tilt (α) and rotation angle (β). Both values were maximal at the end of diastole: 24.41 ± 1.70° (α) and 20.90 ± 2.11° (β) for BAV with stenosis and 31.92 ± 11.51° (α) and 20.84 ± 9.80° (β) for BAV with insufficiency and minimal at peak ejection 23.42 ± 1.65° (α), 20.38 ± 1.61° (β) for stenosis and 26.62 ± 7.86° (α), 19.79 ± 8.45° (β) for insufficiency. In insufficiency, low shear stress of 0-0.08 Pa and moderate pressure (60-80 mmHg) were present. In BAV with stenosis, low shear stress of 0-0.5 Pa and moderate pressure (0-20 mmHg) were present at diastole; at peak ejection high shear stress >2 Pa and elevated pressure of >80 mmHg were present.

CONCLUSIONS

In a BAV with aortic valve stenosis, the haemodynamics are less favourable. The elevated pressure with elevated shear stress may over the long term promote degenerative processes in the leaflets and consequently valve function failure.

Discrete Subaortic Stenosis: Perspective Roadmap to a Complex Disease

Discrete subaortic stenosis (DSS) is a congenital heart disease that results in the formation of a fibro-membranous tissue, causing an increased pressure gradient in the left ventricular outflow tract (LVOT). While surgical resection of the membrane has shown some success in eliminating the obstruction, it poses significant risks associated with anesthesia, sternotomy, and heart bypass, and it remains associated with a high rate of recurrence. Although a genetic etiology had been initially proposed, the association between DSS and left ventricle (LV) geometrical abnormalities has provided more support to a hemodynamic etiology by which congenital or post-surgical LVOT geometric derangements could generate abnormal shear forces on the septal wall, triggering in turn a fibrotic response. Validating this hypothetical etiology and understanding the mechanobiological processes by which altered shear forces induce fibrosis in the LVOT are major knowledge gaps. This perspective paper describes the current state of knowledge of DSS, articulates the research needs to yield mechanistic insights into a significant pathologic process that is poorly understood, and proposes several strategies aimed at elucidating the potential mechanobiological synergies responsible for DSS pathogenesis. The proposed roadmap has the potential to improve DSS management by identifying early targets for prevention of the fibrotic lesion, and may also prove beneficial in other fibrotic cardiovascular diseases associated with altered flow.

Effect of arteriovenous graft flow rate on vascular access hemodynamics in a novel modular anastomotic valve device

PURPOSE

Perturbed vascular access hemodynamics is considered a potential driver of intimal hyperplasia, the leading cause of vascular access failure. To improve vascular access patency, a modular anastomotic valve device has been designed to normalize venous flow between hemodialysis periods while providing normal vascular access during hemodialysis. The objective of this study was to quantify the effects of arteriovenous graft flow rate on modular anastomotic valve device vascular access hemodynamics under realistic hemodialysis conditions.

METHODS

Modular anastomotic valve device inlet and outlet flow conditions and velocity profiles were measured by ultrasound Doppler in a vascular access flow loop replicating arteriovenous graft flow rates of 800, 1000, and 1500 mL/min. Fluid-structure interaction simulations were performed to identify low wall shear stress regions on the vein wall and to characterize them in terms of temporal shear magnitude, oscillatory shear index, and relative residence time. The model was validated with respect to the Doppler measurements.

RESULTS

The low wall shear stress region generated downstream of the anastomosis under low and moderate arteriovenous graft flow rates was eliminated under the highest arteriovenous graft flow rate. Increase in arteriovenous graft flow rate from 800 to 1500 mL/min resulted in a substantial increase in wall shear stress magnitude (27-fold increase in temporal shear magnitude), the elimination of wall shear stress bidirectionality (0.20-point reduction in oscillatory shear index), and a reduction in flow stagnation (98% decrease in relative residence time). While the results suggest the ability of high arteriovenous graft flow rates to protect the venous wall from intimal hyperplasia-prone hemodynamics, they indicate their adverse impact on the degree of venous hemodynamic abnormality.

Wall Shear Stress Directional Abnormalities in BAV Aortas: Toward a New Hemodynamic Predictor of Aortopathy?

The bicuspid aortic valve (BAV) generates wall shear stress (WSS) abnormalities in the ascending aorta (AA) that may be responsible for the high prevalence of aortopathy in BAV patients. While previous studies have analyzed the magnitude and oscillatory characteristics of the total or streamwise WSS in BAV AAs, the assessment of the circumferential component is lacking despite its expected significance in this highly helical flow environment. This gap may have hampered the identification of a robust hemodynamic predictor of BAV aortopathy. The objective of this study was to perform a global and component-specific assessment of WSS magnitude, oscillatory and directional characteristics in BAV AAs. The WSS environments were computed in the proximal and middle convexity of tricuspid aortic valve (TAV) and BAV AAs using our previous valve-aorta fluid-structure interaction (FSI) models. Component-specific WSS characteristics were investigated in terms of temporal shear magnitude (TSM) and oscillatory shear index (OSI). WSS directionality was quantified in terms of mean WSS vector magnitude and angle, and angular dispersion index (D_α). Local WSS magnitude and multidirectionality were captured in a new shear magnitude and directionality index (SMDI) calculated as the product of the mean WSS magnitude and D_α. BAVs subjected the AA to circumferential TSM overloads (2.4-fold increase vs. TAV). TAV and BAV AAs exhibited a unidirectional circumferential WSS (OSI < 0.04) and an increasingly unidirectional longitudinal WSS between the proximal (OSI > 0.21) and middle (OSI < 0.07) sections. BAVs generated mean WSS vectors skewed toward the anterior wall and WSS angular distributions exhibiting decreased uniformity in the proximal AA (0.27-point increase in D_α vs. TAV). SMDI was elevated in all BAV AAs but peaked in the proximal LR-BAV AA (3.6-fold increase vs. TAV) and in the middle RN-BAV AA (1.6-fold increase vs. TAV). This analysis demonstrates the significance of the circumferential WSS component and the existence of substantial WSS directional abnormalities in BAV AAs. SMDI abnormality distributions in BAV AAs follow the morphotype-dependent occurrence of dilation in BAV AAs, suggesting the predictive potential of this metric for BAV aortopathy.

Computational comparison of regional stress and deformation characteristics in tricuspid and bicuspid aortic valve leaflets

The bicuspid aortic valve (BAV) is the most common congenital valvular defect and a major risk factor for secondary calcific aortic valve disease. While hemodynamics is presumed to be a potential contributor to this complication, the validation of this theory has been hampered by the limited knowledge of the mechanical stress abnormalities experienced by BAV leaflets and their dependence on the heterogeneous BAV fusion patterns. The objective of this study was to compare computationally the regional and temporal fluid wall shear stress (WSS) and structural deformation characteristics in tricuspid aortic valve (TAV), type-0, and type-I BAV leaflets. Arbitrary Lagrangian-Eulerian fluid-structure interaction models were designed to simulate the flow and leaflet dynamics in idealized TAV, type-0, and type-I BAV geometries subjected to physiologic transvalvular pressure. The regional leaflet mechanics was quantified in terms of temporal shear magnitude (TSM), oscillatory shear index (OSI), temporal shear gradient (TSG), and stretch. The simulations identified regions of WSS overloads and increased WSS bidirectionality (174% increase in temporal shear magnitude, 0.10 increase in OSI on type-0 leaflets) in BAV leaflets relative to TAV leaflets. BAV leaflets also experienced larger radial deformations than TAV leaflets (4% increase in type-0 BAV leaflets). Type-I BAV leaflets exhibited contrasted WSS environments marked by WSS overloads on the non-coronary leaflet and sub-physiologic WSS levels on the fused leaflet. This study provides important insights into the mechanical characteristics of BAV leaflets, which may further our understanding of the role played by hemodynamic forces in BAV disease. Copyright © 2016 John Wiley & Sons, Ltd.

Morphotype-Dependent Flow Characteristics in Bicuspid Aortic Valve Ascending Aortas: A Benchtop Particle Image Velocimetry Study

The bicuspid aortic valve (BAV) is a major risk factor for secondary aortopathy such as aortic dilation. The heterogeneous BAV morphotypes [left-right-coronary cusp fusion (LR), right-non-coronary cusp fusion (RN), and left-non-coronary cusp fusion (LN)] are associated with different dilation patterns, suggesting a role for hemodynamics in BAV aortopathogenesis. However, assessment of this theory is still hampered by the limited knowledge of the hemodynamic abnormalities generated by the distinct BAV morphotypes. The objective of this study was to compare experimentally the hemodynamics of a normal (i.e., non-dilated) ascending aorta (AA) subjected to tricuspid aortic valve (TAV), LR-BAV, RN-BAV, and NL-BAV flow. Tissue BAVs reconstructed from porcine TAVs were subjected to physiologic pulsatile flow conditions in a left-heart simulator featuring a realistic aortic root and compliant aorta. Phase-locked particle image velocimetry experiments were carried out to characterize the flow in the aortic root and in the tubular AA in terms of jet skewness and displacement, as well as mean velocity, viscous shear stress and Reynolds shear stress fields. While all three BAVs generated skewed and asymmetrical orifice jets (up to 1.7- and 4.0-fold increase in flow angle and displacement, respectively, relative to the TAV at the sinotubular junction), the RN-BAV jet was out of the plane of observation. The LR- and NL-BAV exhibited a 71% increase in peak-systolic orifice jet velocity relative to the TAV, suggesting an inherent degree of stenosis in BAVs. While these two BAV morphotypes subjected the convexity of the aortic wall to viscous shear stress overloads (1.7-fold increase in maximum peak-systolic viscous shear stress relative to the TAV-AA), the affected sites were morphotype-dependent (LR-BAV: proximal AA, NL-BAV: distal AA). Lastly, the LR- and NL-BAV generated high degrees of turbulence in the AA (up to 2.3-fold increase in peak-systolic Reynolds shear stress relative to the TAV) that were sustained from peak systole throughout the deceleration phase. This in vitro study reveals substantial flow abnormalities (increased jet skewness, asymmetry, jet velocity, turbulence, and shear stress overloads) in non-dilated BAV aortas, which differ from those observed in dilated aortas but still coincide with aortic wall regions prone to dilation.

Simulations of morphotype-dependent hemodynamics in non-dilated bicuspid aortic valve aortas

Bicuspid aortic valves (BAVs) generate flow abnormalities that may promote aortopathy. While positive helix fraction (PHF) index, flow angle (θ), flow displacement (d) and wall shear stress (WSS) exhibit abnormalities in dilated BAV aortas, it is unclear whether those anomalies stem from the abnormal valve anatomy or the dilated aorta. Therefore, the objective of this study was to quantify the early impact of different BAV morphotypes on aorta hemodynamics prior to dilation. Fluid-structure interaction models were designed to quantify standard peak-systolic flow metrics and temporal WSS characteristics in a realistic non-dilated aorta connected to functional tricuspid aortic valve (TAV) and type-I BAVs. While BAVs generated increased helicity (PHF>0.68) in the middle ascending aorta (AA), larger systolic flow skewness (θ>11.2°) and displacement (d>6.8mm) relative to the TAV (PHF=0.51; θ<5.5°; d<3.3mm), no distinct pattern was observed between morphotypes. In contrast, WSS magnitude and directionality abnormalities were BAV morphotype- and site-dependent. Type-I BAVs subjected the AA convexity to peak-systolic WSS overloads (up to 1014% difference vs. TAV). While all BAVs increased WSS unidirectionality on the proximal AA relative to the TAV, the most significant abnormality was achieved by the BAV with left-right-coronary cusp fusion on the wall convexity (up to 0.26 decrease in oscillatory shear index vs. TAV). The results indicate the existence of strong hemodynamic abnormalities in non-dilated type-I BAV AAs, their colocalization with sites vulnerable to dilation and the superior specificity of WSS metrics over global hemodynamic metrics to the valve anatomy.

Aortic valve leaflet wall shear stress characterization revisited: impact of coronary flow

Computational characterizations of aortic valve hemodynamics have typically discarded the effects of coronary flow. The objective of this study was to complement our previous fluid-structure interaction aortic valve model with a physiologic coronary circulation model to quantify the impact of coronary flow on aortic sinus hemodynamics and leaflet wall shear stress (WSS). Coronary flow suppressed vortex development in the two coronary sinuses and altered WSS magnitude and directionality on the three leaflets, with the most substantial differences occurring in the belly and tip regions.

Bicuspid aortic valve hemodynamics does not promote remodeling in porcine aortic wall concavity

AIM: To investigate the role of type-I left-right bicuspid aortic valve (LR-BAV) hemodynamic stresses in the remodeling of the thoracic ascending aorta (AA) concavity, in the absence of underlying genetic or structural defects.

METHODS: Transient wall shear stress (WSS) profiles in the concavity of tricuspid aortic valve (TAV) and LR-BAV AAs were obtained computationally. Tissue specimens excised from the concavity of normal (non-dilated) porcine AAs were subjected for 48 h to those stress environments using a shear stress bioreactor. Tissue remodeling was characterized in terms of matrix metalloproteinase (MMP) expression and activity via immunostaining and gelatin zymography.

RESULTS: Immunostaining semi-quantification results indicated no significant difference in MMP-2 and MMP-9 expression between the tissue groups exposed to TAV and LR-BAV AA WSS (P = 0.80 and P = 0.19, respectively). Zymography densitometry revealed no difference in MMP-2 activity (total activity, active form and latent form) between the groups subjected to TAV AA and LR-BAV AA WSS (P = 0.08, P = 0.15 and P = 0.59, respectively).

CONCLUSION: The hemodynamic stress environment present in the concavity of type-I LR-BAV AA does not cause any significant change in proteolytic enzyme expression and activity as compared to that present in the TAV AA.