Experimental Hemodynamics Within the Penn State Fontan Circulatory Assist Device.

For children born with a single functional ventricle, the Fontan operation bypasses the right ventricle by forming a four-way total cavopulmonary connection and adapts the existing ventricle for the systemic circulation. However, upon reaching adulthood, many Fontan patients exhibit low cardiac output and elevated venous pressure, eventually requiring a heart transplantation. Despite efforts in developing a new device or using an existing device for failing Fontan support, there is still no Food and Drug Administration-approved device for subpulmonary support. Penn State University is developing a hydrodynamically levitated Fontan circulatory assist device (FCAD) for bridge-to-transplant or destination therapy. The hemodynamics within the FCAD, at both steady and patient averaged pulsatile conditions for three physiological pump operating conditions, were quantified using particle image velocimetry (PIV) to determine the velocity magnitudes and Reynolds normal and shear stresses within the device. Data were acquired at three planes (0 mm and ±25% of the radius) for the inferior and superior vena cavae inlets and the pulmonary artery outlet. The inlets had a blunt velocity profile that became skewed toward the collecting volute as fluid approached the rotor. At the outlet, regardless of the flow condition, a high-velocity jet exited the volute and moved downstream in a helical pattern. Turbulent stresses observed at the volute exit were influenced by the rotor's rotation. Regardless of inlet conditions, the pump demonstrated advantageous behavior for clinical use with a predictable flow field and a low risk of platelet adhesion and hemolysis based on calculated wall shear rates and turbulent stresses, respectively.

Acquired Von Willebrand Syndrome and Blood Pump Design.

The existence of acquired von Willebrand syndrome (AVWS) in patients with continuous flow left ventricular assist devices (LVADs) is well documented and has been verified by numerous investigators. AVWS has not been observed to occur in pulsatile devices such as the SynCardia total artificial heart (TAH), the HeartMate XVE, and the Thoratec pulsatile ventricular assist device (PVAD) used as a single pump. AVWS can also occur in patients with aortic stenosis, ventricular septal defect, mitral stenosis, and patent ductus arteriosus. It has been experimentally verified that supraphysiologic shear stress that occurs under these conditions can cleave the von Willebrand molecule, but the critical magnitude of stress and duration is unclear. Limited experimental results demonstrate that shear stresses as low as 5 Pa (50 dyne/cm2 ) can cause cleavage. Stresses in current centrifugal pumps can be as high as two orders of magnitude greater than this value. Pulsatile LVADs have stresses almost two orders of magnitude less than continuous flow LVADs. In order to improve continuous flow LVADs, the challenge for designers is to first determine the magnitude and duration of stress that is causing AVWS and then, if possible, design a pump below these stresses.

Characterizing the HeartMate II Left Ventricular Assist Device Outflow Using Particle Image Velocimetry.

Ventricular assist devices (VADs) are implanted in patients with a diseased ventricle to maintain peripheral perfusion as a bridge-to-transplant or as destination therapy. However, some patients with continuous flow VADs (e.g., HeartMate II (HMII)) have experienced gastrointestinal (GI) bleeding, in part caused by the proteolytic cleavage or mechanical destruction of von Willebrand factor (vWF), a clotting glycoprotein. in vitro studies were performed to measure the flow located within the HMII outlet cannula under both steady and physiological conditions using particle image velocimetry (PIV). Under steady flow, a mock flow loop was used with the HMII producing a flow rate of 3.2 L/min. The physiological experiment included a pulsatile pump operated at 105 BPM with a ventricle filling volume of 50 mL and in conjunction with the HMII producing a total flow rate of 5.0 L/min. Velocity fields, Reynolds normal stresses (RNSs), and Reynolds shear stresses (RSSs) were analyzed to quantify the outlet flow's potential contribution to vWF degradation. Under both flow conditions, the HMII generated principal Reynolds stresses that are, at times, orders of magnitude higher than those needed to unfurl vWF, potentially impacting its physiological function. Under steady flow, principal RNSs were calculated to be approximately 500 Pa in the outlet cannula. Elevated Reynolds stresses were observed throughout every phase of the cardiac cycle under physiological flow with principal RNSs approaching 1500 Pa during peak systole. Prolonged exposure to these conditions may lead to acquired von Willebrand syndrome (AvWS), which is accompanied by uncontrollable bleeding episodes.

Platelet adhesion to polyurethane urea under pulsatile flow conditions.

Platelet adhesion to a polyurethane urea surface is a precursor to thrombus formation within blood-contacting cardiovascular devices, and platelets have been found to adhere strongly to polyurethane surfaces below a shear rate of approximately 500 s(-1). The aim of the current work is to determine the properties of platelet adhesion to the polyurethane urea surface as a function of time-varying shear exposure. A rotating disk system was used to study the influence of steady and pulsatile flow conditions (e.g., cardiac inflow and sawtooth waveforms) for platelet adhesion to the biomaterial surface. All experiments were conducted with the same root mean square angular rotation velocity (29.63 rad/s) and waveform period. The disk was rotated in platelet-rich bovine plasma for 2 h, with adhesion quantified by confocal microscopy measurements of immunofluorescently labeled bovine platelets. Platelet adhesion under pulsating flow was found to decay exponentially with increasing shear rate. Adhesion levels were found to depend upon peak platelet flux and shear rate, regardless of rotational waveform. In combination with flow measurements, these results may be useful for predicting regions susceptible to thrombus formation within ventricular assist devices.

In vitro quantification of time dependent thrombus size using magnetic resonance imaging and computational simulations of thrombus surface shear stresses.

Thrombosis and thromboembolization remain large obstacles in the design of cardiovascular devices. In this study, the temporal behavior of thrombus size within a backward-facing step (BFS) model is investigated, as this geometry can mimic the flow separation which has been found to contribute to thrombosis in cardiac devices. Magnetic resonance imaging (MRI) is used to quantify thrombus size and collect topographic data of thrombi formed by circulating bovine blood through a BFS model for times ranging between 10 and 90 min at a constant upstream Reynolds number of 490. Thrombus height, length, exposed surface area, and volume are measured, and asymptotic behavior is observed for each as the blood circulation time is increased. Velocity patterns near, and wall shear stress (WSS) distributions on, the exposed thrombus surfaces are calculated using computational fluid dynamics (CFD). Both the mean and maximum WSS on the exposed thrombus surfaces are much more dependent on thrombus topography than thrombus size, and the best predictors for asymptotic thrombus length and volume are the reattachment length and volume of reversed flow, respectively, from the region of separated flow downstream of the BFS.

Effect of Hematocrit and Elevated Beat Rate on the 12cc Penn State Pediatric Ventricular Assist Device.

Congenital heart disease affects approximately 40,000 infants annually in the United States with 25% requiring invasive treatment. Due to limited number of donor hearts and treatment options available for children, pediatric ventricular assist devices (PVADs) are used as a bridge to transplant. The 12cc pneumatic Penn State PVAD is optimized to prevent platelet adhesion and thrombus formation at patient nominal conditions; however, children demonstrate variable blood hematocrit and elevated heart rates. Therefore, with pediatric patients exhibiting greater variability, particle image velocimetry is used to evaluate the PVAD with three non-Newtonian hematocrit blood analogs (20%, 40%, and 60%) and at two beat rates (75 and 120 bpm) to understand the device's performance. The flow fields demonstrate a strong inlet jet that transitions to a solid body rotation during diastole. During systole, the rotation dissipates and reorganizes into an outlet jet. This flow field is consistent across all hematocrits and beat rates but at a higher velocity magnitude during 120 bpm. There are also minor differences in flow field timing and surface washing due to hematocrit. Therefore, despite patient differences in hematocrit or required pumping output, thorough surface washing can be achieved in the PVAD by altering operating conditions, thus reducing platelet adhesion potential.

An in vitro fluid dynamic study of pediatric cannulae: the value of animal studies to predict human flow.

A challenge to the development of pediatric ventricular assist devices (PVADs) is the use of the aortic cannulae attached to the devices. Cannulae used for pediatric application have small diameters and large pressure drops. Furthermore, during the development of the 12cc Penn State pediatric PVAD, particle image velocimetry (PIV) illustrated that hematocrit levels, through changes in blood viscoelasticity, affected the fluid dynamics. The objective of this study is to compare the fluid dynamics of a pediatric viscoelastic blood analog and a goat viscoelastic blood analog within the PVAD aortic cannula. Two acrylic models were manufactured to model the aortic cannula (6 mm and 8 mm diameters). PIV data was collected to examine the flow at the outlet of the VAD and in the aortic cannula at heart rates of 50 and 75 beats per minute (bpm). Three planes of data were taken, one at the centerline and two 1.5 mm above and below the centerline. Three more planes of data were taken orthogonal to the original planes. While a 75 bpm heart rate was used to represent normal operating conditions, a 50 bpm heart rate represented use of the PVAD during weaning. At 75 bpm, differences were evident between the two different fluids and the two models. Separation zones developed in the plane below the centerline for the higher hematocrit pediatric blood analog. This study raises question to the usefulness of animal testing results in regard to how well they predict the outcome of pediatric patients.

A fluid dynamics study in a 50 cc pulsatile ventricular assist device: influence of heart rate variability.

Although left ventricular assist devices (LVADs) have had success in supporting severe heart failure patients, thrombus formation within these devices still limits their long term use. Research has shown that thrombosis in the Penn State pulsatile LVAD, on a polyurethane blood sac, is largely a function of the underlying fluid mechanics and may be correlated to wall shear rates below 500 s(-1). Given the large range of heart rate and systolic durations employed, in vivo it is useful to study the fluid mechanics of pulsatile LVADs under these conditions. Particle image velocimetry (PIV) was used to capture planar flow in the pump body of a Penn State 50 cubic centimeters (cc) LVAD for heart rates of 75-150 bpm and respective systolic durations of 38-50%. Shear rates were calculated along the lower device wall with attention given to the uncertainty of the shear rate measurement as a function of pixel magnification. Spatial and temporal shear rate changes associated with data collection frequency were also investigated. The accuracy of the shear rate calculation improved by approximately 40% as the resolution increased from 35 to 12 µm/pixel. In addition, data collection in 10 ms, rather than 50 ms, intervals was found to be preferable. Increasing heart rate and systolic duration showed little change in wall shear rate patterns, with wall shear rate magnitude scaling by approximately the kinematic viscosity divided by the square of the average inlet velocity, which is essentially half the friction coefficient. Changes in in vivo operating conditions strongly influence wall shear rates within our device, and likely play a significant role in thrombus deposition. Refinement of PIV techniques at higher magnifications can be useful in moving towards better prediction of thrombosis in LVADs.

Multilaboratory particle image velocimetry analysis of the FDA benchmark nozzle model to support validation of computational fluid dynamics simulations.

This study is part of a FDA-sponsored project to evaluate the use and limitations of computational fluid dynamics (CFD) in assessing blood flow parameters related to medical device safety. In an interlaboratory study, fluid velocities and pressures were measured in a nozzle model to provide experimental validation for a companion round-robin CFD study. The simple benchmark nozzle model, which mimicked the flow fields in several medical devices, consisted of a gradual flow constriction, a narrow throat region, and a sudden expansion region where a fluid jet exited the center of the nozzle with recirculation zones near the model walls. Measurements of mean velocity and turbulent flow quantities were made in the benchmark device at three independent laboratories using particle image velocimetry (PIV). Flow measurements were performed over a range of nozzle throat Reynolds numbers (Re(throat)) from 500 to 6500, covering the laminar, transitional, and turbulent flow regimes. A standard operating procedure was developed for performing experiments under controlled temperature and flow conditions and for minimizing systematic errors during PIV image acquisition and processing. For laminar (Re(throat)=500) and turbulent flow conditions (Re(throat)≥3500), the velocities measured by the three laboratories were similar with an interlaboratory uncertainty of ∼10% at most of the locations. However, for the transitional flow case (Re(throat)=2000), the uncertainty in the size and the velocity of the jet at the nozzle exit increased to ∼60% and was very sensitive to the flow conditions. An error analysis showed that by minimizing the variability in the experimental parameters such as flow rate and fluid viscosity to less than 5% and by matching the inlet turbulence level between the laboratories, the uncertainties in the velocities of the transitional flow case could be reduced to ∼15%. The experimental procedure and flow results from this interlaboratory study (available at http://fdacfd.nci.nih.gov) will be useful in validating CFD simulations of the benchmark nozzle model and in performing PIV studies on other medical device models.

Transcatheter Heart Valve Downstream Fluid Dynamics in an Accelerated Evaluation Environment.

Transcatheter aortic valve replacements (TAVRs) provide minimally invasive delivery of bioprosthetic heart valves (BHVs) for the treatment of aortic valve disease. While surgical BHVs show efficacy for 8-10 years, long-term TAVR durability remains unknown. Pre-clinical testing evaluates BHV durability in an ISO:5840 compliant accelerated wear tester (AWT), yet, the design and development of AWTs and their accuracy in predicting in vivo performance, is unclear. As a result of limited knowledge on AWT environment and BHV loading, durability assessment of candidate valves remains fundamentally empirical. For the first time, high-speed particle image velocimetry quantified an ISO:5840 compliant downstream AWT velocity field, Reynolds stresses, and turbulence intensity. TAVR enface imaging quantified the orifice area and estimated the flow rate. When valve area and flow rate were at their maximum during peak systole (1.49 cm2 and 16.05 L/min, respectively), central jet velocity, Reynolds normal and shear stress, and turbulence intensity grew to 0.50 m/s, 265.1 Pa, 124.6 Pa, and 37.3%, respectively. During diastole, unique AWT recirculation produced retrograde flow and the directional changes created eddies. These novel AWT findings demonstrated a substantially reduced valve fully loaded period and pressure not matching in vivo or in vitro studies, despite the comparable fluid environment and TAVR motion.

The influence of operational protocol on the fluid dynamics in the 12 cc Penn state pulsatile pediatric ventricular assist device: the effect of end-diastolic delay.

The success of adult ventricular assist devices (VADs), coupled with the high transplant waiting list mortality of infants (40%) has prompted Penn State to develop a pediatric version of the clinically successful adult device. Although the primary use of this device will be bridge-to-transplant, there has been sufficient clinical data to demonstrate the efficacy of VADs in a bridge-to-recovery setting. However, removing the patient from the device, a process known as weaning, demands operation of the device at a lower beat rate and concomitant increased risk for thromboembolism. Previous studies have shown that the interrelated flow characteristics necessary for the prevention of thrombosis in a pulsatile VAD are a strong inlet jet, a late diastolic recirculating flow, and a wall shear rate greater than 500/s. In an effort to develop a strong inlet jet and rotational flow pattern at a lower beat and flow rate, we have compressed diastole by altering the end-diastolic delay time (EDD). Particle image velocimetry was used to compare the flow fields and wall shear rates in the chamber of the 12 cc Penn State pulsatile pediatric VAD operated at 50 beats per minute using EDDs of 10, 50, and 100 ms. Although we expected the 100 ms EDD to have the best wall shear profiles, we found that the 50 ms EDD condition was superior to both the 10 and 100 EDD conditions, due to a longer sustained inlet jet.

Hemodynamics of an end-to-side anastomotic graft for a pulsatile pediatric ventricular assist device.

Numerical simulations are performed to investigate the flow within the end-to-side proximal anastomosis of a pulsatile pediatric ventricular assist device (PVAD) to an aorta. The anastomotic model is constructed from a patient-specific pediatric aorta. The three great vessels originating from the aortic arch--brachiocephalic (innominate), left common carotid, and left subclavian arteries--are included. An implicit large eddy simulation method based on a finite volume approach is used to study the resulting turbulent flow. A resistance boundary condition is applied at each branch outlet to study flow splitting. The PVAD anastomosis is found to alter the aortic flow dramatically. More flow is diverted into the great vessels with the PVAD support. Turbulence is found in the jet impingement area at peak systole for 100% bypass, and a maximum principal normal Reynolds stress of 7081 dyn/cm(2) is estimated based on ten flow cycles. This may be high enough to cause hemolysis and platelet activation. Regions prone to intimal hyperplasia are identified by combining the time-averaged wall shear stress and oscillatory shear index. These regions are found to vary, depending on the percentage of the flow bypass.

Impact of design parameters on bileaflet mechanical heart valve flow dynamics.

BACKGROUND AND AIM OF THE STUDY: One significant problem encountered during surgery to implant mechanical heart valve prostheses is the propensity for thrombus formation near the valve leaflet and housing. This may be caused by the high shear stresses present in the leakage jet flows through small gaps between leaflets and the valve housing during the valve closure phase. METHODS: A two-dimensional (2D) study was undertaken to demonstrate that design changes in bileaflet mechanical valves result in notable changes in the flow-induced stresses and prediction of platelet activation. A Cartesian grid technique was used for the 2D simulation of blood flow through two models of bileaflet mechanical valves, and their flow patterns, closure characteristics and platelet activation potential were compared. A local mesh refinement algorithm allowed efficient and fast flow computations with mesh adaptation based on the gradients of the flow field in the gap between the leaflet and housing at the instant of valve closure. Leaflet motion was calculated dynamically, based on the fluid forces acting on it. Platelets were modeled and tracked as point particles by a Lagrangian particle tracking method which incorporated the hemodynamic forces on the particles. RESULTS: A comparison of results showed that the velocity, wall shear stress and simulated platelet activation parameter were lower in the valve model, with a smaller angle of leaflet traverse between the fully open to the fully closed position. The parameters were also affected to a lesser extent by local changes in the leaflet and housing geometry. CONCLUSION: Computational simulations can be used to examine local design changes to help minimize the fluid-induced stresses that may play a key role in thrombus initiation with the implanted mechanical valves.

Numerical study of blood flow at the end-to-side anastomosis of a left ventricular assist device for adult patients.

We use an implicit large eddy simulation (ILES) method based on a finite volume approach to capture the turbulence in the anastomoses of a left ventricular assist device (LVAD) to the aorta. The order-of-accuracy of the numerical schemes is computed using a two-dimensional decaying Taylor-Green vortex. The ILES method is carefully validated by comparing to documented results for a fully developed turbulent channel flow at Re(tau)=395. Two different anastomotic flows (proximal and distal) are simulated for 50% and 100% LVAD supports and the results are compared with a healthy aortic flow. All the analyses are based on a planar aortic model under steady inflow conditions for simplification. Our results reveal that the outflow cannulae induce high exit jet flows in the aorta, resulting in turbulent flow. The distal configuration causes more turbulence in the aorta than the proximal configuration. The turbulence, however, may not cause any hemolysis due to low Reynolds stresses and relatively large Kolmogorov length scales compared with red blood cells. The LVAD support causes an acute increase in flow splitting in the major branch vessels for both anastomotic configurations, although its long-term effect on the flow splitting remains unknown. A large increase in wall shear stress is found near the cannulation sites during the LVAD support. This work builds a foundation for more physiologically realistic simulations under pulsatile flow conditions.

Flow behavior within the 12-cc Penn State pulsatile pediatric ventricular assist device: an experimental study of the initial design.

Planar particle image velocimetry was used to explore the flow behavior of the newly designed 12-cc Penn State pneumatic pediatric assist pump. Wall shear maps complemented the velocity data. Bjork-Shiley Monostrut 17-mm mechanical heart valves were used in the inlet and outlet ports. In comparison with larger Penn State pumps, the 12-cc device is not only smaller but has reduced valve effective orifice areas and more highly angled valve ports. In contrast to results from the larger pumps, the flow field was highly three dimensional during early diastole with poorer penetration by the valve inlet jet. This led to a later start to a "wall washing" rotational pattern. A significant separation region, never before observed, was created upstream of the outlet valve leaflet during late diastole--effectively reducing the area and increasing the pressure drop through the valve. Wall shear maps suggest that regions of low shear might persist throughout the cycle at the bottom of the pump on the outlet side. An attempt to improve the flow field characteristics by exploring different valves, valve orientations and inlet valve angles, systolic/diastolic flow timing, and perhaps a larger outlet valve was planned.

An enlarging pacemaker pocket: A case report of a plasmablastic lymphoma arising as a primary tumor around a cardiac pacemaker and systematic literature review of various malignancies arising at the pacemaker pocket.

To date, there have been limited reports of oncogenesis occurring within pacemaker pockets. We report the case of a 100-year-old male who presented to the emergency department complaining of expansion of his pacemaker pocket over the period of 8 days. Dissection of the pacemaker pocket and pathological analysis of tissue samples revealed plasmablastic lymphoma, a subset of diffuse large B-cell lymphoma, commonly seen in immunocompromised elderly patients. This is the first known reported case of plasmablastic lymphoma occurring within a pacemaker pocket. .

Asynchronous Pumping of a Pulsatile Ventricular Assist Device in a Pediatric Anastomosis Model.

BACKGROUND: Both pulsatile and continuous flow ventricular assist devices are being developed for pediatric congenital heart defect patients. Pulsatile devices are often operated asynchronously with the heart in either an "automatic" or a fixed beat rate mode. However, most studies have only investigated synchronized ejection. METHODS: A previously validated viscoelastic blood solver is used to investigate the parameters of pulsatility, power loss, and graft failure in a pediatric aortic anastomosis model. RESULTS: Pulsatility was highest with synchronized flow and lowest at a 90° phase shift. Power loss decreased at 90° and 180° phase shifts but increased at a 270° phase shift. Similar regions of potential intimal hyperplasia and graft failure were seen in all cases but with phase-shifted ejection leading to higher wall shear stress on the anastomotic floor and oscillatory shear index on the anastomotic toe. CONCLUSION: The ranges of pulsatility and hemodynamics that can result clinically using asynchronous pulsatile devices were investigated in a pediatric anastomosis model. These results, along with the different postoperative benefits of pump modulation, can be used to design an optimal weaning protocol.

Development of a platelet adhesion transport equation for a computational thrombosis model.

Thrombosis is a significant issue for cardiovascular device development and use. While thrombosis models are available, very few are device-related and none have been thoroughly validated experimentally. Here, we introduce a surface adherent platelet transport equation into a continuum model to account for the biomaterial interface/blood interaction. Using a rotating disc system and polyurethane-urea material, we characterize steady and pulsatile flow fields using laser Doppler velocimetry. In vitro measurements of platelet adhesion are used in combination with the LDV data to provide further experimental validation. The rotating disc system is computationally studied using the device-induced thrombosis model with the surface platelet adherent transport equation. The results indicate that the flow field is in excellent agreement to the experimental LDV data and that the platelet adhesion simulations are in good agreement with the in vitro platelet data. These results provide good evidence that this transport equation can be used to express the relationship between blood and a biomaterial if the correct platelet adhesion characteristics are known for the biomaterial. Further validation is necessary with other materials.

FDA Benchmark Medical Device Flow Models for CFD Validation.

Computational fluid dynamics (CFD) is increasingly being used to develop blood-contacting medical devices. However, the lack of standardized methods for validating CFD simulations and blood damage predictions limits its use in the safety evaluation of devices. Through a U.S. Food and Drug Administration (FDA) initiative, two benchmark models of typical device flow geometries (nozzle and centrifugal blood pump) were tested in multiple laboratories to provide experimental velocities, pressures, and hemolysis data to support CFD validation. In addition, computational simulations were performed by more than 20 independent groups to assess current CFD techniques. The primary goal of this article is to summarize the FDA initiative and to report recent findings from the benchmark blood pump model study. Discrepancies between CFD predicted velocities and those measured using particle image velocimetry most often occurred in regions of flow separation (e.g., downstream of the nozzle throat, and in the pump exit diffuser). For the six pump test conditions, 57% of the CFD predictions of pressure head were within one standard deviation of the mean measured values. Notably, only 37% of all CFD submissions contained hemolysis predictions. This project aided in the development of an FDA Guidance Document on factors to consider when reporting computational studies in medical device regulatory submissions. There is an accompanying podcast available for this article. Please visit the journal's Web site (www.asaiojournal.com) to listen.

The 50cc Penn State left ventricular assist device: a parametric study of valve orientation flow dynamics.

We investigated the flow fields associated with the Bjork-Shiley Monostrut mechanical heart valve in the mitral position of the 50 cc Penn State left ventricular assist device. The valve orientation was adjusted whereby flow field data was collected using planar particle image velocimetry. The mitral valve was rotated from 0 to 45 degrees, in 15-degree increments. For each valve orientation, measurements were made in three planes (3, 5, and 8 mm from the front wall) parallel to the device pusher plate. Penetration of the inlet jet was affected by the valve orientation with more intense and longer duration wall washing motion occurring at 45 degrees. As a result, the 45-degree mitral valve orientation is recommended to help prevent areas of thrombus deposition. Valve orientation is an important aspect of assist device design.