Immersogeometric fluid-structure interaction modeling and simulation of transcatheter aortic valve replacement

The role of computational methods in cardiovascular medicine: a narrative review

Background and Objective

Computational models of the cardiovascular system allow for a detailed and quantitative investigation of both physiological and pathological conditions, thanks to their ability to combine clinical-possibly patient-specific-data with physical knowledge of the processes underlying the heart function. These models have been increasingly employed in clinical practice to understand pathological mechanisms and their progression, design medical devices, support clinicians in improving therapies. Hinging upon a long-year experience in cardiovascular modeling, we have recently constructed a computational multi-physics and multi-scale integrated model of the heart for the investigation of its physiological function, the analysis of pathological conditions, and to support clinicians in both diagnosis and treatment planning. This narrative review aims to systematically discuss the role that such model had in addressing specific clinical questions, and how further impact of computational models on clinical practice are envisaged.

Methods

We developed computational models of the physical processes encompassed by the heart function (electrophysiology, electrical activation, force generation, mechanics, blood flow dynamics, valve dynamics, myocardial perfusion) and of their inherently strong coupling. To solve the equations of such models, we devised advanced numerical methods, implemented in a flexible and highly efficient software library. We also developed computational procedures for clinical data post-processing-like the reconstruction of the heart geometry and motion from diagnostic images-and for their integration into computational models.

Key Content and Findings

Our integrated computational model of the heart function provides non-invasive measures of indicators characterizing the heart function and dysfunctions, and sheds light on its underlying processes and their coupling. Moreover, thanks to the close collaboration with several clinical partners, we addressed specific clinical questions on pathological conditions, such as arrhythmias, ventricular dyssynchrony, hypertrophic cardiomyopathy, degeneration of prosthetic valves, and the way coronavirus disease 2019 (COVID-19) infection may affect the cardiac function. In multiple cases, we were also able to provide quantitative indications for treatment.

Conclusions

Computational models provide a quantitative and detailed tool to support clinicians in patient care, which can enhance the assessment of cardiac diseases, the prediction of the development of pathological conditions, and the planning of treatments and follow-up tests.

Artificial Intelligence in Transcatheter Aortic Valve Replacement: Its Current Role and Ongoing Challenges

Transcatheter aortic valve replacement (TAVR) has emerged as a viable alternative to surgical aortic valve replacement, as accumulating clinical evidence has demonstrated its safety and efficacy. TAVR indications have expanded beyond high-risk or inoperable patients to include intermediate and low-risk patients with severe aortic stenosis. Artificial intelligence (AI) is revolutionizing the field of cardiology, aiding in the interpretation of medical imaging and developing risk models for at-risk individuals and those with cardiac disease. This article explores the growing role of AI in TAVR procedures and assesses its potential impact, with particular focus on its ability to improve patient selection, procedural planning, post-implantation monitoring and contribute to optimized patient outcomes. In addition, current challenges and future directions in AI implementation are highlighted.

Current progress toward isogeometric modeling of the heart biophysics

In this paper, we review a powerful methodology to solve complex numerical simulations, known as isogeometric analysis, with a focus on applications to the biophysical modeling of the heart. We focus on the hemodynamics, modeling of the valves, cardiac tissue mechanics, and on the simulation of medical devices and treatments. For every topic, we provide an overview of the methods employed to solve the specific numerical issue entailed by the simulation. We try to cover the complete process, starting from the creation of the geometrical model up to the analysis and post-processing, highlighting the advantages and disadvantages of the methodology.

A contact model based on the coefficient of restitution for simulations of bio-prosthetic heart valves

A new general contact model is proposed for preventing inter-leaflet penetration of bio-prosthetic heart valves (BHV) at the end of the systole, which has the advantage of applying kinematic constraints directly and creating smooth free edges. At the end of each time step, the impenetrability constraints and momentum exchange between the impacting bodies are applied separately based on the coefficient of restitution. The contact method is implemented in a rotation-free, large deformation, and thin shell finite-element (FE) framework based on loop's subdivision surfaces. A nonlinear, anisotropic material model for a BHV is employed which uses Fung-elastic constitutive laws for in-plane and bending responses, respectively. The contact model is verified and validated against several benchmark problems. For a BHV-specific validation, the computed strains on different regions of a BHV under constant pressure are compared with experimentally measured data. Finally, dynamic simulations of BHV under physiological pressure waveform are performed for symmetrical and asymmetrical fiber orientations incorporating the new contact model and compared with the penalty contact method. The proposed contact model provides the coaptation area of a functioning BHV during the closing phase for both of the fiber orientations. Our results show that fiber orientation affects the dynamic of leaflets during the opening and closing phases. A swirling motion for the BHV with asymmetrical fiber orientation is observed, similar to experimental data. To include the fluid effects, fluid-structure interaction (FSI) simulation of the BHV is performed and compared to the dynamic results.

Latest Developments in Adapting Deep Learning for Assessing TAVR Procedures and Outcomes

Aortic valve defects are among the most prevalent clinical conditions. A severely damaged or non-functioning aortic valve is commonly replaced with a bioprosthetic heart valve (BHV) via the transcatheter aortic valve replacement (TAVR) procedure. Accurate pre-operative planning is crucial for a successful TAVR outcome. Assessment of computational fluid dynamics (CFD), finite element analysis (FEA), and fluid-solid interaction (FSI) analysis offer a solution that has been increasingly utilized to evaluate BHV mechanics and dynamics. However, the high computational costs and the complex operation of computational modeling hinder its application. Recent advancements in the deep learning (DL) domain can offer a real-time surrogate that can render hemodynamic parameters in a few seconds, thus guiding clinicians to select the optimal treatment option. Herein, we provide a comprehensive review of classical computational modeling approaches, medical imaging, and DL approaches for planning and outcome assessment of TAVR. Particularly, we focus on DL approaches in previous studies, highlighting the utilized datasets, deployed DL models, and achieved results. We emphasize the critical challenges and recommend several future directions for innovative researchers to tackle. Finally, an end-to-end smart DL framework is outlined for real-time assessment and recommendation of the best BHV design for TAVR. Ultimately, deploying such a framework in future studies will support clinicians in minimizing risks during TAVR therapy planning and will help in improving patient care.

A review of numerical simulation in transcatheter aortic valve replacement decision optimization

BACKGROUND

Recent trials indicated a further expansion of clinical indication of transcatheter aortic valve replacement to younger and low-risk patients. Factors related to longer-term complications are becoming more important for use in these patients. Accumulating evidence indicates that numerical simulation plays a significant role in improving the outcome of transcatheter aortic valve replacement. Understanding mechanical features' magnitude, pattern, and duration is a topic of ongoing relevance.

METHODS

We searched the PubMed database using keywords such as "transcatheter aortic valve replacement" and "numerical simulation" and reviewed and summarized relevant literature.

FINDINGS

This review integrated recently published evidence into three subtopics: 1) prediction of transcatheter aortic valve replacement outcomes through numerical simulation, 2) implications for surgeons, and 3) trends in transcatheter aortic valve replacement numerical simulation.

INTERPRETATIONS

Our study offers a comprehensive overview of the utilization of numerical simulation in the context of transcatheter aortic valve replacement, and highlights the advantages, potential challenges from a clinical standpoint. The convergence of medicine and engineering plays a pivotal role in enhancing the outcomes of transcatheter aortic valve replacement. Numerical simulation has provided evidence of potential utility for tailored treatments.

An optical coherence tomography and endothelial shear stress study of a novel bioresorbable bypass graft

Endothelial shear stress (ESS) plays a key role in the clinical outcomes in native and stented segments; however, their implications in bypass grafts and especially in a synthetic biorestorative coronary artery bypass graft are yet unclear. This report aims to examine the interplay between ESS and the morphological alterations of a biorestorative coronary bypass graft in an animal model. Computational fluid dynamics (CFD) simulation derived from the fusion of angiography and optical coherence tomography (OCT) imaging was used to reconstruct data on the luminal anatomy of a bioresorbable coronary bypass graft with an endoluminal "flap" identified during OCT acquisition. The "flap" compromised the smooth lumen surface and considerably disturbed the local flow, leading to abnormally low ESS and high oscillatory shear stress (OSI) in the vicinity of the "flap". In the presence of the catheter, the flow is more stable (median OSI 0.02384 versus 0.02635, p < 0.0001; maximum OSI 0.4612 versus 0.4837). Conversely, OSI increased as the catheter was withdrawn which can potentially cause back-and-forth motions of the "flap", triggering tissue fatigue failure. CFD analysis in this report provided sophisticated physiological information that complements the anatomic assessment from imaging enabling a complete understanding of biorestorative graft pathophysiology.

Patient-Specific Immersed Finite Element-Difference Model of Transcatheter Aortic Valve Replacement

Transcatheter aortic valve replacement (TAVR) first received FDA approval for high-risk surgical patients in 2011 and has been approved for low-risk surgical patients since 2019. It is now the most common type of aortic valve replacement, and its use continues to accelerate. Computer modeling and simulation (CM&S) is a tool to aid in TAVR device design, regulatory approval, and indication in patient-specific care. This study introduces a computational fluid-structure interaction (FSI) model of TAVR with Medtronic's CoreValve Evolut R device using the immersed finite element-difference (IFED) method. We perform dynamic simulations of crimping and deployment of the Evolut R, as well as device behavior across the cardiac cycle in a patient-specific aortic root anatomy reconstructed from computed tomography (CT) image data. These IFED simulations, which incorporate biomechanics models fit to experimental tensile test data, automatically capture the contact within the device and between the self-expanding stent and native anatomy. Further, we apply realistic driving and loading conditions based on clinical measurements of human ventricular and aortic pressures and flow rates to demonstrate that our Evolut R model supports a physiological diastolic pressure load and provides informative clinical performance predictions.

Biomechanics of Transcatheter Aortic Valve Replacement Complications and Computational Predictive Modeling

Transcatheter aortic valve replacement (TAVR) is a rapidly growing field enabling replacement of diseased aortic valves without the need for open heart surgery. However, due to the nature of the procedure and nonremoval of the diseased tissue, there are rates of complications ranging from tissue rupture and coronary obstruction to paravalvular leak, valve thrombosis, and permanent pacemaker implantation. In recent years, computational modeling has shown a great deal of promise in its capabilities to understand the biomechanical implications of TAVR as well as help preoperatively predict risks inherent to device-patient-specific anatomy biomechanical interaction. This includes intricate replication of stent and leaflet designs and tested and validated simulated deployments with structural and fluid mechanical simulations. This review outlines current biomechanical understanding of device-related complications from TAVR and related predictive strategies using computational modeling. An outlook on future modeling strategies highlighting reduced order modeling which could significantly reduce the high time and cost that are required for computational prediction of TAVR outcomes is presented in this review paper. A summary of current commercial/in-development software is presented in the final section.

Fluid-Structure Interaction Analysis on the Influence of the Aortic Valve Stent Leaflet Structure in Hemodynamics

Transcatheter aortic valve replacement (TAVR) is a minimally invasive surgical treatment for heart valve disease. At present, personalized TAVR valves are not available for some patients. This study adopts the fluid-structure interaction (FSI) model of the research object that has a three-disc leaflet form and structural design in the valve leaflet area. The valve opening shape, orifice area, stress-strain, and distribution of hemodynamic flow and pressure were compared under the condition of equal contact area between valve and blood. The FSI method was used to simulate the complex three dimensional characteristics of the flow field more accurately around the valve after TAVR stent implantation. Three personalized stent systems were established to study the performance of the leaflet design based on computational fluid dynamics. By comparing the different leaflet geometries, the maximum stress on leaflets and stents of model B was relatively reduced, which effectively improved the reliability of the stent design. Such valve design also causes the opening area of the valve leaflet to increase and the low-velocity area of the flow field to decrease during the working process of the valve, thus reducing the possibility of thrombosis. These findings can underpin breakthroughs in product design, and provide important theoretical support and technical guidance for clinical research.

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.

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.

Finite element analysis of transcatheter aortic valve implantation: Insights on the modelling of self-expandable devices

Computational simulations of Transcatheter Aortic Valve Implantation (TAVI) have reached a high level of complexity and accuracy for the prediction of possible implantation scenarios during the decision-making process. However, when focusing on the prosthetic device, currently different devices are available on the market which not only have different geometries, but also different material properties. The present work focuses on the calibration of Nitinol constitutive parameters of four self-expandable devices starting from experimental radial force tests on the prosthetic samples. Beside providing optimal material properties for each specific device, we also perform a patient-specific simulation, comparing the results obtained using both "literature" and calibrated parameters with the aim of investigating the impact of metallic frame parameters choice on simulation results.

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.

Thinner biological tissues induce leaflet flutter in aortic heart valve replacements

Valvular heart disease has recently become an increasing public health concern due to the high prevalence of valve degeneration in aging populations. For patients with severely impacted aortic valves that require replacement, catheter-based bioprosthetic valve deployment offers a minimally invasive treatment option that eliminates many of the risks associated with surgical valve replacement. Although recent percutaneous device advancements have incorporated thinner, more flexible biological tissues to streamline safer deployment through catheters, the impact of such tissues in the complex, mechanically demanding, and highly dynamic valvular system remains poorly understood. The present work utilized a validated computational fluid-structure interaction approach to isolate the behavior of thinner, more compliant aortic valve tissues in a physiologically realistic system. This computational study identified and quantified significant leaflet flutter induced by the use of thinner tissues that initiated blood flow disturbances and oscillatory leaflet strains. The aortic flow and valvular dynamics associated with these thinner valvular tissues have not been previously identified and provide essential information that can significantly advance fundamental knowledge about the cardiac system and support future medical device innovation. Considering the risks associated with such observed flutter phenomena, including blood damage and accelerated leaflet deterioration, this study demonstrates the potentially serious impact of introducing thinner, more flexible tissues into the cardiac system.

A Deep Learning Framework for Design and Analysis of Surgical Bioprosthetic Heart Valves

Bioprosthetic heart valves (BHVs) are commonly used as heart valve replacements but they are prone to fatigue failure; estimating their remaining life directly from medical images is difficult. Analyzing the valve performance can provide better guidance for personalized valve design. However, such analyses are often computationally intensive. In this work, we introduce the concept of deep learning (DL) based finite element analysis (DLFEA) to learn the deformation biomechanics of bioprosthetic aortic valves directly from simulations. The proposed DL framework can eliminate the time-consuming biomechanics simulations, while predicting valve deformations with the same fidelity. We present statistical results that demonstrate the high performance of the DLFEA framework and the applicability of the framework to predict bioprosthetic aortic valve deformations. With further development, such a tool can provide fast decision support for designing surgical bioprosthetic aortic valves. Ultimately, this framework could be extended to other BHVs and improve patient care.