Impact of modeling fluid–structure interaction in the computational analysis of aortic root biomechanics

A novel mono-physics particle-based approach for the simulation of cardiovascular fluid-structure interaction problems

BACKGROUND AND OBJECTIVE

Fluid-structure interaction (FSI) is required in the study of several cardiovascular engineering problems were the mutual interaction between the pulsatile blood flow and the tissue structures is essential to establish the biomechanics of the system. Traditional FSI methods are partitioned approaches where two independent solvers, one for the fluid and one for the structure, are asynchronously coupled. This process results into high computational costs. In this work, a new FSI scheme which avoids the coupling of different solvers is presented in the framework of the truly incompressible smoothed particle hydrodynamics (ISPH) method.

METHODS

In the proposed FSI method, ISPH particles contribute to define both the fluid and structural domains and are solved together in a unified system. Solid particles, geometrically defined at the beginning of the simulation, are linked through spring bounds with elastic constant providing the material Young's modulus. At each iteration, internal elastic forces are calculated to restore the springs resting length. These forces are added in the predictor step of the fractional-step procedure used to solve the momentum and continuity equations for incompressible flows of all particles.

RESULTS

The method was validated with a benchmark test case consisting of a flexible beam immersed in a channel. Results showed good agreement with the system coupling approach of a well-established commercial software, ANSYS®, both in terms of fluid-dynamics and beam deformation. The approach was then applied to model a complex cardiovascular problem, consisting in the aortic valve operating function. The valve dynamics during opening and closing phases were compared qualitatively with literature results, demonstrating good consistency.

CONCLUSIONS

The method is computationally more efficient than traditional FSI strategies, and overcomes some of their main drawbacks, such as the impossibility of simulating the correct valve coaptation during the closing phase. Thanks to the incompressibility scheme, the proposed FSI method is appropriate to model biological soft tissues. The simplicity and flexibility of the approach also makes it suitable to be expanded for the modelling of thromboembolic phenomena.

Structural study of a polymeric aortic valve prosthesis. Analysis for a hyperelastic material

This work presents a structural computational simulation of a polymeric aortic valve prosthesis, made with a hyperelastic material (Styrene-Ethylene/Propylene-Styrene). The valve has a suture ring, three pillars placed at 120° and three leaflets. The analysis is based on a modification over previous designs consisting in a fillet concave surface to avoid stress concentration at the junctions between leaflets and pillars. Three shapes were simulated. The first one was used to validate the computational method by comparison of the results with a recent paper. The second shape was designed to show that a fillet or "rounding" can be beneficial to the stress leaflet reduction. The third shape was also designed to show that the reduction of leaflet thickness and intercommissural distance between leaflets at the pillar junctions improves the valve opening and closure. The use of fillet with a 0.5 mm radius, reduced 26.5% the maximum Von Mises stresses for the second shape and 33.9% for the third shape. Additionally, for the latter, the opening area was not affected for the high stiffness due to fillet. The results -mainly for the third shape-are promising and give rise to future studies: further shape optimization, analysis for other materials and valve simulation under pathological loads.

Modeling Dynamics of the Cardiovascular System Using Fluid-Structure Interaction Methods

Using fluid-structure interaction algorithms to simulate the human circulatory system is an innovative approach that can provide valuable insights into cardiovascular dynamics. Fluid-structure interaction algorithms enable us to couple simulations of blood flow and mechanical responses of the blood vessels while taking into account interactions between fluid dynamics and structural behaviors of vessel walls, heart walls, or valves. In the context of the human circulatory system, these algorithms offer a more comprehensive representation by considering the complex interplay between blood flow and the elasticity of blood vessels. Algorithms that simulate fluid flow dynamics and the resulting forces exerted on vessel walls can capture phenomena such as wall deformation, arterial compliance, and the propagation of pressure waves throughout the cardiovascular system. These models enhance the understanding of vasculature properties in human anatomy. The utilization of fluid-structure interaction methods in combination with medical imaging can generate patient-specific models for individual patients to facilitate the process of devising treatment plans. This review evaluates current applications and implications of fluid-structure interaction algorithms with respect to the vasculature, while considering their potential role as a guidance tool for intervention procedures.

Analysis of the Effect of Thickness on the Performance of Polymeric Heart Valves

Polymeric heart valves (PHVs) are a promising and more affordable alternative to mechanical heart valves (MHVs) and bioprosthetic heart valves (BHVs). Materials with good durability and biocompatibility used for PHVs have always been the research focus in the field of prosthetic heart valves for many years, and leaflet thickness is a major design parameter for PHVs. The study aims to discuss the relationship between material properties and valve thickness, provided that the basic functions of PHVs are qualified. The fluid-structure interaction (FSI) approach was employed to obtain a more reliable solution of the effective orifice area (EOA), regurgitant fraction (RF), and stress and strain distribution of the valves with different thicknesses under three materials: Carbothane PC-3585A, xSIBS and SIBS-CNTs. This study demonstrates that the smaller elastic modulus of Carbothane PC-3585A allowed for a thicker valve (>0.3 mm) to be produced, while for materials with an elastic modulus higher than that of xSIBS (2.8 MPa), a thickness less than 0.2 mm would be a good attempt to meet the RF standard. What is more, when the elastic modulus is higher than 23.9 MPa, the thickness of the PHV is recommended to be 0.l-0.15 mm. Reducing the RF is one of the directions of PHV optimization in the future. Reducing the thickness and improving other design parameters are reliable means to reduce the RF for materials with high and low elastic modulus, respectively.

A Doppler-exclusive non-invasive computational diagnostic framework for personalized transcatheter aortic valve replacement

Given the associated risks with transcatheter aortic valve replacement (TAVR), it is crucial to determine how the implant will affect the valve dynamics and cardiac function, and if TAVR will improve or worsen the outcome of the patient. Effective treatment strategies, indeed, rely heavily on the complete understanding of the valve dynamics. We developed an innovative Doppler-exclusive non-invasive computational framework that can function as a diagnostic tool to assess valve dynamics in patients with aortic stenosis in both pre- and post-TAVR status. Clinical Doppler pressure was reduced by TAVR (52.2 ± 20.4 vs. 17.3 ± 13.8 [mmHg], p < 0.001), but it was not always accompanied by improvements in valve dynamics and left ventricle (LV) hemodynamics metrics. TAVR had no effect on LV workload in 4 patients, and LV workload post-TAVR significantly rose in 4 other patients. Despite the group level improvements in maximum LV pressure (166.4 ± 32.2 vs 131.4 ± 16.9 [mmHg], p < 0.05), only 5 of the 12 patients (41%) had a decrease in LV pressure. Moreover, TAVR did not always improve valve dynamics. TAVR did not necessarily result in a decrease (in 9 out of 12 patients investigated in this study) in major principal stress on the aortic valve leaflets which is one of the main contributors in valve degeneration and, consequently, failure of heart valves. Diastolic stresses increased significantly post-TAVR (34%, 109% and 81%, p < 0.001) for each left, right and non-coronary leaflets respectively. Moreover, we quantified the stiffness and material properties of aortic valve leaflets which correspond with the reduced calcified region average stiffness among leaflets (66%, 74% and 62%; p < 0.001; N = 12). Valve dynamics post-intervention should be quantified and monitored to ensure the improvement of patient conditions and prevent any further complications. Improper evaluation of biomechanical valve features pre-intervention as well as post-intervention may result in harmful effects post-TAVR in patients including paravalvular leaks, valve degeneration, failure of TAVR and heart failure.

An ultrasound-exclusive non-invasive computational diagnostic framework for personalized cardiology of aortic valve stenosis

Aortic stenosis (AS) is an acute and chronic cardiovascular disease and If left untreated, 50% of these patients will die within two years of developing symptoms. AS is characterized as the stiffening of the aortic valve leaflets which restricts their motion and prevents the proper opening under transvalvular pressure. Assessments of the valve dynamics, if available, would provide valuable information about the patient's state of cardiac deterioration as well as heart recovery and can have incredible impacts on patient care, planning interventions and making critical clinical decisions with life-threatening risks. Despite remarkable advancements in medical imaging, there are no clinical tools available to quantify valve dynamics invasively or noninvasively. In this study, we developed a highly innovative ultrasound-based non-invasive computational framework that can function as a diagnostic tool to assess valve dynamics (e.g. transient 3-D distribution of stress and displacement, 3-D deformed shape of leaflets, geometric orifice area and angular positions of leaflets) for patients with AS at no risk to the patients. Such a diagnostic tool considers the local valve dynamics and the global circulatory system to provide a platform for testing the intervention scenarios and evaluating their effects. We used clinical data of 12 patients with AS not only to validate the proposed framework but also to demonstrate its diagnostic abilities by providing novel analyses and interpretations of clinical data in both pre and post intervention states. We used transthoracic echocardiogram (TTE) data for the developments and transesophageal echocardiography (TEE) data for validation.

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.

Importance of Non-Newtonian Computational Fluid Modeling on Severely Calcified Aortic Valve Geometries - Insights from Quasi-Steady State Simulations

The Newtonian model has commonly been used to represent the viscosity of blood in the aorta, despite blood itself being a non-Newtonian fluid. This is justified where shear rates tend to be large. However, we hypothesized that using the Newtonian model to predict the hemodynamics on the aortic valve, particularly in those with severe calcifications, is inaccurate owing to valve leaflet geometry irregularities inducing multiple regions of low shear rates, &lt;100 s-1, where a Newtonian model is invalid. We investigated the use of three fluid models: Newtonian, Carreau, and Quemada on a healthy and calcified valve and compared their ability to capture important hemodynamic parameters of time averaged wall shear stress (TAWSS) and the oscillatory shear index (OSI). Our findings indicate that when the shear rates were large enough, &gt;100 s-1, use of a Newtonian model was justified for the ventricularis side of the leaflet. However, in spatial regions of low shear rates, &lt;100 s-1, specifically on the inner cusps of the fibrosa side of the valve, TAWSS calculations under a Newtonian model were found to be much lower when compared with their non-Newtonian Carreau and Quemada counterparts. We hereby conclude that in order to facilitate more accurate computational flow simulations in severe aortic valve calcification which is subjected to relatively large spatial regions of low shear (&lt;100 s-1), a non-Newtonian model should be applied.

Outcome of transcatheter aortic valve replacement in bicuspid aortic valve stenosis with new-generation devices

OBJECTIVES

To compare device success and paravalvular leak rates of 3 new-generation transcatheter aortic valve replacement devices in patients with bicuspid aortic valve stenosis and to test their biomechanical performance in a computer-based simulation model of aortic root with increasing ellipticity.

METHODS

This retrospective multicentre study included 56 bicuspid aortic valve patients undergoing transcatheter aortic valve replacement with new-generation devices: Lotus/Lotus Edge (N = 15; 27%), Evolut-R (N = 20; 36%) and ACURATE neo (N = 21; 37%). Three virtual simulation models of aortic root with increasing index of eccentricity (0-0.25-0.5) were implemented. Stress distribution, stent-root contact area and paravalvular orifice area were computed.

RESULTS

Device success was achieved in 43/56 patients (77%) with comparable rates among Lotus (87%), Evolut-R (60%) and ACURATE neo (86%; P = 0.085). Moderate paravalvular leak rate was significantly lower in the Lotus group as compared to Evolut-R group (0% vs 30%; P = 0.027) and comparable to the ACURATE neo group (0% vs 10%; P = 0.33). By index of eccentricity = 0.5, Lotus showed a uniform and symmetric pattern of stress distribution with absent paravalvular orifice area, ACURATE neo showed a mild asymmetry with small paravalvular orifice area (1.1 mm2), whereas a severely asymmetric pattern was evident with Evolut-R, resulting in a large paravalvular orifice area (12.0 mm2).

CONCLUSIONS

Transcatheter aortic valve replacement in bicuspid aortic valve patients with new-generation devices showed comparable device success rates. Lotus showed moderate paravalvular leak rate comparable to that of ACURATE neo and significantly lower than Evolut-R. On simulation, Lotus and ACURATE neo showed optimal adaptability to elliptic anatomies as compared to Evolut-R.

Isogeometric finite element-based simulation of the aortic heart valve: Integration of neural network structural material model and structural tensor fiber architecture representations

The functional complexity of native and replacement aortic heart valves (AVs) is well known, incorporating such physical phenomenons as time-varying non-linear anisotropic soft tissue mechanical behavior, geometric non-linearity, complex multi-surface time varying contact, and fluid-structure interactions to name a few. It is thus clear that computational simulations are critical in understanding AV function and for the rational basis for design of their replacements. However, such approaches continued to be limited by ad-hoc approaches for incorporating tissue fibrous structure, high-fidelity material models, and valve geometry. To this end, we developed an integrated tri-leaflet valve pipeline built upon an isogeometric analysis framework. A high-order structural tensor (HOST)-based method was developed for efficient storage and mapping the two-dimensional fiber structural data onto the valvular 3D geometry. We then developed a neural network (NN) material model that learned the responses of a detailed meso-structural model for exogenously cross-linked planar soft tissues. The NN material model not only reproduced the full anisotropic mechanical responses but also demonstrated a considerable efficiency improvement, as it was trained over a range of realizable fibrous structures. Results of parametric simulations were then performed, as well as population-based bicuspid AV fiber structure, that demonstrated the efficiency and robustness of the present approach. In summary, the present approach that integrates HOST and NN material model provides an efficient computational analysis framework with increased physical and functional realism for the simulation of native and replacement tri-leaflet heart valves.

Predicting TAVI paravalvular regurgitation outcomes based on numerical simulation of the aortic annulus eccentricity and perivalvular areas

Trans-catheter aortic valve implantation (TAVI) is an increasingly adopted technique which provides a minimal invasive solution for patients who suffer from severe aortic stenosis. Some complications of the procedure could be annular rupture or paravalvular leakage, both related with adverse outcome. In TAVI with balloon expandable devices, a mismatch between those two factors leads to a conflict situation, where improving one worsens the other. The presented research proposes a methodology that uses numerical simulation to obtain certain TAVI outcomes related with aortic regurgitation due to paravalvular leakage, such as perivalvular area, aortic eccentricity or annular pressure. The application of the methodology for two patients shows the possibility of predicting those quantities. The highest stress values are distributed along the contact area. Results also show that a great deformation on the aortic annulus does not necessarily imply a higher stress; pressure can either be converted into root reshape or into root stretching. Validation of the results was done using scientific publications, clinical guidelines and clinical reports. Numerical simulation provides a suitable tool that could possibly contribute to optimize the planification procedure adjusting the mismatch between size and pressure.

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.

Mechanics of the dicrotic notch: An acceleration hypothesis

The dicrotic notch is a prominent and distinctive feature of the pressure waveform in the central arteries. It is universally used to demarcate the end of systole and the beginning of diastole in these arteries. Despite its importance clinically, no physical mechanism for the formation of the dicrotic notch has been demonstrated convincingly. We first explore a mechanism based on the reflection of a backward wavefront from the aortic valve at the time of closure. This hypothesis is rejected on the basis of experimental evidence from measurements made in dogs. A new hypothesis is presented involving the acceleration of the aortic valve apparatus at the time of valve closure. This hypothesis is supported by new calculations of the acceleration of the aortic valve apparatus during the cardiac cycle based on computed tomography scans in man.

Transient Study of Flow and Cavitation Inside a Bileaflet Mechanical Heart Valve

A mechanical heart valve (MHV) is an effective device to cure heart disease, which has the advantage of long life and high reliability. Due to the hemodynamic characteristics of blood, mechanical heart valves can lead to potential complications such as hemolysis, which have damage to the blood elements and thrombosis. In this paper, flowing features of the blood in the valve are analyzed and the cavitation mechanism in bileaflet mechanical heart valve (BMHV) is studied. Results show that the water hammer effect and the high-speed leakage flow effect are the primary causes of the cavitation in the valve. Compared with the high-speed leakage flow effect, the water hammer has a greater effect on the cavitation strength. The valve goes through four kinds of working condition within one heart beating period, including, fully opening stage, closing stage and fully closing stage. These four stages, respectively, make up 8.5%, 16.1%, 4.7% and 70.7% of the total period. The cavitation occurs on the fully closing stage. When the valve is in closing stage, the high pressure downstream of the valve lasts for about 20 ms and the high-speed leakage flow lasts for about 200 ms. This study systematically analyzes the causes of cavitation emerged in the process of periodic motion, which proposes the method for characterizing the intensity of the cavitation, and can be referred to for the cavitation suppression of the BHMV and similar valves.

A biomechanical model of the pathological aortic valve: simulation of aortic stenosis

Aortic stenosis (AS) disease is a narrowing of the aortic valve (AV) opening which reduces blood flow from the heart causing several health complications. Although a lot of work has been done in AV simulations, most of the efforts have been conducted regarding healthy valves. In this article, a new three-dimensional patient-specific biomechanical model of the valve, based on a parametric formulation of the stenosis that permits the simulation of different degrees of pathology, is presented. The formulation is based on a double approach: the first one is done from the geometric point of view, reducing the effective ejection area of the AV by joining leaflets using a zipper effect to sew them; the second one, in terms of functionality, is based on the modification of AV tissue properties due to the effect of calcifications. Both healthy and stenotic valves were created using patient-specific data and results of the numerical simulation of the valve function are provided. Analysis of the results shows a variation in the first principal stress, geometric orifice area, and blood velocity which were validated against clinical data. Thus, the possibility to create a pipeline which allows the integration of patient-specific data from echocardiographic images and iFR studies to perform finite elements analysis is proved.

Flow-Structure Interaction Simulations of the Aortic Heart Valve at Physiologic Conditions: The Role of Tissue Constitutive Model

The blood flow patterns in the region around the aortic valve depend on the geometry of the aorta and on the complex flow-structure interaction between the pulsatile flow and the valve leaflets. Consequently, the flow depends strongly on the constitutive properties of the tissue, which can be expected to vary between healthy and diseased heart valves or native and prosthetic valves. The main goal of this work is to qualitatively demonstrate that the choice of the constitutive model of the aortic valve is critical in analysis of heart hemodynamics. To accomplish that two different constitutive models were used in curvilinear immersed boundary-finite element-fluid-structure interaction (CURVIB-FE-FSI) method developed by Gilmanov et al. (2015, "A Numerical Approach for Simulating Fluid Structure Interaction of Flexible Thin Shells Undergoing Arbitrarily Large Deformations in Complex Domains," J. Comput. Phys., 300, pp. 814-843.) to simulate an aortic valve in an anatomic aorta at physiologic conditions. The two constitutive models are: (1) the Saint-Venant (StV) model and (2) the modified May-Newman&Yin (MNY) model. The MNY model is more general and includes nonlinear, anisotropic effects. It is appropriate to model the behavior of both prosthetic and biological tissue including native valves. Both models are employed to carry out FSI simulations of the same valve in the same aorta anatomy. The computed results reveal dramatic differences in both the vorticity dynamics in the aortic sinus and the wall shear-stress patterns on the aortic valve leaflets and underscore the importance of tissue constitutive models for clinically relevant simulations of aortic valves.

Comparison of a fixed-grid and arbitrary Lagrangian–Eulerian methods on modelling fluid–structure interaction of the aortic valve

This study was aimed at assessing the robustness of a fixed-grid fluid–structure interaction method (Multi-Material Arbitrary Lagrangian–Eulerian) to modelling the two-dimensional native aortic valve dynamics and comparing it to the Arbitrary Lagrangian–Eulerian method. For the fixed-grid method, the explicit finite element solver LS-DYNA was utilized, where two independent meshes for the fluid and structure were generated and the penalty method was used to handle the coupling between the fluid and structure domains. For the Arbitrary Lagrangian–Eulerian method, the implicit finite element solver ADINA was used where two separate conforming meshes were used for the valve structure and the fluid domains. The comparison demonstrated that both fluid–structure interaction methods predicted accurately the valve dynamics, fluid flow, and stress distribution, implying that fixed-grid methods can be used in situations where the Arbitrary Lagrangian–Eulerian method fails.

Proposed percutaneous aortic valve prosthesis made of cryogel

Transcatheter heart valves are promising for high-risk patients. Generally, their leaflets are made of pericardium stented in a Nitinol basket. Despite their relative success, they are associated with significant complications such as valve migration, implantation risks, stroke, coronary obstruction, myocardial infraction, acute kidney injury (which all are due to the release of detached solid calcific pieces in to the blood stream) and expected issues existing with tissue valves such as leaflet calcification. This study is an attempt to fabricate the first ever polymeric percutaneous valves made of cryogel following the geometry and mechanical properties of porcine aortic valve to address some of the above-mentioned shortcomings. A novel, one-piece, tricuspid percutaneous valve, consisting of leaflets made entirely from the hydrogel, polyvinyl alcohol cryogel reinforced by bacterial cellulose natural nanocomposite, attached to a Nitinol basket was developed and demonstrated. Following the natural geometry of the valve, a novel approach was applied based on the revolution about an axis of a hyperboloid shape. The geometry was modified based on avoiding sharp warpage of leaflets and removal of the central opening orifice area of the valve when valve is fully closed using the finite element analysis. The modified geometry was replaced by a cloud of (control) points and was essentially converted to Bezier surfaces for further adjustment. A cavity mold was then designed and fabricated to form the valve. The fabricated valve was sewn into the Nitinol basket which is covered by Dacron cloth. The models presented in this study merit further development and revisions for both aortic and mitral positions.

Validation and Extension of a Fluid-Structure Interaction Model of the Healthy Aortic Valve

Purpose

The understanding of the optimum function of the healthy aortic valve is essential in interpreting the effect of pathologies in the region, and in devising effective treatments to restore the physiological functions. Still, there is no consensus on the operating mechanism that regulates the valve opening and closing dynamics. The aim of this study is to develop a numerical model that can support a better comprehension of the valve function and serve as a reference to identify the changes produced by specific pathologies and treatments.

Methods

A numerical model was developed and adapted to accurately replicate the conditions of a previous in vitro investigation into aortic valve dynamics, performed by means of particle image velocimetry (PIV). The resulting velocity fields of the two analyses were qualitatively and quantitatively compared to validate the numerical model. In order to simulate more physiological operating conditions, this was then modified to overcome the main limitations of the experimental setup, such as the presence of a supporting stent and the non-physiological properties of the fluid and vessels.

Results

The velocity fields of the initial model resulted in good agreement with those obtained from the PIV, with similar flow structures and about 90% of the computed velocities after valve opening within the standard deviation of the equivalent velocity measurements of the in vitro model. Once the experimental limitations were removed from the model, the valve opening dynamics changed substantially, with the leaflets opening into the sinuses to a much greater extent, enlarging the effective orifice area by 11%, and reducing greatly the vortical structures previously observed in proximity of the Valsalva sinuses wall.

Conclusions

The study suggests a new operating mechanism for the healthy aortic valve leaflets considerably different from what reported in the literature to date and largely more efficient in terms of hydrodynamic performance. This work also confirms the crucial role that numerical approaches, complemented with experimental findings, can play in overcoming some of the limitations inherent in experimental techniques, supporting the full understanding of complex physiological phenomena.

Electronic supplementary material

The online version of this article (doi:10.1007/s13239-018-00391-1) contains supplementary material, which is available to authorized users.

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.

Aortic expansion induces lumen narrrowing in anomalous coronary arteries: a parametric structural finite element analysis

Anomalous aortic origin of coronary arteries (AAOCA) is a congenital disease that can lead to cardiac ischemia during intense physical activity. Although AAOCA is responsible for sudden cardiac death (SCD) among young athletes and soldiers, the mechanisms underlying the coronary occlusion during physical effort still have to be clarified. The present study investigates the correlation between geometric features of the anomaly and coronary lumen narrowing under aortic root dilatations. Idealized parametric computer-aided designed (CAD) models of the aortic root with anomalous and normal coronary are created and static finite element (FE) simulations of increasing aortic root expansions are carried out. Different coronary take-off angles and intramural penetrations are investigated to assess their role on coronary lumen narrowing. Results show that increasing aortic and coronary pressures lead to lumen expansions in normal coronaries, particularly in the proximal tract, while the expansion of anomalous coronary is impaired especially at the ostium. Concerning the geometric features of the anomaly, acute take-off angles cause elongated coronary ostia, with an eccentricity increasing with aortic expansion; the impact of intramural penetration of coronary on its luminal narrowing is limited. The present study provides a proof of concept of the biomechanical reasons underlying the lumen narrowing in AAOCA during aortic expansion, promoting the role of computational simulations as a tool to assess the mechanisms of this pathology.

Fluid-structure interaction simulation of artificial textile reinforced aortic heart valve: Validation with an in-vitro test

Prosthetic heart valves deployed in the left heart (aortic and mitral) are subjected to harsh hemodynamical conditions. Most of the tissue engineered heart valves have been developed for the low pressure pulmonary position because of the difficulties in fabricating a mechanically strong valve, able to withstand the systemic circulation. This necessitates the use of reinforcing scaffolds, resulting in a tissue-engineered textile reinforced tubular aortic heart valve. Therefore, to better design these implants, material behaviour of the composite, valve kinematics and its hemodynamical response need to be evaluated. Experimental assessment can be immensely time consuming and expensive, paving way for numerical studies. In this work, the material properties obtained using the previously proposed multi-scale numerical method for textile composites was evaluated for its accuracy. An in silico immersed boundary (IB) fluid structure interaction (FSI) simulation emulating the in vitro experiment was set-up to evaluate and compare the geometric orifice area and flow rate for one beat cycle. Results from the in silico FSI simulation were found to be in good coherence with the in vitro test during the systolic phase, while mean deviation of approximately 9% was observed during the diastolic phase of a beat cycle. Merits and demerits of the in silico IB-FSI method for the presented case study has been discussed with the advantages outweighing the drawbacks, indicating the potential towards an effective use of this framework in the development and analysis of heart valves.

The Effects of Left Ventricular Assist Device Support Level on the Biomechanical States of Aortic Valve

Background

Although aortic valve disease caused by left ventricular assist device (LVAD) support has attracted more and more attention, the precise biomechanical effects of LVAD support level on the aortic valve are still unclear.

Material/Methods

A structural finite element models study was conducted using an ideal aortic valve geometric model. Four different study conditions were designed, according to the reduction of the open duration of the aortic valve. The isotropic hyperelastic constitutive equation was chosen to reflect the mechanical property of the leaflets. The distribution of the stress, strain, and transient dynamics of the leaflet were calculated.

Results

Along with the increase of LVAD support level, the open duration of the aortic valve was also reduced by the increase of LVAD support (low support level case 0.23 seconds versus middle support level case 0.2 seconds versus high support level case 0.14 seconds). Moreover, along with the increase of support mode of LVAD, the von Mises stress in most leaflet areas was increased from the low stress level (0–0.4 MPa) to the middle region (0.4–0.8 MPa). Once the leaflets were continuously closed, the high stress level (larger than 0.8 MPa) was observed. In contrast, the support level of LVAD only had slight effects on the distribution of von Mises strain. According to the aforementioned results, maintaining the open duration of aortic valve longer than 0.2 seconds could achieve better performance of biomechanical states of leaflets.

Conclusions

This study could provide useful information on the determination of optimal LVAD support strategy.

Computational methods for the aortic heart valve and its replacements

INTRODUCTION

Replacement with a prosthetic device remains a major treatment option for the patients suffering from heart valve disease, with prevalence growing resulting from an ageing population. While the most popular replacement heart valve continues to be the bioprosthetic heart valve (BHV), its durability remains limited. There is thus a continued need to develop a general understanding of the underlying mechanisms limiting BHV durability to facilitate development of a more durable prosthesis. In this regard, computational models can play a pivotal role as they can evaluate our understanding of the underlying mechanisms and be used to optimize designs that may not always be intuitive. Areas covered: This review covers recent progress in computational models for the simulation of BHV, with a focus on aortic valve (AV) replacement. Recent contributions in valve geometry, leaflet material models, novel methods for numerical simulation, and applications to BHV optimization are discussed. This information should serve not only to infer reliable and dependable BHV function, but also to establish guidelines and insight for the design of future prosthetic valves by analyzing the influence of design, hemodynamics and tissue mechanics. Expert commentary: The paradigm of predictive modeling of heart valve prosthesis are becoming a reality which can simultaneously improve clinical outcomes and reduce costs. It can also lead to patient-specific valve design.

Numerical Investigation of Pulse Wave Propagation in Arteries Using Fluid Structure Interaction Capabilities

The aim of this study is to present a reliable computational scheme to serve in pulse wave velocity (PWV) assessment in large arteries. Clinicians considered it as an indication of human blood vessels' stiffness. The simulation of PWV was conducted using a 3D elastic tube representing an artery. The constitutive material model specific for vascular applications was applied to the tube material. The fluid was defined with an equation of state representing the blood material. The onset of a velocity pulse was applied at the tube inlet to produce wave propagation. The Coupled Eulerian-Lagrangian (CEL) modeling technique with fluid structure interaction (FSI) was implemented. The scaling of sound speed and its effect on results and computing time is discussed and concluded that a value of 60 m/s was suitable for simulating vascular biomechanical problems. Two methods were used: foot-to-foot measurement of velocity waveforms and slope of the regression line of the wall radial deflection wave peaks throughout a contour plot. Both methods showed coincident results. Results were approximately 6% less than those calculated from the Moens-Korteweg equation. The proposed method was able to describe the increase in the stiffness of the walls of large human arteries via the PWV estimates.

Fully-coupled fluid-structure interaction simulation of the aortic and mitral valves in a realistic 3D left ventricle model

In this study, we present a fully-coupled fluid-structure interaction (FSI) framework that combines smoothed particle hydrodynamics (SPH) and nonlinear finite element (FE) method to investigate the coupled aortic and mitral valves structural response and the bulk intraventricular hemodynamics in a realistic left ventricle (LV) model during the entire cardiac cycle. The FSI model incorporates valve structures that consider native asymmetric leaflet geometries, anisotropic hyperelastic material models and human material properties. Comparison of FSI results with subject-specific echocardiography data demonstrates that the SPH-FE approach is able to quantitatively predict the opening and closing times of the valves, the mitral leaflet opening and closing angles, and the large-scale intraventricular flow phenomena with a reasonable agreement. Moreover, comparison of FSI results with a LV model without valves reveals substantial differences in the flow field. Peak systolic velocities obtained from the FSI model and the LV model without valves are 2.56 m/s and 1.16 m/s, respectively, compared to the Doppler echo data of 2.17 m/s. The proposed SPH-FE FSI framework represents a further step towards modeling patient-specific coupled LV-valve dynamics, and has the potential to improve our understanding of cardiovascular physiology and to support professionals in clinical decision-making.

A novel approach to the quantification of aortic root in vivo structural mechanics

Understanding aortic root in vivo biomechanics can help in elucidating key mechanisms involved in aortic root pathologies and in the outcome of their surgical treatment. Numerical models can provide useful quantitative information. For this to be reliable, detailed aortic root anatomy should be captured. Also, since the aortic root is never unloaded throughout the cardiac cycle, the modeled geometry should be consistent with the in vivo loads acting on it. Achieving such consistency is still a challenge, which was tackled only by few numerical studies. Here we propose and describe in detail a new approach to the finite element modeling of aortic root in vivo structural mechanics. Our approach exploits the anatomical information yielded by magnetic resonance imaging by reconstructing the 3-dimensional end-diastolic geometry of the aortic root and makes the reconstructed geometry consistent with end-diastolic loading conditions through the estimation of the corresponding prestresses field. We implemented our approach through a semiautomated modeling pipeline, and we applied it to quantify aortic root biomechanics in 4 healthy participants. Computed results highlighted that including prestresses into the model allowed for pressurizing the aortic root to the end-diastolic pressure while matching the image-based ground truth data. Aortic root dynamics, tissues strains, and stresses computed at relevant time points through the cardiac cycle were consistent with a broad set of data from previous computational and in vivo studies, strongly suggesting the potential of the method. Also, results highlighted the major role played by the anatomy in driving aortic root biomechanics.

Experimental Study of Anisotropic Stress/Strain Relationships of Aortic and Pulmonary Artery Homografts and Synthetic Vascular Grafts

Homografts and synthetic grafts are used in surgery for congenital heart disease (CHD). Determining these materials' mechanical properties will aid in understanding tissue behavior when subjected to abnormal CHD hemodynamics. Homograft tissue samples from anterior/posterior aspects, of ascending/descending aorta (AA, DA), innominate artery (IA), left subclavian artery (LScA), left common carotid artery (LCCA), main/left/right pulmonary artery (MPA, LPA, RPA), and synthetic vascular grafts, were obtained in three orientations: circumferential, diagonal (45 deg relative to circumferential direction), and longitudinal. Samples were subjected to uniaxial tensile testing (UTT). True strain-Cauchy stress curves were individually fitted for each orientation to calibrate Fung model. Then, they were used to calibrate anisotropic Holzapfel–Gasser model (R2 > 0.95). Most samples demonstrated a nonlinear hyperelastic strain–stress response to UTT. Stiffness (measured by tangent modulus at different strains) in all orientations were compared and shown as contour plots. For each vessel segment at all strain levels, stiffness was not significantly different among aspects and orientations. For synthetic grafts, stiffness was significantly different among orientations (p < 0.042). Aorta is significantly stiffer than pulmonary artery at 10% strain, comparing all orientations, aspects, and regions (p = 0.0001). Synthetic grafts are significantly stiffer than aortic and pulmonary homografts at all strain levels (p < 0.046). Aortic, pulmonary artery, and synthetic grafts exhibit hyperelastic biomechanical behavior with anisotropic effect. Differences in mechanical properties among vascular grafts may affect native tissue behavior and ventricular/arterial mechanical coupling, and increase the risk of deformation due to abnormal CHD hemodynamics.

Toward the Virtual Benchmarking of Pneumatic Ventricular Assist Devices: Application of a Novel Fluid-Structure Interaction-Based Strategy to the Penn State 12 cc Device

The pediatric use of pneumatic ventricular assist devices (VADs) as a bridge to heart transplant still suffers for short-term major complications such as bleeding and thromboembolism. Although numerical techniques are increasingly exploited to support the process of device optimization, an effective virtual benchmark is still lacking. Focusing on the 12 cc Penn State pneumatic VAD, we developed a novel fluid-structure interaction (FSI) model able to capture the device functioning, reproducing the mechanical interplay between the diaphragm, the blood chamber, and the pneumatic actuation. The FSI model included the diaphragm mechanical response from uniaxial tensile tests, realistic VAD pressure operative conditions from a dedicated mock loop system, and the behavior of VAD valves. Our FSI-based benchmark effectively captured the complexity of the diaphragm dynamics. During diastole, the initial slow diaphragm retraction in the air chamber was followed by a more rapid phase; asymmetries were noticed in the diaphragm configuration during its systolic inflation in the blood chamber. The FSI model also captured the major features of the device fluid dynamics. In particular, during diastole, a rotational wall washing pattern is promoted by the penetrating inlet jet with a low-velocity region located in the center of the device. Our numerical analysis of the 12 cc Penn State VAD points out the potential of the proposed FSI approach well resembling previous experimental evidences; if further tested and validated, it could be exploited as a virtual benchmark to deepen VAD-related complications and to support the ongoing optimization of pediatric devices.

A patient-specific aortic valve model based on moving resistive immersed implicit surfaces

In this paper, we propose a full computational framework to simulate the hemodynamics in the aorta including the valve. Closed and open valve surfaces, as well as the lumen aorta, are reconstructed directly from medical images using new ad hoc algorithms, allowing a patient-specific simulation. The fluid dynamics problem that accounts from the movement of the valve is solved by a new 3D-0D fluid-structure interaction model in which the valve surface is implicitly represented through level set functions, yielding, in the Navier-Stokes equations, a resistive penalization term enforcing the blood to adhere to the valve leaflets. The dynamics of the valve between its closed and open position is modeled using a reduced geometric 0D model. At the discrete level, a finite element formulation is used and the SUPG stabilization is extended to include the resistive term in the Navier-Stokes equations. Then, after time discretization, the 3D fluid and 0D valve models are coupled through a staggered approach. This computational framework, applied to a patient-specific geometry and data, allows to simulate the movement of the valve, the sharp pressure jump occurring across the leaflets, and the blood flow pattern inside the aorta.

Evaluation of an aortic valve prosthesis: Fluid-structure interaction or structural simulation?

Bio-inspired polymeric heart valves (PHVs) are excellent candidates to mimic the structural and the fluid dynamic features of the native valve. PHVs can be implanted as prosthetic alternative to currently clinically used mechanical and biological valves or as potential candidate for a minimally invasive treatment, like the transcatheter aortic valve implantation. Nevertheless, PHVs are not currently used for clinical applications due to their lack of reliability. In order to investigate the main features of this new class of prostheses, pulsatile tests in an in-house pulse duplicator were carried out and reproduced in silico with both structural Finite-Element (FE) and Fluid-Structure interaction (FSI) analyses. Valve kinematics and geometric orifice area (GOA) were evaluated to compare the in vitro and the in silico tests. Numerical results showed better similarity with experiments for the FSI than for the FE simulations. The maximum difference between experimental and FSI GOA at maximum opening time was only 5%, as compared to the 46.5% between experimental and structural FE GOA. The stress distribution on the valve leaflets clearly reflected the difference in valve kinematics. Higher stress values were found in the FSI simulations with respect to those obtained in the FE simulation. This study demonstrates that FSI simulations are more appropriate than FE simulations to describe the actual behaviour of PHVs as they can replicate the valve-fluid interaction while providing realistic fluid dynamic results.

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.

Immersogeometric cardiovascular fluid-structure interaction analysis with divergence-conforming B-splines

This paper uses a divergence-conforming B-spline fluid discretization to address the long-standing issue of poor mass conservation in immersed methods for computational fluid-structure interaction (FSI) that represent the influence of the structure as a forcing term in the fluid subproblem. We focus, in particular, on the immersogeometric method developed in our earlier work, analyze its convergence for linear model problems, then apply it to FSI analysis of heart valves, using divergence-conforming B-splines to discretize the fluid subproblem. Poor mass conservation can manifest as effective leakage of fluid through thin solid barriers. This leakage disrupts the qualitative behavior of FSI systems such as heart valves, which exist specifically to block flow. Divergence-conforming discretizations can enforce mass conservation exactly, avoiding this problem. To demonstrate the practical utility of immersogeometric FSI analysis with divergence-conforming B-splines, we use the methods described in this paper to construct and evaluate a computational model of an in vitro experiment that pumps water through an artificial valve.

Numerical simulation of closure performance for neo-aortic valve for arterial switch operation

Background

Modeling neo-aortic valve for arterial switch surgical planning to simulate the neo-aortic valve closure performance.

Methods

We created five geometrical models of neo-aortic valve, namely model A, model B, model C, model D and model E with different size of sinotubular junction or sinus. The nodes at the ends of aorta and left ventricle duct fixed all the degrees of freedom. Transvalvular pressure of normal diastolic blood pressure of 54 mmHg was applied on the neo-aortic valve cusps. The neo-aortic valve closure performance was investigated by the parameters, such as stress of neo-aortic root, variation of neo-aortic valve ring as well as aortic valve cusps contact force in the cardiac diastole.

Results

The maximum stress of the five neo-aortic valves were 96.29, 98.34, 96.28, 98.26, and 90.60 kPa, respectively. Compared among five neo-aortic valve, aortic valve cusps contact forces were changed by 43.33, −10.00% enlarging or narrowing the sinotubular junction by 20% respectively based on the reference model A. The cusps contact forces were changed by 6.67, −23.33% with sinus diameter varying 1.2 times and 0.8 times respectively.

Conclusions

Comparing with stress of healthy adult subjects, the neo-aortic valve of infant creates lower stress. It is evident that enlarging or narrowing the sinotubular junction within a range of 20% can increase or decrease the maximum stress and aortic valve cusps contact force of neo-aortic valve.

Fluid-Structure Interaction Study of Transcatheter Aortic Valve Dynamics Using Smoothed Particle Hydrodynamics

Computational modeling of heart valve dynamics incorporating both fluid dynamics and valve structural responses has been challenging. In this study, we developed a novel fully-coupled fluid-structure interaction (FSI) model using smoothed particle hydrodynamics (SPH). A previously developed nonlinear finite element (FE) model of transcatheter aortic valves (TAV) was utilized to couple with SPH to simulate valve leaflet dynamics throughout the entire cardiac cycle. Comparative simulations were performed to investigate the impact of using FE-only models vs. FSI models, as well as an isotropic vs. an anisotropic leaflet material model in TAV simulations. From the results, substantial differences in leaflet kinematics between FE-only and FSI models were observed, and the FSI model could capture the realistic leaflet dynamic deformation due to its more accurate spatial and temporal loading conditions imposed on the leaflets. The stress and the strain distributions were similar between the FE and FSI simulations. However, the peak stresses were different due to the water hammer effect induced by the fluid inertia in the FSI model during the closing phase, which led to 13-28% lower peak stresses in the FE-only model compared to that of the FSI model. The simulation results also indicated that tissue anisotropy had a minor impact on hemodynamics of the valve. However, a lower tissue stiffness in the radial direction of the leaflets could reduce the leaflet peak stress caused by the water hammer effect. It is hoped that the developed FSI models can serve as an effective tool to better assess valve dynamics and optimize next generation TAV designs.

Biomechanical Behavior of Bioprosthetic Heart Valve Heterograft Tissues: Characterization, Simulation, and Performance

The use of replacement heart valves continues to grow due to the increased prevalence of valvular heart disease resulting from an ageing population. Since bioprosthetic heart valves (BHVs) continue to be the preferred replacement valve, there continues to be a strong need to develop better and more reliable BHVs through and improved the general understanding of BHV failure mechanisms. The major technological hurdle for the lifespan of the BHV implant continues to be the durability of the constituent leaflet biomaterials, which if improved can lead to substantial clinical impact. In order to develop improved solutions for BHV biomaterials, it is critical to have a better understanding of the inherent biomechanical behaviors of the leaflet biomaterials, including chemical treatment technologies, the impact of repetitive mechanical loading, and the inherent failure modes. This review seeks to provide a comprehensive overview of these issues, with a focus on developing insight on the mechanisms of BHV function and failure. Additionally, this review provides a detailed summary of the computational biomechanical simulations that have been used to inform and develop a higher level of understanding of BHV tissues and their failure modes. Collectively, this information should serve as a tool not only to infer reliable and dependable prosthesis function, but also to instigate and facilitate the design of future bioprosthetic valves and clinically impact cardiology.

Fluid-Structure Interaction Simulation of Prosthetic Aortic Valves: Comparison between Immersed Boundary and Arbitrary Lagrangian-Eulerian Techniques for the Mesh Representation

In recent years the role of FSI (fluid-structure interaction) simulations in the analysis of the fluid-mechanics of heart valves is becoming more and more important, being able to capture the interaction between the blood and both the surrounding biological tissues and the valve itself. When setting up an FSI simulation, several choices have to be made to select the most suitable approach for the case of interest: in particular, to simulate flexible leaflet cardiac valves, the type of discretization of the fluid domain is crucial, which can be described with an ALE (Arbitrary Lagrangian-Eulerian) or an Eulerian formulation. The majority of the reported 3D heart valve FSI simulations are performed with the Eulerian formulation, allowing for large deformations of the domains without compromising the quality of the fluid grid. Nevertheless, it is known that the ALE-FSI approach guarantees more accurate results at the interface between the solid and the fluid. The goal of this paper is to describe the same aortic valve model in the two cases, comparing the performances of an ALE-based FSI solution and an Eulerian-based FSI approach. After a first simplified 2D case, the aortic geometry was considered in a full 3D set-up. The model was kept as similar as possible in the two settings, to better compare the simulations’ outcomes. Although for the 2D case the differences were unsubstantial, in our experience the performance of a full 3D ALE-FSI simulation was significantly limited by the technical problems and requirements inherent to the ALE formulation, mainly related to the mesh motion and deformation of the fluid domain. As a secondary outcome of this work, it is important to point out that the choice of the solver also influenced the reliability of the final results.

Impact of different aortic valve calcification patterns on the outcome of transcatheter aortic valve implantation: A finite element study

Transcatheter aortic valve implantation (TAVI) can treat symptomatic patients with calcific aortic stenosis. However, the severity and distribution of the calcification of valve leaflets can impair the TAVI efficacy. Here we tackle this issue from a biomechanical standpoint, by finite element simulation of a widely adopted balloon-expandable TAVI in three models representing the aortic root with different scenarios of calcific aortic stenosis. We developed a modeling approach realistically accounting for aortic root pressurization and complex anatomy, detailed calcification patterns, and for the actual stent deployment through balloon-expansion. Numerical results highlighted the dependency on the specific calcification pattern of the "dog-boning" of the stent. Also, local stent distortions were associated with leaflet calcifications, and led to localized gaps between the TAVI stent and the aortic tissues, with potential implications in terms of paravalvular leakage. High stresses were found on calcium deposits, which may be a risk factor for stroke; their magnitude and the extent of the affected regions substantially increased for the case of an "arc-shaped" calcification, running from commissure to commissure. Moreover, high stresses due to the interaction between the aortic wall and the leaflet calcifications were computed in the annular region, suggesting an increased risk for annular damage. Our analyses suggest a relation between the alteration of the stresses in the native anatomical components and prosthetic implant with the presence and distribution of relevant calcifications. This alteration is dependent on the patient-specific features of the calcific aortic stenosis and may be a relevant indicator of suboptimal TAVI results.

Hemodynamic and thrombogenic analysis of a trileaflet polymeric valve using a fluid-structure interaction approach

Surgical valve replacement in patients with severe calcific aortic valve disease using either bioprosthetic or mechanical heart valves is still limited by structural valve deterioration for the former and thrombosis risk mandating anticoagulant therapy for the latter. Prosthetic polymeric heart valves have the potential to overcome the inherent material and design limitations of these valves, but their development is still ongoing. The aim of this study was to characterize the hemodynamics and thrombogenic potential of the Polynova polymeric trileaflet valve prototype using a fluid-structure interaction (FSI) approach. The FSI model replicated experimental conditions of the valve as tested in a left heart simulator. Hemodynamic parameters (transvalvular pressure gradient, flow rate, maximum velocity, and effective orifice area) were compared to assess the validity of the FSI model. The thrombogenic footprint of the polymeric valve was evaluated using a Lagrangian approach to calculate the stress accumulation (SA) values along multiple platelet trajectories and their statistical distribution. In the commissural regions, platelets were exposed to the highest SA values because of highest stress levels combined with local reverse flow patterns and vortices. Stress-loading waveforms from representative trajectories in regions of interest were emulated in our hemodynamic shearing device (HSD). Platelet activity was measured using our platelet activation state (PAS) assay and the results confirmed the higher thrombogenic potential of the commissural hotspots. In conclusion, the proposed method provides an in depth analysis of the hemodynamic and thrombogenic performance of the polymer valve prototype and identifies locations for further design optimization.

Three-dimensional macro-scale assessment of regional and temporal wall shear stress characteristics on aortic valve leaflets

The aortic valve (AV) achieves unidirectional blood flow between the left ventricle and the aorta. Although hemodynamic stresses have been shown to regulate valvular biology, the native wall shear stress (WSS) experienced by AV leaflets remains largely unknown. The objective of this study was to quantify computationally the macro-scale leaflet WSS environment using fluid-structure interaction modeling. An arbitrary Lagrangian-Eulerian approach was implemented to predict valvular flow and leaflet dynamics in a three-dimensional AV geometry subjected to physiologic transvalvular pressure. Local WSS characteristics were quantified in terms of temporal shear magnitude (TSM), oscillatory shear index (OSI) and temporal shear gradient (TSG). The dominant radial WSS predicted on the leaflets exhibited high amplitude and unidirectionality on the ventricularis (TSM>7.50 dyn/cm(2), OSI < 0.17, TSG>325.54 dyn/cm(2) s) but low amplitude and bidirectionality on the fibrosa (TSM < 2.73 dyn/cm(2), OSI>0.38, TSG < 191.17 dyn/cm(2) s). The radial WSS component computed in the leaflet base, belly and tip demonstrated strong regional variability (ventricularis TSM: 7.50-22.32 dyn/cm(2), fibrosa TSM: 1.26-2.73 dyn/cm(2)). While the circumferential WSS exhibited similar spatially dependent magnitude (ventricularis TSM: 1.41-3.40 dyn/cm(2), fibrosa TSM: 0.42-0.76 dyn/cm(2)) and side-specific amplitude (ventricularis TSG: 101.73-184.43 dyn/cm(2) s, fibrosa TSG: 41.92-54.10 dyn/cm(2) s), its temporal variations were consistently bidirectional (OSI>0.25). This study provides new insights into the role played by leaflet-blood flow interactions in valvular function and critical hemodynamic stress data for the assessment of the hemodynamic theory of AV disease.