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.     

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.     

A detailed fluid mechanics study of tilting disk mechanical heart valve closure and the implications to blood damage

Hemolysis and thrombosis are among the most detrimental effects associated with mechanical heart valves. The strength and structure of the flows generated by the closure of mechanical heart valves can be correlated with the extent of blood damage. In this in vitro study, a tilting disk mechanical heart valve has been modified to measure the flow created within the valve housing during the closing phase. This is the first study to focus on the region just upstream of the mitral valve occluder during this part of the cardiac cycle, where cavitation is known to occur and blood damage is most severe. Closure of the tilting disk valve was studied in a "single shot" chamber driven by a pneumatic pump. Laser Doppler velocimetry was used to measure all three velocity components over a 30 ms period encompassing the initial valve impact and rebound. An acrylic window placed in the housing enabled us to make flow measurements as close as 200 microm away from the closed occluder. Velocity profiles reveal the development of an atrial vortex on the major orifice side of the valve shed off the tip of the leaflet. The vortex strength makes this region susceptible to cavitation. Mean and maximum axial velocities as high as 7 ms and 20 ms were recorded, respectively. At closure, peak wall shear rates of 80,000 s(-1) were calculated close to the valve tip. The region of the flow examined here has been identified as a likely location of hemolysis and thrombosis in tilting disk valves. The results of this first comprehensive study measuring the flow within the housing of a tilting disk valve may be helpful in minimizing the extent of blood damage through the combined efforts of experimental and computational fluid dynamics to improve mechanical heart valve designs.     

Lumped parameter model for computing the minimum pressure during mechanical heart valve closure

The cavitation inception threshold of mechanical heart valves has been shown to be highly variable. This is in part due to the random distribution of the initial and final conditions that characterize leaflet closure. While numerous hypotheses exist explaining the mechanisms of inception, no consistent scaling laws have been developed to describe this phenomenon due to the complex nature of these dynamic conditions. Thus in order to isolate and assess the impact of these varied conditions and mechanisms on inception, a system of ordinary differential equations is developed to describe each system component and solved numerically to predict the minimum pressure generated during valve closure. In addition, an experiment was conducted in a mock circulatory loop using an optically transparent size 29 bileaflet valve over a range of conditions to calibrate and validate this model under physiological conditions. High-speed video and high-response pressure measurements were obtained simultaneously to characterize the relationship between the valve motion, fluid motion, and negative pressure transients during closure. The simulation model was calibrated using data from a single closure cycle and then compared to other experimental flow conditions and to results found in the literature. The simulation showed good agreement with the closing dynamics and with the minimum pressure trends in the current experiment. Additionally, the simulation suggests that the variability observed experimentally (when using dP/dt alone as the primary measure of cavitation inception) is predictable. Overall, results from the current form of this lumped parameter model indicate that it is a good engineering assessment tool.     

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.     

A numerical simulation of mechanical heart valve closure fluid dynamics

A computational fluid dynamics model for the analysis of the bileaflet mechanical heart valve closure process is presented. The objective of the study is to demonstrate the ability of the numerical model to simulate the leaflet motion during the closing phase in order to investigate the closure fluid dynamics and to evaluate the effect of alterations in the leaflet tip geometry. The model has been applied to six different combinations of the leaflet tip geometry and the gap width between the leaflet tip and the housing. The results show that the negative pressure quickly develops on the atrial side of the leaflet tip. The pressure becomes more negative as the leaflet closure progresses and the lowest pressure is reached before the leaflet comes to a stop in the closed position. The flow dynamics at the instant of valve closure is strongly dependent on the leaflet velocity during the closing phase. Decrease of the tip velocity by a factor of three in the last four degrees of leaflet motion leads to a 50% reduction in the negative pressure magnitude.     

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.     

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.     

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 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.     

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.     

Fundamental mechanics of aortic heart valve closure

Stresses in a prosthetic heart valve at closure are determined by its geometrical and structural characteristics, by the mechanical support environment, and by the momentum of the valve leaflets or occluder and of the blood at the instant of closure. The mass of blood to be arrested is significantly greater than that of the leaflets or occluder, and is therefore likely to dominate the closure impulse. The kinetic energy of the blood must be transduced into potential energy in the structural components (valve leaflets, aortic root and aorta). This paper presents a methodology for computation and parameterisation of the blood momentum associated with a valve in the aortic position. It is suggested that the influence of physiological parameters, such as systolic waveform and systemic impedance, on the closure characteristics can be investigated based on the fluid dynamic implications. Detailed results are presented for a single leaflet mechanical valve (Bjork-Shiley 60 degrees Convexo-Concave). It is demonstrated that a simple analytical method can yield results that might be adequate for the purposes of valve design.     

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.     

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.