A fluid-structure interaction model of the aortic valve with coaptation and compliant aortic root

Contemporary Valvular Mechanisms of Aortic Regurgitation in Tricuspid Aortic Valves: Importance in Repair Versus Replacement Strategy

BACKGROUND

This study was performed to determine cusp causes of aortic regurgitation in patients with tricuspid aortic valves without significant aortic dilatation and define cusp pathologies amenable to surgical repair (aortic valve repair [AVr]) versus aortic valve replacement.

METHODS AND RESULTS

We retrospectively reviewed surgical reports of consecutive adults with tricuspid aortic valves undergoing surgery for clinically significant aortic regurgitation within a prospective registry from January 2005 to September 2019. Valvular mechanisms were determined by systematic in vivo intraoperative quantification methods. Of 516 patients, 287 (56%) underwent repair (AVr; mean±SD age, 59.9±12.4 years; 81% men) and 229 (44%) underwent replacement (aortic valve replacement; mean±SD age, 62.8±13.8 years [P=0.01 compared to AVr]; 67% men). A single valvular mechanism was present in 454 patients (88%), with cusp prolapse (46%), retraction (24%), and perforation (18%) being the most common. Prolapse involved the right cusp in 86% of cases and was more frequent in men (P<0.001). Two-dimensional transesophageal echocardiography accuracy for predicting mechanisms was 73% to 82% for the right cusp, 55% to 61% for the noncoronary cusp, and 0% for the left-coronary cusp. Cusp prolapse, younger age, and larger patient size were associated with successful AVr (all P<0.03), whereas retraction, perforation, older age, and concomitant mitral repair were associated with aortic valve replacement (all P<0.03).

CONCLUSIONS

Right cusp prolapse is the most frequent single valvular mechanism in patients with tricuspid aortic valve aortic regurgitation, followed by cusp retraction and perforation. The accuracy of 2-dimensional transesophageal echocardiography is limited for left and noncoronary cusp mechanistic assessment. Prolapse is associated with successful AVr, whereas retraction and perforation are associated with aortic valve replacement. With systematic intraoperative quantification methods and current surgical techniques, more than half of tricuspid aortic valve aortic regurgitation cases may be successfully repaired.

A parametric study regarding structural design of a bioprosthetic aortic valve by 3D fluid-structure interaction simulations

Since the introduction of transcatheter aortic valve (AV) implantation as a viable option, surgical bioprosthetic AVs have recently started incorporating shorter struts considering future valve-in-valve procedures. However, the effect of leaflet coaptation geometry on the longevity of these valves remains unexplored. To address this gap, we performed a finite element analysis on bioprosthetic AVs with varying strut heights using a two-way fluid-structure interaction method. To establish a baseline, we used a standard height based on a rendered platform image of the CE PERIMOUNT Magna Ease valve from Edward Lifesciences in Irvine, CA. Bovine pericardium properties were assigned to the leaflets, while normal saline properties were used as the recirculating fluid in hemodynamic simulations. The physiological pressure profile of the cardiac cycle was applied between the aorta and left ventricle. We calculated blood flow velocity, effective orifice area (EOA), and mechanical stress on the leaflets. The results reveal that as the strut height increases, the stroke volume increases, leakage volume decreases, and EOA improves. Additionally, the maximum mechanical stress experienced by the leaflet decreases by 62% as the strut height increases to 1.2 times the standard height. This research highlights that a low-strut design in bioprosthetic AVs may negatively affect their durability, which can be useful in design of next-generation bioprosthetic AVs.

A predictive surrogate model for hemodynamics and structural prediction in abdominal aorta for different physiological conditions

BACKGROUND AND OBJECTIVE

This study investigates the application of a Predictive Surrogate Model (PSM) for the prediction of the fluid and solid variables in the abdominal aorta by integrating Proper Orthogonal Decomposition (POD) and Long Short-Term Memory (LSTM) techniques.

METHODS

The Fluid-Structure Interaction (FSI) solver, which serves as the Full-Order Model (FOM), can capture the blood hemodynamics and structural mechanics precisely for a variety of physiological states, namely the rest and exercise conditions.

RESULTS

Detailed analyses have been conducted on velocity components, pressure, Wall Shear Stress (WSS), and Oscillatory Shear Index (OSI) variables. Firstly, the reconstruction error has been derived based on a specific number of POD bases to assess the Reduced Order Model (ROM). Notably, the reconstruction error for velocity components in the rest condition is one order of magnitude higher than that in the exercise condition, yet both remained below 10%. This error for pressure is even more minimal, being less than 1%.

CONCLUSIONS

The PSM is evaluated against rest and exercise conditions, exhibiting promising results despite the inherent complexities of the physiological conditions. Despite the inherent complexities of phenomena in the aorta, the predictive model demonstrates consistent error magnitudes for velocity components and wall-related indices, while solid variables show slightly higher errors.

Research on fatigue optimization simulation of polymeric heart valve based on the iterative sub-regional thickened method

Prosthetic polymeric heart valves (PHVs) have the potential to overcome the inherent material and design limitations of traditional valves in the treatment of valvular heart disease; however, their durability remains limited. Optimal design of the valve structure is necessary to improve their durability. This study aimed to enhance the fatigue resistance of PHVs by improving the stress distribution. Iterative subregional thickening of the leaflets was used, and the mechanical stress distribution and hemodynamics of these polymeric tri-leaflet valves were characterized using a fluid-structure interaction approach. Subregional thickening led to a reduction in stress concentration on the leaflet, with the effective orifice area still meeting ISO 5840-3 and the regurgitant volume achieving a similar value to those in previous studies. The maximum stress in the final iteration was reduced by 28% compared with that of the prototype. The proposed method shows potential for analyzing the stress distribution and hemodynamic performance of subregional thickened valves and can further improve the durability of PHVs.

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.

Fluid-structure interaction modeling of compliant aortic valves using the lattice Boltzmann CFD and FEM methods

The lattice Boltzmann method (LBM) has been increasingly used as a stand-alone CFD solver in various biomechanical applications. This study proposes a new fluid-structure interaction (FSI) co-modeling framework for the hemodynamic-structural analysis of compliant aortic valves. Toward that goal, two commercial software packages are integrated using the lattice Boltzmann (LBM) and finite element (FE) methods. The suitability of the LBM-FE hemodynamic FSI is examined in modeling healthy tricuspid and bicuspid aortic valves (TAV and BAV), respectively. In addition, a multi-scale structural approach that has been employed explicitly recognizes the heterogeneous leaflet tissues and differentiates between the collagen fiber network (CFN) embedded within the elastin matrix of the leaflets. The CFN multi-scale tissue model is inspired by monitoring the distribution of the collagen in 15 porcine leaflets. Different simulations have been examined, and structural stresses and resulting hemodynamics are analyzed. We found that LBM-FE FSI approach can produce good predictions for the flow and structural behaviors of TAV and BAV and correlates well with those reported in the literature. The multi-scale heterogeneous CFN tissue structural model enhances our understanding of the mechanical roles of the CFN and the elastin matrix behaviors. The importance of LBM-FE FSI also emerges in its ability to resolve local hemodynamic and structural behaviors. In particular, the diastolic fluctuating velocity phenomenon near the leaflets is explicitly predicted, providing vital information on the flow transient nature. The full closure of the contacting leaflets in BAV is also demonstrated. Accordingly, good structural kinematics and deformations are captured for the entire cardiac cycle.

A design-based model of the aortic valve for fluid-structure interaction

This paper presents a new method for modeling the mechanics of the aortic valve and simulates its interaction with blood. As much as possible, the model construction is based on first principles, but such that the model is consistent with experimental observations. We require that tension in the leaflets must support a pressure, then derive a system of partial differential equations governing its mechanical equilibrium. The solution to these differential equations is referred to as the predicted loaded configuration; it includes the loaded leaflet geometry, fiber orientations and tensions needed to support the prescribed load. From this configuration, we derive a reference configuration and constitutive law. In fluid-structure interaction simulations with the immersed boundary method, the model seals reliably under physiological pressures and opens freely over multiple cardiac cycles. Further, model closure is robust to extreme hypo- and hypertensive pressures. Then, exploiting the unique features of this model construction, we conduct experiments on reference configurations, constitutive laws and gross morphology. These experiments suggest the following conclusions: (1) The loaded geometry, tensions and tangent moduli primarily determine model function. (2) Alterations to the reference configuration have little effect if the predicted loaded configuration is identical. (3) The leaflets must have sufficiently nonlinear material response to function over a variety of pressures. (4) Valve performance is highly sensitive to free edge length and leaflet height. These conclusions suggest appropriate gross morphology and material properties for the design of prosthetic aortic valves. In future studies, our aortic valve modeling framework can be used with patient-specific models of vascular or cardiac flow.

Model-based analysis of the sensitivities and diagnostic implications of FFR and CFR under various pathological conditions

Although fractional flow reserve (FFR) and coronary flow reserve (CFR) are both frequently used to assess the functional severity of coronary artery stenosis, discordant results of diagnosis between FFR and CFR in some patient cohorts have been reported. In the present study, a computational model was employed to quantify the impacts of various pathophysiological factors on FFR and CFR. In addition, a hyperemic myocardial ischemic index (HMIx) was proposed as a reference for comparing the diagnostic performances of FFR and CFR. Obtained results showed that CFR was more susceptible than FFR to the influence of many pathophysiological factors unrelated to coronary artery stenosis. In particular, the numerical study proved that increasing hyperemic coronary microvascular resistance significantly elevated FFR while reducing CFR despite fixed severity of coronary artery stenosis, whereas introducing aortic valve disease only caused a significant decrease in CFR with little influence on FFR. These results provided theoretical evidence for explaining some clinical observations, such as the increased risk of discordant diagnostic results between FFR and CFR in patients with increased hyperemic microvascular resistance, and significant increase in CFR after surgical relief of severe aortic valve disease. When evaluated with respect to the predictive value for hyperemic myocardial ischemia, the performance of FFR was found to be considerably compromised in the presence of severe coronary vasodilation dysfunction or aortic valve disease, whereas the relationship between CFR and HMIx remained relatively stable, suggesting that CFR may be a more reliable indicator of myocardial ischemia under complex pathophysiological conditions.

Progressive Calcification in Bicuspid Valves: A Coupled Hemodynamics and Multiscale Structural Computations

Bicuspid aortic valve (BAV) is the most common congenital heart disease. Calcific aortic valve disease (CAVD) accounts for the majority of aortic stenosis (AS) cases. Half of the patients diagnosed with AS have a BAV, which has an accelerated progression rate. This study aims to develop a computational modeling approach of both the calcification progression in BAV, and its biomechanical response incorporating fluid-structure interaction (FSI) simulations during the disease progression. The calcification is patient-specifically reconstructed from Micro-CT images of excised calcified BAV leaflets, and processed with a novel reverse calcification technique that predicts prior states of CAVD using a density-based criterion, resulting in a multilayered calcified structure. Four progressive multilayered calcified BAV models were generated: healthy, mild, moderate, and severe, and were modeled by FSI simulations during the full cardiac cycle. A valve apparatus model, composed of the excised calcified BAV leaflets, was tested in an in-vitro pulse duplicator, to validate the severe model. The healthy model was validated against echocardiography scans. Progressive AS was characterized by higher systolic jet flow velocities (2.08, 2.3, 3.37, and 3.85 m s^-1), which induced intense vortices surrounding the jet, coupled with irregular recirculation backflow patterns that elevated viscous shear stresses on the leaflets. This study shed light on the fluid-structure mechanism that drives CAVD progression in BAV patients.

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.

Structural Responses of Integrated Parametric Aortic Valve in an Electro-Mechanical Full Heart Model

The aortic valve (AV) is located between the left ventricle and the aorta and responsible for maintaining an outward unidirectional flow. Many AV hemodynamic and structural aspects of have been extensively studied, however, more sophisticated models are needed to better understand the AV biomechanical behavior. This study deals with integrating a new parametric AV structural model with the electro-mechanical Living Heart Human Model® (LHHM). The LHHM is a finite element model simulating human heart capable of realistic electro-mechanical simulations. Different geometric metrics of AV have been examined. New integrated structural AV model within the LHHM better predict local stresses during the cardiac cycle due to the realistic boundary condition derived from the LHHM. It was found that ellipticity index (EI), calculated as the ratio between the maximal (Max) and minimal (Min) aortic annulus (AA) diameters, well correlates with measured clinical data obtained from patients undergoing computed tomography (CT) while the annular perimeter (Perim) matches the same trend. This increases the confidence in the predicted kinematic behavior, leaflets coaptation, and the overall stresses. From the clinical aspect, the new proposed coupled and integrated AV modeling can serve as a platform for design and implementation of pre-transcatheter aortic valve replacement (TAVR) procedures.

State-of-the art bicuspid aortic valve repair in 2020

Patients with a bicuspid aortic valve (BAV) frequently require surgical intervention for aortic regurgitation (AR) and/or aneurysm. Valve-preserving surgery and repair of regurgitant BAVs have evolved into an increasingly used alternative to replacement. Anatomic predictors of possible repair failures have been identified and solutions developed. Using current techniques most non-calcified BAVs can be preserved or repaired. Excellent repair durability and freedom from valve-related complications can be achieved if all pathologic components of aortic valve and root including annular dilatation are corrected. Anatomic variations must be addressed using tailored approaches.

Hemodynamic effects of support modes of LVADs on the aortic valve

As the alternative treatment for heart failure, left ventricular assist devices (LVADs) have been widely applied to clinical practice. However, the effects of the support modes of LVADs on the biomechanical states of the aortic valve are still poorly understood. Hence, the present study investigates such effects and proposes a novel fluid-structure interaction (FSI) approach that combines the lattice Boltzmann method (LBM) and finite element (FE) method. Two support modes of LVADs, namely constant speed mode and constant flow mode, which have been widely applied to clinical practice, are also designed. Results demonstrate that the support modes of LVADs could significantly affect the biomechanical states of the aortic valve and the blood flow pattern of the ascending aorta. Compared with those in the constant flow mode, the leaflets in the constant speed mode could achieve better dynamic performance and lower stress during the systolic phase. The max radial displacement of the leaflets in the constant speed mode is at 8 mm, whereas that in the constant flow mode is at 0.8 mm. Furthermore, the outflow of LVADs directly impacts the aortic surfaces of the leaflets during the diastolic phase by increasing the level of wall shear stress of the leaflets. The leaflets in the constant speed mode receive less impact than those in the constant flow mode. The condition with such minimal impact is conducive to maintaining the normal structure of leaflets and benefits the reduction of the risk of valvular diseases. In sum, the support modes of LVADs exert a crucial effect on the biomechanical environment of the aortic valve. The constant speed mode is better than the constant flow mode in terms of providing a good hemodynamic environment for the aortic valve.

Biomechanical modeling of transcatheter aortic valve replacement in a stenotic bicuspid aortic valve: deployments and paravalvular leakage

Calcific aortic valve disease (CAVD) is characterized by stiffened aortic valve leaflets. Bicuspid aortic valve (BAV) is the most common congenital heart disease. Transcatheter aortic valve replacement (TAVR) is a treatment approach for CAVD where a stent with mounted bioprosthetic valve is deployed on the stenotic valve. Performing TAVR in calcified BAV patients may be associated with post-procedural complications due to the BAV asymmetrical structure. This study aims to develop refined computational models simulating the deployments of Evolut R and PRO TAVR devices in a representative calcified BAV. The paravalvular leakage (PVL) was also calculated by computational fluid dynamics simulations. Computed tomography scan of severely stenotic BAV patient was acquired. The 3D calcium deposits were generated and embedded inside a parametric model of the BAV. Deployments of the Evolut R and PRO inside the calcified BAV were simulated in five bioprosthesis leaflet orientations. The hypothesis of asymmetric and elliptic stent deployment was confirmed. Positioning the bioprosthesis commissures aligned with the native commissures yielded the lowest PVL (15.7 vs. 29.5 mL/beat). The Evolut PRO reduced the PVL in half compared with the Evolut R (15.7 vs. 28.7 mL/beat). The proposed biomechanical computational model could optimize future TAVR treatment in BAV patients. Graphical abstract.

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.

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.

A New Growth Model for Aortic Valve Calcification

Calcific aortic valve disease (CAVD) is a progressive disease in which minerals accumulate in the tissue of the aortic valve cusps, stiffening them and preventing valve opening and closing. The process of valve calcification was found to be similar to that of bone formation including cell differentiation to osteoblast-like cells. Studies have shown the contribution of high strains to calcification initiation and growth process acceleration. In this paper, a new strain-based calcification growth model is proposed. The model aims to explain the unique shape of the calcification and other disease characteristics. The calcification process was divided into two stages: Calcification initiation and calcification growth. The initiation locations were based on previously published findings and a reverse calcification technique (RCT), which uses computed tomography (CT) scans of patients to reveal the calcification initiation point. The calcification growth process was simulated by a finite element model of one aortic valve cusp loaded with cyclic loading. Similar to Wolff's law, describing bone response to stress, our model uses strains to drive calcification formation. The simulation grows calcification from its initiation point to its full typical stenotic shape. Study results showed that the model was able to reproduce the typical calcification growth pattern and shape, suggesting that strain is the main driving force behind calcification progression. The simulation also sheds light on other disease characteristics, such as calcification growth acceleration as the disease progresses, as well as sensitivity to hypertension.

A framework for designing patient-specific bioprosthetic heart valves using immersogeometric fluid-structure interaction analysis

Numerous studies have suggested that medical image derived computational mechanics models could be developed to reduce mortality and morbidity due to cardiovascular diseases by allowing for patient-specific surgical planning and customized medical device design. In this work, we present a novel framework for designing prosthetic heart valves using a parametric design platform and immersogeometric fluid-structure interaction (FSI) analysis. We parameterize the leaflet geometry using several key design parameters. This allows for generating various perturbations of the leaflet design for the patient-specific aortic root reconstructed from the medical image data. Each design is analyzed using our hybrid arbitrary Lagrangian-Eulerian/immersogeometric FSI methodology, which allows us to efficiently simulate the coupling of the deforming aortic root, the parametrically designed prosthetic valves, and the surrounding blood flow under physiological conditions. A parametric study is performed to investigate the influence of the geometry on heart valve performance, indicated by the effective orifice area and the coaptation area. Finally, the FSI simulation result of a design that balances effective orifice area and coaptation area reasonably well is compared with patient-specific phase contrast magnetic resonance imaging data to demonstrate the qualitative similarity of the flow patterns in the ascending aorta.

INFLUENCE OF LEAFLET’S MATRIX STIFFNESS AND FIBER ORIENTATION ON THE OPENING DYNAMICS OF A PROSTHETIC TRILEAFLET HEART VALVE

Biological valves are employed for aortic valve substitution since a long time but there is a growing effort toward the development of new engineered tissues, in which the complex mechanical response of native leaflets is replicated using composite materials consisting of a soft matrix with embedded reinforcing fibers. The main goal of the present study is to investigate the influence that variations on fiber orientation and matrix stiffness may have on valve dynamics. To this aim, a fluid–structure interaction (FSI) model of a trileaflet valve was implemented in which the opening phase was simulated and leaflet matrix stiffness and fiber orientation were varied in the framework of an anisotropic hyperelastic strain energy function. Results show that both parameters may affect significantly transvalvular pressure gradient and effective orifice area (EOA). For the opening phase of the valve examined, less favorable flow conditions were found when preferred fiber orientation is circumferential, due to lower maximum EOA achievable. Such configuration in combination with stiffer matrix may result in significant degradation of valve performances. Overall fiber orientation can potentially be taylored to optimize valve dynamics, provided also structural aspects that may be prominent in the closure phase, are considered.

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.

The impact of the aortic valve impairment on the distant coronary arteries hemodynamics: a fluid-structure interaction study

Atherosclerosis is still the leading cause of death in the developed world. Although its initiation and progression is a complex multifactorial process, it is well known that blood flow-induced wall shear stress (WSS) is an important factor involved in early atherosclerotic plaque initiation. In recent clinical studies, it was established that the regional pathologies of the aortic valve can be involved in the formation of atherosclerotic plaques. However, the impact of hemodynamic effects is not yet fully elucidated for disease initiation and progression. In this study, our developed 3D global fluid-structure interaction model of the aortic root incorporating coronary arteries is used to investigate the possible interaction between coronary arteries and aortic valve pathologies. The coronary hemodynamics was examined and quantified for different degrees of aortic stenosis varying from nonexistent to severe. For the simulated healthy model, the calculated WSS varied between 0.41 and 1.34 Pa which is in the atheroprotective range. However, for moderate and severe aortic stenoses, wide regions of the coronary structures, especially the proximal sections around the first bifurcation, were exposed to lower values of WSS and therefore they were prone to atherosclerosis even in the case of healthy coronary arteries.

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.

Pledget-Armed Sutures Affect the Haemodynamic Performance of Biologic Aortic Valve Substitutes: A Preliminary Experimental and Computational Study

Surgical aortic valve replacement is the most common procedure of choice for the treatment of severe aortic stenosis. Bioprosthetic valves are traditionally sewed-in the aortic root by means of pledget-armed sutures during open-heart surgery. Recently, novel bioprostheses which include a stent-based anchoring system have been introduced to allow rapid implantation, therefore reducing the duration and invasiveness of the intervention. Different effects on the hemodynamics were clinically reported associated with the two technologies. The aim of this study was therefore to investigate whether the differences in hemodynamic performances are an effect of different anchoring systems. Two commercially available bio-prosthetic aortic valves, one sewed-in with pledget-armed sutures and one rapid-deployment, were thus tested in this study by means of a combined approach of experimental and computational tools. In vitro experiments were performed to evaluate the overall hydrodynamic performance under identical standard conditions; computational fluid dynamics analyses were set-up to explore local flow variations due to different design of the anchoring system. The results showed how the performance of cardiac valve substitutes is negatively affected by the presence of pledget-armed sutures. These are causing flow disturbances, which in turn increase the mean pressure gradient and decrease the effective orifice area. The combined approach of experiments and numerical simulations can be effectively used to quantify the detailed relationship between local fluid-dynamics and overall performances associated with different valve technologies.

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.

The effect of pathologic venous valve on neighboring valves: fluid-structure interactions modeling

Understanding the hemodynamics surrounding the venous valve environment is of a great importance for prosthetic valves design. The present study aims to evaluate the effect of leaflets' stiffening process on the venous valve hemodynamics, valve's failure on the next proximal valve hemodynamics and valve's failure in a secondary daughter vein on the healthy valve hemodynamics in the main vein when both of these valves are distal to a venous junction. Fully coupled, two-way fluid-structure interaction computational models were developed and employed. The sinus pocket region experiences the lowest fluid shear stress, and the base region of the sinus side of the leaflet experiences the highest tissue stress. The leaflets' stiffening increases the tissue stress the valve is experiencing in a very low fluid shear region. A similar effect occurs with the proximal healthy valve as a consequence of the distal valve's failure and with the mother vein valve as a consequence of daughter vein valve's failure. Understanding the described mechanisms may be helpful for elucidating the venous valve stiffness-function relationship in nature, the reasons for a retrograde development of reflux and the relationship between venous valves located near venous junctions, and for designing better prosthetic valves and for improving their positioning.

Effects of Leaflet Design on Transvalvular Gradients of Bioprosthetic Heart Valves

Bioprosthetic aortic valves (BAVs) are becoming the prostheses of choice in heart valve replacement. The objective of this paper is to assess the effects of leaflet geometry on the mechanics and hemodynamics of BAVs in a fluid structure interaction model. The curvature and angle of leaflets were varied in 10 case studies whereby the following design parameters were altered: a circular arch, a line, and a parabola for the radial curvature, and a circular arch, a spline, and a parabola for the circumferential curvature. Six different leaflet angles (representative of the inclination of the leaflets toward the surrounding aortic wall) were analyzed. The 3-dimensional geometry of the models were created using SolidWorks, Pointwise was used for meshing, and Comsol Multiphysics was used for implicit finite element calculations. Realistic loading was enforced by considering the time-dependent strongly-coupled interaction between blood flow and leaflets. Higher mean pressure gradients as well as von Mises stresses were obtained with a parabolic or circular curvature for radial curvature or a parabolic or spline curvature for the circumferential curvature. A smaller leaflet angle was associated with a lower pressure gradient, and, a lower von Mises stress. The leaflet curvature and angle noticeably affected the speed of valve opening, and closing. When a parabola was used for circumferential or radial curvature, leaflets displacements were asymmetric, and they opened and closed more slowly. A circular circumferential leaflet curvature, a linear leaflet radial curvature, and leaflet inclination toward the surrounding aortic wall were associated with superior BAVs mechanics.

Immersed boundary-finite element model of fluid-structure interaction in the aortic root

It has long been recognized that aortic root elasticity helps to ensure efficient aortic valve closure, but our understanding of the functional importance of the elasticity and geometry of the aortic root continues to e-volve as increasingly detailed in vivo imaging data become available. Herein, we describe a fluid-structure interaction model of the aortic root, including the aortic valve leaflets, the sinsuses of Valsalva, the aortic annulus, and the sinotubular junction, that employs a version of Peskin's immersed boundary (IB) method with a finite element (FE) description of the structural elasticity. As in earlier work, we use a fiber-based model of the valve leaflets, but this study extends earlier IB models of the aortic root by employing an incompressible hyperelastic model of the mechanics of the sinuses and ascending aorta using a constitutive law fit to experimental data from human aortic root tissue. In vivo pressure loading is accounted for by a backward displacement method that determines the unloaded configurations of the root model. Our model yields realistic cardiac output at physiological pressures, with low transvalvular pressure differences during forward flow, minimal regurgitation during valve closure, and realistic pressure loads when the valve is closed during diastole. Further, results from high-resolution computations indicate that although the detailed leaflet and root kinematics show some grid sensitivity, our IB model of the aortic root nonetheless produces essentially grid-converged flow rates and pressures at practical grid spacings for the high-Reynolds number flows of the aortic root. These results thereby clarify minimum grid resolutions required by such models when used as stand-alone models of the aortic valve as well as when used to provide models of the outflow valves in models of left ventricular fluid dynamics.

Fluid-structure interaction modeling of calcific aortic valve disease using patient-specific three-dimensional calcification scans

Calcific aortic valve disease (CAVD) is characterized by calcification accumulation and thickening of the aortic valve cusps, leading to stenosis. The importance of fluid flow shear stress in the initiation and regulation of CAVD progression is well known and has been studied recently using fluid-structure interaction (FSI) models. While cusp calcifications are three-dimensional (3D) masses, previously published FSI models have represented them as either stiffened or thickened two-dimensional (2D) cusps. This study investigates the hemodynamic effect of these calcifications employing FSI models using 3D patient-specific calcification masses. A new reverse calcification technique (RCT) is used for modeling different stages of calcification growth based on the spatial distribution of calcification density. The RCT is applied to generate the 3D calcification deposits reconstructed from a patient-specific CT scans. Our results showed that consideration of 3D calcification deposits led to both higher fluid shear stresses and unique fluid shear stress distribution on the aortic side of the cusps that may have an impact on the calcification growth rate. However, the flow did not seem to affect the geometry of the calcification during the growth phase.

Imaging analysis of collagen fiber networks in cusps of porcine aortic valves: effect of their local distribution and alignment on valve functionality

The cusps of native aortic valve (AV) are composed of collagen bundles embedded in soft tissue, creating a heterogenic tissue with asymmetric alignment in each cusp. This study compares native collagen fiber networks (CFNs) with a goal to better understand their influence on stress distribution and valve kinematics. Images of CFNs from five porcine tricuspid AVs are analyzed and fluid-structure interaction models are generated based on them. Although the valves had similar overall kinematics, the CFNs had distinctive influence on local mechanics. The regions with dilute CFN are more prone to damage since they are subjected to higher stress magnitudes.

Biomechanical evaluation of the pathophysiologic developmental mechanisms of mitral valve prolapse: effect of valvular morphologic alteration

Mitral valve prolapse (MVP) refers to an excessive billowing of the mitral valve (MV) leaflets across the mitral annular plane into the left atrium during the systolic portion of the cardiac cycle. The underlying mechanisms for the development of MVP and mitral regurgitation in association with MV tissue remodeling are still unclear. We performed computational MV simulations to investigate the pathophysiologic developmental mechanisms of MVP. A parametric MV geometry model was utilized for this study. Posterior leaflet enlargement and posterior chordal elongation models were created by adjusting the geometry of the posterior leaflet and chordae, respectively. Dynamic finite element simulations of MV function were performed over the complete cardiac cycle. Computational simulations demonstrated that enlarging posterior leaflet area increased large stress concentration in the posterior leaflets and chordae, and posterior chordal elongation decreased leaflet coaptation. When MVP was accompanied by both posterior leaflet enlargement and chordal elongation simultaneously, the posterior leaflet was exposed to extremely large prolapse with a substantial lack of leaflet coaptation. These data indicate that MVP development is closely related to tissue alterations of the leaflets and chordae. This biomechanical evaluation strategy can help us better understand the pathophysiologic developmental mechanisms of MVP.

The Influence of the Aortic Root Geometry on Flow Characteristics of a Prosthetic Heart Valve

In this paper, performance of aortic heart valve prosthesis in different geometries of the aortic root is investigated experimentally. The objective of this investigation is to establish a set of parameters, which are associated with abnormal flow patterns due to the flow through a prosthetic heart valve implanted in the patients that had certain types of valve diseases prior to the valve replacement. Specific valve diseases were classified into two clinical categories and were correlated with the corresponding changes in aortic root geometry while keeping the aortic base diameter fixed. These categories correspond to aortic valve stenosis and aortic valve insufficiency. The control case that corresponds to the aortic root of a patient without valve disease was used as a reference. Experiments were performed at test conditions corresponding to 70 beats/min, 5.5 L/min target cardiac output, and a mean aortic pressure of 100 mmHg. By varying the aortic root geometry, while keeping the diameter of the orifice constant, it was possible to investigate corresponding changes in the levels of Reynolds shear stress and establish the possibility of platelet activation and, as a result of that, the formation of blood clots.

Numerical model of total artificial heart hemodynamics and the effect of its size on stress accumulation

The total artificial heart (TAH) is a bi-ventricular mechanical circulatory support device that replaces the heart in patients with end-stage congestive heart failure. The device acts as blood pump via pneumatic activation of diaphragms altering the volume of the ventricular chambers. Flow in and out of the ventricles is controlled by mechanical heart valves. The aim of this study is to evaluate the flow regime in the TAH and to estimate the thrombogenic potential during systole. Toward that goal, three numerical models of TAHs of differing sizes, that include the deforming diaphragm and the blood flow from the left chamber to the aorta, are introduced. A multiphase model with injection of platelet particles is employed to calculate their trajectories. The shear stress accumulation in the three models are calculated along the platelets trajectories and their probability density functions, which represent the `thrombogenic footprint' of the device are compared. The calculated flow regime successfully captures the mitral regurgitation and the flows that open and close the aortic valve during systole. Physiological velocity magnitudes are found in all three models, with higher velocities and increased stress accumulation predicted for smaller devices.

EFFECT OF THE POSITION OF THE CORONARY SINUS ORIFICE ON AORTIC LEAFLET COAPTATION

Background: The various components of the aortic root maintain a particular geometric relationship to guarantee unobstructed blood flow across the aortic valve and valve competence. Objective: To quantify the effect of the position of the coronary sinus orifice (CSO) on aortic leaflet coaptation. Methods: 2D and 3D finite element models of an aortic valve and root were constructed, with the CSO located on the bottom and middle of the sinuses. ADINA fluid-structure interaction solver was employed to perform computational simulation. Results: The mean sinus pressure with left and right CSO was 1.02E+4 Pa and 1.03E+4 Pa, respectively, and the average leaflet pressure with left and right CSO was 1.06E+4 Pa and 1.05E+4 Pa, respectively, for the model with CSO located in the middle and bottom of the sinus. The leaflet summit displacement differences of the CSO position on the bottom and middle between left and right coronary sinuses and none were 11.56, -107.57, 16.17 and -92.86 μm, respectively. Conclusions: The position of the CSO affects the pressure distribution of the aortic root. The local high pressure results in symmetrical deformation of the three leaflets, and decreases the risk of leaflet mismatch in coaptation.

Leaflet reconstructive techniques for aortic valve repair

BACKGROUND

Refining leaflet reconstruction has become a primary issue in aortic valve repair. This descriptive analysis reviews leaflet pathology, repair techniques, and early results in a prospective regulatory trial of aortic valve repair.

METHODS

Sixty-five patients underwent valve repair for predominant moderate to severe aortic insufficiency (AI). The mean age was 63 ± 13 years, and 69% of the patients were male. Ascending aortic/root replacement was required in 62%. As a first step, ring annuloplasty was performed, and then leaflet repair included leaflet plication for prolapse, nodular unfolding, double pericardial patching of commissural defects or holes, complete pericardial leaflet replacement, leaflet extension, and Gore-Tex reinforcement. Leaflet techniques and causes of adverse outcomes were evaluated.

RESULTS

The follow-up time was 2-years maximal and 0.9 years mean, with a survival of 97%. Eighty percent of patients required repair of leaflet defects: leaflet prolapse (52/65-80%), ruptured commissures (6/65-9%), leaflet holes (4/65-6%), and nodular retraction (6/65-9%). The average preoperative AI grade of 2.9 ± 0.8 fell to 0.7 ± 0.7 (p < 0.0001). Three patients (4.6%) required interval valve replacement because of (1) suture untying, (2) iatrogenic leaflet tear, or (3) diphtheroid endocarditis. Five other patients experienced grade 2 or grade 3 AI: probable suture untying in 1 patient, ineffective leaflet extensions in 2 patients, and unsuccessful Gore-Tex reinforcements in 2 patients. Two patients with single pericardial leaflet replacements and all those with double pericardial reconstructions did well.

CONCLUSIONS

Leaflet defects are common in patients with moderate to severe AI. Leaflet plication, nodular unfolding, and double pericardial patching performed well. Gore-Tex and leaflet extension seemed less satisfactory. Standardization and experience with leaflet reconstruction will be important for optimizing the outcomes of aortic valve repair.