Plug Effect of Erythrocytes in Capillary Blood Vessels

As an idealized problem of the motion of blood in small capillary blood vessels, the low Reynolds number flow of plasma (a newtonian fluid) in a circular cylindrical tube involving a series of circular disks is studied. It is assumed in this study that the suspended disks are equally spaced along the axis of the tube, and that their centers remain on the axis of the tube and that their faces are perpendicular to the tube axis. The inertial force of the fluid due to the convective acceleration is neglected on the basis of the smallness of the Reynolds number. The solution of the problem is derived for a quasi-steady flow involving infinitesimally thin disks. The numerical calculation is carried out for a set of different combinations of the interdisk distance and the ratio of the disk radius to the tube radius. The ratio of the velocity of the disk to the average velocity of the fluid is calculated. The different rates of transport of red blood cells and of plasma in capillary blood vessels are discussed. The average pressure gradient along the axis of the tube is computed, and the dependence of the effective viscosity of the blood on the hematocrit and the diameter of the capillary vessel is discussed.     

Measurement of oscillatory flow pressure gradient in an elastic artery model

In vitro experiments were conducted to measure the oscillatory flow pressure gradient along an elastic tube in order to assess the recent nonlinear theory of Wang and Tarbell. According to this theory, in an elastic tube with oscillatory flow, the mean (time-averaged) pressure gradient cannot be calculated using Poiseuille's law. The effect of wall motion creates a nonlinear convective acceleration, and an induced mean pressure gradient is required to balance the convective acceleration. The induced mean pressure gradient depends on the diameter variation over a cycle, the pulsatility and unsteadiness of the flow, and the phase difference between the pressure wave form and the flow wave form. The amplitude of the pressure gradient also depends on these parameters and may deviate significantly from Womersley's rigid tube theory. A flow loop was constructed to produce oscillatory flow in an elastic tube. Flow wave forms were measured with an ultrasonic flow probe, and ultrasonic diameter crystals were used to measure wall movement. A special device for pressure drop measurement was constructed using Millar catheter tip transducers to obtain both forward and backward pressure drops that were then averaged. This averaging method eliminated the static error of the pressure transducers. The pressure-flow phase angle was varied by clamping a distal elastic section at various locations. Pressure gradients were obtained for a range of phase angles between −55 ° and +49 °. The mean and amplitude of the measured pressure gradient were compared to theoretical values. Both positive and negative induced mean pressure gradients were measured over the range of phase angles. The measured pressure gradient amplitudes were always lower than predicted by Womersley's rigid tube theory. The experimental means and amplitudes are in good agreement with the elastic tube theoretical values. Thus, the experiments verify the theory of Wang and Tarbell.     

On the steady laminar flow of a viscous incompressible fluid in an elastic tube

The conditions are examined under which an approximate relation between the radius of the tube and the distance along the axis, as obtained by N. Rashevsky (1945), is consistent with the assumptions made in the solution. These conditions are reduced to relations between two dimensionless parameters of the system. Second approximations are found for three different ranges of values of these parameters.     

Hydrodynamics in the renal tubule

Explicit solutions of the Navier-Stokes equations are presented for axially symmetric slow flow in an infinite cylinder whose walls reabsorb fluid at a rate which varies exponentially with the longitudinal coordinate. Results similar to those of a previous paper which assumed a constant rate of reabsorption are obtained. When the radius of the tube is small the solutions resemble Poiseuille flow; the longitudinal velocity profile is parabolic, and the drop in mean pressure is proportional to the mean axial flow, the length of tube between reference points, and inversely proportional to the fourth power of the radius. By expressing the tubular reabsorption as a Fourier integral, solutions are obtained for the general case where the rate of reabsorption is an arbitrary function of the longitudinal coordinate.     

Velocity field of pulsatile flow in a porous tube

This paper describes velocity fields for fully developed periodic laminar flow in a rigid tube with a porous wall. We obtained an analytical solution of the flow by the linear approximation of the Navier-Stokes equation. Unlike the previous works with a constant seepage rate along the axis, we used a wall velocity which contained hydraulic permeation constant Lp. The axial velocity profile shows a local maximum velocity near the wall at a large Womersley number alpha. This suggests that concentration polarization in porous tubular membrane may be reduced at high frequencies if a membrane device is operated under pulsatile flow conditions. The magnitude of wall permeation velocity decreases linearly along the tube axis because the damping of the pressure difference between the inside and the outside of the tube is very small.     

Contributions to the mathematical biophysics of the cardiovascular system

The reflection of pressure waves in a fluid enclosed within a tube with an elastic wall is studied for the case of a localized change in diameter of the tube. The concept of impedance is introduced. The relation of the reflection characteristics of the parts of the tube at either side of the change is derived on the basis of the continuity of pressure and mass flow at the site of the change. This relations is used to derive the expression for the ratio of the pressure oscillations measured in front of, and behind, the constriction in terms of the constants of the system. As a result, a method is indicated to locate the coarctation from measurements of the pressures in front of, and behind it.     

A nonlinear axisymmetric model with fluid-wall interactions for steady viscous flow in stenotic elastic tubes.

Arteries with high-grade stenoses may compress under physiologic conditions due to negative transmural pressure caused by high-velocity flow passing through the stenoses. To quantify the compressive conditions near the stenosis, a nonlinear axisymmetric model with fluid-wall interactions is introduced to simulate the viscous flow in a compliant stenotic tube. The nonlinear elastic properties of the tube (tube law) are measured experimentally and used in the model. The model is solved using ADINA (Automatic Dynamic Incremental Nonlinear Analysis), which is a finite element package capable of solving problems with fluid-structure interactions. Our results indicate that severe stenoses cause critical flow conditions such as negative pressure and high and low shear stresses, which may be related to artery compression, plaque cap rupture, platelet activation, and thrombus formation. The pressure filed near a stenosis has a complex pattern not seen in one-dimensional models. Negative transmural pressure as low as -24 mmHg for a 78 percent stenosis by diameter is observed at the throat of the stenosis for a downstream pressure of 30 mmHg. Maximum shear stress as a high as 1860 dyn/cm2 occurs at the throat of the stenoses, while low shear stress with reversed direction is observed right distal to the stenosis. Compressive stresses are observed inside the tube wall. The maximal principal stress and hoop stress in the 78 percent stenosis are 80 percent higher than that from the 50 percent stenosis used in our simulation. Flow rates under different pressure drop conditions are calculated and compared with experimental measurements and reasonable agreement is found for the prebuckling stage.     

Estimated flow resistance increase in a spiral human coronary artery segment

Coronary flow estimates were made for a spiral coronary artery segment (identified from a post-mortem replica casting) by using a modified Dean number based on the approximate coil radius of curvature, as suggested earlier. The estimates were found to correlate experimental pressure drop data for helical coiled tubes. Over a physiological range of mean Reynolds numbers from 100 to 400 for blood flow through main coronary arteries, estimates of the flow resistance increase relative to a straight lumen segment ranged from about 20 to 80 percent, and were of similar magnitude to those found in a flow study in a sinuous coronary vessel segment with no spiral.     

The Influence of Drag on Nonlinear Oscillatory Flow through Concentric Annulus

A mathematical model has been developed to study the effect of particle drag parameter and frequency parameter on velocity and pressure gradient in nonlinear oscillatory two phase flow. The main purpose is to apply the model to study the combined effect of introduction of the catheter and elastic properties of the arterial wall on the pulsatile nature of the blood flow. We model the artery as an isotropic thin walled elastic tube and the catheter as a coaxial flexible tube. Blood is modeled as an incompressible particulate viscous Newtonian fluid. Perturbation technique has been applied to find the approximations for velocity and pressure gradient up to second order. Numerical solutions are investigated with graphical presentations to understand the effects of drag parameter, frequency parameter and phase angle on velocity along radial direction and pressure gradient along axial directions. As the drag parameter increases, mean pressure gradient and mean velocity will be decreased. As frequency parameter increases mean velocity profile bends near the outer wall. Due to elastic nature of artery wall, a thin catheter experience small oscillations and a thick catheter remains stationary inside the artery. Finally, the effect of catheterization on various physiologically important flow rate characteristics—mean velocity, mean pressure gradient are studied for a range of different catheter sizes, particle drag parameter and frequency parameters.     

The velocity of the arterial pulse wave: a viscous-fluid shock wave in an elastic tube

Background

The arterial pulse is a viscous-fluid shock wave that is initiated by blood ejected from the heart. This wave travels away from the heart at a speed termed the pulse wave velocity (PWV). The PWV increases during the course of a number of diseases, and this increase is often attributed to arterial stiffness. As the pulse wave approaches a point in an artery, the pressure rises as does the pressure gradient. This pressure gradient increases the rate of blood flow ahead of the wave. The rate of blood flow ahead of the wave decreases with distance because the pressure gradient also decreases with distance ahead of the wave. Consequently, the amount of blood per unit length in a segment of an artery increases ahead of the wave, and this increase stretches the wall of the artery. As a result, the tension in the wall increases, and this results in an increase in the pressure of blood in the artery.

Methods

An expression for the PWV is derived from an equation describing the flow-pressure coupling (FPC) for a pulse wave in an incompressible, viscous fluid in an elastic tube. The initial increase in force of the fluid in the tube is described by an increasing exponential function of time. The relationship between force gradient and fluid flow is approximated by an expression known to hold for a rigid tube.

Results

For large arteries, the PWV derived by this method agrees with the Korteweg-Moens equation for the PWV in a non-viscous fluid. For small arteries, the PWV is approximately proportional to the Korteweg-Moens velocity divided by the radius of the artery. The PWV in small arteries is also predicted to increase when the specific rate of increase in pressure as a function of time decreases. This rate decreases with increasing myocardial ischemia, suggesting an explanation for the observation that an increase in the PWV is a predictor of future myocardial infarction. The derivation of the equation for the PWV that has been used for more than fifty years is analyzed and shown to yield predictions that do not appear to be correct.

Conclusion

Contrary to the theory used for more than fifty years to predict the PWV, it speeds up as arteries become smaller and smaller. Furthermore, an increase in the PWV in some cases may be due to decreasing force of myocardial contraction rather than arterial stiffness.     

Scaling phloem transport: elasticity and pressure-concentration waves

Earlier theoretical analyses of the rate of propagation of pressure-concentration waves in the phloem were performed without adequate attention to the elastic expansion of sieve tube walls. Here, it is shown that the rate of propagation of pressure-concentration waves in phloem sieve tubes is not significantly impeded by wall elasticity, but rather, as previously implicated, by the ratio of sap osmotic pressure to the axial drop in sap hydrostatic pressure. It is also shown that pressure-concentration waves move equally well in both the upstream and downstream directions. These results permit future models to ignore elastic effects, and lend additional theoretical support to the "osmoregulatory flow" hypothesis, which argues that efficient molecular control of the phloem is permitted by maintaining sieve sap hydrostatic pressure at a value that is spatially nearly constant, which in turn permits changes in sieve tube state to be rapidly transmitted throughout the sieve tube via pressure-concentration waves.     

Numerical simulation of sinusoidal flow in a straight elastic tube: effects of phase angles

Numerical simulations of flow in straight elastic (moving wall) tubes subjected to a sinusoidal pressure gradient were performed for conditions prevailing in large and medium sized arteries. The effects of varying the phase angle between the pressure gradient and the tube radius, the amplitude of wall motion, and the unsteadiness parameter (alpha) on flow rate and wall shear stress were investigated. Mean and peak flow rates and shear stresses were found to be strongly affected by the phase angle between the pressure gradient and the tube radius with greater sensitivity at higher diameter variation and higher alpha. In large artery simulations (alpha = 12), means flow rate was found to be 60% higher and peak flow rate to be 73% higher than corresponding rigid tube values for certain phase angles, while a threefold increase in mean wall shear stress and sevenfold increase in peak wall shear stress were observed in a sensitive phase angle range. Significant reversal in the wall shear stress direction occurred in the sensitive phase angle range even when there was negligible flow rate reversal. All effects were greatly diminished in simulations of medium sized vessels (alpha = 4). Some experimental evidence to support the predictions of a strong effect of phase angle on wall shear stress in large vessels is presented. Finally, physiological implications of the present work are discussed from a basis of aortic input impedance data, and a physical explanation for the extreme sensitivity of the flow field to small amplitude wall motion at high alpha is given.     

Comparison of pressure losses in steady non-Newtonian flow through experimental tapered and cylindrical arterial prostheses

The use of an arterial prosthesis with a tapered lumen has several important advantages; for example, improved stability of flow, increased wall shear and better matching of its size with that of the host vessel. Tapering may, however, lead to increased energy losses, particularly if the angle of taper is large and the flow is high. This study is concerned with the determination of pressure drop for steady and laminar converging flow through rigid wall models of tapered arterial grafts. The angles of taper examined ranged from 0.5° to 1.0°. Aqueous solutions of polyacrylamide, with non-Newtonian viscous properties similar to those of blood, were used. The pressure drops across the tapered tubes were measured and the data were measured and the data were related to the pressure loss in cylindrical tubes of equivalent dimensions. Expressions for the ratio of the pressure drop in a tapered tube to that in a cylindrical tube for steady flow of a power law fluid were derived; there was good agreement between the predicted and the measured pressure drop ratios over a wide range of flows. The results of this study may be applied to the design of tapered arterial grafts. The pressure losses to be expected in tapered bypass grafts having various dimensions can easily be computed.     

Numerical simulation for the propagation of nonlinear pulsatile waves in arteries.

The behavior of nonlinear pulsatile flow of incompressible blood contained in an elastic tube is examined. The theory takes into account the nonlinear convective terms of the Navier-Stokes equations. The motion of the arterial wall is characterized by a set of linearized differential equations. The region bounded by the flexible arterial wall is mapped into a fixed area in which numerical discretization takes place. The finite element method (Galerkin weighted residual approach) is used for the solution of this nonlinear system. The results obtained are pressure distribution, velocity profile, flow rate and wall displacements along the elastic tube (20 cm long).     

Wave Propagation in a Viscous Fluid Contained in an Orthotropic Elastic Tube

The problem of pressure wave propagation through a viscous fluid contained in an orthotropic elastic tube is considered in connection with arterial blood flow. Solutions to the fluid flow and elasticity equations are obtained for the presence of a reflected wave. Numerical results are presented for both isotropic and orthotropic elastic tubes. In particular, the pressure pulse, flow rate, axial fluid velocity, and wall displacements are plotted vs. time at various stations along the ascending aorta of man. The results indicate an increase in the peak value of the pressure pulse and a decrease in the flow rate as the pulse propagates away from the heart. Finally, the velocity of wave propagation depends mainly on the tangential modulus of elasticity of the arterial wall, and anisotropy of the wall accounts in part for the reduction of longitudinal movements and an increase in the hydraulic resistance.     

Flow of a newtonian fluid through a permeable tube: The application to the proximal renal tubule

Creeping flow of a Newtonian fluid through a rigid permeable tube is considered and the transmural seepage is assumed to obey Darcy's law. Closed-form solutions for the pressure and velocity fields are presented and equations describing the axial variation of the mean cross-sectional pressure, the axial volumetric flow and the transmural fluid flux are derived. Approximate solutions for small seepage rates are given and are applied to the flow in the proximal renal tubule. Probable values for the epithelium permeability and the intraluminal hydrostatic pressure drop are obtained.     

Periodic flow of a viscous liquid in a thick-walled elastic tube

Arterial blood flow is analyzed on the basis of a realistic model consisting of a viscous liquid contained in a thick-walled viscoelastic tube. Approximate forms of the Navier-Stokes and continuity equations are derived for this model and solved in conjunction with the equations of motion of an elastic solid. Expressions are found for the displacement of the tube wall, velocity distribution, volume flow rate and phase velocity of the pressure wave. Changes in the shape of the pressure wave caused by damping and dispersion are determined, and the effect of viscoelasticity is assessed. Numerical results are presented which correspond to observed parameters of the circulatory systems of living animals.     

A semi-empirical model for flow of blood and other particulate suspensions through narrow tubes

A semi-empirical model applicable to the flow of blood and other particulate suspensions through narrow tubes has been developed. It envisages a central core of blood surrounded by a wall layer of reduced hematocrit. With the help of this model the wall layer thickness and extent of plug flow may be calculated using pressure drop, flow rate and hematocrit reduction data. It has been found from the available data in the literature that for a given sample of blood the extent of plug flow increases with decreasing tube diameter. Also for a flow through a given tube it increases with hematocrit. The wall layer thickness is found to decrease with increase in blood hematocrit. A comparison between the results of rigid particulate suspensions and blood reveals that the thicker wall layer and smaller plug flow radius in the case of blood may be attributed to the deformability of the erythrocytes.     

A mathematical model for the peristaltic flow of chyme movement in small intestine

A mathematical model based on viscoelastic fluid (fractional Oldroyd-B model) flow is considered for the peristaltic flow of chyme in small intestine, which is assumed to be in the form of an inclined cylindrical tube. The peristaltic flow of chyme is modeled more realistically by assuming that the peristaltic rush wave is a sinusoidal wave, which propagates along the tube. The governing equations are simplified by making the assumptions of long wavelength and low Reynolds number. Analytical approximate solutions of problem are obtained by using homotopy analysis method and convergence of the obtained series solution is properly checked. For the realistic values of the emerging parameters such as fractional parameters, relaxation time, retardation time, Reynolds number, Froude number and inclination of tube, the numerical results for the pressure difference and the frictional force across one wavelength are computed and discussed the roles played by these parameters during the peristaltic flow. On the basis of this study, it is found that the first fractional parameter, relaxation time and Froude number resist the movement of chyme, while, the second fractional parameter, retardation time, Reynolds number and inclination of tube favour the movement of chyme through the small intestine during pumping. It is further revealed that size of trapped bolus reduces with increasing the amplitude ratio whereas it is unaltered with other parameters.     

Wave dissipation in flexible tubes in the time domain: in vitro model of arterial waves

Earlier work of wave dissipation in flexible tubes and arteries has been carried out predominantly in the frequency domain and most of the studies used the measured pressure waveform for presenting the results. In this work we investigate the pattern of wave dissipation in the time domain using the separated forward and backward travelling waves in flexible tubes. We tested four sizes of latex tubes of 2m in length each, where a single semi-sinusoidal in shape, pressure wave, was produced at the inlet of each tube. Simultaneous measurements of pressure and flow waveforms were recorded every 5cm along the tubes and wave speed was determined using the pressure-velocity loop method (PU-loop). The measured data and wave speed were used to separate the pressure waveform and wave intensity, into their forward and backward directions, using wave intensity analysis (WIA). Also, the energy carried by the wave was calculated by integrating the relevant area under the wave intensity curve. The peak of the measured pressure waveform increased downstream, however, the peak of the separated forward pressure waveform decreased exponentially along the tube. Wave intensity and energy also dissipated exponentially along the travelling distance. The peaks of the separated pressure and wave intensity decreased in the forward in a similar exponential way to that in the backward direction in all four tube sizes. Also, the smaller the size of the tube the greater wave dissipation it caused. We conclude that wave separation is useful in studying wave dissipation in elastic tubes, and WIA provides a convenient method for determining the dissipation of the energy carried by the wave along the travelled distance. The separated pressure waveform, wave intensity and wave energy dissipate exponentially with the travelling distance, and wave dissipation varies conversely with the diameter of elastic tubes.