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.     

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.     

Device Thrombogenicity Emulator (DTE) − Design optimization methodology for cardiovascular devices: A study in two bileaflet MHV designs

Patients who receive prosthetic heart valve (PHV) implants require mandatory anticoagulation medication after implantation due to the thrombogenic potential of the valve. Optimization of PHV designs may facilitate reduction of flow-induced thrombogenicity and reduce or eliminate the need for post-implant anticoagulants. We present a methodology entitled Device Thrombogenicty Emulator (DTE) for optimizing the thrombo-resistance performance of PHV by combining numerical and experimental approaches. Two bileaflet mechanical heart valves (MHV) designs, St. Jude Medical (SJM) and ATS, were investigated by studying the effect of distinct flow phases on platelet activation. Transient turbulent and direct numerical simulations (DNS) were conducted, and stress loading histories experienced by the platelets were calculated along flow trajectories. The numerical simulations indicated distinct design dependent differences between the two valves. The stress loading waveforms extracted from the numerical simulations were programmed into a hemodynamic shearing device (HSD), emulating the flow conditions past the valves in distinct ‘hot-spot’ flow regions that are implicated in MHV thrombogenicity. The resultant platelet activity was measured with a modified prothrombinase assay, and was found to be significantly higher in the SJM valve, mostly during the regurgitation phase. The experimental results were in excellent agreement with the calculated platelet activation potential. This establishes the utility of the DTE methodology for serving as a test bed for evaluating design modifications for achieving better thrombogenic performance for such devices.     

Particle Image Velocimetry Study of Pulsatile Flow in Bi-leaflet Mechanical Heart Valves with Image Compensation Method

Particle Image Velocimetry (PIV) is an important technique in studying blood flow in heart valves. Previous PIV studies of flow around prosthetic heart valves had different research concentrations, and thus never provided the physical flow field pictures in a complete heart cycle, which compromised their pertinence for a better understanding of the valvular mechanism. In this study, a digital PIV (DPIV) investigation was carried out with improved accuracy, to analyse the pulsatile flow field around the bi-leaflet mechanical heart valve (MHV) in a complete heart cycle. For this purpose a pulsatile flow test rig was constructed to provide the necessary in vitro test environment, and the flow field around a St. Jude size 29 bi-leaflet MHV and a similar MHV model were studied under a simulated physiological pressure waveform with flow rate of 5.2 l/min and pulse rate at 72 beats/min. A phase-locking method was applied to gate the dynamic process of valve leaflet motions. A special image-processing program was applied to eliminate optical distortion caused by the difference in refractive indexes between the blood analogue fluid and the test section. Results clearly showed that, due to the presence of the two leaflets, the valvular flow conduit was partitioned into three flow channels. In the opening process, flow in the two side channels was first to develop under the presence of the forward pressure gradient. The flow in the central channel was developed much later at about the mid-stage of the opening process. Forward flows in all three channels were observed at the late stage of the opening process. At the early closing process, a backward flow developed first in the central channel. Under the influence of the reverse pressure gradient, the flow in the central channel first appeared to be disturbed, which was then transformed into backward flow. The backward flow in the central channel was found to be the main driving factor for the leaflet rotation in the valve closing process. After the valve was fully closed, local flow activities in the proximity of the valve region persisted for a certain time before slowly dying out. In both the valve opening and closing processes, maximum velocity always appeared near the leaflet trailing edges. The flow field features revealed in the present paper improved our understanding of valve motion mechanism under physiological conditions, and this knowledge is very helpful in designing the new generation of MHVs.     

Free emboli formation in the wake of bi-leaflet mechanical heart valves and the effects of implantation techniques

The high incidence of thromboembolic complications of mechanical heart valves (MHV), primarily due to platelet activation by contact with foreign surfaces and by non-physiological flow patterns past the valve, still limits their success as permanent implants. The latter include elevated shear and turbulent stresses and shed vortices formed in the wake of the valve's leaflets during the deceleration phase, potentially entrapping activated and aggregated platelets. It is hypothesized that these flow patterns induce the formation of free emboli which are the source of cerebrovascular microemboli associated with MHV. Implicit to this hypothesis is that free emboli formation will be affected by the implantation technique employed and the valve orientation, as those will alter the flow characteristics past the valve and the interaction of the platelets with the flow. In this study, numerical simulations of turbulent pulsatile flow past a St. Jude Medical bi-leaflet MHV were conducted. Platelet shear histories were calculated along pertinent turbulent platelet trajectories, and the effect of a misaligned valve on platelet activation was quantified and compared to that of an aligned valve. It demonstrated that the combination of a tilted valve and subannularly sutured pledgets had an explicit detrimental effect on platelet activation, with the following entrapment of the platelets within the shed vortices of the wake leading to a significant increase of the thromboembolic potential of the valve. This numerical model depicted a viable course for free emboli formation, and indicated how the implantation technique may enhance the risk of cardioembolism.     

Vortex shedding as a mechanism for free emboli formation in mechanical heart valves

The high incidence of thromboembolic complications of mechanical heart valves (MHV) limits their success as permanent implants. The thrombogenicity of all MHV is primarily due to platelet activation by contact with foreign surfaces and by nonphysiological flow patterns. The latter include elevated flow stresses and regions of recirculation of blood that are induced by valve design characteristics. A numerical simulation of unsteady turbulent flow through a bileaflet MHV was conducted, using the Wilcox k-omega turbulence model for internal low-Reynolds-number flows, and compared to quantitative flow visualization performed in a pulse duplicator system using Digital Particle Image Velocimetry (DPIV). The wake of the valve leaflet during the deceleration phase revealed an intricate pattern of interacting shed vortices. Particle paths showed that platelets that were exposed to the highest flow stresses around the leaflets were entrapped within the shed vortices. Potentially activated, such platelets may tend to aggregate and form free emboli. Once formed, such free emboli would be convected downstream by the shed vortices, increasing the risk of systemic emboli.     

新型机械心脏瓣膜的研究

目前临床使用的各种机械心脏瓣膜的主要问题是血栓栓塞和与抗凝治疗有关的出血,其缺陷在于瓣膜开启时,碟片和支架将瓣膜的整个血流通道分隔成三至四个较小的血流通道。在这种受阻隔的血流通宫,形成容易诱发血栓的高剪应力区、紊流和滞流区。我们研制的两种机械心脏瓣膜在瓣膜开启时,没有任何支架和碟片分隔瓣膜的血流通道,使血流与天然心脏瓣膜中的相类似,可减少对血液的危害,从而可减少换瓣病人对抗凝治疗的依赖程度。     

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.     

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.     

Turbulent shear stress measurements in the vicinity of aortic heart valve prostheses

A two dimensional laser Doppler anemometer system has been used to measure the turbulent shear fields in the immediate downstream vicinity of a variety of mechanical and bioprosthetic aortic heart valves. The measurements revealed that all the mechanical valves studied, created regions of elevated levels of turbulent shear stress during the major portion of systole. The tissue bioprostheses also created elevated levels of turbulence, but they were confined to narrow regions in the bulk of the flow field. The newer generation of bioprostheses create turbulent shear stresses which are considerably lower than those created by the older generation tissue valve designs. All the aortic valves studied (mechanical and tissue) create turbulent shear stress levels which are capable of causing sub-lethal and/or lethal damage to blood elements.     

A novel mathematical model of activation and sensitization of platelets subjected to dynamic stress histories

Blood recirculating devices, such as ventricular assist devices and prosthetic heart valves, are burdened by thromboembolic complications requiring complex and lifelong anticoagulant therapy with its inherent hemorrhagic risks. Pathologic flow patterns occurring in such devices chronically activate platelets, and the optimization of their thrombogenic performance requires the development of flow-induced platelet activation models. However, existing models are based on empirical correlations using the well-established power law paradigm of constant levels of shear stress during certain exposure times as factors for mechanical platelet activation. These models are limited by their range of application and do not account for other relevant phenomena, such as loading rate dependence and platelet sensitization to high stress conditions, which characterize the dynamic flow conditions in devices. These limitations were addressed by developing a new class of phenomenological stress-induced platelet activation models that specifies the rate of platelet activation as a function of the entire stress history and results in a differential equation that can be directly integrated to calculate the cumulative levels of activation. The proposed model reverts to the power law under constant shear stress conditions and is able to describe experimental results in response to a diverse range of highly dynamic stress conditions found in blood recirculating devices. The model was tested in vitro under emulated device flow conditions and correlates well with experimental results. This new model provides a reliable and robust mathematical tool that can be incorporated into computational fluid dynamic studies in order to optimize design, with the goal of improving the thrombogenic performance of blood recirculating devices.     

Validation of a Fluid–Structure Interaction Model of a Heart Valve using the Dynamic Mesh Method in Fluent

Simulations of coupled problems such as fluid–structure interaction (FSI) are becoming more and more important for engineering purposes. This is particularly true when modeling the aortic valve, where the FSI between the blood and the valve determines the valve movement and the valvular hemodynamics. Nevertheless only a few studies are focusing on the opening and closing behavior during the ejection phase (systole). In this paper, we present the validation of a FSI model using the dynamic mesh method of Fluent for the two-dimensional (2D) simulation of mechanical heart valves during the ejection phase of the cardiac cycle. The FSI model is successfully validated by comparing simulation results to experimental data obtained from in vitro studies using a CCD camera.     

Fluid-structure interaction simulation of artificial textile reinforced aortic heart valve: Validation with an in-vitro test

Prosthetic heart valves deployed in the left heart (aortic and mitral) are subjected to harsh hemodynamical conditions. Most of the tissue engineered heart valves have been developed for the low pressure pulmonary position because of the difficulties in fabricating a mechanically strong valve, able to withstand the systemic circulation. This necessitates the use of reinforcing scaffolds, resulting in a tissue-engineered textile reinforced tubular aortic heart valve. Therefore, to better design these implants, material behaviour of the composite, valve kinematics and its hemodynamical response need to be evaluated. Experimental assessment can be immensely time consuming and expensive, paving way for numerical studies. In this work, the material properties obtained using the previously proposed multi-scale numerical method for textile composites was evaluated for its accuracy. An in silico immersed boundary (IB) fluid structure interaction (FSI) simulation emulating the in vitro experiment was set-up to evaluate and compare the geometric orifice area and flow rate for one beat cycle. Results from the in silico FSI simulation were found to be in good coherence with the in vitro test during the systolic phase, while mean deviation of approximately 9% was observed during the diastolic phase of a beat cycle. Merits and demerits of the in silico IB-FSI method for the presented case study has been discussed with the advantages outweighing the drawbacks, indicating the potential towards an effective use of this framework in the development and analysis of heart valves.     

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.     

Transplantation of Pulmonary Valve Using a Mouse Model of Heterotopic Heart Transplantation

Tissue engineered heart valves, especially decellularized valves, are starting to gain momentum in clinical use of reconstructive surgery with mixed results. However, the cellular and molecular mechanisms of the neotissue development, valve thickening, and stenosis development are not researched extensively. To answer the above questions, we developed a murine heterotopic heart valve transplantation model. A heart valve was harvested from a valve donor mouse and transplanted to a heart donor mouse. The heart with a new valve was transplanted heterotopically to a recipient mouse. The transplanted heart showed its own heartbeat, independent of the recipient’s heartbeat. The blood flow was quantified using a high frequency ultrasound system with a pulsed wave Doppler. The flow through the implanted pulmonary valve showed forward flow with minimal regurgitation and the peak flow was close to 100 mm/sec. This murine model of heart valve transplantation is highly versatile, so it can be modified and adapted to provide different hemodynamic environments and/or can be used with various transgenic mice to study neotissue development in a tissue engineered heart valve.     

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.     

Bioprosthetic heart valve leaflet motion monitored by dual camera stereo photogrammetry

Dual camera stereo photogrammetry (DCSP) was applied to investigate the leaflet motion of bioprosthetic heart valves (BHVs) in a physiologic pulse flow loop (PFL). A 25-mm bovine pericardial valve was installed in the aortic valve position of the PFL, which was operated at a pulse rate of 70 beats/min and a cardiac output of 5 l/min. The systolic/diastolic aortic pressure was maintained at 120/80 mmHg to mimic the physiologic load experienced by the aortic valve. The leaflet of the test valve was marked with 80 India ink dots to form a fan-shaped matrix. From the acquired image sequences, 3-D coordinates of the marker matrix were derived and hence the surface contour, local mean and Gaussian curvatures at each opening and closing phase during one cardiac cycle were reconstructed. It is generally believed that the long-term failure rate of BHV is related to the uneven distribution of mechanical stresses occurring in the leaflet material during opening and closing. Unfortunately, a quantitative analysis of the leaflet motion under physiological conditions has not been reported. The newly developed technique permits frame-by-frame mapping of the leaflet surface, which is essential for dynamic analysis of stress-strain behavior in BHV.     

Two-dimensional color-mapping of turbulent shear stress distribution downstream of two aortic bioprosthetic valves in vitro.

Since artificial heart valve related complications such as thrombus formation, hemolysis and calcification are considered related to flow disturbances caused by the inserted valve, a thorough hemodynamic characterization of heart valve prostheses is essential. In a pulsatile flow model, fluid velocities were measured one diameter downstream of a Hancock Porcine (HAPO) and a Ionescu-Shiley Pericardial Standard (ISPS) aortic valve. Hot-film anemometry (HFA) was used for velocity measurements at 41 points in the cross-sectional area of the ascending aorta. Three-dimensional visualization of the velocity profiles, at 100 different instants during one mean pump cycle, was performed. Turbulence analysis was performed as a function of time by calculating the axial turbulence energy within 50 ms overlapping time windows during the systole. The turbulent shear stresses were estimated by using the correlation equation between Reynolds normal stress and turbulent (Reynolds) shear stress. The turbulent shear stress distribution was visualized by two-dimensional color-mapping at different instants during one mean pump cycle. Based on the velocity profiles and the turbulent shear stress distribution, a relative blood damage index (RBDI) was calculated. It has the feature of combining the magnitude and exposure time of the estimated shear stresses in one index, covering the entire cross-sectional area. The HAPO valve showed a skewed jet-type velocity profile with the highest velocities towards the left posterior aortic wall. The ISPS valve revealed a more parabolic-shaped velocity profile during systole. The turbulent shear stresses were highest in areas of high or rapidly changing velocity gradients. For the HAPO valve the maximum estimated turbulent shear stress was 194 N m-2 and for the ISPS valve 154 Nm-2. The RBDI was the same for the two valves. The turbulent shear stresses had magnitudes and exposure times that might cause endothelial damage and sublethal or lethal damage to blood corpuscules. The RBDI makes comparison between different heart valves easier and may prove important when making correlation with clinical observations.     

Simulated bioprosthetic heart valve deformation under quasi-static loading

For more than 40 years, the replacement of diseased natural heart valves with prosthetic devices has dramatically extended the quality and length of the lives of millions of patients worldwide. However, bioprosthetic heart valves (BHV) continue to fail due to structural failure resulting from poor tissue durability and faulty design. Clearly, an in-depth understanding of the biomechanical behavior of BHV at both the tissue and functional prosthesis levels is essential to improving BHV design and to reduce rates of failure. In this study, we simulated quasi-static BHV leaflet deformation under 40, 80, and 120 mm Hg quasi-static transvalvular pressures. A Fung-elastic material model was used that incorporated material parameters and axes derived from actual leaflet biaxial tests and measured leaflet collagen fiber structure. Rigorous experimental validation of predicted leaflet strain field was used to validate the model results. An overall maximum discrepancy of 2.36% strain between the finite element (FE) results and experiment measurements was obtained, indicating good agreement between computed and measured major principal strains. Parametric studies utilizing the material parameter set from one leaflet for all three leaflets resulted in substantial variations in leaflet stress and strain distributions. This result suggests that utilization of actual leaflet material properties is essential for accurate BHV FE simulations. The present study also underscores the need for rigorous experimentation and accurate constitutive models in simulating BHV function and design.