Thinner biological tissues induce leaflet flutter in aortic heart valve replacements

Simulated Effects of Acute Left Ventricular Myocardial Infarction on Mitral Regurgitation in an Ovine Model

Ischemic mitral regurgitation (IMR) occurs from incomplete coaptation of the mitral valve (MV) after myocardial infarction (MI), typically worsened by continued remodeling of the left ventricular (LV). The importance of LV remodeling is clear as IMR is induced by the post-MI dual mechanisms of mitral annular dilation and leaflet tethering from papillary muscle (PM) distension via the MV chordae tendineae (MVCT). However, the detailed etiology of IMR remains poorly understood, in large part due to the complex interactions of the MV and the post-MI LV remodeling processes. Given the patient-specific anatomical complexities of the IMR disease processes, simulation-based approaches represent an ideal approach to improve our understanding of this deadly disease. However, development of patient-specific models of left ventricle-mitral valve (LV-MV) interactions in IMR are complicated by the substantial variability and complexity of the MR etiology itself, making it difficult to extract underlying mechanisms from clinical data alone. To address these shortcomings, we developed a detailed ovine LV-MV finite element (FE) model based on extant comprehensive ovine experimental data. First, an extant ovine LV FE model (Sci. Rep. 2021 Jun 29;11(1):13466) was extended to incorporate the MV using a high fidelity ovine in vivo derived MV leaflet geometry. As it is not currently possible to image the MVCT in vivo, a functionally equivalent MVCT network was developed to create the final LV-MV model. Interestingly, in pilot studies, the MV leaflet strains did not agree well with known in vivo MV leaflet strain fields. We then incorporated previously reported MV leaflet prestrains (J. Biomech. Eng. 2023 Nov 1;145(11):111002) in the simulations. The resulting LV-MV model produced excellent agreement with the known in vivo ovine MV leaflet strains and deformed shapes in the normal state. We then simulated the effects of regional acute infarctions of varying sizes and anatomical locations by shutting down the local myocardial contractility. The remaining healthy (noninfarcted) myocardium mechanical behaviors were maintained, but allowed to adjust their active contractile patterns to maintain the prescribed pressure-volume loop behaviors in the acute post-MI state. For all cases studied, the LV-MV simulation demonstrated excellent agreement with known LV and MV in vivo strains and MV regurgitation orifice areas. Infarct location was shown to play a critical role in resultant MV leaflet strain fields. Specifically, extensional deformations of the posterior leaflets occurred in the posterobasal and laterobasal infarcts, while compressive deformations of the anterior leaflet were observed in the anterobasal infarct. Moreover, the simulated posterobasal infarct induced the largest MV regurgitation orifice area, consistent with experimental observations. The present study is the first detailed LV-MV simulation that reveals the important role of MV leaflet prestrain and functionally equivalent MVCT for accurate predictions of LV-MV interactions. Importantly, the current study further underscored simulation-based methods in understanding MV function as an integral part of the LV.

On the role of aortic valve architecture for physiological hemodynamics and valve replacement, Part I: Flow configuration and vortex dynamics

Aortic valve replacement has become an increasing concern due to the rising prevalence of aortic stenosis in an ageing population. Existing replacement options have limitations, necessitating the development of improved prosthetic aortic valves. In this study, flow characteristics during systole in a stenotic aortic valve case are compared with those downstream of two newly designed surgical bioprosthetic aortic valves (BioAVs). To do so, advanced three-dimensional fluid-structure interaction simulations are conducted and dedicated analysis methods to investigate jet flow configuration and vortex dynamics are developed. Our findings reveal that the stenotic case maintains a high jet flow eccentricity due to a fixed orifice geometry, resulting in flow separation and increased vortex stretching and tilting in the commissural low-flow regions. One BioAV design introduces non-axisymmetric leaflet motion, which reduces the maximum jet velocity and forms more vortical structures. The other BioAV design produces a fixed symmetric triangular jet shape due to non-moving leaflets and exhibits favourable vorticity attenuation, revealed by negative temporally and spatially averaged projected vortex stretching values, and significantly reduced drag. Therefore, this study highlights the benefits of custom-designed aortic valves in the context of their replacement through comprehensive and novel flow analyses. The results emphasise the importance of analysing jet flow, vortical structures, momentum balance and vorticity transport for thoroughly evaluating aortic valve performance.

In Vitro Proof of Concept of a First-Generation Growth-Accommodating Heart Valved Conduit for Pediatric Use

Currently available heart valve prostheses have no growth potential, requiring children with heart valve diseases to endure multiple valve replacement surgeries with compounding risks. This study demonstrates the in vitro proof of concept of a biostable polymeric trileaflet valved conduit designed for surgical implantation and subsequent expansion via transcatheter balloon dilation to accommodate the growth of pediatric patients and delay or avoid repeated open-heart surgeries. The valved conduit is formed via dip molding using a polydimethylsiloxane-based polyurethane, a biocompatible material shown here to be capable of permanent stretching under mechanical loading. The valve leaflets are designed with an increased coaptation area to preserve valve competence at expanded diameters. Four 22 mm diameter valved conduits are tested in vitro for hydrodynamics, balloon dilated to new permanent diameters of 23.26 ± 0.38 mm, and then tested again. Upon further dilation, two valved conduits sustain leaflet tears, while the two surviving devices reach final diameters of 24.38 ± 0.19 mm. After each successful dilation, the valved conduits show increased effective orifice areas and decreased transvalvular pressure differentials while maintaining low regurgitation. These results demonstrate concept feasibility and motivate further development of a polymeric balloon-expandable device to replace valves in children and avoid reoperations.

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.

Biomechanical Engineering Analysis of Pulmonary Valve Leaflet Hemodynamics and Kinematics in the Ross Procedure

The Ross procedure using the inclusion technique with anticommissural plication (ACP) is associated with excellent valve hemodynamics and favorable leaflet kinematics. The objective was to evaluate individual pulmonary cusp's biomechanics and fluttering by including coronary flow in the Ross procedure using an ex vivo three-dimensional-printed heart simulator. Ten porcine and five human pulmonary autografts were harvested from a meat abattoir and heart transplant patients. Five porcine autografts without reinforcement served as controls. The other autografts were prepared using the inclusion technique with and without ACP (ACP and NACP). Hemodynamic and high-speed videography data were measured using the ex vivo heart simulator. Although porcine autografts showed similar leaflet rapid opening and closing mean velocities, human ACP compared to NACP autografts demonstrated lower leaflet rapid opening mean velocity in the right (p = 0.02) and left coronary cusps (p = 0.003). The porcine and human autograft leaflet rapid opening and closing mean velocities were similar in all three cusps. Porcine autografts showed similar leaflet flutter frequencies in the left (p = 0.3) and noncoronary cusps (p = 0.4), but porcine NACP autografts versus controls demonstrated higher leaflet flutter frequency in the right coronary cusp (p = 0.05). The human NACP versus ACP autografts showed higher flutter frequency in the noncoronary cusp (p = 0.02). The leaflet flutter amplitudes were similar in all three cusps in both porcine and human autografts. The ACP compared to NACP autografts in the Ross procedure was associated with more favorable leaflet kinematics. These results may translate to the improved long-term durability of the pulmonary autografts.

Structural analysis of regional transcatheter aortic valve underexpansion and its implications for subclinical leaflet thrombosis

Subclinical leaflet thrombosis has been increasingly recognized following transcatheter aortic valve replacement (TAVR). Determining the risk factors is vital in preventing clinical leaflet thrombosis and ensuring long-term value durability. Clinical data have indicated that regional stent under-expansion of transcatheter aortic valves (TAVs), particularly self-expanding devices, may be associated with an increased risk of subclinical leaflet thrombosis. This study aimed to determine the effects of regional TAV frame under-expansion on leaflet kinematics, leaflet structural characteristics, and explore its impact on the likelihood of leaflet thrombosis. In this study, mild and moderate regional frame under-expansion of a 26-mm CoreValve were examined using experimental testing and computational simulations. The results indicated that regional TAV frame under-expansion impairs leaflet kinematics and reduces the range of motion in leaflets with an angle less than 120°. The reduced range of motion can increase blood stasis on the surface of the TAV leaflets. The results also demonstrated that regional frame under-expansion induced localized high-stress regions in the leaflets close to the fixed boundary edge. The increased mechanical stress can lead to accelerated tissue degeneration. The study improves our understanding of the effects of regional stent under-expansion in TAVR. Post-procedural balloon dilatation of self-expanding TAVs can potentially be advantageous in reducing leaflet distortion and normalizing leaflet stress distribution. Large-scale, prospective, and well-controlled studies are needed to further investigate regional TAV frame under-expansion effects on subclinical leaflet thrombosis and long-term valve durability.

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.

Computational Methods for Fluid-Structure Interaction Simulation of Heart Valves in Patient-Specific Left Heart Anatomies

Given the complexity of human left heart anatomy and valvular structures, the fluid–structure interaction (FSI) simulation of native and prosthetic valves poses a significant challenge for numerical methods. In this review, recent numerical advancements for both fluid and structural solvers for heart valves in patient-specific left hearts are systematically considered, emphasizing the numerical treatments of blood flow and valve surfaces, which are the most critical aspects for accurate simulations. Numerical methods for hemodynamics are considered under both the continuum and discrete (particle) approaches. The numerical treatments for the structural dynamics of aortic/mitral valves and FSI coupling methods between the solid Ωs and fluid domain Ωf are also reviewed. Future work toward more advanced patient-specific simulations is also discussed, including the fusion of high-fidelity simulation within vivo measurements and physics-based digital twining based on data analytics and machine learning techniques.

New results in the theory of plane strain flexure of incompressible isotropic hyperelastic materials

New results on the classical problem of bending by end moments for incompressible isotropic hyperelastic materials within the framework of nonlinear elasticity are investigated and presented in this paper. The particular results of concern here include (i) the adaptation of Rivlin's standard analysis to the case where one end of the beam is fixed and the other end is subjected to a bending moment; and (ii) results on the finite bending of (infinitesimally) thin isotropic hyperelastic plates which are valid for large deformations , extending the classical results from the linear elasticity theory which are restricted to small deformations. An interesting feature observed in this context is that a flexed thin plate develops an oscillatory surface along the circular arc near the free end, due to local (small) deviations of the radius of curvature. A potential application to the bending of a biological soft tissue, namely the aortic valve leaflet, is briefly described by way of an example. Finally, some new results are obtained for finite bending of hyperelastic materials that exhibit limiting chain extensibility at the molecular level and involve constraints on the deformation. The amount of bending that such materials can sustain is limited by the constraint. On using a limiting chain extensibility model, closed-form solutions for the Cauchy stress components, the bending moment and the normal out-of-plane force required to sustain the bending deformation are derived.

Poly-2-methyl-2-oxazoline–modified bioprosthetic heart valve leaflets have enhanced biocompatibility and resist structural degeneration

Bioprosthetic heart valves (BHV) fabricated from glutaraldehyde-fixed heterograft tissue, such as bovine pericardium (BP), are widely used for treating heart valve disease, a group of disorders that affects millions. Structural valve degeneration (SVD) of BHV due to both calcification and the accumulation of advanced glycation end products (AGE) with associated serum proteins limits durability. We hypothesized that BP modified with poly-2-methyl-2-oxazoline (POZ) to inhibit protein entry would demonstrate reduced accumulation of AGE and serum proteins, mitigating SVD. In vitro studies of POZ-modified BP demonstrated reduced accumulation of serum albumin and AGE. BP-POZ in vitro maintained collagen microarchitecture per two-photon microscopy despite AGE incubation, and in cell culture studies was associated with no change in tumor necrosis factor-α after exposure to AGE and activated macrophages. Comparing POZ and polyethylene glycol (PEG)–modified BP in vitro, BP-POZ was minimally affected by oxidative conditions, whereas BP-PEG was susceptible to oxidative deterioration. In juvenile rat subdermal implants, BP-POZ demonstrated reduced AGE formation and serum albumin infiltration, while calcification was not inhibited. However, BP-POZ rat subdermal implants with ethanol pretreatment demonstrated inhibition of both AGE accumulation and calcification. Ex vivo laminar flow studies with human blood demonstrated BP-POZ enhanced thromboresistance with reduced white blood cell accumulation. We conclude that SVD associated with AGE and serum protein accumulation can be mitigated through POZ functionalization that both enhances biocompatibility and facilitates ethanol pretreatment inhibition of BP calcification.

Improving transcatheter aortic valve interventional predictability via fluid-structure interaction modelling using patient-specific anatomy

Transcatheter aortic valve replacement (TAVR) is now a standard treatment for high-surgical-risk patients with severe aortic valve stenosis. TAVR is being explored for broader indications including degenerated bioprosthetic valves, bicuspid valves and for aortic valve (AV) insufficiency. It is, however, challenging to predict whether the chosen valve size, design or its orientation would produce the most-optimal haemodynamics in the patient. Here, we present a novel patient-specific evaluation framework to realistically predict the patient's AV performance with a high-fidelity fluid-structure interaction analysis that included the patient's left ventricle and ascending aorta (AAo). We retrospectively evaluated the pre- and post-TAVR dynamics of a patient who underwent a 23 mm TAVR and evaluated against the patient's virtually de-calcified AV serving as a hypothetical benchmark. Our model predictions were consistent with clinical data. Stenosed AV produced a turbulent flow during peak-systole, while aortic flow with TAVR and de-calcified AV were both in the laminar-to-turbulent transitional regime with an estimated fivefold reduction in viscous dissipation. For TAVR, dissipation was highest during early systole when valve deformation was the greatest, suggesting that an efficient valve opening may reduce energy loss. Our study demonstrates that such patient-specific modelling frameworks can be used to improve predictability and in the planning of AV interventions.

Effects of membrane and flexural stiffnesses on aortic valve dynamics: identifying the mechanics of leaflet flutter in thinner biological tissues

Valvular pathologies that induce deterioration in the aortic valve are a common cause of heart disease among aging populations. Although there are numerous available technologies to treat valvular conditions and replicate normal aortic function by replacing the diseased valve with a bioprosthetic implant, many of these devices face challenges in terms of long-term durability. One such phenomenon that may exacerbate valve deterioration and induce undesirable hemodynamic effects in the aorta is leaflet flutter, which is characterized by oscillatory motion in the biological tissues. While this behavior has been observed for thinner bioprosthetic valves, the specific underlying mechanics that lead to leaflet flutter have not previously been identified. This work proposes a computational approach to isolate the fundamental mechanics that induce leaflet flutter in thinner biological tissues during the cardiac cycle. The simulations in this work identify reduced flexural stiffness as the primary factor that contributes to increased leaflet flutter in thinner biological tissues, while decreased membrane stiffness and mass of the thinner tissues do not directly induce flutter in these valves. The results of this study provide an improved understanding of the mechanical tissue properties that contribute to flutter and offer significant insights into possible developments in the design of bioprosthetic tissues to account for and reduce the incidence of flutter.

Anisotropic elastic behavior of a hydrogel-coated electrospun polyurethane: Suitability for heart valve leaflets

Although xenograft biomaterials have been used for decades in replacement heart valves, they continue to face multiple limitations, including limited durability, mineralization, and restricted design space due to their biological origins. These issues necessitate the need for novel replacement heart valve biomaterials that are durable, non-thrombogenic, and compatible with transcatheter aortic valve replacement devices. In this study, we explored the suitability of an electrospun poly(carbonate urethane) (ES-PCU) mesh coated with a poly(ethylene glycol) diacrylate (PEGDA) hydrogel as a synthetic biomaterial for replacement heart valve leaflets. In this material design, the mesh provides the mechanical support, while the hydrogel provides the required surface hemocompatibility. We conducted a comprehensive study to characterize the structural and mechanical properties of the uncoated mesh as well as the hydrogel-coated mesh (composite biomaterial) over the estimated operational range. We found that the composite biomaterial was functionally robust with reproducible stress-strain behavior within and beyond the functional ranges for replacement heart valves, and was able to withstand the rigors of mechanical evaluation without any observable damage. In addition, the composite biomaterial displayed a wide range of mechanical anisotropic responses, which were governed by fiber orientation of the mesh, which in turn, was controlled with the fabrication process. Finally, we developed a novel constitutive modeling approach to predict the mechanical behavior of the composite biomaterial under in-plane extension and shear deformation modes. This model identified the existence of fiber-fiber mechanical interactions in the mesh that have not previously been reported. Interestingly, there was no evidence of fiber-hydrogel mechanical interactions. This important finding suggests that the hydrogel coating can be optimized for hemocompatibility independent of the structural mechanical responses required by the leaflet. This initial study indicated that the composite biomaterial has mechanical properties well-suited for replacement heart valve applications and that the electrospun mesh microarchitecture and hydrogel biological properties can be optimized independently. It also reveals that the structural mechanisms contributing to the mechanical response are more complicated than what was previously established and paves the pathway for more detailed future studies.

Leaflet Stresses During Full Device Simulation of Crimping to 6 mm in Transcatheter Aortic Valve Implantation, TAVI

BACKGROUND

With continuing growth in transcatheter aortic valve implantation for the treatment of a failing aortic valve, there is increasing interest in prosthetic valve durability and the potential damage caused to leaflets by stress. Whilst most available research into the computational prediction of leaflet stresses using finite element analysis, FEA, has focussed on variations during dynamic loading, very little appears to have been reported for the impact of crimping, even though awareness of this effect is widespread. Potentially, this has been due to the difficulty of performing full model simulations of crimping to clinically meaningful diameters.

METHOD

A full model comprising a self-expanding frame, skirt and leaflets has been developed and crimped to a final diameter of 6 mm. A detailed description is provided of the FEA setup, emphasising the importance of the skirt definition needed to successfully crimp to this small diameter. Then, an analysis of leaflet folding and stresses is presented, particularly with respect to the differences produced between leaflet thicknesses of 0.20, 0.25 and 0.30 mm and for bioprosthetic and polymeric leaflet material models.

RESULTS

In all cases, peak stresses occurred close to the modelled suture lines joining the leaflets and the skirt and high stresses were also present along axially aligned folds in the leaflets. Stresses were lower for the polymeric leaflets.

CONCLUSION

Successful simulation of crimping requires a finely resolved skirt mesh. Leaflet stresses during crimping are dependent on leaflet thickness, material properties and the ratio of leaflet volume to the available volume inside the crimped valve.

Effect of Leaflet Type and Leaflet-Stent Attachment Height on Transcatheter Aortic Valve Leaflet Thrombosis Potential

Transcatheter aortic valve replacement devices vary in leaflet material and in the height for which leaflets attach to the stented valve frame. Combinations of these features can influence leaflet dynamics, neo-sinus geometries, and fluid dynamics, thereby reducing or exacerbating the potential for blood flow stasis and leaflet thrombosis. To investigate these interconnected relationships, this study evaluated the effects of transcatheter valve leaflet type (porcine vs. bovine pericardium) and the leaflet-stent attachment height (low, mid, and high) on flow stasis and potential for leaflet thrombosis. Transcatheter valve models were manufactured and tested within an aortic simulator under pulsatile left heart hemodynamic conditions. Transvalvular hemodynamics, leaflet kinematics, and flow structures were evaluated by direct measurement, high-speed imaging, and two differing techniques of particle image velocimetry. Transcatheter valves with porcine pericardial leaflets were observed to be less stiff, exhibit a lesser resistance to flow, were associated with reduced regions of neo-sinus flow stasis, and superior sinus washout times. More elevated attachments of the leaflets were associated with less neo-sinus flow stasis. These initial results and observations suggest combinations of leaflet type and stent attachment height may reduce transcatheter aortic valve flow stasis and the potential for leaflet thrombosis.

Simulating the time evolving geometry, mechanical properties, and fibrous structure of bioprosthetic heart valve leaflets under cyclic loading

Currently, the most common replacement heart valve design is the 'bioprosthetic' heart valve (BHV), which has important advantages in that it does not require permanent anti-coagulation therapy, operates noiselessly, and has blood flow characteristics similar to the native valve. BHVs are typically fabricated from glutaraldehyde-crosslinked pericardial xenograft tissue biomaterials (XTBs) attached to a rigid, semi-flexible, or fully collapsible stent in the case of the increasingly popular transcutaneous aortic valve replacement (TAVR). While current TAVR assessments are positive, clinical results to date are generally limited to <2 years. Since TAVR leaflets are constructed using thinner XTBs, their mechanical demands are substantially greater than surgical BHV due to the increased stresses during in vivo operation, potentially resulting in decreased durability. Given the functional complexity of heart valve operation, in-silico predictive simulations clearly have potential to greatly improve the TAVR development process. As such simulations must start with accurate material models, we have developed a novel time-evolving constitutive model for pericardial xenograft tissue biomaterials (XTB) utilized in BHV (doi: 10.1016/j.jmbbm.2017.07.013). This model was able to simulate the observed tissue plasticity effects that occur in approximately in the first two years of in vivo function (50 million cycles). In the present work, we implemented this model into a complete simulation pipeline to predict the BHV time evolving geometry to 50 million cycles. The pipeline was implemented within an isogeometric finite element formulation that directly integrated our established BHV NURBS-based geometry (doi: 10.1007/s00466-015-1166-x). Simulations of successive loading cycles indicated continual changes in leaflet shape, as indicated by spatially varying increases in leaflet curvature. While the simulation model assumed an initial uniform fiber orientation distribution, anisotropic regional changes in leaflet tissue plastic strain induced a complex changes in regional fiber orientation. We have previously noted in our time-evolving constitutive model that the increases in collagen fiber recruitment with cyclic loading placed an upper bound on plastic strain levels. This effect was manifested by restricting further changes in leaflet geometry past 50 million cycles. Such phenomena was accurately captured in the valve-level simulations due to the use of a tissue-level structural-based modeling approach. Changes in basic leaflet dimensions agreed well with extant experimental studies. As a whole, the results of the present study indicate the complexity of BHV responses to cyclic loading, including changes in leaflet shape and internal fibrous structure. It should be noted that the later effect also influences changes in local mechanical behavior (i.e. changes in leaflet anisotropic tissue stress-strain relationship) due to internal fibrous structure resulting from plastic strains. Such mechanism-based simulations can help pave the way towards the application of sophisticated simulation technologies in the development of replacement heart valve technology.

Parameterization, geometric modeling, and isogeometric analysis of tricuspid valves

Approximately 1.6 million patients in the United States are affected by tricuspid valve regurgitation, which occurs when the tricuspid valve does not close properly to prevent backward blood flow into the right atrium. Despite its critical role in proper cardiac function, the tricuspid valve has received limited research attention compared to the mitral and aortic valves on the left side of the heart. As a result, proper valvular function and the pathologies that may cause dysfunction remain poorly understood. To promote further investigations of the biomechanical behavior and response of the tricuspid valve, this work establishes a parameter-based approach that provides a template for tricuspid valve modeling and simulation. The proposed tricuspid valve parameterization presents a comprehensive description of the leaflets and the complex chordae tendineae for capturing the typical three-cusp structural deformation observed from medical data. This simulation framework develops a practical procedure for modeling tricuspid valves and offers a robust, flexible approach to analyze the performance and effectiveness of various valve configurations using isogeometric analysis. The proposed methods also establish a baseline to examine the tricuspid valve's structural deformation, perform future investigations of native valve configurations under healthy and disease conditions, and optimize prosthetic valve designs.

Bioprosthetic aortic valve diameter and thickness are directly related to leaflet fluttering: Results from a combined experimental and computational modeling study

OBJECTIVE

Bioprosthetic heart valves (BHVs) are commonly used in surgical and percutaneous valve replacement. The durability of percutaneous valve replacement is unknown, but surgical valves have been shown to require reintervention after 10 to 15 years. Further, smaller-diameter surgical BHVs generally experience higher rates of prosthesis-patient mismatch, which leads to higher rates of failure. Bioprosthetic aortic valves can flutter in systole, and fluttering is associated with fatigue and failure in flexible structures. The determinants of flutter in BHVs have not been well characterized, despite their potential to influence durability.

METHODS

We use an experimental pulse duplicator and a computational fluid-structure interaction model of this system to study the role of device geometry on BHV dynamics. The experimental system mimics physiological conditions, and the computational model enables precise control of leaflet biomechanics and flow conditions to isolate the effects of variations in BHV geometry on leaflet dynamics.

RESULTS

Both experimental and computational models demonstrate that smaller-diameter BHVs yield markedly higher leaflet fluttering frequencies across a range of conditions. The computational model also predicts that fluttering frequency is directly related to leaflet thickness. A scaling model is introduced that rationalizes these findings.

CONCLUSIONS

We systematically characterize the influence of BHV diameter and leaflet thickness on fluttering dynamics. Although this study does not determine how flutter influences device durability, increased flutter in smaller-diameter BHVs may explain how prosthesis-patient mismatch could induce BHV leaflet fatigue and failure. Ultimately, understanding the effects of device geometry on leaflet kinematics may lead to more durable valve replacements.

Computational investigation of left ventricular hemodynamics following bioprosthetic aortic and mitral valve replacement

The left ventricle of the heart is a fundamental structure in the human cardiac system that pumps oxygenated blood into the systemic circulation. Several valvular conditions can cause the aortic and mitral valves associated with the left ventricle to become severely diseased and require replacement. However, the clinical outcomes of such operations, specifically the postoperative ventricular hemodynamics of replacing both valves, are not well understood. This work uses computational fluid-structure interaction (FSI) to develop an improved understanding of this effect by modeling a left ventricle with the aortic and mitral valves replaced with bioprostheses. We use a hybrid Arbitrary Lagrangian-Eulerian/immersogeometric framework to accommodate the analysis of cardiac hemodynamics and heart valve structural mechanics in a moving fluid domain. The motion of the endocardium is obtained from a cardiac biomechanics simulation and provided as an input to the proposed numerical framework. The results from the simulations in this work indicate that the replacement of the native mitral valve with a tri-radially symmetric bioprosthesis dramatically changes the ventricular hemodynamics. Most significantly, the vortical motion in the left ventricle is found to reverse direction after mitral valve replacement. This study demonstrates that the proposed computational FSI framework is capable of simulating complex multiphysics problems and can provide an in-depth understanding of the cardiac mechanics.