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.     

Passive flow control of bileaflet mechanical heart valve leakage flow

Blood damage and platelet activation are inherent problems with present day mechanical heart valve designs. We investigate the approach of passive flow control applied to bileaflet mechanical heart valve (BMHV) flows as a means of optimizing leakage flow hemodynamics at length scales relevant to blood damage and platelet activation. Rectangular and hemispherical vortex generator (VG) arrays were mounted on the downstream surfaces of a 25 mm St. Jude Medical valve adjacent to the b-datum leaflet edge (central line where the two leaflets touch in closed position). The effect of VGs on the flow structure emanating from the b-datum line under both pulsatile and steady flow conditions was measured using high resolution particle image velocimetry technique. The VGs were seen to spatially disperse and dissipate the coherent leakage jet structure emanating from the b-datum line. This resulted in a significant diminution of turbulence stresses, particularly with the rectangular VG configuration. This study shows that passive flow control techniques deployed on BHMVs is potentially beneficial as significant control of flow at small length scales may be achieved without altering large scale designs of the valve.     

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.     

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.     

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.     

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.     

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.     

Observation and quantification of gas bubble formation on a mechanical heart valve

Clinical studies using transcranial Doppler ultrasonography in patients with mechanical heart valves (MHV) have detected gaseous emboli. The relationship of gaseous emboli release and cavitation on MHV has been a subject of debate in the literature. To study the influence of cavitation and gas content on the formation and growth of stable gas bubbles, a mock circulatory loop, which employed a Medtronic-Hall pyrolytic carbon disk valve in the mitral position, was used. A high-speed video camera allowed observation of cavitation and gas bubble release on the inflow valve surfaces as a function of cavitation intensity and carbon dioxide (CO2) concentration, while an ultrasonic monitoring system scanned the aortic outflow tract to quantify gas bubble production by calculating the gray scale levels of the images. In the absence of cavitation, no stable gas bubbles were formed. When gas bubbles were formed, they were first seen a few milliseconds after and in the vicinity of cavitation collapse. The volume of the gas bubbles detected in the aortic track increased with both increased CO2 and increased cavitation intensity. No correlation was observed between O2 concentration and bubble volume. We conclude that cavitation is an essential precursor to stable gas bubble formation, and CO2, the most soluble blood gas, is the major component of stable gas bubbles.     

The relation between collagen fibril kinematics and mechanical properties in the mitral valve anterior leaflet

We have recently demonstrated that the mitral valve anterior leaflet (MVAL) exhibited minimal hysteresis, no strain rate sensitivity, stress relaxation but not creep (Grashow et al., 2006, Ann Biomed Eng., 34(2), pp. 315-325; Grashow et al., 2006, Ann Biomed. Eng., 34(10), pp. 1509-1518). However, the underlying structural basis for this unique quasi-elastic mechanical behavior is presently unknown. As collagen is the major structural component of the MVAL, we investigated the relation between collagen fibril kinematics (rotation and stretch) and tissue-level mechanical properties in the MVAL under biaxial loading using small angle X-ray scattering. A novel device was developed and utilized to perform simultaneous measurements of tissue level forces and strain under a planar biaxial loading state. Collagen fibril D-period strain (epsilonD) and the fibrillar angular distribution were measured under equibiaxial tension, creep, and stress relaxation to a peak tension of 90 N/m. Results indicated that, under equibiaxial tension, collagen fibril straining did not initiate until the end of the nonlinear region of the tissue-level stress-strain curve. At higher tissue tension levels, epsilonD increased linearly with increasing tension. Changes in the angular distribution of the collagen fibrils mainly occurred in the tissue toe region. Using epsilonD, the tangent modulus of collagen fibrils was estimated to be 95.5+/-25.5 MPa, which was approximately 27 times higher than the tissue tensile tangent modulus of 3.58+/-1.83 MPa. In creep tests performed at 90 N/m equibiaxial tension for 60 min, both tissue strain and epsilonD remained constant with no observable changes over the test length. In contrast, in stress relaxation tests performed for 90 min epsilonD was found to rapidly decrease in the first 10 min followed by a slower decay rate for the remainder of the test. Using a single exponential model, the time constant for the reduction in collagen fibril strain was 8.3 min, which was smaller than the tissue-level stress relaxation time constants of 22.0 and 16.9 min in the circumferential and radial directions, respectively. Moreover, there was no change in the fibril angular distribution under both creep and stress relaxation over the test period. Our results suggest that (1) the MVAL collagen fibrils do not exhibit intrinsic viscoelastic behavior, (2) tissue relaxation results from the removal of stress from the fibrils, possibly by a slipping mechanism modulated by noncollagenous components (e.g. proteoglycans), and (3) the lack of creep but the occurrence of stress relaxation suggests a "load-locking" behavior under maintained loading conditions. These unique mechanical characteristics are likely necessary for normal valvular function.     

Comparison of the hemodynamic and thrombogenic performance of two bileaflet mechanical heart valves using a CFD/FSI model

The hemodynamic and the thrombogenic performance of two commercially available bileaflet mechanical heart valves (MHVs)--the ATS Open Pivot Valve (ATS) and the St. Jude Regent Valve (SJM), was compared using a state of the art computational fluid dynamics-fluid structure interaction (CFD-FSI) methodology. A transient simulation of the ATS and SJM valves was conducted in a three-dimensional model geometry of a straight conduit with sudden expansion distal the valves, including the valve housing and detailed hinge geometry. An aortic flow waveform (60 beats/min, cardiac output 4 l/min) was applied at the inlet. The FSI formulation utilized a fully implicit coupling procedure using a separate solver for the fluid problem (FLUENT) and for the structural problem. Valve leaflet excursion and pressure differences were calculated, as well as shear stress on the leaflets and accumulated shear stress on particles released during both forward and backward flow phases through the open and closed valve, respectively. In contrast to the SJM, the ATS valve opened to less than maximal opening angle. Nevertheless, maximal and mean pressure gradients and velocity patterns through the valve orifices were comparable. Platelet stress accumulation during forward flow indicated that no platelets experienced a stress accumulation higher than 35 dyne x s/cm2, the threshold for platelet activation (Hellums criterion). However, during the regurgitation flow phase, 0.81% of the platelets in the SJM valve experienced a stress accumulation higher than 35 dyne x s/cm2, compared with 0.63% for the ATS valve. The numerical results indicate that the designs of the ATS and SJM valves, which differ mostly in their hinge mechanism, lead to different potential for platelet activation, especially during the regurgitation phase. This numerical methodology can be used to assess the effects of design parameters on the flow induced thrombogenic potential of blood recirculating devices.     

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.     

Comparison of hinge microflow fields of bileaflet mechanical heart valves implanted in different sinus shape and downstream geometry

The characterization of the bileaflet mechanical heart valves (BMHVs) hinge microflow fields is a crucial step in heart valve engineering. Earlier in vitro studies of BMHV hinge flow at the aorta position in idealized straight pipes have shown that the aortic sinus shapes and sizes may have a direct impact on hinge microflow fields. In this paper, we used a numerical study to look at how different aortic sinus shapes, the downstream aortic arch geometry, and the location of the hinge recess can influence the flow fields in the hinge regions. Two geometric models for sinus were investigated: a simplified axisymmetric sinus and an idealized three-sinus aortic root model, with two different downstream geometries: a straight pipe and a simplified curved aortic arch. The flow fields of a 29-mm St Jude Medical BMHV with its four hinges were investigated. The simulations were performed throughout the entire cardiac cycle. At peak systole, recirculating flows were observed in curved downsteam aortic arch unlike in straight downstream pipe. Highly complex three-dimensional leakage flow through the hinge gap was observed in the simulation results during early diastole with the highest velocity at 4.7 m/s, whose intensity decreased toward late diastole. Also, elevated wall shear stresses were observed in the ventricular regions of the hinge recess with the highest recorded at 1.65 kPa. Different flow patterns were observed between the hinge regions in straight pipe and curved aortic arch models. We compared the four hinge regions at peak systole in an aortic arch downstream model and found that each individual hinge did not vary much in terms of the leakage flow rate through the valves.     

Wavelet transforms in the analysis of mechanical heart valve cavitation

Cavitation is known to cause blood element damage and may introduce gaseous emboli into the cerebral circulation, increasing the patient's risk of stroke. Discovering methods to reduce the intensity of cavitation induced by mechanical heart valves (MHVs) has long been an area of interest. A novel approach for analyzing MHV cavitation is presented. A wavelet denoising method is explored because currently used analytical techniques fail to suitably unmask the cavitation signal from other valve closing sounds and noise detected with a hydrophone. Wavelet functions are used to denoise the cavitation signal during MHV closure and rebound. The wavelet technique is applied to the signal produced by closure of a 29-mm Medtronic-Hall MHV in degassed water with a gas content of 5 ppm. Valve closing dynamics are investigated under loading conditions of 500, 2500, and 4500 mm Hg/s. The results display a marked improvement in the quantity and quality of information that can be extracted from acoustic cavitation signals using the wavelet technique compared to conventional analytical techniques. Time and frequency data indicate the likelihood and characteristics of cavitation formation under specified conditions. Using this wavelet technique we observe an improved signal-to-noise ratio, an enhanced time-dependent aspect, and the potential to minimize valve closing sounds, which disguise individual cavitation events. The overall goal of this work is to eventually link specific valves with characteristic waveforms or distinct types of cavitation, thus promoting improved valve designs.     

The impact of imperfect frame deployment and rotational orientation on stress within the prosthetic leaflets during transcatheter aortic valve implantation

TAVI devices are manufactured with cylindrical frames. However, the frames are rarely cylindrical post-deployment since deformation due to localised under expansion can be induced by calcified material on the native valve leaflets exerting irregular forces upon the frame. Consequently, the leaflets within a deformed TAVI device may undergo elevated stress during operation, which may lead to premature device failure.Using computational analysis a complete TAVI device model was simulated undergoing deployment into an aortic root model derived from CT data for a patient with severe calcific aortic stenosis, followed by a pressure simulated cardiac cycle. The complete analysis was performed eight times, each with the device at a different rotational orientation relative to the native valve, with an increment spacing of 15°.The TAVI device frames consistently featured significant distortions associated with bulky calcified material at the base of the non-coronary sinus. It was found that the average von Mises stress in the prosthetic valves was only increased in one of the cases relative to an idealised device. However, the maximum von Mises stress in the prosthetic valves was elevated in the majority of the cases.Furthermore, it was found that there were preferable orientations to deploy the prosthetic device, in this case, when the prosthetic leaflets were aligned with the native leaflets. As device orientation deviated from this orientation, the stresses in the valve increased because the distance between the prosthetic commissures decreased. This potentially could represent a sufficient increase in stress to induce variation in device lifespan.     

Dynamic analysis of flutter in disk type mechanical heart valve prostheses

Parametric study of the low frequency oscillations occasionally observed in certain types of disc type prosthetic heart valves (PHV) are carried out using a finite element technique. The analysis is performed to determine the frequencies of the dynamic fluttering with the help of the 'ANSYS' computer program. The results show that the frequencies of the dynamic fluttering for both the circular occluders and the semi-circular occluders are at least two orders of magnitude higher than that observed in vivo. It is thus concluded that the clinically observed leaflet oscillations should not be a dynamic flutter phenomenon. Rather, the vortex shedding has been assumed to be the cause of these oscillations.