A review of state-of-the-art numerical methods for simulating flow through mechanical heart valves

Biomechanical Analysis of Age-Dependent Changes in Fontan Power Loss

The hemodynamics in Fontan patients with single ventricles rely on favorable flow and energetics, especially in the absence of a subpulmonary ventricle. Age-related changes in energetics for extracardiac and lateral tunnel Fontan procedures are not well understood. Vorticity (VOR) and viscous dissipation rate (VDR) are two descriptors that can provide insights into flow dynamics and dissipative areas in Fontan pathways, potentially contributing to power loss. This study examined power loss and its correlation with spatio-temporal flow descriptors (vorticity and VDR). Data from 414 Fontan patients were used to establish a relationship between the superior vena cava (SVC) to inferior vena cava (IVC) flow ratio and age. Computational flow modeling was conducted for both extracardiac conduits (ECC, n = 16) and lateral tunnels (LT, n = 25) at different caval inflow ratios of 2, 1, and 0.5 that corresponded with ages 3, 8, and 15+. In both cohorts, vorticity and VDR correlated well with PL, but ECC cohort exhibited a slightly stronger correlation for PL-VOR (>0.83) and PL-VDR (>0.89) than that for LT cohort (>0.76 and > 0.77, respectively) at all ages. Our data also suggested that absolute and indexed PL increase (p < 0.02) non-linearly as caval inflow changes with age and are highly patient-specific. Comparison of indexed power loss between our ECC and LT cohort showed that while ECC had a slightly higher median PL for all 3 caval inflow ratio examined (3.3, 8.3, 15.3) as opposed to (2.7, 7.6, 14.8), these differences were statistically non-significant. Lastly, there was a consistent rise in pressure gradient across the TCPC with age-related increase in IVC flows for both ECC and LT Fontan patient cohort. Our study provided hemodynamic insights into Fontan energetics and how they are impacted by age-dependent change in caval inflow. This workflow may help assess the long-term sustainability of the Fontan circulation and inform the design of more efficient Fontan conduits.

Fluid-structure interaction simulation of mechanical aortic valves: a narrative review exploring its role in total product life cycle

Over the last years computer modelling and simulation has emerged as an effective tool to support the total product life cycle of cardiovascular devices, particularly in the device preclinical evaluation and post-market assessment. Computational modelling is particularly relevant for heart valve prostheses, which require an extensive assessment of their hydrodynamic performance and of risks of hemolysis and thromboembolic complications associated with mechanically-induced blood damage. These biomechanical aspects are typically evaluated through a fluid-structure interaction (FSI) approach, which enables valve fluid dynamics evaluation accounting for leaflets movement. In this context, the present narrative review focuses on the computational modelling of bileaflet mechanical aortic valves through FSI approach, aiming to foster and guide the use of simulations in device total product life cycle. The state of the art of FSI simulation of heart valve prostheses is reviewed to highlight the variety of modelling strategies adopted in the literature. Furthermore, the integration of FSI simulations in the total product life cycle of bileaflet aortic valves is discussed, with particular emphasis on the role of simulations in complementing and potentially replacing the experimental tests suggested by international standards. Simulations credibility assessment is also discussed in the light of recently published guidelines, thus paving the way for a broader inclusion of in silico evidence in regulatory submissions. The present narrative review highlights that FSI simulations can be successfully framed within the total product life cycle of bileaflet mechanical aortic valves, emphasizing that credible in silico models evaluating the performance of implantable devices can (at least) partially replace preclinical in vitro experimentation and support post-market biomechanical evaluation, leading to a reduction in both time and cost required for device development.

Simulation of Mechanical Heart Valve Dysfunction and the Non-Newtonian Blood Model Approach

The mechanical heart valve (MHV) is commonly used for the treatment of cardiovascular diseases. Nonphysiological hemodynamic in the MHV may cause hemolysis, platelet activation, and an increased risk of thromboembolism. Thromboembolism may cause severe complications and valve dysfunction. This paper thoroughly reviewed the simulation of physical quantities (velocity distribution, vortex formation, and shear stress) in healthy and dysfunctional MHV and reviewed the non-Newtonian blood flow characteristics in MHV. In the MHV numerical study, the dysfunction will affect the simulation results, increase the pressure gradient and shear stress, and change the blood flow patterns, increasing the risks of hemolysis and platelet activation. The blood flow passes downstream and has obvious recirculation and stagnation region with the increased dysfunction severity. Due to the complex structure of the MHV, the non-Newtonian shear-thinning viscosity blood characteristics become apparent in MHV simulations. The comparative study between Newtonian and non-Newtonian always shows the difference. The shear-thinning blood viscosity model is the basics to build the blood, also the blood exhibiting viscoelastic properties. More details are needed to establish a complete and more realistic simulation.

Image-Based Computational Hemodynamics Analysis of Systolic Obstruction in Hypertrophic Cardiomyopathy

Hypertrophic Cardiomyopathy (HCM) is a pathological condition characterized by an abnormal thickening of the myocardium. When affecting the medio-basal portion of the septum, it is named Hypertrophic Obstructive Cardiomyopathy (HOCM) because it induces a flow obstruction in the left ventricular outflow tract. In any type of HCM, the myocardial function can become compromised, possibly resulting in cardiac death. In this study, we investigated with computational analysis the hemodynamics of patients with different types of HCM. The aim was quantifying the effects of this pathology on the intraventricular blood flow and pressure gradients, and providing information potentially useful to guide the indication and the modality of the surgical treatment (septal myectomy). We employed an image-based computational approach, integrating fluid dynamics simulations with geometric and functional data, reconstructed from standard cardiac cine-MRI acquisitions. We showed that with our approach we can better understand the patho-physiological behavior of intraventricular blood flow dynamics due to the abnormal morphological and functional aspect of the left ventricle. The main results of our investigation are: (a) a detailed patient-specific analysis of the blood velocity, pressure and stress distribution associated to HCM; (b) a computation-based classification of patients affected by HCM that can complement the current clinical guidelines for the diagnosis and treatment of HOCM.

Computational Analysis of Wall Shear Stress Patterns on Calcified and Bicuspid Aortic Valves: Focus on Radial and Coaptation Patterns

Calcification and bicuspid valve formation are important aortic valve disorders that disturb the hemodynamics and the valve function. The detailed analysis of aortic valve hemodynamics would lead to a better understanding of the disease’s etiology. We computationally modeled the aortic valve using simplified three-dimensional geometry and inlet velocity conditions obtained via echocardiography. We examined various calcification severities and bicuspid valve formation. Fluid-structure interaction (FSI) analyses were adapted using ANSYS Workbench to incorporate both flow dynamics and leaflet deformation accurately. Simulation results were validated by comparing leaflet movements in B-mode echo recordings. Results indicate that the biomechanical environment is significantly changed for calcified and bicuspid valves. High flow jet velocities are observed in the calcified valves which results in high transvalvular pressure difference (TPG). Wall shear stresses (WSS) increased with the calcification on both fibrosa (aorta side) and ventricularis (left ventricle side) surfaces of the leaflet. The WSS distribution is regular on the ventricularis, as the WSS values proportionally increase from the base to the tip of the leaflet. However, WSS patterns are spatially complex on the fibrosa side. Low WSS levels and spatially complex WSS patterns on the fibrosa side are considered as promoting factors for further calcification and valvular diseases.

Application of physics-based flow models in cardiovascular medicine: Current practices and challenges

Personalized physics-based flow models are becoming increasingly important in cardiovascular medicine. They are a powerful complement to traditional methods of clinical decision-making and offer a wealth of physiological information beyond conventional anatomic viewing using medical imaging data. These models have been used to identify key hemodynamic biomarkers, such as pressure gradient and wall shear stress, which are associated with determining the functional severity of cardiovascular diseases. Importantly, simulation-driven diagnostics can help researchers understand the complex interplay between geometric and fluid dynamic parameters, which can ultimately improve patient outcomes and treatment planning. The possibility to compute and predict diagnostic variables and hemodynamics biomarkers can therefore play a pivotal role in reducing adverse treatment outcomes and accelerate development of novel strategies for cardiovascular disease management.

The Apex bileaflet mechanical heart valve

Mechanical Heart Valves (MHVs) are known for their excellent lifespan and feasibly are the most reliable and stable valves amongst all prosthetic valves. Successful bileaflet MHVs such as the St. Jude Medical (SJM) are known for providing central blood flow and minimal pressure drop across the valve. However, due to their non-physiological flow conditions, they still suffer from hemodynamic complications, that is, red blood cell (RBC) lysis and/or thrombogenicity, to date. Our hypothesis is that the design of MHVs can be improved so that their hemodynamics can be comparable to those of tissue valves. In this study, a new concept for the design of MHVs is proposed. To accomplish this, we identified the major design limitations of bileaflet MHVs, such as the gold standard SJM valve as well as the believed contributing factors to their thrombogenicity. We developed a novel design architecture for bileaflet MHVs that addressed these limitations, and from it, the Apex Valve (AV). Our experimental assessment of the AV found that its hemodynamics were closer to that of a bioprosthetic valve than of a bileaflet MHV. This design has been filed as a US Provisional Patent.

Immersed boundary method for multiphase transport phenomena

Multiphase flows with momentum, heat, and mass transfer exist widely in a variety of industrial applications. With the rapid development of numerical algorithms and computer capacity, advanced numerical simulation has become a promising tool in investigating multiphase transport problems. Immersed boundary (IB) method has recently emerged as such a popular interface capturing method for efficient simulations of multiphase flows, and significant achievements have been obtained. In this review, we attempt to give an overview of recent progresses on IB method for multiphase transport phenomena. Firstly, the governing equations, the basic ideas, and different boundary conditions for the IB methods are introduced. This is followed by numerical strategies, from which the IB methods are classified into two types, namely the artificial boundary method and the authentic boundary method. Discussions on the implementation of various boundary conditions at the interphase surface with momentum, heat, and mass transfer for different IB methods are then presented, together with a summary. Then, the state-of-the-art applications of IB methods to multiphase flows, including the isothermal flows, the heat transfer flows, and the mass transfer problems are outlined, with particular emphasis on the latter two topics. Finally, the conclusions and future challenges are identified.

Next-generation tissue-engineered heart valves with repair, remodelling and regeneration capacity

Valvular heart disease is a major cause of morbidity and mortality worldwide. Surgical valve repair or replacement has been the standard of care for patients with valvular heart disease for many decades, but transcatheter heart valve therapy has revolutionized the field in the past 15 years. However, despite the tremendous technical evolution of transcatheter heart valves, to date, the clinically available heart valve prostheses for surgical and transcatheter replacement have considerable limitations. The design of next-generation tissue-engineered heart valves (TEHVs) with repair, remodelling and regenerative capacity can address these limitations, and TEHVs could become a promising therapeutic alternative for patients with valvular disease. In this Review, we present a comprehensive overview of current clinically adopted heart valve replacement options, with a focus on transcatheter prostheses. We discuss the various concepts of heart valve tissue engineering underlying the design of next-generation TEHVs, focusing on off-the-shelf technologies. We also summarize the latest preclinical and clinical evidence for the use of these TEHVs and describe the current scientific, regulatory and clinical challenges associated with the safe and broad clinical translation of this technology.

Computational fluid dynamics: a primer for congenital heart disease clinicians

Computational fluid dynamics has become an important tool for studying blood flow dynamics. As an in-silico collection of methods, computational fluid dynamics is noninvasive and provides numerical values for the most important parameters of blood flow, such as velocity and pressure that are crucial in hemodynamic studies. In this primer, we briefly explain the basic theory and workflow of the two most commonly applied computational fluid dynamics techniques used in the congenital heart disease literature: the finite element method and the finite volume method. We define important terminology and include specific examples of how using these methods can answer important clinical questions in congenital cardiac surgery planning and perioperative patient management.

An image-based computational hemodynamics study of the Systolic Anterior Motion of the mitral valve

Systolic Anterior Motion (SAM) of the mitral valve - often associated with Hypertrophic Obstructive Cardiomyopathy (HOCM) - is a cardiac pathology in which a functional subaortic stenosis is induced during systole by the mitral leaflets partially obstructing the outflow tract of the left ventricle. Its assessment by diagnostic tests is often difficult, possibly underestimating its severity and thus increasing the risk of heart failure. In this paper, we propose a new computational pipeline, based on cardiac cine Magnetic Resonance Imaging (cine-MRI) data, for the assessment of SAM. The pipeline encompasses image processing of the left ventricle and the mitral valve, and numerical investigation of cardiac hemodynamics by means of Computational Fluid Dynamics (CFD) in a moving domain with image-based prescribed displacement. Patient-specific geometry and motion of the left ventricle are considered in view of an Arbitrary Lagrangian-Eulerian approach for CFD, while the reconstructed mitral valve is immersed in the computational domain by means of a resistive method. We assess clinically relevant flow and pressure indicators in a parametric study for different degrees of SAM severity, in order to provide a better quantitative evaluation of the pathological condition. Moreover, we provide specific indications for its possible surgical treatment, i.e. septal myectomy.

Antibacterial Coatings for Improving the Performance of Biomaterials

Biomedical devices have become essential in the health care. Every day, an enormous number of these devices are used or implanted in humans. In this context, the bacterial contamination that could be developed in implanted devices is critical since it is estimated that infections kill more people than other medical causes. Commonly, these infections are treated with antibiotics, but the biofilm formation on implant surfaces could significantly reduce the effectiveness of these antibiotics since bacteria inside the biofilm is protected from the drug. In some cases, a complete removal of the implant is necessary in order to overcome the infection. In this context, antibacterial coatings are considered an excellent strategy to avoid biofilm formation and, therefore, mitigate the derived complications. In this review, the main biomaterials used in biomedical devices, the mechanism of biofilm formation, and the main strategies for the development of antibacterial coatings, are reviewed. Finally, the main polymer-based strategies to develop antibacterial coatings are summarized, with the aim of these coatings being to avoid the bacteria proliferation by controlling the antibacterial mechanisms involved and enhancing long-term stability.

Effect of Aortic Wall Deformation with Healthy and Calcified Annulus on Hemodynamic Performance of Implanted On-X Valve

INTRODUCTION

In this research, the hemodynamic performance of a 23-mm On-X bileaflet mechanical heart valve (BMHV) was investigated with the realistic geometry model of the valve and the deformable aorta in accelerating systole. In addition, the effect of ascending aorta flexibility and aortic annulus calcification on the complex blood flow characteristics were investigated.

METHODS

The geometry of the aorta is derived from the medical images, and the Ogden model has been utilized for the mechanical behavior of the ascending aorta. The 3D numerical simulation by a two-way Fluid-Structure Interaction (FSI) analysis using the Arbitrary Lagrangian-Eulerian (ALE) method was performed throughout the accelerating systolic phase.

RESULTS

The dynamics of the leaflets are investigated, and blood flow characteristics such as velocities, vorticities as well as viscous and turbulent shear stress were precisely captured in the flow domain specifically in the hinge region. Streamline results are in accordance with the previously reported data, which show that the flared On-X valves inlet yields a more uniform flow in accelerating systole. Simulations show that aorta flexibility or valve annulus calcification causes variations up to 7% in maximum fluid velocity and 20% in Turbulence Kinetic Energy (TKE).

CONCLUSIONS

In this study, the complex flow field characteristics in the new generation of BMHVs considering aorta flexibility with healthy and calcified annulus were investigated. It was found that the blood flow around the hinges region is in the danger of hemolysis and platelet activation and subsequently thromboembolism. Furthermore, the results show that similar to vessel wall deformation, considering the probable annulus calcification after valve replacement is also essential.

Some Effects of Different Constitutive Laws on FSI Simulation for the Mitral Valve

In this paper, three different constitutive laws for mitral leaflets and two laws for chordae tendineae are selected to study their effects on mitral valve dynamics with fluid-structure interaction. We first fit these three mitral leaflet constitutive laws and two chordae tendineae laws with experimental data. The fluid-structure interaction is implemented in an immersed boundary framework with finite element extension for solid, that is the hybrid immersed boundary/finite element(IB/FE) method. We specifically compare the fluid-structure results of different constitutive laws since fluid-structure interaction is the physiological loading environment. This allows us to look at the peak jet velocity, the closure regurgitation volume, and the orifice area. Our numerical results show that different constitutive laws can affect mitral valve dynamics, such as the transvalvular flow rate, closure regurgitation and the orifice area, while the differences in fiber strain and stress are insignificant because all leaflet constitutive laws are fitted to the same set of experimental data. In addition, when an exponential constitutive law of chordae tendineae is used, a lower closure regurgitation flow is observed compared to that of a linear material model. In conclusion, combining numerical dynamic simulations and static experimental tests, we are able to identify suitable constitutive laws for dynamic behaviour of mitral leaflets and chordae under physiological conditions.

PROPOSED OVAL ST. JUDE MEDICAL VALVE: EFFECT OF HEART RATE

In this study, the hemodynamic performance of the conventional St. Jude Medical (SJM) valve and our proposed design known as the oval SJM valve are studied and compared. These studies are based on a wide range of physiological heart rates, i.e., 70–130[Formula: see text]bpm, in the opening phase. We designed and developed a precise computational platform to assess the hemodynamics of bileaflet mechanical heart valves for laminar and turbulent regimes. Also, as one of the fundamental changes applied to the conventional SJM vales, the housing is considered oval similar to oval shape of annulus. Results clearly indicate hemodynamic improvements in the proposed design over the SJM valve. The improvements are characterized by lower shear stress and wall shear stress distributions around the valve and leaflets, and lower valve pressure drop compared to that of the conventional SJM model. The proposed design shows potential and merits additional development.

Modelling mitral valvular dynamics-current trend and future directions

Dysfunction of mitral valve causes morbidity and premature mortality and remains a leading medical problem worldwide. Computational modelling aims to understand the biomechanics of human mitral valve and could lead to the development of new treatment, prevention and diagnosis of mitral valve diseases. Compared with the aortic valve, the mitral valve has been much less studied owing to its highly complex structure and strong interaction with the blood flow and the ventricles. However, the interest in mitral valve modelling is growing, and the sophistication level is increasing with the advanced development of computational technology and imaging tools. This review summarises the state-of-the-art modelling of the mitral valve, including static and dynamics models, models with fluid-structure interaction, and models with the left ventricle interaction. Challenges and future directions are also discussed.

Non-Newtonian Blood Flow Simulation of Diastolic Phase in Bileaflet Mechanical Heart Valve Implanted in a Realistic Aortic Root Containing Coronary Arteries

Coronary arteries, which are branched from the sinuses, have tangible effects on the hemodynamic performance of the bileaflet mechanical heart valve (BMHV), especially in the diastolic phase. To better understand this issue, a computer model of ascending aorta including realistic sinus shapes and coronary arteries has been generated in this study in order to investigate the BMHV performance during diastole. Three-dimensional transient numerical analysis is conducted to simulate the diastolic blood flow through the hinges and in coronary arteries under the assumption of non-Newtonian behavior. Results indicate that as blood flows to the coronary arteries mainly during diastole, leakage flow from the hinge and other gaps will change considering the influence of coronary arteries. In addition, BMHV in the case of aortic replacement will increase blood flow rate into the coronary arteries about 100% as the mechanical valve resistance is higher than a native heart valve. Also, it will change the wall shear stress (WSS) distribution and increase coronary artery disease (CAD) potential. It is found out that although less leakage flow reduces the velocity magnitudes through the gaps, the shear stress acting on blood elements with non-Newtonian assumption will be detrimental in the hinge corner at the ventricular side. High WSS of 1800 Pa is observed at beginning of diastole at this region.

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

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

Micro and nanotechnologies in heart valve tissue engineering

Due to the increased morbidity and mortality resulting from heart valve diseases, there is a growing demand for off-the-shelf implantable tissue engineered heart valves (TEHVs). Despite the significant progress in recent years in improving the design and performance of TEHV constructs, viable and functional human implantable TEHV constructs have remained elusive. The recent advances in micro and nanoscale technologies including the microfabrication, nano-microfiber based scaffolds preparation, 3D cell encapsulated hydrogels preparation, microfluidic, micro-bioreactors, nano-microscale biosensors as well as the computational methods and models for simulation of biological tissues have increased the potential for realizing viable, functional and implantable TEHV constructs. In this review, we aim to present an overview of the importance and recent advances in micro and nano-scale technologies for the development of TEHV constructs.

On-X Heart Valve Prosthesis: Numerical Simulation of Hemodynamic Performance in Accelerating Systole

Numerical simulation of the bileaflet mechanical heart valves (BMHVs) has been of interest for many researchers due to its capability of predicting hemodynamic performance. A lot of studies have tried to simulate this three-dimensional complex flow in order to analyze the effect of different valve designs on the blood flow pattern. However, simplified models and prescribed motion for the leaflets were utilized. In this paper, transient complex blood flow in the location of ascending aorta has been investigated in a realistic model by fully coupled simulation. Geometry model for the aorta and the replaced valve is constructed based on the medical images and extracted point clouds. A 23-mm On-X Medical BMHV as the new generation design has been selected for the flow field analysis. The two-way coupling simulation is conducted throughout the accelerating phase in order to obtain valve dynamics in the opening process. The complex flow field in the hinge recess is captured precisely for all leaflet positions and recirculating zones and elevated shear stress areas have been observed. Results indicate that On-X valve yields relatively less transvalvular pressure gradient which would lower cardiac external work. Furthermore, converging inlet leads to a more uniform flow and consequently less turbulent eddies. However, the leaflets cannot open fully due to middle diffuser-shaped orifice. In addition, asymmetric butterfly-shaped hinge design and converging orifice leads to better hemodynamic performance. With the help of two-way fluid solid interaction simulation, leaflet angle follows the experimental trends more precisely rather than the prescribed motion in previous 3D simulations.

Transport physics and biorheology in the setting of hemostasis and thrombosis

The biophysics of blood flow can dictate the function of molecules and cells in the vasculature with consequent effects on hemostasis, thrombosis, embolism, and fibrinolysis. Flow and transport dynamics are distinct for (i) hemostasis vs. thrombosis and (ii) venous vs. arterial episodes. Intraclot transport changes dramatically the moment hemostasis is achieved or the moment a thrombus becomes fully occlusive. With platelet concentrations that are 50- to 200-fold greater than platelet-rich plasma, clots formed under flow have a different composition and structure compared with blood clotted statically in a tube. The platelet-rich, core/shell architecture is a prominent feature of self-limiting hemostatic clots formed under flow. Importantly, a critical threshold concentration of surface tissue factor is required for fibrin generation under flow. Once initiated by wall-derived tissue factor, thrombin generation and its spatial propagation within a clot can be modulated by γ'-fibrinogen incorporated into fibrin, engageability of activated factor (FIXa)/activated FVIIIa tenase within the clot, platelet-derived polyphosphate, transclot permeation, and reduction of porosity via platelet retraction. Fibrin imparts tremendous strength to a thrombus to resist embolism up to wall shear stresses of 2400 dyne cm(-2) . Extreme flows, as found in severe vessel stenosis or in mechanical assist devices, can cause von Willebrand factor self-association into massive fibers along with shear-induced platelet activation. Pathological von Willebrand factor fibers are A Disintegrin And Metalloprotease with ThromboSpondin-1 domain 13 resistant but are a substrate for fibrin generation due to FXIIa capture. Recently, microfluidic technologies have enhanced the ability to interrogate blood in the context of stenotic flows, acquired von Willebrand disease, hemophilia, traumatic bleeding, and drug action.

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.

Effect of heart rate on the hemodynamics of bileaflet mechanical heart valves’ prostheses (St. Jude Medical) in the aortic position and in the opening phase: A computational study

To date, to the best of the authors’ knowledge, in almost all of the studies performed around the hemodynamics of bileaflet mechanical heart valves, a heart rate of 70–72 beats/min has been considered. In fact, the heart rate of ~72 beats/min does not represent the entire normal physiological conditions under which the aortic or prosthetic valves function. The heart rates of 120 or 50 beats/min may lead to hemodynamic complications, such as plaque formation and/or thromboembolism in patients. In this study, the hemodynamic performance of the bileaflet mechanical heart valves in a wide range of normal and physiological heart rates, that is, 60–150 beats/min, was studied in the opening phase. The model considered in this study was a St. Jude Medical bileaflet mechanical heart valve with the inner diameter of 27 mm in the aortic position. The hemodynamics of the native valve and the St. Jude Medical valve were studied in a variety of heart rates in the opening phase and the results were carefully compared. The results indicate that peak values of the velocity profile downstream of the valve increase as heart rate increases, as well as the location of the maximum velocity changes with heart rate in the St. Jude Medical valve model. Also, the maximum values of shear stress and wall shear stresses downstream of the valve are proportional to heart rate in both models. Interestingly, the maximum shear stress and wall shear stress values in both models are in the same range when heart rate is <90 beats/min; however, these values significantly increase in the St. Jude Medical valve model when heart rate is >90 beats/min (up to ~40% growth compared to that of the native valve). The findings of this study may be of importance in the hemodynamic performance of bileaflet mechanical heart valves. They may also play an important role in design improvement of conventional prosthetic heart valves and the design of the next generation of prosthetic valves, such as percutaneous valves.

Hemodynamic study of the elliptic St. Jude Medical valve: A computational study

Despite successful implantation of St. Jude Medical bileaflet mechanical heart valves, red blood cell lysis and thrombogenic complications associated with these types of valves are yet to be addressed. In our previous study, we proposed an elliptic housing where 10% ovality was applied to the housing of St. Jude Medical valves. Our preliminary results suggested that the overall hemodynamic performance of St. Jude Medical valves improved in both the closing and opening phases. In this study, we evaluated the hemodynamics around the leaflets in the opening phase using a more sophisticated computational platform, computational fluid dynamics. Results suggested both lower shear stress and wall shear stress values and an overall improved hemodynamic performance in the proposed design. This improvement is characterized by lower values of shear stress and wall shear stress in the regions downstream of the leaflets, lower pressure drop across the valve and smaller recirculation zones in the sinuses areas. The proposed design may open a new chapter in the concept of design and hemodynamic improvement of the next generation of mechanical heart valves.

Computational fluid dynamics modelling in cardiovascular medicine

This paper reviews the methods, benefits and challenges associated with the adoption and translation of computational fluid dynamics (CFD) modelling within cardiovascular medicine. CFD, a specialist area of mathematics and a branch of fluid mechanics, is used routinely in a diverse range of safety-critical engineering systems, which increasingly is being applied to the cardiovascular system. By facilitating rapid, economical, low-risk prototyping, CFD modelling has already revolutionised research and development of devices such as stents, valve prostheses, and ventricular assist devices. Combined with cardiovascular imaging, CFD simulation enables detailed characterisation of complex physiological pressure and flow fields and the computation of metrics which cannot be directly measured, for example, wall shear stress. CFD models are now being translated into clinical tools for physicians to use across the spectrum of coronary, valvular, congenital, myocardial and peripheral vascular diseases. CFD modelling is apposite for minimally-invasive patient assessment. Patient-specific (incorporating data unique to the individual) and multi-scale (combining models of different length- and time-scales) modelling enables individualised risk prediction and virtual treatment planning. This represents a significant departure from traditional dependence upon registry-based, population-averaged data. Model integration is progressively moving towards ‘digital patient’ or ‘virtual physiological human’ representations. When combined with population-scale numerical models, these models have the potential to reduce the cost, time and risk associated with clinical trials. The adoption of CFD modelling signals a new era in cardiovascular medicine. While potentially highly beneficial, a number of academic and commercial groups are addressing the associated methodological, regulatory, education- and service-related challenges.

The Role of Shape and Heart Rate on the Performance of the Left Ventricle

The left ventricle function is to pump the oxygenated blood through the circulatory system. Ejection fraction is the main noninvasive parameter for detecting heart disease (healthy >55%), and it is thought to be the main parameter affecting efficiency. However, the effects of other parameters on efficiency have yet to be investigated. We investigate the effect of heart rate and left ventricle shape by carrying out 3D numerical simulations of a left ventricle at different heart rates and perturbed geometries under constant, normal ejection fraction. The simulation using the immersed boundary method provide the 3D flow and pressure fields, which enable direct calculation of a new hemodynamic efficiency (H-efficiency) parameter, which does not depend on any reference pressure. The H-efficiency is defined as the ratio of flux of kinetic energy (useful power) to the total cardiac power into the left ventricle control volume. Our simulations show that H-efficiency is not that sensitive to heart rate but is maximized at around normal heart rate (72 bpm). Nevertheless, it is more sensitive to the shape of the left ventricle, which affects the H-efficiency by as much as 15% under constant ejection fraction.

Numerical Modeling of Intraventricular Flow during Diastole after Implantation of BMHV

This work presents a numerical simulation of intraventricular flow after the implantation of a bileaflet mechanical heart valve at the mitral position. The left ventricle was simplified conceptually as a truncated prolate spheroid and its motion was prescribed based on that of a healthy subject. The rigid leaflet rotation was driven by the transmitral flow and hence the leaflet dynamics were solved using fluid-structure interaction approach. The simulation results showed that the bileaflet mechanical heart valve at the mitral position behaved similarly to that at the aortic position. Sudden area expansion near the aortic root initiated a clockwise anterior vortex, and the continuous injection of flow through the orifice resulted in further growth of the anterior vortex during diastole, which dominated the intraventricular flow. This flow feature is beneficial to preserving the flow momentum and redirecting the blood flow towards the aortic valve. To the best of our knowledge, this is the first attempt to numerically model intraventricular flow with the mechanical heart valve incorporated at the mitral position using a fluid-structure interaction approach. This study facilitates future patient-specific studies.

Micro-scale blood particulate dynamics using a non-uniform rational B-spline-based isogeometric analysis

The current research presents a novel method in which blood particulates - biconcave red blood cells (RBCs) and spherical cells are modeled using isogeometric analysis, specifically Non-Uniform Rational B-Splines (NURBS) in 3-D. The use of NURBS ensures that even with a coarse representation, the geometry of the blood particulates maintains an accurate description when subjected to large deformations. The fundamental advantage of this method is the coupling of the geometrical description and the stress analysis of the cell membrane into a single, unified framework. Details on the modeling approach, implementation of boundary conditions and the membrane mechanics analysis using isogeometric modeling are presented, along with validation cases for spherical and biconcave cells. Using NURBS - based isogeometric analysis, the behavior of individual cells in fluid flow is presented and analyzed in different flow regimes using as few as 176 elements for a spherical cell and 220 elements for a biconcave RBC. This work provides a framework for modeling a large number of 3-D deformable biological cells, each with its own geometric description and membrane properties. To the best knowledge of the authors, this is the first application of the NURBS - based isogeometric analysis to model and simulate blood particulates in flow in 3D.

Validation of a numerical FSI simulation of an aortic BMHV by in vitro PIV experiments

In this paper, a validation of a recently developed fluid-structure interaction (FSI) coupling algorithm to simulate numerically the dynamics of an aortic bileaflet mechanical heart valve (BMHV) is performed. This validation is done by comparing the numerical simulation results with in vitro experiments. For the in vitro experiments, the leaflet kinematics and flow fields are obtained via the particle image velocimetry (PIV) technique. Subsequently, the same case is numerically simulated by the coupling algorithm and the resulting leaflet kinematics and flow fields are obtained. Finally, the results are compared, revealing great similarity in leaflet motion and flow fields between the numerical simulation and the experimental test. Therefore, it is concluded that the developed algorithm is able to capture very accurately all the major leaflet kinematics and dynamics and can be used to study and optimize the design of BMHVs.

Numerical simulation of patient-specific left ventricular model with both mitral and aortic valves by FSI approach

Intraventricular flow is important in understanding left ventricular function; however, relevant numerical simulations are limited, especially when heart valve function is taken into account. In this study, intraventricular flow in a patient-specific left ventricle has been modelled in two-dimension (2D) with both mitral and aortic valves integrated. The arbitrary Lagrangian-Eulerian (ALE) approach was employed to handle the large mesh deformation induced by the beating ventricular wall and moving leaflets. Ventricular wall deformation was predefined based on MRI data, while leaflet dynamics were predicted numerically by fluid-structure interaction (FSI). Comparisons of simulation results with in vitro and in vivo measurements reported in the literature demonstrated that numerical method in combination with MRI was able to predict qualitatively the patient-specific intraventricular flow. To the best of our knowledge, we are the first to simulate patient-specific ventricular flow taking into account both mitral and aortic valves.

Fluid-structure interaction of an aortic heart valve prosthesis driven by an animated anatomic left ventricle

We develop a novel large-scale kinematic model for animating the left ventricle (LV) wall and use this model to drive the fluid-structure interaction (FSI) between the ensuing blood flow and a mechanical heart valve prosthesis implanted in the aortic position of an anatomic LV/aorta configuration. The kinematic model is of lumped type and employs a cell-based, FitzHugh-Nagumo framework to simulate the motion of the LV wall in response to an excitation wavefront propagating along the heart wall. The emerging large-scale LV wall motion exhibits complex contractile mechanisms that include contraction (twist) and expansion (untwist). The kinematic model is shown to yield global LV motion parameters that are well within the physiologic range throughout the cardiac cycle. The FSI between the leaflets of the mechanical heart valve and the blood flow driven by the dynamic LV wall motion and mitral inflow is simulated using the curvilinear immersed boundary (CURVIB) method [1, 2] implemented in conjunction with a domain decomposition approach. The computed results show that the simulated flow patterns are in good qualitative agreement with in vivo observations. The simulations also reveal complex kinematics of the valve leaflets, thus, underscoring the need for patient-specific simulations of heart valve prosthesis and other cardiac devices.

Left ventricular flow analysis: recent advances in numerical methods and applications in cardiac ultrasound

The left ventricle (LV) pumps oxygenated blood from the lungs to the rest of the body through systemic circulation. The efficiency of such a pumping function is dependent on blood flow within the LV chamber. It is therefore crucial to accurately characterize LV hemodynamics. Improved understanding of LV hemodynamics is expected to provide important clinical diagnostic and prognostic information. We review the recent advances in numerical and experimental methods for characterizing LV flows and focus on analysis of intraventricular flow fields by echocardiographic particle image velocimetry (echo-PIV), due to its potential for broad and practical utility. Future research directions to advance patient-specific LV simulations include development of methods capable of resolving heart valves, higher temporal resolution, automated generation of three-dimensional (3D) geometry, and incorporating actual flow measurements into the numerical solution of the 3D cardiovascular fluid dynamics.

Toward patient-specific simulations of cardiac valves: state-of-the-art and future directions

Recent computational methods enabling patient-specific simulations of native and prosthetic heart valves are reviewed. Emphasis is placed on two critical components of such methods: (1) anatomically realistic finite element models for simulating the structural dynamics of heart valves; and (2) fluid structure interaction methods for simulating the performance of heart valves in a patient-specific beating left ventricle. It is shown that the significant progress achieved in both fronts paves the way toward clinically relevant computational models that can simulate the performance of a range of heart valves, native and prosthetic, in a patient-specific left heart environment. The significant algorithmic and model validation challenges that need to be tackled in the future to realize this goal are also discussed.

Evaluation of shear stress accumulation on blood components in normal and dysfunctional bileaflet mechanical heart valves using smoothed particle hydrodynamics

Evaluating shear induced hemodynamic complications is one of the major concerns in design of the mechanical heart valves (MHVs). The monitoring of these events relies on both numerical simulations and experimental measurements. Currently, numerical approaches are mainly based on a combined Eulerian-Lagrangian approach. A more straightforward evaluation can be based on the Lagrangian analysis of the whole blood. As a consequence, Lagrangian meshfree methods are more adapted to such evaluation. In this study, smoothed particle hydrodynamics (SPH), a fully meshfree particle method originated to simulate compressible astrophysical flows, is applied to study the flow through a normal and a dysfunctional bileaflet mechanical heart valves (BMHVs). The SPH results are compared with the reference data. The accumulation of shear stress patterns on blood components illustrates the important role played by non-physiological flow patterns and mainly vortical structures in this issue. The statistical distribution of particles with respect to shear stress loading history provides important information regarding the relative number of blood components that can be damaged. This can be used as a measure of the response of blood components to the presence of the valve implant or any implantable medical device. This work presents the first attempt to simulate pulsatile flow through BMHVs using SPH method.

A robust and efficient valve model based on resistive immersed surfaces

A procedure for modeling the heart valves is presented. Instead of modeling complete leaflet motion, leaflets are modeled in open and closed configurations. The geometry of each configuration can be defined, for example, from in vivo image data. This method enables significant computational savings compared with complete fluid-structure interaction and contact modeling, while maintaining realistic three-dimensional velocity and pressure distributions near the valve, which is not possible from lumped parameter modeling. Leaflets are modeled as immersed, fixed surfaces over which a resistance to flow is assigned. On the basis of local flow conditions, the resistance values assigned for each configuration are changed to switch the valve between open and closed states. This formulation allows for the pressure to be discontinuous across the valve. To illustrate the versatility of the model, realistic and patient-specific simulations are presented, as well as comparison with complete fluid-structure interaction simulation.

A simple, versatile valve model for use in lumped parameter and one-dimensional cardiovascular models

Lumped parameter and one-dimensional models of the cardiovascular system generally employ ideal cardiac and/or venous valves that open and close instantaneously. However, under normal or pathological conditions, valves can exhibit complex motions that are mainly determined by the instantaneous difference between upstream and downstream pressures. We present a simple valve model that predicts valve motion on the basis of this pressure difference, and can be used to investigate not only valve pathology, but a wide range of cardiac and vascular factors that are likely to influence valve motion.

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

While aortic valve root compliance and leaflet coaptation have significant influence on valve closure, their implications have not yet been fully evaluated. The present study developed a full fluid-structure interaction (FSI) model that is able to cope with arbitrary coaptation between the leaflets of the aortic valve during the closing phase. Two simplifications were also evaluated for the simulation of the closing phase only. One employs an FSI model with a rigid root and the other uses a "dry" (without flow) model. Numerical tests were performed to verify the model. New metrics were defined to process the results in terms of leaflet coaptation area and contact pressure. The axial displacement of the leaflets, closure time and coaptation parameters were similar in the two FSI models, whereas the dry model, with imposed uniform load on the leaflets, produced larger coaptation area and contact pressure, larger axial displacement and faster closure time compared with the FSI model. The differences were up to 30% in the coaptation area, 55% in the contact pressure and 170% in the closure time. Consequently, an FSI model should be used to accurately resolve the kinematics of the aortic valve and leaflet coaptation details during the end-closing stage.

Design and validation of a novel bioreactor to subject aortic valve leaflets to side-specific shear stress

Hemodynamic stresses are presumed to play an important role in the development of calcific aortic valve disease (CAVD). The elucidation of the shear stress mechanisms involved in the pathogenesis of CAVD has been hampered by the complexity of the native unsteady and side-specific valvular flow environment. To address this gap, this article describes the design and validation of a novel device to expose leaflet samples to time-dependent side-specific shear stress. The device built on a double cone-and-plate geometry was dimensioned based on our previous single-sided shear stress device that minimizes secondary flow effects inherent to this geometry. A fluid-structure interaction (FSI) model was designed to predict the actual shear stress produced on a tissue sample mounted in the new device. Staining was performed on porcine leaflets conditioned in the new bioreactor to assess endothelial integrity and cellular apoptosis. The FSI results demonstrated good agreement between the target (native) and the actual side-specific shear stress produced on a tissue sample. No significant difference in endothelial integrity and cellular apoptosis was detected between samples conditioned for 96 h and fresh controls. This new device will enable the investigation of valvular response to normal and pathologic hemodynamics and the potential mechano-etiology of CAVD.

Dynamics of Flow in a Mechanical Heart Valve: The Role of Leaflet Inertia and Leaflet Compliance

In this work, we examine the dynamics of fluid flow in a mechanical heart valve when the solid inertia and leaflet compliance are important. The fluid is incompressible and Newtonian, and the leaflet is an incompressible neo-Hookean material. In the case of an inertialess leaflet, we find that the maximum valve opening angle and the time that the valve remains closed increase as the shear modulus of the leaflet decreases. More importantly, the regurgitant volume decreases with decreasing shear modulus. When we examined the forces exerted on the leaflet, we found that the downward motion of the leaflet is initiated by a vertical force exerted on its right side and, later on, by a vertical force exerted on the top side of the leaflet. In the case of solid inertia, we find that the maximum valve opening angle and the regurgitant volume are larger than in the case of an inertialess leaflet. These results highlight the importance of solid compliance in the dynamics of blood flow in a mechanical heart valve. More importantly, they indicate that mechanical heart valves with compliant leaflets may have smaller regurgitant volumes and smaller shear stresses than the ones with rigid leaflets.

The effect of implantation orientation of a bileaflet mechanical heart valve on kinematics and hemodynamics in an anatomic aorta

We carry out three-dimensional high-resolution numerical simulations of a bileaflet mechanical heart valve under physiologic pulsatile flow conditions implanted at different orientations in an anatomic aorta obtained from magnetic resonance imaging (MRI) of a volunteer. We use the extensively validated for heart valve flow curvilinear-immersed boundary (CURVIB) fluid-structure interaction (FSI) solver in which the empty aorta is discretized with a curvilinear, aorta-conforming grid while the valve is handled as an immersed boundary. The motion of the valve leaflets are calculated through a strongly coupled FSI algorithm implemented in conjunction with the Aitken convergence acceleration technique. We perform simulations for three valve orientations, which differ from each other by 45 deg and compare the results in terms of leaflet motion and flow field. We show that the valve implanted symmetrically relative to the symmetry plane of the ascending aorta curvature exhibits the smallest overall asymmetry in the motion of its two leaflets and lowest rebound during closure. Consequently, we hypothesize that this orientation is beneficial to reduce the chance of intermittent regurgitation. Furthermore, we find that the valve orientation does not significantly affect the shear stress distribution in the aortic lumen, which is in agreement with previous studies.