Biomedical doctoral students' research practices when facing dilemmas: two vignette-based randomized control trials

Our aim was to describe the research practices of doctoral students facing a dilemma to research integrity and to assess the impact of inappropriate research environments, i.e. exposure to (a) a post-doctoral researcher who committed a Detrimental Research Practice (DRP) in a similar situation and (b) a supervisor who did not oppose the DRP. We conducted two 2-arm, parallel-group randomized controlled trials. We created 10 vignettes describing a realistic dilemma with two alternative courses of action (good practice versus DRP). 630 PhD students were randomized through an online system to a vignette (a) with (n = 151) or without (n = 164) exposure to a post-doctoral researcher; (b) with (n = 155) or without (n = 160) exposure to a supervisor. The primary outcome was a score from - 5 to + 5, where positive scores indicated the choice of DRP and negative scores indicated good practice. Overall, 37% of unexposed participants chose to commit DRP with important variation across vignettes (minimum 10%; maximum 66%). The mean difference [95%CI] was 0.17 [- 0.65 to 0.99;], p = 0.65 when exposed to the post-doctoral researcher, and 0.79 [- 0.38; 1.94], p = 0.16, when exposed to the supervisor. In conclusion, we did not find evidence of an impact of postdoctoral researchers and supervisors on student research practices.Trial registration: NCT04263805, NCT04263506 (registration date 11 February 2020).

Practical implications of the erroneous treatment of exposure time in the Eulerian hemolysis power law model

BACKGROUND

Eulerian and Lagrangian power-law formulations are both widely used for computational fluid dynamics (CFD) to predict flow-induced hemolysis in blood-contacting medical devices. Both are based on the same empirical power-law correlation between hemolysis and the shear stress and exposure time. In the Lagrangian approach, blood damage is predicted by tracking both the stress and exposure time along a finite number of pathlines in the domain. In the Eulerian approach, a scalar transport equation is solved for a time-linearized damage index within the entire domain. Previous analytical work has demonstrated that there is a fundamental problem with the treatment of exposure time in the Eulerian model formulation such that the only condition under which the model correctly represents the true exposure time is in a flow field with no streamwise velocity variation. However, the practical implications of this limitation have yet to be thoroughly investigated.

METHODS

In this study, we demonstrate the inaccuracy of Eulerian hemolysis power-law model predictions due to the erroneous treatment of exposure time by systematically considering four benchmark test cases with increasing degrees of flow acceleration: Poiseuille flow through a straight tube, inclined Couette flow, and flow through a converging tube with two different convergence ratios.

RESULTS

Compared with Lagrangian predictions, we show that, as flow acceleration becomes more pronounced, the resultant inaccuracy in the Eulerian predictions increases. For the inclined Couette flow case, there is a small degree of flow acceleration that yields a discrepancy in the range of 10% between Lagrangian and Eulerian predictions. For flows with a larger degree of acceleration, as occurs in the converging tube flow cases, the discrepancy is considerably larger (up to 257%).

CONCLUSION

The inaccuracy of hemolysis predictions due to the erroneous treatment of exposure time in the Eulerian power-law model can be significant when there is streamwise velocity variation in the flow field. These results may partially explain the extremely large variability in CFD hemolysis predictions reported in the literature between Lagrangian and Eulerian models.

Pulsatile cerebral paraarterial flow by peristalsis, pressure and directional resistance

BACKGROUND

A glymphatic system has been proposed that comprises flow that enters along cerebral paraarterial channels between the artery wall and the surrounding glial layer, continues through the parenchyma, and then exits along similar paravenous channels. The mechanism driving flow through this system is unclear. The pulsatile (oscillatory plus mean) flow measured in the space surrounding the middle cerebral artery (MCA) suggests that peristalsis created by intravascular blood pressure pulses is a candidate for the paraarterial flow in the subarachnoid spaces. However, peristalsis is ineffective in driving significant mean flow when the amplitude of channel wall motion is small, as has been observed in the MCA artery wall. In this paper, peristalsis in combination with two additional mechanisms, a longitudinal pressure gradient and directional flow resistance, is evaluated to match the measured MCA paraarterial oscillatory and mean flows.

METHODS

Two analytical models are used that simplify the paraarterial branched network to a long continuous channel with a traveling wave in order to maximize the potential effect of peristalsis on the mean flow. The two models have parallel-plate and annulus geometries, respectively, with and without an added longitudinal pressure gradient. The effect of directional flow resistors was also evaluated for the parallel-plate geometry.

RESULTS

For these models, the measured amplitude of arterial wall motion is too large to cause the small measured amplitude of oscillatory velocity, indicating that the outer wall must also move. At a combined motion matching the measured oscillatory velocity, peristalsis is incapable of driving enough mean flow. Directional flow resistance elements augment the mean flow, but not enough to provide a match. With a steady longitudinal pressure gradient, both oscillatory and mean flows can be matched to the measurements.

CONCLUSIONS

These results suggest that peristalsis drives the oscillatory flow in the subarachnoid paraarterial space, but is incapable of driving the mean flow. The effect of directional flow resistors is insufficient to produce a match, but a small longitudinal pressure gradient is capable of creating the mean flow. Additional experiments are needed to confirm whether the outer wall also moves, as well as to validate the pressure gradient.

Mechanisms of tracer transport in cerebral perivascular spaces

Tracers infused into the brain appear to be transported along channels surrounding cerebral blood vessels. Bulk fluid flow has been hypothesized in paravascular "glymphatic" channels (outer space between the pial membrane and astrocyte endfeet), as well as in the periarterial space (inner space between smooth muscle cells). The plausibility of net flow in these channels due to steady and oscillatory pressures is reviewed, as is that of transport by oscillatory shear-enhanced dispersion in the absence of net flow. Models including 1D branching networks of annular channels and an expanded compartmental model for humans both predict that flow driven by physiologic steady pressure differences is unlikely in both periarterial and paraarterial spaces, whether the spaces are open or filled with porous media. One exception is that a small additional steady pressure difference could drive paraarterial flow if the space is open. The potential that the tracer injection itself could present such a pressure difference is outlined. Oscillatory (peristaltic) wall motion alone has been found to be insufficient to drive significant forward flow. However, a number of hypothesized mechanisms that have yet to be experimentally verified in the brain may create directional flow in combination with wall motion. Shear-augmented dispersion due to oscillatory pressure in channels with a range of porosity has been modeled analytically. Enhancement of axial dispersion is small for periarterial channels. In open paraarterial channels, dispersion enhancement with optimal lateral mixing is large enough that it may explain observed tracer transport without net forward fluid flow.

On the Discretization of the Power-Law Hemolysis Model

Flow-induced hemolysis remains a concern for blood-contacting devices, and computer-based prediction of hemolysis could facilitate faster and more economical refinement of such devices. While evaluation of convergence of velocity fields obtained by computational fluid dynamics (CFD) simulations has become conventional, convergence of hemolysis calculations is also essential. In this paper, convergence of the power-law hemolysis model is compared for simple flows, including pathlines with exponentially increasing and decreasing stress, in gradually expanding and contracting Couette flows, in a sudden radial expansion and in the Food and Drug Administration (FDA) channel. In the exponential cases, convergence along a pathline required from one to tens of thousands of timesteps, depending on the exponent. Greater timesteps were required for rapidly increasing (large exponent) stress and for rapidly decreasing (small exponent) stress. Example pathlines in the Couette flows could be fit with exponential curves, and convergence behavior followed the trends identified from the exponential cases. More complex flows, such as in the radial expansion and the FDA channel, increase the likelihood of encountering problematic pathlines. For the exponential cases, comparison of converged hemolysis values with analytical solutions demonstrated that the error of the converged solution may exceed 10% for both rapidly decreasing and rapidly increasing stress.

In vitro and numerical simulation of blood removal from cerebrospinal fluid: comparison of lumbar drain to Neurapheresis therapy

BACKGROUND

Blood removal from cerebrospinal fluid (CSF) in post-subarachnoid hemorrhage patients may reduce the risk of related secondary brain injury. We formulated a computational fluid dynamics (CFD) model to investigate the impact of a dual-lumen catheter-based CSF filtration system, called Neurapheresis™ therapy, on blood removal from CSF compared to lumbar drain.

METHODS

A subject-specific multiphase CFD model of CSF system-wide solute transport was constructed based on MRI measurements. The Neurapheresis catheter geometry was added to the model within the spinal subarachnoid space (SAS). Neurapheresis flow aspiration and return rate was 2.0 and 1.8 mL/min, versus 0.2 mL/min drainage for lumbar drain. Blood was modeled as a bulk fluid phase within CSF with a 10% initial tracer concentration and identical viscosity and density as CSF. Subject-specific oscillatory CSF flow was applied at the model inlet. The dura and spinal cord geometry were considered to be stationary. Spatial-temporal tracer concentration was quantified based on time-average steady-streaming velocities throughout the domain under Neurapheresis therapy and lumbar drain. To help verify CFD results, an optically clear in vitro CSF model was constructed with fluorescein used as a blood surrogate. Quantitative comparison of numerical and in vitro results was performed by linear regression of spatial-temporal tracer concentration over 24-h.

RESULTS

After 24-h, tracer concentration was reduced to 4.9% under Neurapheresis therapy compared to 6.5% under lumbar drain. Tracer clearance was most rapid between the catheter aspiration and return ports. Neurapheresis therapy was found to have a greater impact on steady-streaming compared to lumbar drain. Steady-streaming in the cranial SAS was ~ 50× smaller than in the spinal SAS for both cases. CFD results were strongly correlated with the in vitro spatial-temporal tracer concentration under Neurapheresis therapy (R² = 0.89 with + 2.13% and - 1.93% tracer concentration confidence interval).

CONCLUSION

A subject-specific CFD model of CSF system-wide solute transport was used to investigate the impact of Neurapheresis therapy on tracer removal from CSF compared to lumbar drain over a 24-h period. Neurapheresis therapy was found to substantially increase tracer clearance compared to lumbar drain. The multiphase CFD results were verified by in vitro fluorescein tracer experiments.

Impact of Neurapheresis System on Intrathecal Cerebrospinal Fluid Dynamics: A Computational Fluid Dynamics Study

It has been hypothesized that early and rapid filtration of blood from cerebrospinal fluid (CSF) in postsubarachnoid hemorrhage patients may reduce hospital stay and related adverse events. In this study, we formulated a subject-specific computational fluid dynamics (CFD) model to parametrically investigate the impact of a novel dual-lumen catheter-based CSF filtration system, the Neurapheresis™ system (Minnetronix Neuro, Inc., St. Paul, MN), on intrathecal CSF dynamics. The operating principle of this system is to remove CSF from one location along the spine (aspiration port), externally filter the CSF routing the retentate to a waste bag, and return permeate (uncontaminated CSF) to another location along the spine (return port). The CFD model allowed parametric simulation of how the Neurapheresis system impacts intrathecal CSF velocities and steady-steady streaming under various Neurapheresis flow settings ranging from 0.5 to 2.0 ml/min and with a constant retentate removal rate of 0.2 ml/min simulation of the Neurapheresis system were compared to a lumbar drain simulation with a typical CSF removal rate setting of 0.2 ml/min. Results showed that the Neurapheresis system at a maximum flow of 2.0 ml/min increased average steady streaming CSF velocity 2× in comparison to lumbar drain (0.190 ± 0.133 versus 0.093 ± 0.107 mm/s, respectively). This affect was localized to the region within the Neurapheresis flow loop. The mean velocities introduced by the flow loop were relatively small in comparison to normal cardiac-induced CSF velocities.

Anthropomorphic Model of Intrathecal Cerebrospinal Fluid Dynamics Within the Spinal Subarachnoid Space: Spinal Cord Nerve Roots Increase Steady-Streaming

Cerebrospinal fluid (CSF) dynamics are thought to play a vital role in central nervous system (CNS) physiology. The objective of this study was to investigate the impact of spinal cord (SC) nerve roots (NR) on CSF dynamics. A subject-specific computational fluid dynamics (CFD) model of the complete spinal subarachnoid space (SSS) with and without anatomically realistic NR and nonuniform moving dura wall deformation was constructed. This CFD model allowed detailed investigation of the impact of NR on CSF velocities that is not possible in vivo using magnetic resonance imaging (MRI) or other noninvasive imaging methods. Results showed that NR altered CSF dynamics in terms of velocity field, steady-streaming, and vortical structures. Vortices occurred in the cervical spine around NR during CSF flow reversal. The magnitude of steady-streaming CSF flow increased with NR, in particular within the cervical spine. This increase was located axially upstream and downstream of NR due to the interface of adjacent vortices that formed around NR.

Deformation of human red blood cells in extensional flow through a hyperbolic contraction

Flow-induced damage to red blood cells has been an issue of considerable importance since the introduction of the first cardiovascular devices. Early blood damage prediction models were based on measurements of damage by shear stress only. Subsequently, these models were extrapolated to include other components of the fluid stress tensor. However, the expanded models were not validated by measurements of damage in response to the added types of stress. Recent investigations have proposed that extensional stress might be more damaging to red cells than shear stress. In this study, experiments were conducted to compare human red cell deformation under laminar extensional stress versus laminar shear stress. It was found that the deformation caused by shear stress is matched by that produced by an extensional stress that is approximately 34 times smaller. Assuming that blood damage scales directly with cell deformation, this result indicates that mechanistic blood damage prediction models should weigh extensional stress more than shear stress.

Dispersion in porous media in oscillatory flow between flat plates: applications to intrathecal, periarterial and paraarterial solute transport in the central nervous system

Background

As an alternative to advection, solute transport by shear-augmented dispersion within oscillatory cerebrospinal fluid flow was investigated in small channels representing the basement membranes located between cerebral arterial smooth muscle cells, the paraarterial space surrounding the vessel wall and in large channels modeling the spinal subarachnoid space (SSS).

Methods

Geometries were modeled as two-dimensional. Fully developed flows in the channels were modeled by the Darcy–Brinkman momentum equation and dispersion by the passive transport equation. Scaling of the enhancement of axial dispersion relative to molecular diffusion was developed for regimes of flow including quasi-steady, porous and unsteady, and for regimes of dispersion including diffusive and unsteady.

Results

Maximum enhancement occurs when the characteristic time for lateral dispersion is matched to the cycle period. The Darcy–Brinkman model represents the porous media as a continuous flow resistance, and also imposes no-slip boundary conditions at the walls of the channel. Consequently, predicted dispersion is always reduced relative to that of a channel without porous media, except when the flow and dispersion are both unsteady.

Discussion/conclusions

In the basement membranes, flow and dispersion are both quasi-steady and enhancement of dispersion is small even if lateral dispersion is reduced by the porous media to achieve maximum enhancement. In the paraarterial space, maximum enhancement R_max = 73,200 has the potential to be significant. In the SSS, the dispersion is unsteady and the flow is in the transition zone between porous and unsteady. Enhancement is 5.8 times that of molecular diffusion, and grows to a maximum of 1.6E+6 when lateral dispersion is increased. The maximum enhancement produces rostral transport time in agreement with experiments.

Electronic supplementary material

The online version of this article (10.1186/s12987-019-0132-y) contains supplementary material, which is available to authorized users.

On Eulerian versus Lagrangian models of mechanical blood damage and the linearized damage function

Two limitations have been discovered in the derivation of the Eulerian method of hemolysis prediction using a linearized blood damage function. First is that in the transformation from the Lagrangian material volume of the original power-law model to a fixed Eulerian control volume, the spatial dependence of duration of exposure to fluid stress was neglected. This omission has the implication that the Eulerian method as reported is valid only for steady, uniaxial flow in which velocity is constant along streamlines. The second issue is related to linearization, which involves distributing an exponent across an integral. This operation is valid only for limited conditions that include the exponent being unity (which is not the case for any power-law hemolysis models) or the blood damage function being constant throughout the flow regime. These constraints severely restrict the applicability of the Eulerian method. An example problem is presented that demonstrates that the source term of the Eulerian method as reported does not account for differences in velocity between 2 similar flows. Correcting the source term to match the hemolysis prediction to that of the original, unlinearized method requires an analytical description of the flow field that may not be easily obtained for the complex flows in some cardiovascular devices.

Is bulk flow plausible in perivascular, paravascular and paravenous channels?

BACKGROUND

Transport of solutes has been observed in the spaces surrounding cerebral arteries and veins. Indeed, transport has been found in opposite directions in two different spaces around arteries. These findings have motivated hypotheses of bulk flow within these spaces. The glymphatic circulation hypothesis involves flow of cerebrospinal fluid from the cortical subarachnoid space to the parenchyma along the paraarterial (extramural, Virchow-Robin) space around arteries, and return flow to the cerebrospinal fluid (CSF) space via paravenous channels. The second hypothesis involves flow of interstitial fluid from the parenchyma to lymphatic vessels along basement membranes between arterial smooth muscle cells.

METHODS

This article evaluates the plausibility of steady, pressure-driven flow in these channels with one-dimensional branching models.

RESULTS

According to the models, the hydraulic resistance of arterial basement membranes is too large to accommodate estimated interstitial perfusion of the brain, unless the flow empties to lymphatic ducts after only several generations (still within the parenchyma). The estimated pressure drops required to drive paraarterial and paravenous flows of the same magnitude are not large, but paravenous flow back to the CSF space means that the total pressure difference driving both flows is limited to local pressure differences among the different CSF compartments, which are estimated to be small.

CONCLUSIONS

Periarterial flow and glymphatic circulation driven by steady pressure are both found to be implausible, given current estimates of anatomical and fluid dynamic parameters.

Characterization of erythrocyte membrane tension for hemolysis prediction in complex flows

Hemolysis is a persistent issue with blood-contacting devices. Many experimental and theoretical contributions over the last few decades have increased insight into the mechanisms of hemolysis in both laminar and turbulent flows, with the ultimate goal of developing a comprehensive, mechanistic hemolysis model. Many models assume that hemolysis scales with a resultant, scalar stress representing all components of the fluid stress tensor. This study critically evaluates this scalar stress hypothesis by calculating the response of the red blood cell membrane to different types of fluid stress (laminar shear and extension, and three turbulent shear and extension cases), each with the same scalar stress. It was found that even though the scalar stress is the same for all cases, membrane tension varied by up to three orders of magnitude. In addition, extensional flow causes constant tension, while tank-treading in shear flow causes periodic tension, with tank-treading frequency varying by three orders of magnitude among the cases. For turbulent flow, tension also depends on eddy size. It is concluded, therefore, that scalar stress alone is inadequate for scaling hemolysis. Fundamental investigations are needed to establish a new index of the fluid stress tensor that provides reliable hemolysis prediction across the wide range of complex flows that occur in cardiovascular 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.

Extending the Power-Law Hemolysis Model to Complex Flows

Hemolysis (damage to red blood cells) is a long-standing problem in blood contacting devices, and its prediction has been the goal of considerable research. The most popular model relating hemolysis to fluid stresses is the power-law model, which was developed from experiments in pure shear only. In the absence of better data, this model has been extended to more complex flows by replacing the shear stress in the power-law equation with a von Mises-like scalar stress. While the validity of the scalar stress also remains to be confirmed, inconsistencies exist in its application, in particular, two forms that vary by a factor of 2 have been used. This article will clarify the proper extension of the power law to complex flows in a way that maintains correct results in the limit of pure shear.

In vitro validation of flow measurement with phase contrast MRI at 3 tesla using stereoscopic particle image velocimetry and stereoscopic particle image velocimetry-based computational fluid dynamics

PURPOSE

To validate conventional phase-contrast MRI (PC-MRI) measurements of steady and pulsatile flows through stenotic phantoms with various degrees of narrowing at Reynolds numbers mimicking flows in the human iliac artery using stereoscopic particle image velocimetry (SPIV) as gold standard.

MATERIALS AND METHODS

A series of detailed experiments are reported for validation of MR measurements of steady and pulsatile flows with SPIV and CFD on three different stenotic models with 50%, 74%, and 87% area occlusions at three sites: two diameters proximal to the stenosis, at the throat, and two diameters distal to the stenosis.

RESULTS

Agreement between conventional spin-warp PC-MRI with Cartesian read-out and SPIV was demonstrated for both steady and pulsatile flows with mean Reynolds numbers of 130, 160, and 190 at the inlet by evaluating the linear regression between the two methods. The analysis revealed a correlation coefficient of > 0.99 and > 0.96 for steady and pulsatile flows, respectively. Additionally, it was found that the most accurate measures of flow by the sequence were at the throat of the stenosis (error < 5% for both steady and pulsatile mean flows). The flow rate error distal to the stenosis was primarily found to be a function of narrowing severity including dependence on proper Venc selection.

CONCLUSION

SPIV and CFD provide excellent approaches to in vitro validation of new or existing PC-MRI flow measurement techniques.

Space physiology IV: mathematical modeling of the cardiovascular system in space exploration

Mathematical modeling represents an important tool for analyzing cardiovascular function during spaceflight. This review describes how modeling of the cardiovascular system can contribute to space life science research and illustrates this process via modeling efforts to study postflight orthostatic intolerance (POI), a key issue for spaceflight. Examining this application also provides a context for considering broader applications of modeling techniques to the challenges of bioastronautics. POI, which affects a large fraction of astronauts in stand tests upon return to Earth, presents as dizziness, fainting and other symptoms, which can diminish crew performance and cause safety hazards. POI on the Moon or Mars could be more critical. In the field of bioastronautics, POI has been the dominant application of cardiovascular modeling for more than a decade, and a number of mechanisms for POI have been investigated. Modeling approaches include computational models with a range of incorporated factors and hemodynamic sophistication, and also physical models tested in parabolic and orbital flight. Mathematical methods such as parameter sensitivity analysis can help identify key system mechanisms. In the case of POI, this could lead to more effective countermeasures. Validation is a persistent issue in modeling efforts, and key considerations and needs for experimental data to synergistically improve understanding of cardiovascular responses are outlined. Future directions in cardiovascular modeling include subject-specific assessment of system status, as well as research on integrated physiological responses, leading, for instance, to assessment of subject-specific susceptibility to POI or effects of cardiovascular alterations on muscular, vision and cognitive function.

Maximizing information from space data resources: a case for expanding integration across research disciplines

Regulatory systems are affected in space by exposure to weightlessness, high-energy radiation or other spaceflight-induced changes. The impact of spaceflight occurs across multiple scales and systems. Exploring such interactions and interdependencies via an integrative approach provides new opportunities for elucidating these complex responses. This paper argues the case for increased emphasis on integration, systematically archiving, and the coordination of past, present and future space and ground-based analogue experiments. We also discuss possible mechanisms for such integration across disciplines and missions. This article then introduces several discipline-specific reviews that show how such integration can be implemented. Areas explored include: adaptation of the central nervous system to space; cerebral autoregulation and weightlessness; modelling of the cardiovascular system in space exploration; human metabolic response to spaceflight; and exercise, artificial gravity, and physiologic countermeasures for spaceflight. In summary, spaceflight physiology research needs a conceptual framework that extends problem solving beyond disciplinary barriers. Administrative commitment and a high degree of cooperation among investigators are needed to further such a process. Well-designed interdisciplinary research can expand opportunities for broad interpretation of results across multiple physiological systems, which may have applications on Earth.

Percutaneous Transtracheal Ventilation: Resuscitation Bags Do Not Provide Adequate Ventilation

Introduction:

Percutaneous, transtracheal jet ventilation (percutaneous transtracheal jet ventilation) is an effective way to ventilate both adults and children. However, some authors suggest that a resuscitation bag can be utilized to ventilate through a cannula placed into the trachea.

Hypothesis:

Percutaneous transtracheal ventilation (percutaneous transtracheal ventilation) through a 14-gauge catheter is ineffective when attempted using a resuscitation bag.

Methods:

Eight insufflation methods were studied. A 14-gauge intravenous catheter was attached to an adult resuscitation bag, a pediatric resuscitation bag, wall-source (wall) oxygen, portable-tank oxygen with a regulator, and a jet ventilator (JV) at two flow rates. The resuscitation bags were connected to the 14-gauge catheter using a 7 mm adult endotracheal tube adaptor connected to a 3 cc syringe barrel. The wall and tank oxygen were connected to he 14-gauge catheter using a three-way stopcock. The wall oxygen was tested with the regulator set at 15 liters per minute (LPM) and with the regulator wide open. The tank was tested with the regulator set at 15 and 25 LPM. The JV was connected directly to the 14-gauge catheter using JV tubing supplied by the manufacturer. Flow was measured using an Ohmeda 5420 Volume Monitor. A total of 30 measurements were taken, each during four seconds of insufflation, and the results averaged (milliliters (ml) per second (sec)) for each device.

Results:

Flow rates obtained using both resuscitation bags, tank oxygen, and regulated wall oxygen were extremely low (adult 215 ±20 ml/sec; pediatric 195 ±19 ml/sec; tank 358 ±13 ml/sec; wall at 15 l/min 346 ±20 ml/sec). Flow rates of 1,394 ±13 ml were obtained using wall oxygen with the regulator wide open. Using the JV with the regulator set at 50 pounds per square inch (psi), a flow rate of 1,759 ±40 was obtained.These were the only two methods that produced flow rates high enough to provide an adequate tidal volume to an adult.

Conclusions:

Resuscitation bags should not be used to ventilate adult patients through a 14-gauge, transtracheal catheter. Jet ventilation is needed when percutaneous transtracheal ventilation is attempted. If jet ventilation is attempted using oxygen supply tubing, it must be connected to an unregulated oxygen source of at least 50 psi.

Effects of biaxial oscillatory shear stress on endothelial cell proliferation and morphology

Wall shear stress (WSS) on anchored cells affects their responses, including cell proliferation and morphology. In this study, the effects of the directionality of pulsatile WSS on endothelial cell proliferation and morphology were investigated for cells grown in a Petri dish orbiting on a shaker platform. Time and location dependent WSS was determined by computational fluid dynamics (CFD). At low orbital speed (50 rpm), WSS was shown to be uniform (0-1 dyne/cm(2)) across the bottom of the dish, while at higher orbital speed (100 and 150 rpm), WSS remained fairly uniform near the center and fluctuated significantly (0-9 dyne/cm(2)) near the side walls of the dish. Since WSS on the bottom of the dish is two-dimensional, a new directional oscillatory shear index (DOSI) was developed to quantify the directionality of oscillating shear. DOSI approached zero for biaxial oscillatory shear of equal magnitudes near the center and approached one for uniaxial pulsatile shear near the wall, where large tangential WSS dominated a much smaller radial component. Near the center (low DOSI), more, smaller and less elongated cells grew, whereas larger cells with greater elongation were observed in the more uniaxial oscillatory shear (high DOSI) near the periphery of the dish. Further, cells aligned with the direction of the largest component of shear but were randomly oriented in low magnitude biaxial shear. Statistical analyses of the individual and interacting effects of multiple factors (DOSI, shear magnitudes and orbital speeds) showed that DOSI significantly affected all the responses, indicating that directionality is an important determinant of cellular responses.

Modelling of cardiovascular response to graded orthostatic stress: role of capillary filtration

BACKGROUND

  Of the many possible factors that may contribute to orthostatic intolerance, the loss of circulating blood because of capillary filtration is one of the few that can explain the gradual decline of arterial pressure during stand tests. This study used a computer model to investigate the relative importance of haemodynamic parameters, including capillary filtration, as potential contributors to orthostatic intolerance. Simulated orthostatic tolerance times were compared to previous experiments combining head-up tilt and lower body negative pressure graded orthostatic stress, which provided haemodynamic data, in particular haematocrit measurements that allowed subject-specific modelling of capillary transport.

MATERIALS AND METHODS

  The cardiovascular system was simulated using a seven-compartment model with measured heart rate, stroke volume, total peripheral resistance, mean arterial pressure and haematocrit data for 12 subjects. Simulations were controlled by decreasing the total blood volume at the measured rates of capillary filtration until cerebral pressure dropped below a threshold for consciousness. Predicted times to syncope were compared to actual times to presyncope, and sensitivity of arterial pressure and cardiac output to independent system parameters were determined.

RESULTS

  There was no statistical difference in modelled times to syncope and actual times to presyncope. Both arterial pressure and cardiac output were most sensitive to total blood volume and least sensitive to caudal compliance parameters.

CONCLUSIONS

  The feasibility of subject-specific simulations of cardiovascular response to orthostatic stress was demonstrated, providing stronger evidence that capillary filtration is a prominent mechanism in causing orthostatic intolerance. These results may have clinical and spaceflight applications.

Evaluation of Mechanisms of Postflight Orthostatic Intolerance with a Simple Cardiovascular System Model

A significant fraction of astronauts experience postflight orthostatic intolerance (POI) during 10-min stand tests conducted on landing day. The average time that nonfinishers can stand is about 7 min. This phenomenon, including the delay in occurrence of presyncope, was studied with a five-compartment model of the cardiovascular system incorporating compartments for the heart/lungs, systemic arteries and cephalic, central, and caudal veins. The model included 28 independent parameters, including factors characterizing cardiac performance, vascular resistance, intrathoracic pressure, nonlinear venous compliance and circulating blood volume, and 13 dependent parameters, including cardiac output and cardiac and vascular compartment pressures and volumes. First, a sensitivity analysis of hemodynamic indicators of presyncope to independent parameters was performed. Results demonstrated that both cardiac output and arterial pressure were most sensitive to volume-related parameters, particularly total blood volume, and less sensitive to peripheral resistance. Next, a simulated postflight stand test confirmed that fluid loss due to capillary filtration, particularly from the caudal region where transmural pressure is high during standing, is a plausible mechanism of POI that also explains the delayed onset of symptoms in most astronauts. An accumulated drop in arterial pressure sufficient to compromise cerebral perfusion and, therefore, cause syncope was reached in about 7 min with a fluid loss of 280 mL. Finally, additional simulations showed that a 75% increase in peripheral resistance, similar to finishers of stand tests, was insufficient to overcome the loss of circulating fluid associated with capillary filtration, and extended the time that the modeled astronaut could stand by only about 1 min. It is therefore concluded that capillary filtration may play a key role in producing POI and that development of countermeasures should perhaps focus on reducing postflight capillary permeability or on stimulating volume-compensating mechanisms.

The Effect of Gravitational Acceleration on Cardiac Diastolic Function: A Biofluid Mechanical Perspective with Initial Results

Echocardiographic measurements of astronaut cardiac function have documented an initial increase, followed by a progressive reduction in both left ventricular end-diastolic volume index and stroke volume with entry into microgravity (micro-G). The investigators hypothesize that the observed reduction in cardiac filling may, in part, be due to the absence of a gravitational acceleration dependent, intraventricular hydrostatic pressure difference in micro-G that exists in the ventricle in normal gravity (1-G) due to its size and anatomic orientation. This acceleration-dependent pressure difference, DeltaP(LV), between the base and the apex of the heart for the upright posture can be estimated to be 6660 dynes/cm(2) ( approximately 5 mm Hg) on Earth. DeltaP(LV) promotes cardiac diastolic filling on Earth, but is absent in micro-G. If the proposed hypothesis is correct, cardiac pumping performance would be diminished in micro-G. To test this hypothesis, ventricular function experiments were conducted in the 1-G environment using an artificial ventricle pumping on a mock circulation system with the longitudinal axis anatomically oriented for the upright posture at 45 degrees to the horizon. Additional measurements were made with the ventricle horizontally oriented to null DeltaP(LV)along the apex-base axis of the heart as would be the case for the supine posture, but resulting in a lesser hydrostatic pressure difference along the minor (anterior-posterior) axis. Comparative experiments were also conducted in the micro-G environment of orbital space flight on board the Space Shuttle. This paper reviews the use of an automated cardiovascular simulator flown on STS-85 and STS-95 as a Get Away Special payload to test this hypothesis. The simulator consisted of a pneumatically actuated, artificial ventricle connected to a closed-loop, fluid circuit with adjustable compliance and resistance elements to create physiologic pressure and flow conditions. Ventricular instrumentation included pressure transducers in the apex and base as well as immediately upstream of the inflow valve and downstream of the outflow valve, and a flow probe downstream of the outflow valve. By varying the circulating fluid volume, ventricular function could be determined for varying preload pressures at a regulated, mean afterload pressure of 95 mm Hg. This variation in preload condition permitted the construction of a ventricular function curve for the micro-G environment for comparison to the same curve for the 1-G environment. Data were collected from both missions at the upper end of the ventricular function curve. Experiment operation in the 1-G, supine orientation or in the micro-G environment eliminated the DeltaP(LV) observed in the 1-G, upright orientation. Consistent with the hypothesis, additional atrial pressure was required in micro-G to obtain stroke volumes and flow rates similar to those measured in 1-G for the upright posture. The necessary increase in atrial pressure was approximately 5 mm Hg in these experiments. In the same range of flow rates and stroke volumes, similar flows were observed in the 1-G supine posture for atrial pressures intermediate to the 1-G upright and micro-G values, also consistent with the hypothesis. Additional experiments on board the Space Shuttle are in preparation to gather data across the rest of the normal physiologic range of the ventricular function curve.

Hemolysis in needleless connectors for phlebotomy

Needleless connectors have been developed recently as a means of reducing transmission of AIDS and other blood-borne diseases by accidental needle sticks. However, the potential for hemolysis induced by fluid stresses within the connector remained to be determined. The influence of needleless connectors on hemolysis was evaluated in simulated clinical blood draws. Blood from five volunteers was drawn with vacuum tubes and syringes through 16, 18, and 22 gauge needles and PosiFlow (Becton Dickinson, Franklin Lakes, NJ) and Clave (ICU Medical, San Clemente, CA) needleless connectors. Hemolysis was measured in the samples by a spectrophotometric technique. Results showed that hemolysis increased when needleless connectors were used. The PosiFlow connector produced more hemolysis than the Clave device. Curiously, erythrocyte damage was greatest for connectors used with 18 gauge needles. Hemolysis was larger for samples drawn with vacuum tubes than with syringes. However, no combinations of connector, needle size, or blood draw device resulted in mean hemolysis values large enough to interfere with clinical assays.

Numerical simulation of the influence of gravity and posture on cardiac performance

A numerical model of the cardiovascular system was used to quantify the influences on cardiac function of intrathoracic pressure and intravascular and intraventricular hydrostatic pressure, which are fundamental biomechanical stimuli for orthostatic response. The model included a detailed arterial circulation with lumped parameter models of the atria, ventricles, pulmonary circulation, and venous circulation. The venous circulation was divided into cranial, central, and caudal regions with nonlinear compliance. Changes in intrathoracic pressure and the effects of hydrostatic pressure were simulated in supine, launch, sitting, and standing postures for 0, 1, and 1.8 G. Increasing intrathoracic pressure experienced with increasing gravity caused 12% and 14% decreases in cardiac output for 1 and 1.8 G supine, respectively, compared to 0 G. Similar results were obtained for launch posture, in which the effects of changing intrathoracic pressure dominated those of hydrostatic pressure. Compared to 0 G, cardiac output decreased 0.9% for 1 G launch and 15% for 1.8 G launch. In sitting and standing, the position of the heart above the hydrostatic indifference level caused the effects of changing hydrostatic pressure to dominate those of intrathoracic pressure. Compared to 0 G, cardiac output decreased 13% for 1 G sitting and 23% for 1.8 G sitting, and decreased 17% for 1 G standing and 31% for 1.8 G standing. For a posture change from supine to standing in 1 G, cardiac output decreased, consistent with the trend necessary to explain orthostatic intolerance in some astronauts during postflight stand tests. Simulated lower body negative pressure (LBNP) in 0 G reduced cardiac output and mean aortic pressure similar to I G standing, suggesting that LBNP provides at least some cardiovascular stimuli that may be useful in preventing postflight orthostatic intolerance. A unifying concept, consistent with the Frank-Starling mechanism of the heart, was that cardiac output was proportional to cardiac diastolic transmural pressure for all postures and gravitational accelerations.

Stability of elliptical cylinders in two-dimensional channel flow

The flow around rigid cylinders of elliptical cross section placed transverse to Poiseuille flow between parallel plates was simulated to investigate issues related to the tumbling of red blood cells and other particles of moderate aspect ratio in the similar flow in a Field Flow Fractionation (FFF) channel. The torque and transverse force on the cylinder were calculated with the cylinder freely translating, but prevented from rotating, in the flow. The aspect ratios (long axis to short axis) of the elliptical cylinders were 2, 3, 4, and 5. The cylinder was placed transversely at locations of y0/H = 0.1, 0.2, 0.3, and 0.4, where y0 is the distance from the bottom of the channel and H is the height of the channel, and the orientation of the cylinder was varied from 0 to 10 deg with respect to the axis of the channel for a channel Reynolds number of 20. The results showed that equilibrium orientations (indicated by a zero net torque on the cylinder) were possible for high-aspect-ratio cylinders at transverse locations y0/H < 0.2. Otherwise, the net torque on the cylinder was positive, indicating that the cylinder would rotate. For cylinders with a stable orientation, however, a transverse lift forced existed up to about y0/H = 0.25. Thus, a cylinder of neutral or low buoyancy might be lifted with a stable orientation from an initial position near the wall until it reached y0/H < 0.2, whereupon it would begin to tumble or oscillate. The dependence of lift and torque on cylinder orientation suggested that neutral or low-buoyancy cylinders may oscillate in both transverse location and angular velocity. Cylinders more dense than the carrier fluid could be in equilibrium both in terms of orientation and transverse location if their sedimentation force matched their lift force for a location y0/H < 0.2.

Aortic input impedance in infants and children

Flow and pressure measurements were performed in the ascending aortas of six pediatric patients ranging in age from 1 to 4 yr and in weight from 7.2 to 16.4 kg. From these measurements, input impedance was calculated. It was found that total vascular resistance decreased with increasing patient weight and was approximately one to three times higher than those of adults. Conductance per unit weight was relatively constant but was approximately three times higher than for adults. Strong inertial character was observed in the impedance of four of the six patients. Among a three-element and two four-element lumped-parameter models, the model with characteristic aortic resistor (R(c)) and inertance in series followed by parallel peripheral resistor (R(p)) and compliance fitted the data best. R(p) decreased with increasing patient weight and was one to three times higher than in adults, and R(c) decreased with increasing patient weight and was 2 to 15 times higher. The R(p)-to-R(c) ratio differed significantly between infants and children vs. adults. The results suggested that R(p) developed more rapidly with patient weight than did R(c). Compliance values increased with increasing patient weight and were 3 to 16 times lower than adult values.

Feedback control of mean aortic pressure in a dynamic model of the cardiovascular system

Orbital measurements of the cardiac function of Space Shuttle crew members have shown an initial increase in cardiac stroke volume upon entry into weightlessness, followed by a gradual reduction in stroke volume to a level approximately 15% less than preflight values. In an effort to explain this response, it was hypothesized that gravity plays a role in cardiac filling. A mock circulatory system was designed to investigate this effect. Preliminary studies carried out with this system on the NASA KC-135 aircraft, which provides brief periods of weightlessness, showed a strong correlation between cardiac filling, stroke volume, and the presence or absence of gravity. The need for extended periods of high quality zero gravity was identified to verify this observation. To accomplish this, the aircraft version of the experiment was reduced in size and fully automated for eventual integration into a Get Away Special canister to conduct an orbital version of the experiment. This article describes the automated system, as well as the development and implementation of a control algorithm for the servoregulation of the mean aortic pressure in the orbital experiment. Three nonlinearities that influence the ability of the apparatus to regulate to a mean aortic pressure of 95 mm Hg were identified and minimized. In preparation for a Space Shuttle flight, the successful function of the servoregulatory scheme was demonstrated during ground tests and additional test flights aboard the KC-135. The control algorithm was successful in carrying out the experimental protocol, including regulation of mean aortic pressure. The algorithm could also be used for the automated operation of long-term tests of circulatory support systems, which may require a scheduled cycling of the pumping conditions on a daily basis.

Development of a hydraulic model of the human systemic circulation

Physical and numeric models of the human circulation are constructed for a number of objectives, including studies and training in physiologic control, interpretation of clinical observations, and testing of prosthetic cardiovascular devices. For many of these purposes it is important to quantitatively validate the dynamic response of the models in terms of the input impedance (Z = oscillatory pressure/oscillatory flow). To address this need, the authors developed an improved physical model. Using a computer study, the authors first identified the configuration of lumped parameter elements in a model of the systemic circulation; the result was a good match with human aortic input impedance with a minimum number of elements. Design, construction, and testing of a hydraulic model analogous to the computer model followed. Numeric results showed that a three element model with two resistors and one compliance produced reasonable matching without undue complication. The subsequent analogous hydraulic model included adjustable resistors incorporating a sliding plate to vary the flow area through a porous material and an adjustable compliance consisting of a variable-volume air chamber. The response of the hydraulic model compared favorably with other circulation models.

Finite element analysis of the lift on a slightly deformable and freely rotating and translating cylinder in two-dimensional channel flow

Motivated by the lateral migration phenomena of fresh and glutaraldehyde-fixed red blood cells in a field flow fractionation (FFF) separation system, we studied the transverse hydrodynamic lift on a slightly flexible cylinder in a two-dimensional channel flow. The finite element method was used to analyze the flow field with the cylinder at different transverse locations in the channel. The shape of the cylinder was determined by the pressure on the surface of the cylinder from the flow field solution and by the internal elastic stress. The cylinder deformation and the flow field were solved simultaneously. The transverse lift exerted on the cylinder was then calculated. The axial and angular speed of the cylinder were iterated such that the drag and torque on the cylinder were nulled to represent a freely translating and rotating state. The results showed that the transverse lift on a deformable cylinder increased greatly and the equilibrium position moved closer to the center of the channel compared to a rigid cylinder. Also, with the same elastic modulus but a higher flow rate, a larger deformation and higher equilibrium location were found. The maximum deformation of the cylinder occurred when the cylinder was closest to the wall where a larger shear rate existed. The numerical results and experimental studies are discussed.

Influence of gravity on cardiac performance

Results obtained by the investigators in ground-based experiments and in two parabolic flight series of tests aboard the NASA KC-135 aircraft with a hydraulic simulator of the human systemic circulation have confirmed that a simple lack of hydrostatic pressure within an artificial ventricle causes a decrease in stroke volume of 20%-50%. A corresponding drop in stroke volume (SV) and cardiac output (CO) was observed over a range of atrial pressures (AP), representing a rightward shift of the classic CO versus AP cardiac function curve. These results are in agreement with echocardiographic experiments performed on space shuttle flights, where an average decrease in SV of 15% was measured following a three-day period of adaptation to weightlessness. The similarity of behavior of the hydraulic model to the human system suggests that the simple physical effects of the lack of hydrostatic pressure may be an important mechanism for the observed changes in cardiac performance in astronauts during the weightlessness of space flight.

Scaling of hemolysis in needles and catheters

Hemolysis in clinical blood samples leads to inaccurate assay results and often to the need for repeated blood draws. In vitro experiments were conducted to determine the influence on hemolysis in phlebotomy needles and catheters of pressure difference, cannula diameter, and cannula material. Fresh blood from five human volunteers was forced from a syringe inside a pressurized chamber through 14, 18, and 22 gauge 304 stainless steel needles and polyurethane and Teflon catheters, all 40 mm long. Hemolysis was measured in the samples by a spectrophotometer. It was found that hemolysis increased with increases in pressure difference and cannula diameter and no consistent trend could be identified with regard to cannula material. The pressure differences required for significant hemolysis were above those typical of clinical venipuncture blood draws. While there was substantial variability among individuals, the hemolysis values scaled with exponent S = (t/t0)[(tau/tau0)-1]2, where t is the characteristic duration of shear, t0 is a time constant, tau is the wall shear stress, and tau0 is the wall shear stress threshold below which no hemolysis occurs. A hemolysis threshold including both time and shear stress was also defined for S = constant. The threshold implies that a threshold shear stress exists below which erythrocytes are not damaged for any length of exposure time, but that red cells may be damaged by an arbitrarily short period of exposure to sufficiently large shear stress.

Compact compliance chamber design for the study of cardiac performance in microgravity

A need was identified for a Mock Circulation System (MCS) of small size and weight that could function in a microgravity environment for the investigation of cardiovascular response to the weightlessness of space flight. Part of the MCS development involved the redesign of the compliance chamber from a Penn State MCS using a coil spring instead of the leaf spring system employed in the Penn State system. The new compliance chambers achieve a weight reduction of 47% and a volume reduction of 64% over the original Penn State design. Testing showed the coil spring compliance chambers retained physiologic characteristics and adjustability by using coil springs of various stiffness, and functioned equivalently to the original Penn State design.

A blood analog for laser-induced photochemical anemometry

A transparent viscoelastic blood analog fluid was developed for use with Laser-Induced Photochemical Anemometry. To provide solubility of the photochemical tracer, 1', 3', 3'-trimethyl-6-nitroindoline-2-benzopyran (TNSB dye, Kodak Chemicals), the analog solvent needed to be nonpolar, thus currently available aqueous blood analogs were not suitable. An analog consisting of 0.04% ethylhydroxyethylcellulose dissolved in gamma-butyrolactone produced a pseudoplastic steady shear response with low elasticity in unsteady shear, while being compatible with the photochemical tracer.

The effect of blood viscoelasticity on pulsatile flow in stationary and axially moving tubes

An analytical solution for pulsatile flow of a generalized Maxwell fluid in straight rigid tubes, with and without axial vessel motion, has been used to calculate the effect of blood viscoelasticity on velocity profiles and shear stress in flows representative of those in the large arteries. Measured bulk flow rate Q waveforms were used as starting points in the calculations for the aorta and femoral arteries, from which axial pressure gradient delta P waves were derived that would reproduce the starting Q waves for viscoelastic flow. The delta P waves were then used to calculate velocity profiles for both viscoelastic and purely viscous flow. For the coronary artery, published delta P and axial vessel acceleration waveforms were used in a similar procedure to determine the separate and combined influences of viscoelasticity and vessel motion. Differences in local velocities, comparing viscous flow to viscoelastic flow, were in all cases less than about 2% of the peak local velocity. Differences in peak wall shear stress were less than about 3%. In the coronary artery, wall shear stress differences between viscous and viscoelastic flow were small, regardless of whether axial vessel motion was included. The shape of the wall shear stress waveform and its difference, however, changed dramatically between the stationary and moving vessel cases. The peaks in wall shear stress difference corresponded with large temporal gradients in the combined driving force for the flow.

An orbiting scroll blood pump without valves or rotating seals

Valves in blood pumps are expensive and provide modes of failure. Rotating seals offer sites of thrombus formation and infection. In this study, a prototype pump incorporating no valves or rotating seals was constructed and tested. In this device, fluid is pumped by the orbital action of a spiral shaped scroll relative to an identical stationary scroll whose starting axis is rotated 180 degrees with respect to the orbiting scroll. The two scrolls, which are machined integral with scroll plates, form pockets that are filled from the outside and then ejected in the center as the orbiting scroll completes each cycle. The orbiting scroll is driven by a crank mechanism connected to a motor. Fluid is contained in the space around the scrolls by a flexible collar and does not contact the driving mechanism. The prototype pump is approximately 7.6 cm in diameter and 2.5 cm thick and has an orbiting radius of 5.1 mm. The output of the pump was very sensitive to the clearance between the scroll tip and the base of the opposite scroll plate. For a clearance of 51 microns, pressure differences as high as 400 mmHg and flows as high as 7.7 l/min (of water) were produced at 260 rpm. At 450 rpm with a 330 microns clearance, pressure differences as high as 185 mmHg and flows as high as 7.3 l/min resulted. The relationships between pressure difference and flow were very linear in all cases. Volumetric efficiency was as high as 70% and increased with speed.

Shear-augmented dispersion in non-Newtonian fluids

The rate of spread of a passive species is modified by the superposition of a velocity gradient on the concentration field. Taylor (18) solved for the rate of axial dispersion in fully developed steady Newtonian flow in a straight pipe under the conditions that the dispersion be relatively steady and that longitudinal transport be controlled by convection rather than diffusion. He found that the resulting effective axial diffusivity was proportional to the square of the Peclet number Pec and inversely proportional to the molecular diffusivity. This article shows that under similar conditions in Casson and power law fluids, both simplified models for blood, and in Bingham fluids the same proportionalities are found. Solutions are presented for fully developed steady flow in a straight tube and between flat plates. The proportionality factor, however, is dependent upon the specific rheology of the fluid. For Bingham and Casson fluids, the controlling parameter is the radius of the constant-velocity core in which the shear stress does not exceed the yield stress of the fluid. For a core radius of one-tenth the radius of the tube, the effective axial diffusivity in Casson fluids is reduced to approximately 0.78 times that in a Newtonian fluid at the same flow. Using average flow conditions, it is found that the core radius/tube radius ratio is 0 (10(-2)) to 0 (10(-1)) in canine arteries and veins. Even at these small values, the effective diffusivity is diminished by 5% to 18%. For power law fluids, Pec2 dependence is again found, but with a proportionality constant dependent upon the power law exponent n.(ABSTRACT TRUNCATED AT 250 WORDS)

Influence of viscosity and pressure on prosthetic valve regurgitation

Blood viscosity varies during the course of artificial heart implants and is affected by pathological conditions. To gauge the potential effect of changing viscosity on valve performance, leakage rates were measured across a closed Medtronic-Hall valve with water, water/glycerol and fresh whole bovine blood for aortic and pulmonary pressure ranges. As might be expected from the low Reynolds numbers (< 140), losses across the valve were found to be primarily viscous. For the two Newtonian fluids, leakage was slightly less than linearly proportional to pressure. This is comparable with empirical data for orifice flow, which predicts three fifths power dependence on pressure. For blood, however, the greater than linear dependence on pressure found suggests that the pseudoplasticity (shear-thinning behavior) of blood is important. These data provide evidence that the viscous and non-Newtonian properties of blood must be taken into account in modelling prosthetic valve performance and may affect the test methods and flow regulation strategies for prosthetic blood pumps.

Sensitivity of the artificial heart to changes in vascular resistance

For conditions near those for normal operation, the cardiac output of the healthy natural heart is inherently sensitive to systemic venous resistance but is relatively insensitive to arterial resistance. A mathematical comparison was undertaken to discover the differences in the sensitivity of two configurations of artificial hearts to these resistances. It was found that one design incorporating independently pumping ventricles can be tailored to passively mimic the sensitivity of the natural heart. However, the other design with volumetrically coupled pumping is incapable of exhibiting similar sensitivity through passive means.