Numerical analysis on the hemodynamics and leaflet dynamics in a bileaflet mechanical heart valve using a fluid-structure interaction method

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.

Overset meshing in combination with novel blended weak-strong fluid-structure interactions for simulations of a translating valve in series with a second valve

Mechanical circulatory support (MCS) devices can bridge the gap to transplant whilst awaiting a viable donor heart. The Realheart Total Artificial Heart is a novel positive-displacement MCS that generates pulsatile flow via bileaflet mechanical valves. This study developed a combined computational fluid dynamics and fluid-structure interaction (FSI) methodology for simulating positive displacement bileaflet valves. Overset meshing discretised the fluid domain, and a blended weak-strong coupling FSI algorithm was combined with variable time-stepping. Four operating conditions of relevant stroke lengths and rates were assessed. The results demonstrated this modelling strategy is stable and efficient for modelling positive-displacement artificial hearts.

Automatic segmentation and tracking of biological prosthetic heart valves

Purpose: Prosthetic heart valve designs must be rigorously tested using cardiovascular equipment. The valve orifice area over time constitutes a key quality metric which is typically assessed manually, thus a tedious and error-prone task. From a computer vision viewpoint, a major unsolved issue lies in the orifice being partly occluded by the leaflets' inner side or inaccurately depicted due to its transparency. Here, we address this issue, which allows us to focus on the accurate and automatic computation of valve orifice areas. Approach: We propose a segmentation approach based on the detection of the leaflets' free edges. Using video frames recorded with a high-speed digital camera during in vitro simulations, an initial estimation of the orifice area is first obtained via active contouring and thresholding and then refined to capture the leaflet free edges via a curve transformation mechanism. Results: Experiments on video data from pulsatile flow testing demonstrate the effectiveness of our approach: a root-mean-square error (RMSE) on the temporal extracted orifice areas between 0.8% and 1.2%, an average Jaccard similarity coefficient between 0.933 and 0.956, and an average Hausdorff distance between 7.2 and 11.9 pixels. Conclusions: Our approach significantly outperformed a state-of-the-art algorithm in terms of evaluation metrics related to valve design (RMSE) and computer vision (accuracy of the orifice shape). It can also cope with lower quality videos and is better at processing frames showing an almost closed valve, a crucial quality for assessing valve design malfunctions related to their improper closing.

Cavitation Suppression of Bileaflet Mechanical Heart Valves

PURPOSE

Mechanical heart valves (MHVs) are widely used to replace diseased heart valves, but it may suffer from cavitation due to the rapid closing velocity of the leaflets, resulting in the damage of red blood cells and platelets. The aim of this study is to apply computational fluid dynamics (CFD) method to investigate the cavitation in bileaflets mechanical heart valves (BMHVs) and discuss the effects of the conduit and leaflet geometries on cavitation intensity.

METHODS

Firstly, CFD method together with moving-grid technology were applied and validated by comparing with experimental results obtained from other literature. Then the leaflets movement and the flow rate of BMHVs with different conduit geometries and leaflet geometries are compared. At last, the duration time of the saturated vapor pressure and the closing velocity of leaflets at the instant of valve closure were used to represent the cavitation intensity.

RESULTS

Larger closing velocity of leaflets at the instant of valve closure means higher cavitation intensity. For BMHVs with different conduit geometries, the conduit with Valsalva sinuses has the maximum cavitation intensity and the straight conduit has the minimum cavitation intensity, but the leaflets cannot reach the fully opened state in a straight conduit. For BMHVs with different leaflet geometries, in order to minimize the cavitation intensity, the leaflets are better to have a large thickness and a small rotational radius.

CONCLUSION

CFD method is a promising method to deal with cavitation in BMHVs, and the closing velocity of leaflets has the same trend with the cavitation intensity. By using CFD method, the effects of the conduit geometry and the leaflet geometry on cavitaion in BMHVs are found out.

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

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

3D Fluid-Structure Interaction Simulation of Aortic Valves Using a Unified Continuum ALE FEM Model

Due to advances in medical imaging, computational fluid dynamics algorithms and high performance computing, computer simulation is developing into an important tool for understanding the relationship between cardiovascular diseases and intraventricular blood flow. The field of cardiac flow simulation is challenging and highly interdisciplinary. We apply a computational framework for automated solutions of partial differential equations using Finite Element Methods where any mathematical description directly can be translated to code. This allows us to develop a cardiac model where specific properties of the heart such as fluid-structure interaction of the aortic valve can be added in a modular way without extensive efforts. In previous work, we simulated the blood flow in the left ventricle of the heart. In this paper, we extend this model by placing prototypes of both a native and a mechanical aortic valve in the outflow region of the left ventricle. Numerical simulation of the blood flow in the vicinity of the valve offers the possibility to improve the treatment of aortic valve diseases as aortic stenosis (narrowing of the valve opening) or regurgitation (leaking) and to optimize the design of prosthetic heart valves in a controlled and specific way. The fluid-structure interaction and contact problem are formulated in a unified continuum model using the conservation laws for mass and momentum and a phase function. The discretization is based on an Arbitrary Lagrangian-Eulerian space-time finite element method with streamline diffusion stabilization, and it is implemented in the open source software Unicorn which shows near optimal scaling up to thousands of cores. Computational results are presented to demonstrate the capability of our framework.

Effects of Leaflet Design on Transvalvular Gradients of Bioprosthetic Heart Valves

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

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.

A novel computational model for the hemodynamics of bileaflet mechanical valves in the opening phase

A powerful alternative means to study the hemodynamics of bileaflet mechanical heart valves is the computational fluid dynamics method. It is well recognized that computational fluid dynamics allows reliable physiological blood flow simulation and measurements of flow parameters. To date, in almost all of the modeling studies on the hemodynamics of bileaflet mechanical heart valves, a velocity (mass flow)-based boundary condition and an axisymmetric geometry for the aortic root have been assigned, which, to some extent, are erroneous. Also, there have been contradictory reports of the profile of velocity in downstream of leaflets, that is, in some studies, it is suggested that the maximum blood velocity occurs in the lateral orifice, and in some other studies, it is postulated that the maximum velocities in the main and lateral orifices are identical. The reported values for the peak velocities range from 1 to 3 m/s, which highly depend on the model assumptions. The objective of this study is to demonstrate the importance of the exact anatomical model of the aortic root and the realistic boundary conditions in the hemodynamics of the bileaflet mechanical heart valves. The model considered in this study is based on the St Jude Medical valve in a novel modeling platform. Through a more realistic geometrical model for the aortic root and the St Jude Medical valve, we have developed a new set of boundary conditions in order to be used for the assessment of the hemodynamics of aortic bileaflet mechanical heart valves. The results of this study are significant for the design improvement of conventional bileaflet mechanical heart valves and for the design of the next generation of prosthetic valves.

FSI simulation of asymmetric mitral valve dynamics during diastolic filling

In this article, we present a fluid-structure interaction algorithm accounting for the mutual interaction between two rigid bodies. The algorithm was used to perform a numerical simulation of mitral valve (MV) dynamics during diastolic filling. In numerical simulations of intraventricular flow and MV motion, the asymmetry of the leaflets is often neglected. In this study the MV was rendered as two rigid, asymmetric leaflets. The 2D simulations incorporated the dynamic interaction of blood flow and leaflet motion and an imposed subject-specific, transient left ventricular wall movement obtained from ultrasound recordings. By including the full Jacobian matrix in the algorithm, the speed of the simulation was enhanced by more than 20% compared to using a diagonal Jacobian matrix. Furthermore, our results indicate that important features of the flow field may not be predicted by the use of symmetric leaflets or in the absence of an adequate model for the left atrium.

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.

Role of inertia in the interaction between oscillatory flow and a wall-mounted leaflet

The recognition of phenomena emerging from fluid-structure interaction represents a central problem in applied fluid mechanics, especially relevant in those biological systems when the tissue properties are not measurable with accuracy. The problem of an oscillatory flow over a wall-mounted leaflet is studied here with the objective to better define the role of inertia. A numerical model based on time-dependent conformal transformation is introduced to ensure an accurate solution and an intrinsic fluid-solid coupling. The results evidence the development of long wave modulation whose amplitude and frequency depend on the leaflet inertia. When inertia is above a certain threshold, the modulation becomes large enough that the leaflet ends impacting on the wall. The potential relevance in terms of mechanical valves design, or long-term testing of natural systems subjected to period forcing, is briefly discussed.