The blood flow dynamics in a stenosed, subject-specific carotid bifurcation is numerically simulated using Direct Numerical Simulation (DNS) and Reynolds-averaged Navier--Stokes (RANS) equations closed with turbulence models. The former is meant to provide a term of comparison for the RANS calculations, that include classic two-equations models ($k-\epsilon$ and $k-\omega$) as well as a transitional three-equations eddy-viscosity model ($k_T-k_L-\omega$). Pulsatile inlet conditions based on in-vivo ultrasound measurements of blood velocity are used. The blood is modelled as a Newtonian fluid, and the vessel walls are rigid. The main purpose of this work is to highlight the problems related to the use of classic RANS models in the numerical simulation of such flows. The time-averaged DNS results, interpreted in view of their finite-time averaging error, are used to demonstrate the superiority of the transitional RANS model, which is found to provide results closer to DNS than those of conventional models. The transitional model is shown to possess better predictive capabilities in terms of turbulence intensity, temporal evolution of the pressure along the cardiac cycle, and the oscillatory shear index (OSI). Indeed, DNS brings to light the locally transitional or weakly turbulent state of the blood flow, which presents velocity and pressure fluctuations only in the post-stenotic region of the internal carotid artery during systole, while the flow is laminar during diastole.

The purpose of this study is to quantify the impact of enzyme activity on the vitreous humor structure over time to understand the mechanical characteristics of the vitreous humor gel. Changes in the mechanical behavior of the vitreous occur due to many reasons including aging, which may lead to many vitreoretinal diseases. The degeneration process of the vitreous has been studied; however, in situ experimental procedures to validate the existing hypotheses are limited. We examined thirty-eight porcine eyes using in situ rheological creep tests to measure the mechanical properties of the vitreous humor of the eyes prior to, 1 h and 24 h after the intravitreal injection. Eyes in one group were injected with collagenase type II solution and eyes in the control group were injected with phosphate buffered saline solution (PBS) with calcium and magnesium chloride. Prior to the injection, viscosity and creep compliance intercept values between both groups were not statistically different. At 1 h and 24 h after the injection, vitreous properties in the eyes from the first group showed a statistically significant increase in the J intercept values (representing the inverse of elasticity) compared to the control group. In addition, 1 h and 24 h after the injection, vitreous viscosity was lower in the eyes from the first group than in the eyes from the control group. These findings are a foundation for future studies on the effectiveness of intravitreal drugs that modify the mechanical properties of the vitreous humor.

This work aims to provide a contribution to determine a proper model for the study of fluid film lubrication for the reduction of knee prostheses failure due to polyethylene wear. The Phan-Thien and Tanner (PTT) rheological law and the elastic deformation of the articular surfaces were considered in this modeling. The governing equations were solved numerically for different geometries and different Weissenberg numbers. The lubrication approximation applied to the PTT rheological law leads to an expression for the apparent viscosity similar to the Cross model. The results attest the importance of considering the non-Newtonian behavior of the synovial fluid, the elastic deformation, and the geometrical features of the prostheses to obtain quantitative information.

Cytoplasmic viscosity-dependent margination of red blood cells (RBC) for flow inside microchannels was studied using numerical simulations, and the results were verified with microfluidic experiments. Wide range of suspension volume fractions or hematocrits was considered in this study. Lattice Boltzmann method for fluid-phase coupled with spectrin-link method for RBC membrane deformation was used for accurate analysis of cell margination. RBC margination behavior shows strong dependence on the internal viscosity of the RBCs. At equilibrium, RBCs with higher internal viscosity marginate closer to the channel wall and the RBCs with normal internal viscosity migrate to the central core of the channel. Same margination pattern has been verified through experiments conducted with straight channel microfluidic devices. Segregation between RBCs of different internal viscosity is enhanced as the shear rate and the hematocrit increases. Stronger separation between normal RBCs and RBCs with high internal viscosity is obtained as the width of a high aspect ratio channel is reduced. Overall, the margination behavior of RBCs with different internal viscosities resembles with the margination behavior of RBCs with different levels of deformability. Observations from this work will be useful in designing microfluidic devices for separating the subpopulations of RBCs with different levels of deformability that appear in many hematologic diseases such as sickle cell disease (SCD), malaria, or cancer.

Soft connective tissues sustain large strains of viscoelastic nature. The rate-independent component is frequently modeled by means of anisotropic hyperelastic models. The rate-dependent component is usually modeled through linear rheological models or quasi-linear viscoelastic (QLV) models. These viscoelastic models are unable, in general, to capture the strain-level dependency of the viscoelastic properties present in many viscoelastic tissues. In linear viscoelastic models, strain-level dependency is frequently accounted for by including the dependence of multipliers of Prony series on strains through additional evolution laws, but the determination of the material parameters is a difficult task and the obtained accuracy is usually not sufficient. In this work, we introduce a model for fully nonlinear viscoelasticity in which the instantaneous and quasi-static behaviors are exactly captured and the relaxation curves are predicted to a high accuracy. The model is based on a fully nonlinear standard rheological model and does not necessitate optimization algorithms to obtain material parameters. Furthermore, in contrast to most models used in modeling the viscoelastic behavior of soft tissues, it is valid for the large deviations from thermodynamic equilibrium typically observed in soft tissues.

The Pacinian corpuscle (PC) is a cutaneous mechanoreceptor that senses low-amplitude, high-frequency vibrations. The PC contains a nerve fiber surrounded by alternating layers of solid lamellae and interlamellar fluid, and this structure is hypothesized to contribute to the PC's role as a band-pass filter for vibrations. In this study, we sought to evaluate the relationship between the PC's material and geometric parameters and its response to vibration. We used a spherical finite element mechanical model based on shell theory and lubrication theory to model the PC's outer core. Specifically, we analyzed the effect of the following structural properties on the PC's frequency sensitivity: lamellar modulus (E), lamellar thickness (h), fluid viscosity ( μ ), PC outer radius (R o ), and number of lamellae (N). The frequency of peak strain amplification (henceforth “peak frequency”) and frequency range over which strain amplification occurred (henceforth “bandwidth”) increased with lamellar modulus or lamellar thickness and decreased with an increase in fluid viscosity or radius. All five structural parameters were combined into expressions for the relationship between the parameters and peak frequency, ω peak = 1.605 × 10 − 6 N 3.475 ( E h / μ R o ) , or bandwidth, B = 1.747 × 10 − 6 N 3.951 ( E h / μ R o ) . Although further work is needed to understand how mechanical variability contributes to functional variability in PCs and how factors such as PC eccentricity also affect PC behavior, this study provides two simple expressions that can be used to predict the impact of structural or material changes with aging or disease on the frequency response of the PC.

Given the current emphasis on accurate computational fluid dynamics (CFD) modeling of cardiovascular flows, which incorporates realistic blood vessel geometries and cardiac waveforms, it is necessary to revisit the conventional wisdom regarding the influences of non-Newtonian effects. In this study, patient-specific reconstructed 3D geometries, whole blood viscosity data, and venous pulses postdialysis access surgery are used as the basis for the hemodynamic simulations of renal failure patients with native fistula access. Rheological analysis of the viscometry data initially suggested that the correct choice of constitutive relations to capture the non-Newtonian behavior of blood is important because the end-stage renal disease (ESRD) patient cohort under observation experience drastic variations in hematocrit (Hct) levels and whole blood viscosity throughout the hemodialysis treatment. For this purpose, various constitutive relations have been tested and implemented in CFD practice, namely Quemada and Casson. Because of the specific interest in neointimal hyperplasia and the onset of stenosis in this study, particular attention is placed on differences in nonhomeostatic wall shear stress (WSS) as that drives the venous adaptation process that leads to venous geometric evolution over time in ESRD patients. Surprisingly, the CFD results exhibit no major differences in the flow field and general flow characteristics of a non-Newtonian simulation and a corresponding identical Newtonian counterpart. It is found that the vein's geometric features and the dialysis-induced flow rate have far greater influence on the WSS distribution within the numerical domain.

The purpose of this work is to predict the effect of impaired red blood cells (RBCs) deformability on blood flow conditions in human carotid artery bifurcation. First, a blood viscosity model is developed that predicts the steady-state blood viscosity as a function of shear rate, plasma viscosity, and mechanical (and geometrical) properties of RBC's. Viscosity model is developed by modifying the well-known Krieger and Dougherty equation for monodisperse suspensions by using the dimensional analysis approach. With the approach, we manage to account for the microscopic properties of RBC's, such as their deformability, in the macroscopic behavior of blood via blood viscosity. In the second part of the paper, the deduced viscosity model is used to numerically predict blood flow conditions in human carotid artery bifurcation. Simulations are performed for different values of RBC's deformability and analyzed by investigating parameters, such as the temporal mean wall shear stress (WSS), oscillatory shear index (OSI), and mean temporal gradient of WSS. The analyses show that the decrease of RBC's deformability decrease the regions of low WSS (i.e., sites known to be prevalent at atherosclerosis-prone regions); increase, in average, the value of WSS along the artery; and decrease the areas of high OSI. These observations provide an insight into the influence of blood's microscopic properties, such as the deformability of RBC's, on hemodynamics in larger arteries and their influence on parameters that are known to play a role in the initiation and progression of atherosclerosis.

This work deals with the viscoelasticity of the arterial wall and its influence on the pulse waves. We describe the viscoelasticity by a nonlinear Kelvin–Voigt model in which the coefficients are fitted using experimental time series of pressure and radius measured on a sheep's arterial network. We obtained a good agreement between the results of the nonlinear Kelvin–Voigt model and the experimental measurements. We found that the viscoelastic relaxation time—defined by the ratio between the viscoelastic coefficient and the Young's modulus—is nearly constant throughout the network. Therefore, as it is well known that smaller arteries are stiffer, the viscoelastic coefficient rises when approaching the peripheral sites to compensate the rise of the Young's modulus, resulting in a higher damping effect. We incorporated the fitted viscoelastic coefficients in a nonlinear 1D fluid model to compute the pulse waves in the network. The damping effect of viscoelasticity on the high-frequency waves is clear especially at the peripheral sites.

To assess how posture affects the transmission of mechanical energy up the spinal column during vibration, 18 porcine functional spinal units (FSUs) were exposed to a sinusoidal force (1500 ± 1200 N) at 5 Hz for 120 min in either a flexed, extended, or neutral posture. Force and FSU height were measured continuously throughout the collection. From these data, specimen height loss, dynamic stiffness, hysteresis, and parameters from a standard linear solid (SLS) model were determined and analyzed for differences between postures. Posture had an influence on all of these parameters. In extension, the FSU had higher dynamic stiffness values than when neutral or flexed (p < 0.0001). In flexion, the FSU had higher hysteresis than both an extended or neutral posture (p < 0.0001). Height loss was greatest in a flexed posture and smallest in an extended posture (p < 0.0001). In extension, the series spring element in the SLS model had a stiffness value higher than both flexed and neutral posture conditions, whereas the stiffness in the parallel spring was the same between extension and neutral (p < 0.01), both higher than in flexion. Viscosity coefficients were highest in extension compared to both flexed and neutral (p < 0.01). Based on these results, it was determined that posture had a significant influence in determining the mechanical properties of the spine when exposed to cyclic compressive loading.

Blood is a complex fluid that, among other things, has been established to behave as a shear thinning, non-Newtonian fluid when exposed to low shear rates (SR). Many hemodynamic investigations use a Newtonian fluid to represent blood when the flow field of study has relatively high SR (>200 s −1 ). Shear thinning fluids have been shown to exhibit differences in transition to turbulence (TT) compared to that of Newtonian fluids. Incorrect prediction of the transition point in a simulation could result in erroneous hemodynamic force predictions. The goal of the present study was to compare velocity profiles near TT of whole blood and Newtonian blood analogs in a straight rigid pipe with a diameter 6.35 mm under steady flow conditions. Rheology was measured for six samples of whole porcine blood and three samples of a Newtonian fluid, and the results show blood acts as a shear thinning non-Newtonian fluid. Measurements also revealed that blood viscosity at SR = 200 s −1 is significantly larger than at SR = 1000 s −1 (13.8%, p < 0.001). Doppler ultrasound (DUS) was used to measure velocity profiles for blood and Newtonian samples at different flow rates to produce Reynolds numbers (Re) ranging from 1000 to 3300 (based on viscosity at SR = 1000 s −1 ). Two mathematically defined methods, based on the velocity profile shape change and turbulent kinetic energy (TKE), were used to detect TT. Results show similar parabolic velocity profiles for both blood and the Newtonian fluid for Re < 2200. However, differences were observed between blood and Newtonian fluid velocity profiles for larger Re. The Newtonian fluid had blunt-like velocity profiles starting at Re = 2403 ± 8 which indicated transition. In contrast, blood did not show this velocity profile change until Re = 2871 ± 104. The Newtonian fluid had large velocity fluctuations (root mean square (RMS) > 20%) with a maximum TKE near the pipe center at Re = 2316 ± 34 which indicated transition. In contrast, blood results showed the maximum TKE at Re = 2806 ± 109. Overall, the critical Re was delayed by ∼20% (p < 0.001) for blood compared to the Newtonian fluid. Thus, a Newtonian assumption for blood at flow conditions near transition could lead to large errors in velocity prediction for steady flow in a straight pipe. However, these results are specific to this pipe diameter and not generalizable since SR is highly dependent on pipe diameter. Further research is necessary to understand this relation in different pipe sizes, more complex geometries, and under pulsatile flow conditions.

This paper reviews and further develops pore-scale computational flow modeling techniques used for creeping flow through orthotropic fiber bundles used in blood oxygenators. Porous model significantly reduces geometrical complexity by taking a homogenization approach to model the fiber bundles. This significantly simplifies meshing and can avoid large time-consuming simulations. Analytical relationships between permeability and porosity exist for Newtonian flow through regular arrangements of fibers and are commonly used in macroscale porous models by introducing a Darcy viscous term in the flow momentum equations. To this extent, verification of analytical Newtonian permeability–porosity relationships has been conducted for parallel and transverse flow through square and staggered arrangements of fibers. Similar procedures are then used to determine the permeability–porosity relationship for non-Newtonian blood. The results demonstrate that modeling non-Newtonian shear-thinning fluids in porous media can be performed via a generalized Darcy equation with a porous medium viscosity decomposed into a constant term and a directional expression through least squares fitting. This concept is then investigated for various non-Newtonian blood viscosity models. The proposed methodology is conducted with two different porous model approaches, homogeneous and heterogeneous, and validated against a high-fidelity model. The results of the heterogeneous porous model approach yield improved pressure and velocity distribution which highlights the importance of wall effects.

Activated neutrophils have been reported to affect peripheral resistance, for example, by plugging capillaries or adhering to the microvasculature. In vivo and ex vivo data indicate that activated neutrophils circulating in the blood also influence peripheral resistance. We used viscometry and microvascular mimics for in vitro corroboration. The rheological impact of differentiated neutrophil-like HL-60 promyelocytes (dHL60s) or human neutrophil suspensions stimulated with 10 nM fMet-Leu-Phe (fMLP) was quantified using a cone-plate rheometer (450 s−1 shear rate). To evaluate their impact on microscale flow resistance, we used 10-μm Isopore® membranes to model capillaries as well as single 200 × 50 μm microchannels and networks of twenty 20 × 50 μm microfluidic channels to mimic noncapillary microvasculature. Stimulation of dHL60 and neutrophil populations significantly altered their flow behavior as evidenced by their impact on suspension viscosity. Notably, hematocrit abrogated the impact of leukocyte activation on blood cell suspension viscosity. In micropore filters, activated cell suspensions enhanced flow resistance. This effect was further enhanced by the presence of erythrocytes. The resistance of our noncapillary microvascular mimics to flow of activated neutrophil suspensions was significantly increased only with hematocrit. Notably, it was elevated to a higher extent within the micronetwork chambers compared to the single-channel chambers. Collectively, our findings provide supportive evidence that activated neutrophils passing through the microcirculation may alter hemodynamic resistance due to their altered rheology in the noncapillary microvasculature. This effect is another way neutrophil activation due to chronic inflammation may, at least in part, contribute to the elevated hemodynamic resistance associated with cardiovascular diseases (e.g., hypertension and hypercholesterolemia).

This study reformulates Murray's well-known principle of minimum work as applied to the cardiovascular system to include the effects of the shear-thinning rheology of blood. The viscous behavior is described using the extended modified power law (EMPL), which is a time-independent, but shear-thinning rheological constitutive equation. The resulting minimization problem is solved numerically for typical parameter ranges. The non-Newtonian analysis still predicts the classical cubic diameter dependence of the volume flow rate and the cubic branching law. The current analysis also predicts a constant wall shear stress throughout the vascular tree, albeit with a numerical value about 15–25% higher than the Newtonian analysis. Thus, experimentally observed deviations from the cubic branching law or the predicted constant wall shear stress in the vasculature cannot likely be attributed to blood's shear-thinning behavior. Further differences between the predictions of the non-Newtonian and the Newtonian analyses are highlighted, and the limitations of the Newtonian analysis are discussed. Finally, the range and limits of applicability of the current results as applied to the human arterial tree are also discussed.

Bone adapts to habitual loading through mechanobiological signaling. Osteocytes are the primary mechanical sensors in bone, upregulating osteogenic factors and downregulating osteoinhibitors, and recruiting osteoclasts to resorb bone in response to microdamage accumulation. However, most of the cell populations of the bone marrow niche, which are intimately involved with bone remodeling as the source of bone osteoblast and osteoclast progenitors, are also mechanosensitive. We hypothesized that the deformation of trabecular bone would impart mechanical stress within the entrapped bone marrow consistent with mechanostimulation of the constituent cells. Detailed fluid-structure interaction models of porcine femoral trabecular bone and bone marrow were created using tetrahedral finite element meshes. The marrow was allowed to flow freely within the bone pores, while the bone was compressed to 2000 or 3000 microstrain at the apparent level. Marrow properties were parametrically varied from a constant 400 mPa·s to a power-law rule exceeding 85 Pa·s. Deformation generated almost no shear stress or pressure in the marrow for the low viscosity fluid, but exceeded 5 Pa when the higher viscosity models were used. The shear stress was higher when the strain rate increased and in higher volume fraction bone. The results demonstrate that cells within the trabecular bone marrow could be mechanically stimulated by bone deformation, depending on deformation rate, bone porosity, and bone marrow properties. Since the marrow contains many mechanosensitive cells, changes in the stimulatory levels may explain the alterations in bone marrow morphology with aging and disease, which may in turn affect the trabecular bone mechanobiology and adaptation.

It is well known that blood has non-Newtonian properties, but it is generally accepted that blood behaves as a Newtonian fluid at shear rates above 100 s−1. However, in transient conditions, there are times and locations where the shear rate is well below 100 s−1, and it is reasonable to infer that non-Newtonian effects could become important. In this study, purely viscous non-Newtonian (generalized Newtonian) properties of blood are incorporated into the simulation-based framework for cardiovascular surgery planning developed by Taylor et al. (1999, “Predictive Medicine: Computational Techniques in Therapeutic Decision Making,” Comput. Aided Surg., 4, pp. 231–247; 1998, “Finite Element Modeling of Blood Flow in Arteries,” Comput. Methods Appl. Mech. Eng., 158, pp. 155–196). Equations describing blood flow are solved in a patient-based abdominal aortic aneurysm model under steady and physiological flow conditions. Direct numerical simulation (DNS) is used, and the complex flow is found to be constantly transitioning between laminar and turbulent in both the spatial and temporal sense. It is found for the case simulated that using the non-Newtonian viscosity modifies the solution in subtle ways that yield a mesh-independent solution with fewer degrees of freedom than the Newtonian counterpart. It appears that in regions of separated flow, the lower shear rate produces higher viscosity with the non-Newtonian model, which reduces the associated resolution needs. When considering the real case of pulsatile flow, high shear layers lead to greater unsteadiness in the Newtonian case relative to the non-Newtonian case. This, in turn, results in a tendency for the non-Newtonian model to need fewer computational resources even though it has to perform additional calculations for the viscosity. It is also shown that both viscosity models predict comparable wall shear stress distribution. This work suggests that the use of a non-Newtonian viscosity models may be attractive to solve cardiovascular flows since it can provide simulation results that are presumably physically more realistic with at least comparable computational effort for a given level of accuracy.

Steady flow simulations of blood flow in an axisymmetric stenosed artery, subjected to a static magnetic field, are performed to investigate the influence of artery size, magnetic field strength, and non-Newtonian behavior on artery wall shear stress and pressure drop in the stenosed section. It is found that wall shear stress and pressure drop increase by decreasing artery size, assuming non-Newtonian fluid, and increasing magnetic field strength. In the computations, the shear thinning behavior of blood is accounted for by the Carreau–Yasuda model. Computational results are compared and found to be inline with available experimental data.

Understanding structure-function relationships in the temporomandibular joint (TMJ) disc is a critical first step toward creating functional tissue replacements for the large population of patients suffering from TMJ disc disorders. While many of these relationships have been identified for the collagenous fraction of the disc, this same understanding is lacking for the next most abundant extracellular matrix component, sulfated glycosaminoglycans (GAGs). Though GAGs are known to play a major role in maintaining compressive integrity in GAG-rich tissues such as articular cartilage, their role in fibrocartilaginous tissues in which GAGs are much less abundant is not clearly defined. Therefore, this study investigates the contribution of GAGs to the regional viscoelastic compressive properties of the temporomandibular joint (TMJ) disc. Chondroitinase ABC (C-ABC) was used to deplete GAGs in five different disc regions, and the time course for >95% GAG removal was defined. The compressive properties of GAG depleted regional specimens were then compared to non-treated controls using an unconfined compression stress-relaxation test. Additionally, treated and non-treated specimens were assayed biochemically and histologically to confirm GAG removal. Compared to untreated controls, the only regions affected by GAG removal in terms of biomechanical properties were in the intermediate zone, the most GAG-rich portion of the disc. Without GAGs, all intermediate zone regions showed decreased tissue viscosity, and the intermediate zone lateral region also showed a 12.5% decrease in modulus of relaxation. However, in the anterior and posterior band regions, no change in compressive properties was observed following GAG depletion, though these regions showed the highest compressive properties overall. Although GAGs are not the major extracellular matrix molecule of the TMJ disc, they are responsible for some of the viscoelastic compressive properties of the tissue. Furthermore, the mechanical role of sulfated GAGs in the disc varies regionally in the tissue, and GAG abundance does not always correlate with higher compressive properties. Overall, this study found that sulfated GAGs are important to TMJ disc mechanics in the intermediate zone, an important finding for establishing design characteristics for future tissue engineering efforts.

Pulse wave evaluation is an effective method for arteriosclerosis screening. In a previous study, we verified that pulse waveforms change markedly due to arterial stiffness. However, a pulse wave consists of two components, the incident wave and multireflected waves. Clarification of the complicated propagation of these waves is necessary to gain an understanding of the nature of pulse waves in vivo. In this study, we built a one-dimensional theoretical model of a pressure wave propagating in a flexible tube. To evaluate the applicability of the model, we compared theoretical estimations with measured data obtained from basic tube models and a simple arterial model. We constructed different viscoelastic tube set-ups: two straight tubes; one tube connected to two tubes of different elasticity; a single bifurcation tube; and a simple arterial network with four bifurcations. Soft polyurethane tubes were used and the configuration was based on a realistic human arterial network. The tensile modulus of the material was similar to the elasticity of arteries. A pulsatile flow with ejection time 0.3 s was applied using a controlled pump. Inner pressure waves and flow velocity were then measured using a pressure sensor and an ultrasonic diagnostic system. We formulated a 1D model derived from the Navier-Stokes equations and a continuity equation to characterize pressure propagation in flexible tubes. The theoretical model includes nonlinearity and attenuation terms due to the tube wall, and flow viscosity derived from a steady Hagen-Poiseuille profile. Under the same configuration as for experiments, the governing equations were computed using the MacCormack scheme. The theoretical pressure waves for each case showed a good fit to the experimental waves. The square sum of residuals (difference between theoretical and experimental wave-forms) for each case was <10.0%. A possible explanation for the increase in the square sum of residuals is the approximation error for flow viscosity. However, the comparatively small values prove the validity of the approach and indicate the usefulness of the model for understanding pressure propagation in the human arterial network.

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.