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

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.     

The 12 cc Penn State pulsatile pediatric ventricular assist device: fluid dynamics associated with valve selection

The mortality rate for infants awaiting a heart transplant is 40% because of the extremely limited number of donor organs. Ventricular assist devices (VADs), a common bridge-to-transplant solution in adults, are becoming a viable option for pediatric patients. A major obstacle faced by VAD designers is 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(-1). Particle image velocimetry was used to compare the flow fields in the chamber of the 12 cc Penn State pediatric pulsatile VAD using two mechanical heart valves: Bjork-Shiley monostrut (BSM) tilting disk valves and CarboMedics (CM) bileaflet valves. In conjunction with the flow evaluation, wall shear data were calculated and analyzed to help quantify wall washing. The major orifice inlet jet of the device containing BSM valves was more intense, which led to better recirculation and wall washing than the three jets produced by the CM valves. Regurgitation through the CM valve served as a significant hindrance to the development of the rotational flow.     

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 novel mouse-driven ex vivo flow chamber for the study of leukocyte and platelet function

Various in vitro and in vivo techniques exist for study of the microcirculation. Whereas in vivo systems impress with their physiological fidelity, in vitro systems excel in the amount of reduction that can be achieved. Here we introduce the autoperfused ex vivo flow chamber designed to study murine leukocytes and platelets under well-defined hemodynamic conditions. In our model, the murine heart continuously drives the blood flow through the chamber, providing a wide range of physiological shear rates. We used a balance of force approach to quantify the prevailing forces at the chamber walls. Numerical simulations show the flow characteristics in the chamber based on a shear-thinning fluid model. We demonstrate specific rolling of wild-type leukocytes on immobilized P-selectin, abolished by a blocking MAb. When uncoated, the surfaces having a constant shear rate supported individual platelet rolling, whereas on areas showing a rapid drop in shear platelets interacted in previously unreported grapelike conglomerates, suggesting an influence of shear rate on the type of platelet interaction. In summary, the ex vivo chamber amounts to an external vessel connecting the arterial and venous systems of a live mouse. This method combines the strengths of existing in vivo and in vitro systems in the study of leukocyte and platelet function. autoperfused flow chamber; intravital microscopy; inflammation; thrombus formation     

PIV measurements of flow in a centrifugal blood pump: steady flow

Magnetically suspended left ventricular assist devices have only one moving part, the impeller. The impeller has absolutely no contact with any of the fixed parts, thus greatly reducing the regions of stagnant or high shear stress that surround a mechanical or fluid bearing. Measurements of the mean flow patterns as well as viscous and turbulent stresses were made in a shaft-driven prototype of a magnetically suspended centrifugal blood pump at several constant flow rates (3-9 L/min) using particle image velocimetry (PIV). The chosen range of flow rates is representative of the range over which the pump may operate while implanted. Measurements on a three-dimensional measurement grid within several regions of the pump, including the inlet, blade passage, exit volute, and diffuser are reported. The measurements are used to identify regions of potential blood damage due to high shear stress and/or stagnation of the blood, both of which have been associated with blood damage within artificial heart valves and diaphragm-type pumps. Levels of turbulence intensity and Reynolds stresses that are comparable to those in artificial heart valves are reported. At the design flow rate (6 L/min), the flow is generally well behaved (no recirculation or stagnant flow) and stress levels are below levels that would be expected to contribute to hemolysis or thrombosis. The flow at both high (9 L/min) and low (3 L/min) flow rates introduces anomalies into the flow, such as recirculation, stagnation, and high stress regions. Levels of viscous and Reynolds shear stresses everywhere within the pump are below reported threshold values for damage to red cells over the entire range of flow rates investigated; however, at both high and low flow rate conditions, the flow field may promote activation of the clotting cascade due to regions of elevated shear stress adjacent to separated or stagnant flow.     

Performance analysis of a cardiac assist device in counterpulsation.

Performance of a cardiac assist device pumping chamber in counterpulsation was evaluated using numerical simulations of the unsteady, three-dimensional flow inside the chamber and an analytical model of the force required to eject and fill the chamber. The wall shear stress within the device was similarly computed and modeled. The analytical model was scaled to match the numerical results and then used to predict performance at physiological operating conditions. According to these models for a stroke volume of 70 ml, between 0.4 and 1.0 W is required for counterpulsation at a frequency of 1.33 Hz against a restorative spring, depending on the spring constant chosen. The power and the maximum force calculated are within the ranges a trained skeletal muscle is capable of providing. Shear stress predictions show that platelet activation in the absence of surface effects and hemolysis due to high shear are unlikely to occur with this design. Furthermore, vortices that develop in the chamber during filling are predicted to increase blood mixing and provide favorable washing of the chamber walls. A computational-analytical approach such as this may have potential to aid rapid performance evaluation of new devices and streamline the design optimization process.     

Study of local hydrodynamic environment in cell-substrate adhesion using side-view μPIV technology

Tumor cell adhesion to the endothelium under shear flow conditions is a critical step that results in circulation-mediated tumor metastasis. This study presents experimental and computational techniques for studying the local hydrodynamic environment around adherent cells and how local shear conditions affect cell-cell interactions on the endothelium in tumor cell adhesion. To study the local hydrodynamic profile around heterotypic adherent cells, a side-view flow chamber assay coupled with micro particle imaging velocimetry (μPIV) technique was developed, where interactions between leukocytes and tumor cells in the near-endothelial wall region and the local shear flow environment were characterized. Computational fluid dynamics (CFD) simulations were also used to obtain quantitative flow properties around those adherent cells. Results showed that cell dimension and relative cell-cell positions had strong influence on local shear rates. The velocity profile above leukocytes and tumor cells displayed very different patterns. Larger cell deformations led to less disturbance to the flow. Local shear rates above smaller cells were observed to be more affected by relative positions between two cells.     

Design of a side-view particle imaging velocimetry flow system for cell-substrate adhesion studies

Experimental models that mimic the flow conditions in microcapillaries have suggested that the local shear stresses and shear rates can mediate tumor cell and leukocyte arrest on the endothelium and subsequent sustained adhesion. However, further investigation has been limited by the lack of experimental models that allow quantitative measurement of the hydrodynamic environment over adherent cells. The purpose of this study was to develop a system capable of acquiring quantitative flow profiles over adherent cells. By combining the techniques of side-view imaging and particle image velocimetry (PIV), an in vitro model was constructed that is capable of obtaining quantitative flow data over cells adhering to the endothelium. The velocity over an adherent leukocyte was measured and the shear rate was calculated under low and high upstream wall shear. The microcapillary channel was modeled using computational fluid dynamics (CFD) and the calculated velocity profiles over cells under the low and high shear rates were compared to experimental results. The drag force applied to each cell by the fluid was then computed. This system provides a means for future study of the forces underlying adhesion by permitting characterization of the local hydrodynamic conditions over adherent cells.     

Experimental measurement of dynamic fluid shear stress on the aortic surface of the aortic valve leaflet

Aortic valve (AV) calcification is a highly prevalent disease with serious impact on mortality and morbidity. Although exact causes and mechanisms of AV calcification are unclear, previous studies suggest that mechanical forces play a role. Since calcium deposits occur almost exclusively on the aortic surfaces of AV leaflets, it has been hypothesized that adverse patterns of fluid shear stress on the aortic surface of AV leaflets promote calcification. The current study characterizes AV leaflet aortic surface fluid shear stresses using Laser Doppler velocimetry and an in vitro pulsatile flow loop. The valve model used was a native porcine valve mounted on a suturing ring and preserved using 0.15% glutaraldehyde solution. This valve model was inserted in a mounting chamber with sinus geometries, which is made of clear acrylic to provide optical access for measurements. To understand the effects of hemodynamics on fluid shear stress, shear stress was measured across a range of conditions: varying stroke volumes at the same heart rate and varying heart rates at the same stroke volume. Systolic shear stress magnitude was found to be much higher than diastolic shear stress magnitude due to the stronger flow in the sinuses during systole, reaching up to 20 dyn/cm² at mid-systole. Upon increasing stroke volume, fluid shear stresses increased due to stronger sinus fluid motion. Upon increasing heart rate, fluid shear stresses decreased due to reduced systolic duration that restricted the formation of strong sinus flow. Significant changes in the shear stress waveform were observed at 90 beats/min, most likely due to altered leaflet dynamics at this higher heart rate. Overall, this study represents the most well-resolved shear stress measurements to date across a range of conditions on the aortic side of the AV. The data presented can be used for further investigation to understand AV biological response to shear stresses.     

Fluid Dynamics of Ventricular Filling in the Embryonic Heart

The vertebrate embryonic heart first forms as a valveless tube that pumps blood using waves of contraction. As the heart develops, the atrium and ventricle bulge out from the heart tube, and valves begin to form through the expansion of the endocardial cushions. As a result of changes in geometry, conduction velocities, and material properties of the heart wall, the fluid dynamics and resulting spatial patterns of shear stress and transmural pressure change dramatically. Recent work suggests that these transitions are significant because fluid forces acting on the cardiac walls, as well as the activity of myocardial cells that drive the flow, are necessary for correct chamber and valve morphogenesis. In this article, computational fluid dynamics was used to explore how spatial distributions of the normal forces acting on the heart wall change as the endocardial cushions grow and as the cardiac wall increases in stiffness. The immersed boundary method was used to simulate the fluid-moving boundary problem of the cardiac wall driving the motion of the blood in a simplified model of a two-dimensional heart. The normal forces acting on the heart walls increased during the period of one atrial contraction because inertial forces are negligible and the ventricular walls must be stretched during filling. Furthermore, the force required to fill the ventricle increased as the stiffness of the ventricular wall was increased. Increased endocardial cushion height also drastically increased the force necessary to contract the ventricle. Finally, flow in the moving boundary model was compared to flow through immobile rigid chambers, and the forces acting normal to the walls were substantially different.     

PIV-measured versus CFD-predicted flow dynamics in anatomically realistic cerebral aneurysm models

Computational fluid dynamics (CFD) modeling of nominally patient-specific cerebral aneurysms is increasingly being used as a research tool to further understand the development, prognosis, and treatment of brain aneurysms. We have previously developed virtual angiography to indirectly validate CFD-predicted gross flow dynamics against the routinely acquired digital subtraction angiograms. Toward a more direct validation, here we compare detailed, CFD-predicted velocity fields against those measured using particle imaging velocimetry (PIV). Two anatomically realistic flow-through phantoms, one a giant internal carotid artery (ICA) aneurysm and the other a basilar artery (BA) tip aneurysm, were constructed of a clear silicone elastomer. The phantoms were placed within a computer-controlled flow loop, programed with representative flow rate waveforms. PIV images were collected on several anterior-posterior (AP) and lateral (LAT) planes. CFD simulations were then carried out using a well-validated, in-house solver, based on micro-CT reconstructions of the geometries of the flow-through phantoms and inlet/outlet boundary conditions derived from flow rates measured during the PIV experiments. PIV and CFD results from the central AP plane of the ICA aneurysm showed a large stable vortex throughout the cardiac cycle. Complex vortex dynamics, captured by PIV and CFD, persisted throughout the cardiac cycle on the central LAT plane. Velocity vector fields showed good overall agreement. For the BA, aneurysm agreement was more compelling, with both PIV and CFD similarly resolving the dynamics of counter-rotating vortices on both AP and LAT planes. Despite the imposition of periodic flow boundary conditions for the CFD simulations, cycle-to-cycle fluctuations were evident in the BA aneurysm simulations, which agreed well, in terms of both amplitudes and spatial distributions, with cycle-to-cycle fluctuations measured by PIV in the same geometry. The overall good agreement between PIV and CFD suggests that CFD can reliably predict the details of the intra-aneurysmal flow dynamics observed in anatomically realistic in vitro models. Nevertheless, given the various modeling assumptions, this does not prove that they are mimicking the actual in vivo hemodynamics, and so validations against in vivo data are encouraged whenever possible.     

Fluid flow and plaque formation in an aortic bifurcation

Considering steady laminar flow in a two-dimensional symmetric branching channel with local occlusions, a finite element model has been developed to study velocity fields including reverse flow regions, pressure profiles and wall shear stress distributions for different Reynolds numbers, bifurcation angles and lumen reductions. The flow analysis has been extended to include a new submodel for the pseudo-transient formation of plaque at sites and deposition rates defined by the physical characteristics of the flow. Specifically, simulating the onset of atherosclerotic lesions, sinusoidal plaque layers have been placed in areas of critically low wall shear stresses, and simulating the growth of particle depositions, plaque layers have been added in a stepwise fashion in regions of critically high and low shear. Thus two somewhat conflicting hypothetical correlations between critical wall shear stress levels and atheroma have been tested and a solution has been postulated. The validated computer simulation model is a predictive tool for analyzing the effects of local changes in wall curvature due to surgical reconstruction and/or atherosclerotic lesions, and for investigating the design of aortic bifurcations which mitigate plaque formation.     

Shear flow induced changes in apolipoprotein C-II conformation and amyloid fibril formation

The misfolding and self-assembly of proteins into amyloid fibrils that occur in several debilitating diseases are affected by a variety of environmental factors, including mechanical factors associated with shear flow. We examined the effects of shear flow on amyloid fibril formation by human apolipoprotein C-II (apoC-II). Shear fields (150, 300, and 500 s(-1)) accelerated the rate of apoC-II fibril formation (1 mg/mL) approximately 5-10-fold. Fibrils produced at shear rates of 150 and 300 s(-1) were similar to the twisted ribbon fibrils formed in the absence of shear, while at 500 s(-1), tangled ropelike structures were observed. The mechanism of the shear-induced acceleration of amyloid fibril formation was investigated at low apoC-II concentrations (50 μg/mL) where fibril formation does not occur. Circular dichroism and tryptophan fluorescence indicated that shear induced an irreversible change in apoC-II secondary structure. Fluorescence resonance energy transfer experiments using the single tryptophan residue in apoC-II as the donor and covalently attached acceptors showed that shear flow increased the distance between the donor and acceptor molecules. Shear-induced higher-order oligomeric species were identified by sedimentation velocity experiments using fluorescence detection, while fibril seeding experiments showed that species formed during shear flow are on the fibril formation pathway. These studies suggest that physiological shear flow conditions and conditions experienced during protein manufacturing can exert significant effects on protein conformation, leading to protein misfolding, aggregation, and amyloid fibril formation.     

Theory and design of disposable clinical blood viscometer

A disposable clinical whole blood viscometer which can produce viscosity measures over a wide range of shear rates in a single rapid determination has been developed and is currently under test. The design is based upon the time varying flow of blood through a capillary. The flow is driven by the pressure in a fixed volume air chamber and transmitted to the sample through a compliant membrane. The time varying pressure in the air chamber is measured by a suitable transducer. The instantaneous shear stress of the blood in the capillary is proportional to the air pressure, while the instantaneous shear rate is proportional to the pressure-time derivative. Proper design ensures that the system operates as a first order dynamic system with flow resistance entirely determined by the nonlinear sample viscosity. By constructing the air chamber in two parts coupled by a quick disconnect fitting the design can allow for the blood-containing part of the instrument to be discarded, eliminating handling and cleaning of blood contacted components. The entire determination is completed in less than a minute, so that anticoagulants are not necessary. Tests on a prototype show that the instrument gives results in excellent agreement with those obtained on a cone-plate rheogoniometer.     

A dual role of extracellular DNA during biofilm formation of Neisseria meningitidis

Major pathogenic clonal complexes (cc) of Neisseria meningitidis differ substantially in their point prevalence among healthy carriers. We show that frequently carried pathogenic cc (e.g. sequence type ST‐41/44 cc and ST‐32 cc) depend on extracellular DNA (eDNA) to initiate in vitro biofilm formation, whereas biofilm formation of cc with low point prevalence (ST‐8 cc and ST‐11 cc) was eDNA‐independent. For initial biofilm formation, a ST‐32 cc type strain, but not a ST‐11 type strain, utilized eDNA. The release of eDNA was mediated by lytic transglycosylase and cytoplasmic N‐acetylmuramyl‐l ‐alanine amidase genes. In late biofilms, outer membrane phospholipase A‐dependent autolysis, which was observed in most cc, but not in ST‐8 and ST‐11 strains, was required for shear force resistance of microcolonies. Taken together, N. meningitidis evolved two different biofilm formation strategies, an eDNA‐dependent one yielding shear force resistant microcolonies, and an eDNA‐independent one. Based on the experimental findings and previous epidemiological observations, we hypothesize that most meningococcal cc display a settler phenotype, which is eDNA‐dependent and results in a stable interaction with the host. On the contrary, spreaders (ST‐11 and ST‐8 cc) are unable to use eDNA for biofilm formation and might compensate for poor colonization properties by high transmission rates.     

The effect of shape on the margination dynamics of non-neutrally buoyant particles in two-dimensional shear flows

The margination dynamics of microparticles with different shapes has been analyzed within a laminar flow mimicking the hydrodynamic conditions in the microcirculation. Silica spherical particles, quasi-hemispherical and discoidal silicon particles have been perfused in a parallel plate flow chamber. The effect of the shape and density on their margination propensity has been investigated at different physiologically relevant shear rates S. Simple scaling laws have been derived showing that the number n of marginating particles scales as S(-0.63) for the spheres; S(-0.85) for discoidal and S(-1) for quasi-hemispherical particles, regardless of their density and size. Within the range considered for the shear rate, discoidal particles marginate in a larger number compared to quasi-hemispherical and spherical particles. These results may be of interest in drug delivery and bio-imaging applications, where particles are expected to drift towards and interact with the walls of the blood vessels.     

Modeling and design of optimal flow perfusion bioreactors for tissue engineering applications

Perfusion bioreactors have been used in different tissue engineering applications because of their consistent distribution of nutrients and flow-induced shear stress within the tissue-engineering scaffold. A widely used configuration uses a scaffold with a circular cross-section enclosed within a cylindrical chamber and inlet and outlet pipes which are connected to the chamber on either side through which media is continuously circulated. However, fluid-flow experiments and simulations have shown that the majority of the flow perfuses through the center. This pattern creates stagnant zones in the peripheral regions as well as in those of high flow rate near the inlet and outlet. This non-uniformity of flow and shear stress, owing to a circular design, results in limited cell proliferation and differentiation in these areas. The focus of this communication is to design an optimized perfusion system using computational fluid dynamics as a mathematical tool to overcome the time-consuming trial and error experimental method. We compared the flow within a circular and a rectangular bioreactor system. Flow simulations within the rectangular bioreactor are shown to overcome the limitations in the circular design. This communication challenges the circular cross-section bioreactor configuration paradigm and provides proof of the advantages of the new design over the existing one.     

Velocity and turbulence measurements past mitral valve prostheses in a model left ventricle

Thrombogenesis and hemolysis have both been linked to the flow dynamics past heart valve prostheses. To learn more about the particular flow dynamics past mitral valve prostheses in the left ventricle under controlled experimental conditions, an in vitro study was performed. The experimental methods included velocity and turbulent shear stress measurements past caged-ball, tilting disc, bileaflet, and polyurethane trileaflet mitral valves in an acrylic rigid model of the left ventricle using laser Doppler anemometry. The results indicate that all four prosthetic heart valves studied create at least mildly disturbed flow fields. The effect of the left ventricular geometry on the flow development is to produce a stabilizing vortex which engulfs the entire left ventricular cavity, depending on the orientation of the valve. The measured turbulent shear stress magnitudes for all four valves did not exceed the reported value for hemolytic damage. However, the measured turbulent shear stresses were near or exceeded the critical shear stress reported in the literature for platelet lysis, a known precursor to thrombus formation.