Near valve flows and potential blood damage during closure of a bileaflet mechanical heart valve

Blood damage and thrombosis are major complications that are commonly seen in patients with implanted mechanical heart valves. For this in vitro study, we isolated the closing phase of a bileaflet mechanical heart valve to study near valve fluid velocities and stresses. By manipulating the valve housing, we gained optical access to a previously inaccessible region of the flow. Laser Doppler velocimetry and particle image velocimetry were used to characterize the flow regime and help to identify the key design characteristics responsible for high shear and rotational flow. Impact of the closing mechanical leaflet with its rigid housing produced the highest fluid stresses observed during the cardiac cycle. Mean velocities as high as 2.4 m/s were observed at the initial valve impact. The velocities measured at the leaflet tip resulted in sustained shear rates in the range of 1500-3500 s(-1), with peak values on the order of 11,000-23,000 s(-1). Using velocity maps, we identified regurgitation zones near the valve tip and through the central orifice of the valve. Entrained flow from the transvalvular jets and flow shed off the leaflet tip during closure combined to generate a dominant vortex posterior to both leaflets after each valve closing cycle. The strength of the peripheral vortex peaked within 2 ms of the initial impact of the leaflet with the housing and rapidly dissipated thereafter, whereas the vortex near the central orifice continued to grow during the rebound phase of the valve. Rebound of the leaflets played a secondary role in sustaining closure-induced vortices.     

Computational simulations of flow dynamics and blood damage through a bileaflet mechanical heart valve scaled to pediatric size and flow

Despite pressing needs, there are currently no FDA approved prosthetic valves available for use in the pediatric population. This study is performed for predictive assessment of blood damage in bileaflet mechanical heart valves (BMHVs) with pediatric sizing and flow conditions. A model of an adult-sized 23 mm St. Jude Medical (SJM) Regent^™ valve is selected for use in simulations, which is scaled in size for a 5-year old child and 6-month old infant. A previously validated lattice-Boltzmann method (LBM) is used to simulate pulsatile flow with thousands of suspended platelets for cases of adult, child, and infant BMHV flows. Adult BMHV flows demonstrate more disorganized small-scale flow features, but pediatric flows are associated with higher fluid shear stresses. Platelet damage in the pediatric cases is higher than in adult flow, highlighting thrombus complication dangers of pediatric BMHV flows. This does not necessarily suggest clinically important differences in thromboembolic potential. Highly damaged platelets in pediatric flows are primarily found far downstream of the valve, as there is less flow recirculation in pediatric flows. In addition, damage levels are well below expected thresholds for platelet activation. The extent of differences here documented between the pediatric and adult cases is of concern, demanding particular attention when pediatric valves are designed and manufactured. However, the differences between the pediatric and adult cases are not such that development of pediatric sized valves is untenable. This study may push for eventual approval of prosthetic valves resized for the pediatric population. Further studies will be necessary to determine the validity and potential thrombotic and clinical implications of these findings.     

An analysis of turbulent shear stresses in leakage flow through a bileaflet mechanical prostheses

In this work, estimates of turbulence were made from pulsatile flow laser Doppler velocimetry measurements using traditional phase averaging and averaging after the removal of cyclic variation. These estimates were compared with estimates obtained from steady leakage flow LDV measurements and an analytical method. The results of these studies indicate that leakage jets which are free and planar in shape may be more unstable than other leakage jets, and that cyclic variation does not cause a gross overestimation of the Reynolds stresses at large distances from the leakage jet orifice.     

Flow in a mechanical bileaflet heart valve at laminar and near-peak systole flow rates: CFD simulations and experiments

Time-accurate, fully 3D numerical simulations and particle image velocity laboratory experiments are carried out for flow through a fully open bileaflet mechanical heart valve under steady (nonpulsatile) inflow conditions. Flows at two different Reynolds numbers, one in the laminar regime and the other turbulent (near-peak systole flow rate), are investigated. A direct numerical simulation is carried out for the laminar flow case while the turbulent flow is investigated with two different unsteady statistical turbulence modeling approaches, unsteady Reynolds-averaged Navier-Stokes (URANS) and detached-eddy simulation (DES) approach. For both the laminar and turbulent cases the computed mean velocity profiles are in good overall agreement with the measurements. For the turbulent simulations, however, the comparisons with the measurements demonstrate clearly the superiority of the DES approach and underscore its potential as a powerful modeling tool of cardiovascular flows at physiological conditions. The study reveals numerous previously unknown features of the flow.     

Flow and thrombosis at orifices simulating mechanical heart valve leakage regions

BACKGROUND: While it is established that mechanical heart valves (MHVs) damage blood elements during leakage and forward flow, the role in thrombus formation of platelet activation by high shear flow geometries remains unclear. In this study, continuously recalcified blood was used to measure the effects of blood flow through orifices, which model MHVs, on the generation of procoagulant thrombin and the resulting formation of thrombus. The contribution of platelets to this process was also assessed. METHOD OF APPROACH: 200, 400, 800, and 1200 microm orifices simulated the hinge region of bileaflet MHVs, and 200, 400, and 800 microm wide slits modeled the centerline where the two leaflets meet when the MHV is closed. To assess activation of coagulation during blood recirculation, samples were withdrawn over 0-47 min and the plasmas assayed for thrombin-antithrombin-llI (TAT) levels. Model geometries were also inspected visually. RESULTS: The 200 and 400 microm round orifices induced significant TAT generation and thrombosis over the study interval. In contrast, thrombin generation by the slit orifices, and by the 800 and 1200 microm round orifices, was negligible. In additional experiments with nonrecalcified or platelet-depleted blood, TAT levels were markedly reduced versus the studies with fully anticoagulated whole blood (p < 0.05). CONCLUSIONS: Using the present method, a significant increase in TAT concentration was found for 200 and 400 microm orifices, but not 800 and 1200 microm orifices, indicating that these flow geometries exhibit a critical threshold for activation of coagulation and resulting formation of thrombus. Markedly lower TAT levels were produced in studies with platelet-depleted blood, documenting a key role for platelets in the thrombotic process.     

Fluid-structure interaction of deformable aortic prostheses with a bileaflet mechanical valve

Two different aortic prostheses can be used for performing the Bentall procedure: a standard straight graft and the Valsalva graft that better reproduces the aortic root anatomy. The aim of the present work is to study the effect of the graft geometry on the blood flow when a bileaflet mechanical heart valve is used, as well as to evaluate the stress concentration near the suture line where the coronary arteries are connected to graft. An accurate three-dimensional numerical method is proposed, based on the immersed boundary technique. The method accounts for the interactions between the flow and the motion of the rigid leaflets and of the deformable aortic root, under physiological pulsatile conditions. The results show that the graft geometry only slightly influences the leaflets dynamics, while using the Valsalva graft the stress level near the coronary-root anastomoses is about half that obtained using the standard straight graft.     

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.     

Time-resolved DPIV analysis of vortex dynamics in a left ventricular model through bileaflet mechanical and porcine heart valve prostheses

The performance of the heart after a mitral valve replacement operation greatly depends on the flow character downstream of the valve. The design and implanting orientation of valves may considerably affect the flow development. A study of the hemodynamics of two orientations, anatomical and anti-anatomical, of the St. Jude Medical (SJM) bileaflet valve are presented and compared with those of the SJM Biocor porcine valve, which served also to represent the natural valve. We document the velocity field in a flexible, transparent (LV) using time-resolved digital particle image velocimetry (TRDPIV). Vortex formation and vortex interaction are two important physical phenomena that dominate the filling and emptying of the ventricle. For the three configurations, the following effects were examined: mitral valve inlet jet asymmetry, survival of vortical structures upstream of the aortic valve, vortex-induced velocities and redirection of theflow in abidance of the Biot-Savart law, domain segmentation, resonant times of vortical structures, and regions of stagnantflow. The presence of three distinct flow patterns, for the three configurations, was identified by the location of vortical structures and level of coherence corresponding to a significant variation in the turbulence level distribution inside the LV. The adverse effect of these observations could potentially compromise the efficiency of the LV and result in flow patterns that deviate from those in the natural heart.     

Numerical and experimental investigations of pulsatile blood flow pattern through a dysfunctional mechanical heart valve

Around 250,000 heart valve replacements are performed every year around the world. Due their higher durability, approximately 2/3 of these replacements use mechanical prosthetic heart valves (mainly bileaflet valves). Although very efficient, these valves can be subject to valve leaflet malfunctions. These malfunctions are usually the consequence of pannus ingrowth and/or thrombus formation and represent serious and potentially fatal complications. Hence, it is important to investigate the flow field downstream of a dysfunctional mechanical heart valve to better understand its impact on blood components (red blood cells, platelets and coagulation factors) and to improve the current diagnosis techniques. Therefore, the objective of this study will be to numerically and experimentally investigate the pulsatile turbulent flow downstream of a dysfunctional bileaflet mechanical heart valve in terms of velocity field, vortex formation and potential negative effect on blood components. The results show that the flow downstream of a dysfunctional valve was characterized by abnormally elevated velocities and shear stresses as well as large scale vortices. These characteristics can predispose to blood components damage. Furthermore, valve malfunction led to an underestimation of maximal transvalvular pressure gradient, using Doppler echocardiography, when compared to numerical results. This could be explained by the shifting of the maximal velocity towards the normally functioning leaflet. As a consequence, clinicians should try, when possible, to check the maximal velocity position not only at the central orifice but also through the lateral orifices. Finding the maximal velocity in the lateral orifice could be an indication of valve dysfunction.     

Turbulence characteristics downstream of bileaflet aortic valve prostheses

This study was focused on a series of in vitro tests on the turbulent flow characteristics of three bileaflet aortic valves: St. Jude Medical (SJM), CarboMedics (CM), and Edwards Tekna (modified Duromedics, DM). The flow fields of the valves were measured in a pulsatile flow model with a laser-Doppler anemometer (LDA) at the aortic sinus area downstream of the valves. The heart rate was set at 70 beats per minute, the cardiac output was maintained at 5 liters per minute, and the aortic pressure wave forms were kept within the physiological range. Cycle-resolved analysis was applied to obtain turbulence data, including mean velocity, Reynolds stresses, autocorrelation coefficients, energy spectral density functions, and turbulence scales. The Reynolds shear stresses of all three valves induced only minor damage to red blood cells, but directly damaged the platelets, increasing the possibility of thrombosis. The smallest turbulence length scale, which offers a more reliable estimate of the effects of turbulence on blood cell damage, was three times the size of red blood cells and five times the size of platelets. This suggests that there is more direct interaction with the blood cells, thus causing more damage.     

Numerical simulation of the dynamics of a bileaflet prosthetic heart valve using a fluid-structure interaction approach

The main purpose of this study is to reproduce in silico the dynamics of a bileaflet mechanical heart valve (MHV; St Jude Hemodynamic Plus, 27mm characteristic size) by means of a fully implicit fluid-structure interaction (FSI) method, and experimentally validate the results using an ultrafast cinematographic technique. The computational model was constructed to realistically reproduce the boundary condition (72 beats per minute (bpm), cardiac output 4.5l/min) and the geometry of the experimental setup, including the valve housing and the hinge configuration. The simulation was carried out coupling a commercial computational fluid dynamics (CFD) package based on finite-volume method with user-defined code for solving the structural domain, and exploiting the parallel performance of the whole numerical setup. Outputs are leaflets excursion from opening to closure and the fluid dynamics through the valve. Results put in evidence a favorable comparison between the computed and the experimental data: the model captures the main features of the leaflet motion during the systole. The use of parallel computing drastically limited the computational costs, showing a linear scaling on 16 processors (despite the massive use of user-defined subroutines to manage the FSI process). The favorable agreement obtained between in vitro and in silico results of the leaflet displacements confirms the consistency of the numerical method used, and candidates the application of FSI models to become a major tool to optimize the MHV design and eventually provides useful information to surgeons.     

Cavitation caused by mechanical heart valve prostheses--a review

Heart valve dysfunction often necessitates surgical implantation of a mechanical heart valve (MHV). Although implantation of a MHV is a life-saving procedure, the patient still faces potentially complications such as thromboembolic events and material failure. These complications may be caused by cavitation, which can occur during valve closure. Cavitation is an erosive phenomenon that can be generated in fluids when the pressure locally drops below the vapor pressure. This paper reviews the literature on cavitation and MHVs and particular features of the valve and closing conditions that potentially increase the intensity of cavitation. Techniques for detecting cavitation will be discussed. Of these, an acoustic approach will be emphasized, since it is currently the only technique able to detect and quantify cavitation in vivo.     

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.     

Comparison of first and second heart sounds after mechanical heart valve replacement

In this article, the spectral features of first heart sounds (S1) and second heart sounds (S2), which comprise the mechanical heart valve sounds obtained after aortic valve replacement (AVR) and mitral valve replacement (MVR), are compared to find out the effect of mechanical heart valve replacement and recording area on S1 and S2. For this aim, the Welch method and the autoregressive (AR) method are applied on the S1 and S2 taken from 66 recordings of 8 patients with AVR and 98 recordings from 11 patients with MVR, thereby yielding power spectrum of the heart sounds. Three features relating to frequency of heart sounds and three features relating to energy of heart sounds are obtained. Results show that in comparison to natural heart valves, mechanical heart valves contain higher frequency components and energy, and energy and frequency components do not show common behaviour for either AVR or MVR depending on the recording areas. Aside from the frequency content and energy of the sound generated by mechanical heart valves being affected by the structure of the lungs–thorax and the recording areas, the pressure across the valve incurred during AVR or MVR is a significant factor in determining the frequency and energy levels of the valve sound produced. Though studies on native heart sounds as a non-invasive diagnostic method has been done for many years, it is observed that studies on mechanical heart valves sounds are limited. The results of this paper will contribute to other studies on using a non-invasive method for assessing the mechanical heart valve sounds.     

Model studies at mechanical aortic heart valve prostheses--Part I: Steady-state flow fields and pressure loss coefficients

In a 3:1 scaled model of the human aorta models of the Omniscience, Bj?rk-Shiley Convexo-Concave, Bj?rk-Shiley Monostrut, Medtronic-Hall, Duromedics (Hemex) and the Saint Jude Medical heart valve prostheses are studied in steady flow representing the systolic peak flow phase. Detailed flow visualization experiments show flow separations at all inner ring surfaces as well as at most of the occluders. The resulting stagnation areas increase the risk of thrombus accumulation. Flow separations also stimulate vortex formation and turbulent mixing at the downstream jet boundaries and thus may intensify blood damage by turbulent shear stresses. The different influences of struts and occluder guides on the flow around the occluders are discussed. The effects of the individual valve components on the flow fields are analyzed and correlated with the resulting pressure losses.     

The influence of the leaflets' curvature on the flow field in two bileaflet prosthetic heart valves

A successful mechanical prosthetic heart valve design is the bileaflet valve, which has been implanted for the first time more than 20 years ago. A key feature of bileaflet valves is the geometry of the two leaflets, which can be very important in determining the flow field. Laser Doppler anemometry (LDA) was used to perform an accurate study of the velocity and turbulence shear stress peak values (TSS(max)) fields at four distances from the valve plane. TSS(max) is a relevant parameter to assess the risk of hemolysis and platelet activation associated to the implantation of a prosthetic device, continuously interacting with blood. Two bileaflet valves were tested: the St. Jude HP and the Sorin Bicarbon, of the same nominal size (19mm). The former has flat leaflets, whereas the latter's leaflets have a cylindrical surface. A high regime (CO: 6l/min) was imposed, in order to test the two valves at maximum Reynolds number and consequent turbulence generation. The flat-leaflet design of the St. Jude generates a TSS field constant with distance; on the contrary, the Bicarbon's shear stress field undergoes an evident development, with an unexpected central peak at a distance comparable to the valve's dimensions (21mm). The two bileaflet valves tested, although very similar in design, behave very differently as for their turbulence properties. In particular, the concept of curved wake leads to conclude that the curvature of the leaflets' surface must be identified as an important parameter, which deserves careful attention in PHV design and development.     

Time-dependent biaxial mechanical behavior of the aortic heart valve leaflet

Despite continued progress in the treatment of aortic valve (AV) disease, current treatments continue to be challenged to consistently restore AV function for extended durations. Improved approaches for AV repair and replacement rests upon our ability to more fully comprehend and simulate AV function. While the elastic behavior the AV leaflet (AVL) has been previously investigated, time-dependent behaviors under physiological biaxial loading states have yet to be quantified. In the current study, we performed strain rate, creep, and stress-relaxation experiments using porcine AVL under planar biaxial stretch and loaded to physiological levels (60 N/m equi-biaxial tension), with strain rates ranging from quasi-static to physiologic. The resulting stress-strain responses were found to be independent of strain rate, as was the observed low level of hysteresis ( approximately 17%). Stress relaxation and creep results indicated that while the AVL exhibited significant stress relaxation, it exhibited negligible creep over the 3h test duration. These results are all in accordance with our previous findings for the mitral valve anterior leaflet (MVAL) [Grashow, J.S., Sacks, M.S., Liao, J., Yoganathan, A.P., 2006a. Planar biaxial creep and stress relaxatin of the mitral valve anterior leaflet. Annals of Biomedical Engineering 34 (10), 1509-1518; Grashow, J.S., Yoganathan, A.P., Sacks, M.S., 2006b. Biaxial stress-stretch behavior of the mitral valve anterior leaflet at physiologic strain rates. Annals of Biomedical Engineering 34 (2), 315-325], and support our observations that valvular tissues are functionally anisotropic, quasi-elastic biological materials. These results appear to be unique to valvular tissues, and indicate an ability to withstand loading without time-dependent effects under physiologic loading conditions. Based on a recent study that suggested valvular collagen fibrils are not intrinsically viscoelastic [Liao, J., Yang, L., Grashow, J., Sacks, M.S., 2007. The relation between collagen fibril kinematics and mechanical properties in the mitral valve anterior leaflet. Journal of Biomechanical Engineering 129 (1), 78-87], we speculate that the mechanisms underlying this quasi-elastic behavior may be attributed to inter-fibrillar structures unique to valvular tissues. These mechanisms are an important functional aspect of native valvular tissues, and are likely critical to improve our understanding of valvular disease and help guide the development of valvular tissue engineering and surgical repair.     

A mechanical analysis of the closed Hancock heart valve prosthesis

In order to obtain mechanical specifications for the design of an artificial leaflet valve prosthesis, a geometrically non-linear numerical model is developed of a closed Hancock leaflet valve prosthesis. In this model, the fibre reinforcement of the leaflet and the viscoelastic properties of frame and leaflets are incorporated. The calculations are primarily restricted to 1/6 part of the valve and a time varying pressure load is applied. The calculations are verified experimentally by measuring the commissure displacements and leaflet centre displacement of a Hancock valve. The numerically obtained commissure displacements are found to be linearly dependent on the pressure load, and the slope of the curves is hardly dependent on loading type and loading velocity. Experimentally a difference is found between the three commissure displacements, which is also predicted numerically using a simplified asymmetric total valve model. Besides, experimentally a clear dependency of commissure displacements on frame size is found. For the leaflet centre displacement, a qualitative agreement exists between numerical prediction and experimental result, although the numerical predicted values are systematically higher. The numerically obtained stress distributions revealed that the maximum von Mises intensity in the membranes occurs in the vicinity of the commissure in the free leaflet area (0.2 N mm-2). Wrinkling of the membranes may occur in the coaptation area near the leaflet suspension. The maximum fibre stress is found near the aortic ring in the fibres which form the boundaries of the coaptation area (0.64 N mm-2). These locations seem to correlate with some common regions of tissue valve failure.