Fluid dynamics of a pediatric ventricular assist device

Hemocompatibility and biophysical interface of left ventricular assist devices and total artificial hearts

ABSTRACT

Over the past 2 decades, there has been a significant increase in the utilization of long-term mechanical circulatory support (MCS) for the treatment of cardiac failure. Left ventricular assist devices (LVADs) and total artificial hearts (TAHs) have been developed in parallel to serve as bridge-to-transplant and destination therapy solutions. Despite the distinct hemodynamic characteristics introduced by LVADs and TAHs, a comparative evaluation of these devices regarding potential complications in supported patients, has not been undertaken. Such a study could provide valuable insights into the complications associated with these devices. Although MCS has shown substantial clinical benefits, significant complications related to hemocompatibility persist, including thrombosis, recurrent bleeding, and cerebrovascular accidents. This review focuses on the current understanding of hemostasis, specifically thrombotic and bleeding complications, and explores the influence of different shear stress regimens in long-term MCS. Furthermore, the role of endothelial cells in protecting against hemocompatibility-related complications of MCS is discussed. We also compared the diverse mechanisms contributing to the occurrence of hemocompatibility-related complications in currently used LVADs and TAHs. By applying the existing knowledge, we present, for the first time, a comprehensive comparison between long-term MCS options.

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.

Effect of the bileaflet inlet valve angle on the flow of a pediatric ventricular assist device: experimental analysis

BACKGROUND

Mechanical heart valves (MHV) and its fluid dynamics inside a pulsatile pediatric ventricular assist device (PVAD) can be associated with blood degradation. In this paper, flow structures are analyzed and compared by an experimental investigation on the effect of bileaflet MHV positioned at varying angles in the inlet port orifice of a PVAD.

METHODS

Time-resolved particle image velocimetry (TR-PIV) was applied to characterize the internal flow of the device. St Jude Medical bileaftlet valves were used on the inlet orifice and positioned at 0°, 15°, 30°, 45°, 60°, and 90° in relation to the centerline of the device. Three planes with bidimensional velocity magnitude fields were considered in the analysis with visualization of diastolic jets, device wall washing patterns and flow circulation during emptying or systole of the pump. Also, the washing vortex area, and vertical velocity probabilities of regurgitant flows in the inlet valve were evaluated.

RESULTS

The results show that a variation in the angle of the MHV at the inlet port produced distinct velocities, fluid structures, and regurgitant flow probabilities within the device. MHV positioned at an angle of 0° generated the strongest inlet jet, larger vortex area during filling, more prominent outgoing flow, and less regurgitation compared to the angles studied. The presence of unfavorable fluid structures, such as small vortices, and/or sudden flow structure interruption, and/or regurgitation, were identified at 45° and 90° angles.

CONCLUSIONS

The 0° inlet angle had better outcomes than other angles due to its consistency in the multiple parameters analyzed.

Novel use of cangrelor in pediatrics: A pilot cohort study demonstrating use in ventricular assist devices

Thromboembolic events and bleeding are major sources of morbidity among pediatric patients supported on a ventricular assist device (VAD). Pharmacokinetics and pharmacodynamics of enteral antiplatelet agents are affected and variable due to erratic enteral absorption in end-stage heart failure and VAD circulation. Additionally, 20%-40% of the population are poor metabolizers of clopidogrel, a prodrug, making cangrelor an alternative when antiplatelet therapy is crucial. Cangrelor has been used effectively and safely for short durations in adults during percutaneous coronary interventions, but the use of cangrelor is still under investigation in pediatrics. This case series utilized cangrelor, a novel short-acting, reversible, intravenous P2Y₁₂ platelet inhibitor in managing pediatric patients supported with a VAD. We performed a retrospective, single-center review of patients admitted to a tertiary medical center with end-stage heart failure requiring mechanical circulatory support and concomitant cangrelor administration between January 2019 and March 2020. Platelet function testing, cangrelor dose, bleeding complications, thromboembolic events, and frequency of circuit interventions during the use of cangrelor were recorded. Optimal platelet reactivity, defined as P2Y₁₂ < 180 platelet reaction units (PRU), was measured with serial point-of-care testing (VerifyNow). Seven patients, median age of 4.9 years, met the above criteria. Three patients had a diagnosis of complex congenital heart disease. Four patients had dilated or restrictive cardiomyopathy. All patients were on continuous flow VADs. The median VAD duration was 84.5 days (IQR 61.5-103). The median duration on cangrelor was 43 days (IQR 8-70). The median cangrelor dose to reach the therapeutic threshold was 0.75 μg/kg/min with the mean P2Y₁₂ , while on cangrelor of 164.75 PRU. Bleeding complications included mild gastrointestinal bleeding and hematuria. There was one patient with pump thrombosis requiring intervention. There were no cerebrovascular events while on cangrelor. We report the first successful long-term use of cangrelor in pediatric patients. The reversibility and short half-life of cangrelor make it a feasible antiplatelet agent in selected patients. This data supports the use of cangrelor in children as a viable antiplatelet option; with minimal bleeding complications and no cerebrovascular events demonstrated in this cohort.

A 4D flow MRI evaluation of the impact of shear-dependent fluid viscosity on in vitro Fontan circulation flow

The Fontan procedure for univentricular heart defects creates a nonphysiologic circulation where systemic venous blood drains directly into the pulmonary arteries, leading to multiorgan dysfunction secondary to chronic low-shear nonpulsatile pulmonary blood flow and central venous hypertension. Although blood viscosity increases exponentially in this low-shear environment, the role of shear-dependent ("non-Newtonian") blood viscosity in this pathophysiology is unclear. We studied three-dimensional (3D)-printed Fontan models in an in vitro flow loop with a Philips 3-T magnetic resonance imaging (MRI) scanner. A 4D flow phase-contrast sequence was used to acquire a time-varying 3D velocity field for each experimental condition. On the basis of blood viscosity of a cohort of patients who had undergone the Fontan procedure, it was decided to use 0.04% xanthan gum as a non-Newtonian blood analog; 45% glycerol was used as a Newtonian control fluid. MRI data were analyzed using GTFlow and MATLAB software. The primary outcome, power loss, was significantly higher with the Newtonian fluid [14.8 (13.3, 16.4) vs. 8.1 (6.4, 9.8)%, medians with 95% confidence interval, P < 0.0001]. The Newtonian fluid also demonstrated marginally higher right pulmonary artery flow, marginally lower shear stress, and a trend toward higher caval flow mixing. Outcomes were modulated by Fontan model complexity, cardiac output, and caval flow ratio. Vortexes, helical flow, and stagnant flow were more prevalent with the non-Newtonian fluid. Our data demonstrate that shear-dependent viscosity significantly alters qualitative flow patterns, power loss, pulmonary flow distribution, shear stress, and caval flow mixing in synthetic models of the Fontan circulation. Potential clinical implications include effects on exercise capacity, ventilation-perfusion matching, risk of pulmonary arteriovenous malformations, and risk of thromboembolism.NEW & NOTEWORTHY Although blood viscosity increases exponentially in low-shear environments, the role of shear-dependent ("non-Newtonian") blood viscosity in the pathophysiology of the low-shear Fontan circulation is unclear. We demonstrate that shear-dependent viscosity significantly alters qualitative flow patterns, power loss, pulmonary flow distribution, shear stress, and caval flow mixing in synthetic models of the Fontan circulation. Potential clinical implications include effects on exercise capacity, ventilation-perfusion matching, risk of pulmonary arteriovenous malformations, and risk of thromboembolism.

Experimental Investigation of the Effect of Non-Newtonian Behavior of Blood Flow in the Fontan Circulation

The Fontan procedure for univentricular heart defects creates a unique circulation where all pulmonary blood flow is passively supplied directly from systemic veins. Computational simulations, aimed at optimizing the surgery, have assumed blood to be a Newtonian fluid without evaluating the potential error introduced by this assumption. We compared flow behavior between a non-Newtonian blood analog (0.04% xanthan gum) and a control Newtonian fluid (45% glycerol) in a simplified model of the Fontan circulation. Particle image velocimetry was used to examine flow behavior at two different cardiac outputs and two caval blood flow distributions. Pressure and flow rates were measured at each inlet and outlet. Velocity, shear strain, and shear stress maps were derived from velocity data. Power loss was calculated from pressure, flow, and velocity data. Power loss was increased in all test conditions with xanthan gum vs. glycerol (mean 10±2.9% vs. 5.6±1.3%, p=0.032). Pulmonary blood flow distribution differed in all conditions, more so at low cardiac output. Caval blood flow mixing patterns and shear stress were also qualitatively different between the solutions in all conditions. We conclude that assuming blood to be a Newtonian fluid introduces considerable error into simulations of the Fontan circulation, where low-shear flow predominates.

Toward the Virtual Benchmarking of Pneumatic Ventricular Assist Devices: Application of a Novel Fluid-Structure Interaction-Based Strategy to the Penn State 12 cc Device

The pediatric use of pneumatic ventricular assist devices (VADs) as a bridge to heart transplant still suffers for short-term major complications such as bleeding and thromboembolism. Although numerical techniques are increasingly exploited to support the process of device optimization, an effective virtual benchmark is still lacking. Focusing on the 12 cc Penn State pneumatic VAD, we developed a novel fluid-structure interaction (FSI) model able to capture the device functioning, reproducing the mechanical interplay between the diaphragm, the blood chamber, and the pneumatic actuation. The FSI model included the diaphragm mechanical response from uniaxial tensile tests, realistic VAD pressure operative conditions from a dedicated mock loop system, and the behavior of VAD valves. Our FSI-based benchmark effectively captured the complexity of the diaphragm dynamics. During diastole, the initial slow diaphragm retraction in the air chamber was followed by a more rapid phase; asymmetries were noticed in the diaphragm configuration during its systolic inflation in the blood chamber. The FSI model also captured the major features of the device fluid dynamics. In particular, during diastole, a rotational wall washing pattern is promoted by the penetrating inlet jet with a low-velocity region located in the center of the device. Our numerical analysis of the 12 cc Penn State VAD points out the potential of the proposed FSI approach well resembling previous experimental evidences; if further tested and validated, it could be exploited as a virtual benchmark to deepen VAD-related complications and to support the ongoing optimization of pediatric devices.

Comparison of Blood Viscoelasticity in Pediatric and Adult Cardiac Patients

Evidence is accumulating that blood flow patterns in the cardiovascular system and in cardiovascular devices do, in some instances, depend on blood viscoelasticity. Thus, to better understand the challenges to providing circulatory support and surgical therapies for pediatric and adult patients, viscous and elastic components of complex blood viscoelasticity of 31 pediatric patients were compared to those of 29 adult patients with a Vilastic-3 rheometer. A random effects model with categorical age covariates found statistically significant differences between pediatric and adult patients for log viscosity (p = 0.005). Log strain (p < 0.0001) and hematocrit (p < 0.0001) effects were also significant, as were the hematocrit-by-log-strain (p = 0.0006) and age-by-log strain (p = 0.001) interactions. The hematocrit-by-age interaction was not significant. For log elasticity, age differences were insignificant (p = 0.39). The model for log elasticity had significant log strain (p < 0.0001), log strain squared (p < 0.0001) and hematocrit (p < 0.0001) effects, as well as hematocrit-by-log-strain and hematocrit-by-log-strain-squared interactions (p = 0.014). A model for log viscosity with continuous age was also fit to the data, which can be used to refine cardiovascular device design and operation to the age of the patient. We conclude that there are distinct differences between pediatric and adult blood viscosity, as well as substantial variation within the pediatric population, that may impact the performance of devices and procedures.

Hemodynamics in a Pediatric Ascending Aorta Using a Viscoelastic Pediatric Blood Model

Congenital heart disease is the leading cause of infant death in the United States with over 36,000 newborns affected each year. Despite this growing problem there are few mechanical circulatory support devices designed specifically for pediatric and neonate patients. Previous research has been done investigating pediatric ventricular assist devices (PVADs) assuming blood to be a Newtonian fluid in computational fluid dynamics (CFD) simulations, ignoring its viscoelastic and shear-thinning properties. In contrast to adult VADs, PVADs may be more susceptible to hemolysis and thrombosis due to altered flow into the aorta, and therefore, a more accurate blood model should be used. A CFD solver that incorporates a modified Oldroyd-B model designed specifically for pediatric blood is used to investigate important hemodynamic parameters in a pediatric aortic model under pulsatile flow conditions. These results are compared to Newtonian blood simulations at three physiological pediatric hematocrits. Minor differences are seen in both velocity and wall shear stress (WSS) during early stages of the cardiac cycle between the Newtonian and viscoelastic models. During diastole, significant differences are seen in the velocities in the descending aorta (up to 12%) and in the aortic branches (up to 30%) between the two models. Additionally, peak WSS differences are seen between the models throughout the cardiac cycle. At the onset of diastole, peak WSS differences of 43% are seen between the Newtonian and viscoelastic model and between the 20 and 60% hematocrit viscoelastic models at peak systole of 41%.

Assessment of turbulence models for pulsatile flow inside a heart pump

Computational fluid dynamics (CFD) is applied to study the unsteady flow inside a pulsatile pump left ventricular assist device, in order to assess the sensitivity to a range of commonly used turbulence models. Levels of strain and wall shear stress are directly relevant to the evaluation of risk from haemolysis and thrombosis, and thus understanding the sensitivity to these turbulence models is important in the assessment of uncertainty in CFD predictions. The study focuses on a positive displacement or pulsatile pump, and the CFD model includes valves and moving pusher plate. An unstructured dynamic layering method was employed to capture this cyclic motion, and valves were simulated in their fully open position to mimic the natural scenario, with in/outflow triggered at control planes away from the valves. Six turbulence models have been used, comprising three relevant to the low Reynolds number nature of this flow and three more intended to investigate different transport effects. In the first group, we consider the shear stress transport (SST) [Formula: see text] model in both its standard and transition-sensitive forms, and the 'laminar' model in which no turbulence model is used. In the second group, we compare the one equation Spalart-Almaras model, the standard two equation [Formula: see text] and the full Reynolds stress model (RSM). Following evaluation of spatial and temporal resolution requirements, results are compared with available experimental data. The model was operated at a systolic duration of 40% of the pumping cycle and a pumping rate of 86 BPM (beats per minute). Contrary to reasonable preconception, the 'transition' model, calibrated to incorporate additional physical modelling specifically for these flow conditions, was not noticeably superior to the standard form of the model. Indeed, observations of turbulent viscosity ratio reveal that the transition model initiates a premature increase of turbulence in this flow, when compared with both experimental and higher order numerical results previously reported in the literature. Furthermore, the RSM is indicated to provide the most accurate prediction over much of the flow, due to its ability to more correctly account for three-dimensional effects. Finally, the clinical relevance of the results is reported along with a discussion on the impact of such modelling uncertainties.

The influence of device position on the flow within the Penn State 12 cc pediatric ventricular assist device

Ventricular assist devices are a commonly used heart failure therapy for adult patients as bridge-to-transplant or bridge-to-recovery tools. The application of adult ventricular assist devices in pediatric patients has led to increased thrombotic events. Therefore, we have been developing a pediatric ventricular assist device (PVAD), the Penn State 12 cc PVAD. It is designed for patients with a body weight of 5-15 kg and has a stroke volume of 12 cc. Clot formation is the major concern. It is correlated to the coagulability of blood, the blood contacting materials and the fluid dynamics within the system. The intent is for the PVAD to be a long term therapy. Therefore, the system may be oriented in different positions according to the patient's behavior. This study evaluates for the first time the impact of position on the flow patterns within the Penn State 12 cc PVAD, which may help to improve the PVAD design concerning chamber and ports geometries. The fluid dynamics are visualized by particle image velocimetry. The evaluation is based on inlet jet behavior and calculated wall shear rates. Vertical and horizontal model orientations are compared, both with a beat rate of 75, outlet pressures of 90/60 mm Hg and a flow rate of 1.3 l/min. The results show a significant change of the inlet jet behavior and the development of a rotational flow pattern. Vertically, the inlet jet is strong along the wall. It initiates a rotational flow pattern with a wandering axis of rotation. In contrast, the horizontal model orientation results show a weaker inlet jet along the wall with a nearly constant center of rotation location, which can be correlated to a higher risk of thrombotic events. In addition, high speed videography illustrates differences in the diaphragm motion during diastole. Diaphragm opening trajectories measurements determine no significant impact of the density of the blood analog fluids. Hence, the results correlate to human blood.

A parametric study of valve orientation on the flow patterns of the Penn State pulsatile pediatric ventricular assist device

Because of the shortage of organs for transplant in pediatric patients with end-stage heart failure, Penn State is developing a pneumatically driven 12 cc pulsatile pediatric ventricular assist device (PVAD). A major concern is the flow field changes related to the volume decrease and its effect on device thrombogenicity. Previous studies of similar devices have shown that changes in the orientation of the inlet valve can lead to improvement in the flow field. Herein, the fluid dynamic effects of orientation changes at both the inlet and outlet valves are studied. Using two-dimensional particle image velocimetry, we examine the flow field in vitro using an acrylic model of the PVAD in a mock circulatory loop. Regardless of valve orientation, the overall flow pattern inside the PVAD remains similar, but important differences were seen locally in the wall shear rates, which is notable because shear rates >500 s may prevent thrombus formation. As the inlet valve was rotated toward the fluid side of the PVAD, we observed an increase in inlet jet velocity and wall shear rates along the inlet port wall. A corresponding rotation of the outlet valve increases the wall shear rate along the outer wall near the device outlet. Wall shear rates were all higher when both valves were rotated toward the fluid side of the device, with the best rates found at orientations of +15 degrees for both the inlet and outlet valves. Overall, orientations of +15 degrees or +30 degrees of both the inlet and outlet valve resulted in an increase in wall shear rates and could aid in the reduction of thrombus formation inside the PVAD.

Validation of a CFD methodology for positive displacement LVAD analysis using PIV data

Computational fluid dynamics (CFD) is used to asses the hydrodynamic performance of a positive displacement left ventricular assist device. The computational model uses implicit large eddy simulation direct resolution of the chamber compression and modeled valve closure to reproduce the in vitro results. The computations are validated through comparisons with experimental particle image velocimetry (PIV) data. Qualitative comparisons of flow patterns, velocity fields, and wall-shear rates demonstrate a high level of agreement between the computations and experiments. Quantitatively, the PIV and CFD show similar probed velocity histories, closely matching jet velocities and comparable wall-strain rates. Overall, it has been shown that CFD can provide detailed flow field and wall-strain rate data, which is important in evaluating blood pump performance.

Flow visualization of three-dimensionality inside the 12 cc Penn State pulsatile pediatric ventricular assist device

In order to aid the ongoing concern of limited organ availability for pediatric heart transplants, Penn State has continued development of a pulsatile Pediatric Ventricular Assist Device (PVAD). Initial studies of the PVAD observed an increase in thrombus formation due to differences in flow field physics when compared to adult sized devices, which included a higher degree of three-dimensionality. This unique flow field brings into question the use of 2D planar particle image velocimetry (PIV) as a flow visualization technique, however the small size and high curvature of the PVAD make other tools such as stereoscopic PIV impractical. In order to test the reliability of the 2D results, we perform a pseudo-3D PIV study using planes both parallel and normal to the diaphragm employing a mock circulatory loop containing a viscoelastic fluid that mimics 40% hematocrit blood. We find that while the third component of velocity is extremely helpful to a physical understanding of the flow, particularly of the diastolic jet and the development of a desired rotational pattern, the flow data taken parallel to the diaphragm is sufficient to describe the wall shear rates, a critical aspect to the study of thrombosis and design of such pumps.

Towards non-thrombogenic performance of blood recirculating devices

Implantable blood recirculating devices have provided life saving solutions to patients with severe cardiovascular diseases. However, common problems of hemolysis and thromboembolism remain an impediment to these devices. In this article, we present a brief review of the work by several groups in the field that has led to the development of new methodologies that may facilitate achieving the daunting goal of optimizing the thrombogenic performance of blood recirculating devices. The aim is to describe work which pertains to the interaction between flow-induced stresses and the blood constituents, and that supports the hypothesis that thromboembolism in prosthetic blood recirculating devices is initiated and maintained primarily by the non-physiological flow patterns and stresses that activate and enhance the aggregation of blood platelets, increasing the risk of thromboembolism and cardioembolic stroke. Such work includes state-of-the-art numerical and experimental tools used to elucidate flow-induced mechanisms leading to thromboembolism in prosthetic devices. Following the review, the paper describes several efforts conducted by some of the groups active in the field, and points to several directions that should be pursued in the future in order to achieve the goal for blood recirculating prosthetic devices becoming more effective as destination therapy in the future.

Flow visualization techniques in a mock ventricle supported by a nonpulsatile left ventricular assist device

Little is known about flow patterns in ventricles supported by continuous flow left ventricular assist devices (LVADs), and valuable information can be obtained with simple flow visualization experiments. We describe the application of several experimental techniques for the in vitro study of ventricular flow patterns (e.g., unsteadiness, vortical motions, stagnation regions) in the presence of a continuous flow LVAD. We used dye streaks, particle paths, and hydrogen bubble techniques to visualize fluid flow in an idealized, static, transparent mock ventricle attached to a Jarvik 2000 continuous flow LVAD. We recorded ventricular flow behavior at various pump speeds while independently adjusting pump flow (by varying the afterload) to emulate in vivo pump flow at various phases of the cardiac cycle. Changes in ventricular flow behavior at different pump flow rates may be of clinical relevance, because continuous flow pumps are extremely sensitive to inflow and outflow pressures and instantaneous pump flow varies significantly at different points throughout the cardiac cycle. Further work is needed to quantitatively compare the flow behavior of different continuous flow devices in a pulsatile ventricular model.

Micro-flow visualization of red blood cell-enhanced platelet concentration at sudden expansion

Microscopic steps and crevices are inevitable features within prosthetic blood-contacting devices. This study aimed to elucidate the thrombogenicity of the associated microscopic flow features by studying the transport of fluorescent platelet-sized particles in a suspension of red blood cells (RBCs) flowing through a 100 microm:200 microm sudden expansion. Micro-flow visualization revealed a strong influence of hematocrit upon the path of RBCs and spatial concentration of particles. At all flow rates studied (Re = 8.3-41.7) and hematocrit 20% and lower, RBC streamlines were found to detach from the microchannel wall creating an RBC-depleted zone inside the step that was much larger than the cells themselves. However, the observed distribution of particles was relatively homogeneous. By contrast, the RBC streamlines of samples with hematocrit equal to or greater than 30% more closely followed the contour of the microchannel, yet exhibited enhanced concentration of particles within the corner. The corresponding size of the cell depletion layer was comparable with the size of the cells. This study implies that local platelet concentration in blood within the physiological range of hematocrit can be elevated within the flow separation region of a sudden expansion and implicates the role of RBCs in causing this effect.

Assessment of hydraulic performance and biocompatibility of a MagLev centrifugal pump system designed for pediatric cardiac or cardiopulmonary support

The treatment of children with life-threatening cardiac and cardiopulmonary failure is a large and underappreciated public health concern. We have previously shown that the CentriMag is a magnetically levitated centrifugal pump system, having the utility for treating adults and large children (1,500 utilized worldwide). We present here the PediVAS, a pump system whose design was modified from the CentriMag to meet the physiological requirements of young pediatric and neonatal patients. The PediVAS is comprised of a single-use centrifugal blood pump, reusable motor, and console, and is suitable for right ventricular assist device (RVAD), left ventricular assist device (LVAD), biventricular assist device (BVAD), or extracorporeal membrane oxygenator (ECMO) applications. It is designed to operate without bearings, seals and valves, and without regions of blood stasis, friction, or wear. The PediVAS pump is compatible with the CentriMag hardware, although the priming volume was reduced from 31 to 14 ml, and the port size reduced from 3/8 to (1/4) in. For the expected range of pediatric flow (0.3-3.0 L/min), the PediVAS exhibited superior hydraulic efficiency compared with the CentriMag. The PediVAS was evaluated in 14 pediatric animals for up to 30 days, demonstrating acceptable hydraulic function and hemocompatibility. The current results substantiate the performance and biocompatibility of the PediVAS cardiac assist system and are likely to support initiation of a US clinical trial in the future.

Development and validation of a computational fluid dynamics methodology for simulation of pulsatile left ventricular assist devices

An unsteady computational fluid dynamic methodology was developed so that design analyses could be undertaken for devices such as the 50cc Penn State positive-displacement left ventricular assist device (LVAD). The piston motion observed in vitro was modeled, yielding the physiologic flow waveform observed during pulsatile experiments. Valve closure was modeled numerically by locally increasing fluid viscosity during the closed phase. Computational geometry contained Bjork-Shiley Monostrut mechanical heart valves in mitral and aortic positions. Cases for computational analysis included LVAD operation under steady-flow and pulsatile-flow conditions. Computations were validated by comparing simulation results with previously obtained in vitro particle image velocimetry (PIV) measurements. The steady portion of the analysis studied effects of mitral valve orientation, comparing the computational results with in vitro data obtained from mock circulatory loop experiments. The velocity field showed good qualitative agreement with the in vitro PIV data. The pulsatile flow simulations modeled the unsteady flow phenomena associated with a positive-displacement LVAD operating through several beat cycles. Flow velocity gradients allowed computation of the scalar wall strain rate, an important factor for determining hemodynamics of the device. Velocity magnitude contours compared well with PIV data throughout the cycle. Computational wall shear rates over the pulsatile cycle were found to be in the same range as wall shear rates observed in vitro.

Models of Flow-Induced Loading on Blood Cells in Laminar and Turbulent Flow, with Application to Cardiovascular Device Flow

Viscous shear stress and Reynolds stress are often used to predict hemolysis and thrombosis due to flow-induced stress on blood elements in cardiovascular devices. These macroscopic stresses are distinct from the true stress on an individual cell, which is determined by the local microscale flow field. In this paper the flow-induced stress on blood cells is calculated for laminar and turbulent flow, using simplified models for cells and for turbulent eddies. The model is applied to estimate shear stress on red blood cells in flow through a prosthetic heart valve, using the energy spectral density measured by Liu et al. [J. Biomech. Eng. 122:118-124, 2000]. Results show that in laminar flow, the maximum stress on a cell is approximately equal to the macroscopic viscous shear stress. In turbulent flow through a prosthetic heart valve, the estimated root mean square of flow-induced stress on a cell is at least an order of magnitude less than the Reynolds stress. The results support the hypothesis that smaller turbulent eddies cause higher stress on cells. However, the stress due to an eddy depends on the velocity scale of the eddy as well as its length scale. For the heart valve flow investigated, turbulence contributes to flow-induced stress on cells almost equally across a broad range of the frequency spectrum. The model suggests that Reynolds stress alone is not an adequate predictor of cell damage in turbulent flow, and highlights the importance of the energy spectral density.

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.

Multiscale analysis of surface thrombosis in vivo in a left ventricular assist system

Thrombosis limits the success of ventricular assist devices as the demand for alternatives to heart transplants is increasing. This study mapped the occurrence of thrombosis in a left ventricular assist system (LVAS) to better understand the biologic response to these devices. Nine calves divided into two groups were implanted with LVAS for 28 to 30 days. One group was anticoagulated, whereas the second group received no long-term anticoagulation. The blood-contacting poly(urethane urea) surfaces of blood sacs in the LVAS were examined for macroscopic thrombi upon retrieval. The sac was partitioned into eight sections and imaged for thrombi by scanning electron microscopy. No difference in thrombosis was observed macroscopically between the groups. Anticoagulation appeared to result in reduction of platelet-like structures, but the presence of fibrin-like structures remained similar between groups. Regional differences correlating with high and low shear stress regions were observed. At the macroscale, fewer thrombi were recorded in the high shear stress ports. At the microscale, features resembling fibrin were observed primarily in the ports and platelet-like features were common in lower shear stress regions. These variations in thrombosis with anticoagulation and location are likely due to varied fluid dynamics within the LVAS blood sac.

Feasibility of a Miniature Centrifugal Rotary Blood Pump for Low-Flow Circulation in Children and Infants

In this study, a seal-less, tiny centrifugal rotary blood pump was designed for low-flow circulatory support in children and infants. The design was targeted to yield a compact and priming volume of 5 ml with a flow rate of 0.5-4 l/min against a head pressure of 40-100 mm Hg. To meet the design requirements, the first prototype had an impeller diameter of 30 mm with six straight vanes. The impeller was supported with a needle-type hydrodynamic bearing and was driven with a six-pole radial magnetic driver. The external pump dimensions included a pump head height of 20 mm, diameter of 49 mm, and priming volume of 5 ml. The weight was 150 g, including the motor driver. In the mock circulatory loop, using fresh porcine blood, the pump yielded a flow of 0.5-4.0 l/min against a head pressure of 40-100 mm Hg at a rotational speed of 1800-4000 rpm using 1/4" inflow and outflow conduits. The maximum flow and head pressure of 5.25 l/min and 244 mm Hg, respectively, were obtained at a rotational speed of 4400 rpm. The maximum electrical-to-hydraulic efficiency occurred at a flow rate of 1.5-3.5 l/min and at a rotational speed of 2000-4400 rpm. The normalized index of hemolysis, which was evaluated using fresh porcine blood, was 0.0076 g/100 l with the impeller in the down-mode and a bearing clearance of 0.1 mm. Further refinement in the bearing and magnetic coupler are required to improve the hemolytic performance of the pump. The durability of the needle-type hydrodynamic bearing and antithrombotic performance of the pump will be performed before clinical applications. The tiny centrifugal blood pump meets the flow requirements necessary to support the circulation of pediatric patients.

Pulsatile Pediatric Ventricular Assist Devices

A pulsatile pediatric ventricular assist device (PVAD) is being developed at The Pennsylvania State University to provide mechanical circulatory support in infants and children. The PVAD is based on the design of the adult-sized Pierce-Donachy VAD (Thoratec VAD). The infant-sized PVAD has a dynamic stroke volume of approximately 13 ml. A larger 25-ml size for children is also planned. The stroke volumes, beat rates, and pulsatility are comparable to normal physiologic values. The expected maximum duration of use is 6 months. The PVAD is intended to be placed paracorporeally, although the pump may be implanted for bridge-to-transplant applications. The pneumatically actuated PVAD uses a full-to-empty control mode, which allows maximum flexibility in its application for left, right, or biventricular assistance. The device may be used with atrial or ventricular inlet cannulation, with blood return to the aorta or pulmonary artery. In vitro testing is underway to measure hydrodynamic performance, hemolysis, and flow velocities using particle image velocimetry. In vivo implantation studies will be performed in juvenile goats or sheep after the completion of baseline studies to assess hematology and surgical fit.

Development of Novel Submicron Textured Polyether(Urethane Urea) for Decreasing Platelet Adhesion

Ventricular assist devices have proven to be a useful clinical option for providing circulatory support as a bridge to transplantation and a mode of destination therapy. Thromboembolism is prevented by designing devices that use blood interfaces that either encourage biological material deposition and strong adhesion, or discourage deposition via surface chemistry, surface finish, and fluid flow fields. Minimum continuous or periodic wall shear forces and maximum time at reduced shear are important, and sometimes difficult-to-satisfy, design constraints. We present an approach to reducing platelet adhesion via surface topography, reducing surface area for platelet-material interaction. Large areas of polyether(urethane urea) were textured with two different sizes of ordered pillar arrays via two-stage replication molding without affecting surface chemistry. Pillars had subplatelet dimensions designed to reduce the surface area a platelet may contact. Platelet adhesion was assessed in a physiologically relevant shear stress range from 0-10 dyn/cm2 using a rotating disk and compared to smooth control. Adhesion was highest from 0-5 dyn/cm2. Surface texturing reduced platelet adhesion without increasing platelet activation in bulk suspension. This study demonstrates that material surface texture is an additional variable that may be used to reduce platelet adhesion under low shear stresses potentially reducing thromboembolism.

Viscoelasticity of Pediatric Blood and its Implications for the Testing of a Pulsatile Pediatric Blood Pump

Red blood cell hematocrit, aggregation and deformability, and plasma protein concentration influence the viscosity and elasticity of whole blood. These parameters affect the flow properties, especially at low shear rates (< 50 s(-1)). In particular, we have previously shown that the viscoelasticity of fluid affects the inlet filling characteristics and regions of flow separation in small pulsatile blood pumps. Although the viscosity of pediatric blood has been thoroughly studied, its elasticity has not been previously measured. Here we present the viscosity and elasticity of pediatric blood against shear rate for hematocrits from 19-56, measured using an oscillatory rheometer. There is little effect of patient age on blood viscoelasticity. A statistical analysis showed that when compared at constant hematocrit, blood from adult and pediatric patients had similar viscoelastic properties. We present blood analog solutions, as a function of hematocrit, constructed on the basis of the pediatric measurements. Flow field results for viscoelastic analogs of 20, 40 and 60% hematocrit and a Newtonian analog will be compared in the initial, in vitro testing of the Penn State pediatric blood pump, to determine the importance of incorporating a viscoelastic analog into the desigh interaction.

Off-design Considerations of the 50cc Penn State Ventricular Assist Device

Particle image velocimetry is used to study the flow field associated with off-design conditions of the 50cc Penn State Ventricular Assist Device. Phase average data, for 300 consecutive beats, was acquired for 60, 75, and 90 beats per minute (bpm) at 35/65 systolic duration and at 75 bpm for 50/50 systolic duration. Diastolic flow is dominated by the inflow jet and shows more penetration and better wall washing at higher beat rates and shorter diastolic time. A flow separation near the minor orifice of the aortic valve, however, is larger under these conditions. Turbulent energy, which peaks during diastole, is found to be strongest at the edges of the jet and near the valves and scales well with the maximum inlet jet velocity. A region of potential separation on the right lateral wall is identified and will be relieved by geometric modifications of the evolving design.

Fluid Dynamic Analysis of the 50 cc Penn State Artificial Heart Under Physiological Operating Conditions Using Particle Image Velocimetry

In order to bridge the gap of existing artificial heart technology to the diverse needs of the patient population, we have been investigating the viability of a scaled-down design of the current 70 cc Penn State artificial heart. The issues of clot formation and hemolysis may become magnified within a 50 cc chamber compared to the existing 70 cc one. Particle image velocimetry (PIV) was employed to map the entire 50 cc Penn State artificial heart chamber. Flow fields constructed from PIV data indicate a rotational flow pattern that provides washout during diastole. In addition, shear rate maps were constructed for the inner walls of the heart chamber. The lateral walls of the mitral and aortic ports experience high shear rates while the upper and bottom walls undergo low shear rates, with sufficiently long exposure times to potentially induce platelet activation or thrombus formation. In this study, we have demonstrated that PIV may adequately map the flow fields accurately in a reasonable amount of time. Therefore, the potential exists of employing PIV as a design tool.

Wall Shear-Rate Estimation Within the 50cc Penn State Artificial Heart Using Particle Image Velocimetry

Particle image velocimetry (PIV) has been gaining acceptance as a routine tool to evaluate the flow fields associated with fluid mechanical devices. We have developed algorithms to investigate the wall shear-rates within the 50cc Penn State artificial heart using low magnification, conventional particle image velocimetry (PIV). Wall shear has been implicated in clot formation, a major post-implant problem with artificial hearts. To address the issues of wall scattering and incomplete measurement volumes, associated with near wall measurements, we have introduced a zero masking and a fluid centroid shifting technique. Simulations using different velocity fields were conducted with the techniques to assess their viability. Subsequently, the techniques were applied to the experimental data collected. The results indicate that the size of the interrogation region should be chosen to be as small as possible to maximize resolution while large enough to ensure an adequate number of particles per region. In the current study, a 16×16 interrogation window performed well with good spatial resolution and particle density for the estimation of wall shear rate. The techniques developed with PIV allow wall shear-rate estimates to be obtained from a large number of sites at one time. Because a planar image of a flow field can be determined relatively rapidly, PIV may prove useful in any preliminary design procedure.