Interaction of peristaltic motion with Poiseuille flow

The effect of Poiseuille flow on peristaltic transport has been investigated in a two-dimensional mathematical model of peristalsis for the case when the wall of the channel executes a sinusoidal motion of small amplitude. Closed-form solutions have been obtained for limiting values of Reynolds number and the Poiseuille flow parameter, while the method of Frobenius series solution has been used for the general case. It is found that the mean flow reversal is strongly dependent on the Poiseuille flow. The position of flow reversal may change drastically from the center of the channel to the boundaries. Numberical results are reported for various values of the physical parameters of interest.     

A model of a radially expanding and contracting lymphangion

The lymphatic system is an extensive vascular network featuring valves and contractile walls that pump interstitial fluid and plasma proteins back to the main circulation. Immune function also relies on the lymphatic system's ability to transport white blood cells. Failure to drain and pump this excess fluid results in edema characterized by fluid retention and swelling of limbs. It is, therefore, important to understand the mechanisms of fluid transport and pumping of lymphatic vessels. Unfortunately, there are very few studies in this area, most of which assume Poiseuille flow conditions. In vivo observations reveal that these vessels contract strongly, with diameter changes of the order of magnitude of the diameter itself over a cycle that lasts typically 2-3s. The radial velocity of the contracting vessel is on the order of the axial fluid velocity, suggesting that modeling flow in these vessels with a Poiseuille model is inappropriate. In this paper, we describe a model of a radially expanding and contracting lymphatic vessel and investigate the validity of assuming Poiseuille flow to estimate wall shear stress, which is presumably important for lymphatic endothelial cell mechanotransduction. Three different wall motions, periodic sinusoidal, skewed sinusoidal and physiologic wall motions, were investigated with steady and unsteady parabolic inlet velocities. Despite high radial velocities resulting from the wall motion, wall shear stress values were within 4% of quasi-static Poiseuille values. Therefore, Poiseuille flow is valid for the estimation of wall shear stress for the majority of the lymphangion contractile cycle.     

Pressure flow patterns in a cylinder with reabsorbing walls

The problem of fluid motion in renal tubules, in contrast to ordinary flow through cylinders with impermeable walls, is complicated by the existence of radial velocities generated by reabsorption processes. As a first approach to this problem, the Navier Stokes equations for axially symmetric, slow flow in an infinite cylinder whose walls reabsorb fluid are integrated. If the rate of reabsorption is constant, the solutions resemble the conventional Poiseuille flow, i.e., the longitudinal velocity profile is parabolic. In addition 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. If the rate of reabsorption is a linear function of the distance from the origin, the presence of an additive term alters these relations. If, for example, the gradient in reabsorption is positive, the axial velocity profile tends to flatten and when the gradient is sufficiently large, the maximum velocity moves from the center of the stream toward the periphery, leaving a relative minimum at the center. In passing from the center of the tube to the walls, the radial velocity passes through a miximum, regardless of the reabsorption properties of the wall.     

Effects of elastic property of the wall on flow characteristics through arterial stenoses

Hemodynamic characteristics of blood flow through arterial stenoses are numerically investigated in this work. The blood is assumed as a Newtonian fluid and the pulsatile nature of flow is modeled by using measured values of the flowrate and pressure for the canine femoral artery. An isotropic elastic and incompressible material is assumed for the wall at each axial section, but a non-uniform distribution of the shear modulus in axial direction is used to model the high stiffness of the wall at the stenosis location. Full Navier equations for a thick wall are used as the governing equations for the wall displacements. A continuous grid extending over the flow field and the wall is considered and governing equations are transformed for use in the computational domain. Discretized forms of the transformed wall and flow equations, which are coupled through the boundary conditions at their interface, are obtained by control volume method and simultaneously solved using the well-known SIMPLER algorithm. To study the effects of wall deformability, solutions are obtained for both rigid and elastic walls. The results indicate that deformability of the wall causes an increase in the time average of pressure drop, but a decrease in the maximum wall shear stress. Displacement and stress distributions in the wall are presented.     

Wave Propagation through a Newtonian Fluid Contained within a Thick-Walled,Viscoelastic Tube

The propagation of harmonic pressure waves through a Newtonian fluid contained within a thick-walled, viscoelastic tube is considered as a model of arterial blood flow. The fluid is assumed to be homogeneous and Newtonian, and its motion to be laminar and axisymmetric. The wall is assumed to be isotropic, incompressible, linear, and viscoelastic. It is also assumed that the motion is such that the convective acceleration is negligible. The motion of the fluid is described by the linearized form of the Navier-Stokes equations and the motion of the wall by classical elasticity theory. The frequency dependence of the wall mechanical properties are represented by a three parameter, relaxation-type model. Using boundary conditions describing the continuity of stress and velocity components in the fluid and the wall, explicit solutions for the system of equations of the model have been obtained. The longitudinal fluid impedance has been expressed in terms of frequency and the system parameters. The frequency equation has been solved and the propagation constant also expressed in terms of frequency and system parameters. The results indicate that the fluid impedance is smaller than predicted by the rigid tube model or by Womersley''s constrained elastic tube model. Also, the velocity of propagation is generally slower and the transmission per wavelength less than predicted by Womersley''s elastic tube model. The propagation constant is very sensitive to changes in the degree of wall viscoelasticity.     

Simulation of two-dimensional fully developed laminar flow for a magneto-hydrodynamic (MHD) pump

MHD micro-pumps circumvent the wear and fatigue caused by high pressure-drop across the check valves of mechanical micro-pumps in micro-fluidic systems. Early analyses of the fluid flow for MHD micro-pumps were mostly made possible by the Poiseuille flow theory; however, this conventional laminar approach cannot illustrate the effects of various channel sizes and shapes. This paper, therefore, presents a simplified MHD flow model based upon steady state, incompressible and fully developed laminar flow theory to investigate the characteristics of a MHD pump. Inside the pump, flowing along the channel is the electrically conducting fluid flowing driven by the Lorentz forces in the direction perpendicular to both dc magnetic field and applied electric currents. The Lorentz forces were converted into a hydrostatic pressure gradient in the momentum equations of the MHD channel flow model. The numerical simulations conducted with the explicit finite difference method show that the channel dimensions and the induced Lorentz forces have significant influences on the flow velocity profile. Furthermore, the simulation results agree well with the experimental results published by other researchers.     

Subcritical flutter in collapsible tube flow: a model of expiratory flow in the trachea

Using an axisymmetric geometry that retains certain qualitative features of the trachea, we extend one-dimensional modeling of flow in collapsible tubes to include both curved shell effects and, for untethered tubes, wall inertia. A systematic scaling of the finite deformation membrane equations leads to an approximate set which is consistent with the one-dimensional fluid model; axial and normal wall variables are coupled elastically, but only axial inertia is retained. Transverse curvature causes elastic coupling that can give rise to axial wall motion and a flutter instability. The source of instability is the product of a nonzero reference axial curvature with axial tension variation due to axial stretching. The numerical results suggest that this mechanism may be significant even in processes which cannot be assumed one-dimensional.     

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.     

Numerical simulation of peristaltic flow of a biorheological fluid with shear-dependent viscosity in a curved channel

Peristaltic motion of a non-Newtonian Carreau fluid is analyzed in a curved channel under the long wavelength and low Reynolds number assumptions, as a simulation of digestive transport. The flow regime is shown to be governed by a dimensionless fourth-order, nonlinear, ordinary differential equation subject to no-slip wall boundary conditions. A well-tested finite difference method based on an iterative scheme is employed for the solution of the boundary value problem. The important phenomena of pumping and trapping associated with the peristaltic motion are investigated for various values of rheological parameters of Carreau fluid and curvature of the channel. An increase in Weissenberg number is found to generate a small eddy in the vicinity of the lower wall of the channel, which is enhanced with further increase in Weissenberg number. For shear-thinning bio-fluids (power-law rheological index, n < 1) greater Weissenberg number displaces the maximum velocity toward the upper wall. For shear-thickening bio-fluids, the velocity amplitude is enhanced markedly with increasing Weissenberg number.     

Effects of velocity jumps on blood flow in large arteries

When a human being experiences a sudden velocity change, the blood flow is disturbed. A theoretical analysis to predict the effects of sudden velocity changes on blood flow in large arteries is presented. The situations is modelled as a one-dimensional flow problem in a viscoelastic tube where the fluid viscosity convective term in the equation of motion and nonlinearity in the elastic modulus of the tube wall are neglected. The governing equations of the model are solved by Laplace transformation. The computed results show that relatively high blood pressures, capable of harming circulation, are produced even by relatively moderate velocity jumps.     

The fluid shear stress distribution on the membrane of leukocytes in the microcirculation

Recent in-vivo and in-vitro evidence indicates that fluid shear stress on the membrane of leukocytes has a powerful control over several aspects of their cell function. This evidence raises a question about the magnitude of the fluid shear stress on leukocytes in the circulation. The flow of plasma on the surface of a leukocyte at a very low Reynolds number is governed by the Stokes equation for the motion of a Newtonian fluid. We numerically estimated the distribution of fluid shear stress on a leukocyte membrane in a microvessel for the cases when the leukocyte is freely suspended, as well as rolling along or attached to a microvessel wall. The results indicate that the fluid shear stress distribution on the leukocyte membrane is nonuniform with a sharp increase when the leukocyte makes membrane attachment to the microvessel wall. In a microvessel (10 microns diameter), the fluid shear stress on the membrane of a freely suspended leukocyte (8 microns diameter) is estimated to be several times larger than the wall shear stress exerted by the undisturbed Poiseuille flow, and increases on an adherent leukocyte up to ten times. High temporal stress gradients are present in freely suspended leukocytes in shear flow due to cell rotation, which are proportional to the local shear rate. In comparison, the temporal stress gradients are reduced on the membrane of leukocytes that are rolling or firmly adhered to the endothelium. High temporal gradients of shear stress are also present on the endothelial wall. At a plasma viscosity of 1 cPoise, the peak shear stresses for suspended and adherent leukocytes are of the order of 10 dyn/cm2 and 100 dyn/cm2, respectively.     

A stochastic model of cell aggregation under planar flow in the dilute regime

Models of the adhesion of a population of cells in a plane flow are developed, considering the dilute regime. Cells considered as rigid punctual entities are virtually injected at regular times within a plane channel limited by two fixed planes. The pressure profile is supposed to be triangular (constant gradient), in accordance with the assumptions of a Poiseuille flow. The cell adherence to the channel wall is governed by the balance of forces, accounting for gravity, non-specific physical interactions, such as electrostatic effects (repulsive) and Van der Waals forces (attractive), specific adhesive forces representing the ligand–receptor interactions, and friction between cells and the fluid in the vicinity of the endothelium wall. The spatial distribution of the adhesion molecules along the wall is supposed to be a random event, accounted for by a stochastic spatial variability of the dipolar moments of those molecules, according to a Gaussian process. Experimental trends reported for the rate of aggregation of L-selectin mediated leukocytes under shear flow are in qualitative accordance with the evolution versus time of adhering cells obtained by the present simulations. The effect of the maximal injection pressure on those kinetics is assessed.     

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

Effects of vessel compliance on flow pattern in porcine epicardial right coronary arterial tree

The compliance of the vessel wall affects hemodynamic parameters which may alter the permeability of the vessel wall. Based on experimental measurements, the present study established a finite element (FE) model in the proximal elastic vessel segments of epicardial right coronary arterial (RCA) tree obtained from computed tomography. The motion of elastic vessel wall was measured by an impedance catheter and the inlet boundary condition was measured by an ultrasound flow probe. The Galerkin FE method was used to solve the Navier–Stokes and Continuity equations, where the convective term in the Navier–Stokes equation was changed in the arbitrary Lagrangian–Eulerian (ALE) framework to incorporate the motion due to vessel compliance. Various hemodynamic parameters (e.g., wall shear stress—WSS, WSS spatial gradient—WSSG, oscillatory shear index—OSI) were analyzed in the model. The motion due to vessel compliance affects the time-averaged WSSG more strongly than WSS at bifurcations. The decrease of WSSG at flow divider in elastic bifurcations, as compared to rigid bifurcations, implies that the vessel compliance decreases the permeability of vessel wall and may be atheroprotective. The model can be used to predict coronary flow pattern in subject-specific anatomy as determined by noninvasive imaging.     

Intracellular microrheology of motile Amoeba proteus

The motility of Amoeba proteus was examined using the technique of passive particle tracking microrheology, with the aid of newly developed particle tracking software, a fast digital camera, and an optical microscope. We tracked large numbers of endogeneous particles in the amoebae, which displayed subdiffusive motion at short timescales, corresponding to thermal motion in a viscoelastic medium, and superdiffusive motion at long timescales due to the convection of the cytoplasm. Subdiffusive motion was characterized by a rheological scaling exponent of 3/4 in the cortex, indicative of the semiflexible dynamics of the actin fibers. We observed shear-thinning in the flowing endoplasm, where exponents increased with increasing flow rate; i.e., the endoplasm became more fluid-like. The rheology of the cortex is found to be isotropic, reflecting an isotropic actin gel. A clear difference was seen between cortical and endoplasmic layers in terms of both viscoelasticity and flow velocity, where the profile of the latter is close to a Poiseuille flow for a Newtonian fluid.     

Elastic expulsion from a long tube

The paper presents a theoretical analysis of elastic expulsion from a long pressurized tube following instantaneous severance. The governing equations of motion and continuity are used, together with an equation representing radial equilibrium of the tube wall, to construct a differential equation for the inner radius of the tube as a function of axial position and time. Similarity solutions for this equation are presented for a Newtonian fluid. The effect of the relaxed/fully-extended radius ratio on wall dilatation and expulsion rate is discussed.     

Flutter in collapsible tubes: a theoretical model of wheezes

A mathematical analysis of flow through a flexible channel is examined as a model of flow-induced flutter oscillations that pertain to the production of wheezing breath sounds. The model provides predictions for the critical fluid speed that will initiate flutter waves of the wall, as well as their frequency and wavelength. The mathematical results are separated into linear theory (small oscillations) and nonlinear theory (larger oscillations). Linear theory determines the onset of the flutter, whereas nonlinear theory determines the relationships between the fluid speed and both the wave amplitudes and frequencies. The linear theory predictions correlate well with data taken at the onset of flutter and flow limitation during experiments of airflow in thick-walled collapsible tubes. The nonlinear theory predictions correlate well with data taken as these flows are forced to higher velocities while keeping the flow rate constant. Particular ranges of the parameters are selected to investigate and discuss the applications to airway flows. According to this theory, the mechanism of generation of wheezes is based in the interactions of fluid forces and friction and wall elastic-restoring forces and damping. In particular, a phase delay between the fluid pressure and wall motion is necessary. The wave speed theory of flow limitation is discussed with respect to the specific data and the flutter model.     

Theoretical analysis of pressure pulse propagation in arterial vessels.

An original mathematical model of viscous fluid motion in a tapered and distensible tube is presented. The model equations are deduced by assuming a two-dimensional flow and taking into account the nonlinear terms in the fluid motion equations, as well as the nonlinear deformation of the tube wall. One distinctive feature of the model is the formal integration with respect to the radial coordinate of the Navier-Stokes equations by power series expansion. The consequent computational frame allows an easy, accurate evaluation of the effects produced by changing the values of all physical and geometrical tube parameters. The model is employed to study the propagation along an arterial vessel of a pressure pulse produced by a single flow pulse applied at the proximal vessel extremity. In particular, the effects of the natural taper angle of the arterial wall on pulse propagation are investigated. The simulation results show that tapering considerably influences wave attenuation but not wave velocity. The substantially different behavior of pulse propagation, depending upon whether it travels towards the distal extremity or in the opposite direction, is observed: natural tapering causes a continuous increase in the pulse amplitude as it moves towards the distal extremity; on the contrary, the reflected pulse, running in the opposite direction, is greatly damped. For a vessel with physical and geometrical properties similar to those of a canine femoral artery and 0.1 degree taper angle, the forward amplification is about 0.9 m-1 and the backward attenuation is 1.4 m-1, so that the overall tapering effect gives a remarkably damped pressure response. For a natural taper angle of 0.14 degrees the perturbation is almost extinct when the pulse wave returns to the proximal extremity.     

Syringomyelia hydrodynamics: an in vitro study based on in vivo measurements

A simplified in vitro model of the spinal canal, based on in vivo magnetic resonance imaging, was used to examine the hydrodynamics of the human spinal cord and subarachnoid space with syringomyelia. In vivo magnetic resonance imaging (MRI) measurements of subarachnoid (SAS) geometry and cerebrospinal fluid velocity were acquired in a patient with syringomyelia and used to aid in the in vitro model design and experiment. The in vitro model contained a fluid-filled coaxial elastic tube to represent a syrinx. A computer controlled pulsatile pump was used to subject the in vitro model to a CSF flow waveform representative of that measured in vivo. Fluid velocity was measured at three axial locations within the in vitro model using the same MRI scanner as the patient study. Pressure and syrinx wall motion measurements were conducted external to the MR scanner using the same model and flow input. Transducers measured unsteady pressure both in the SAS and intra-syrinx at four axial locations in the model A laser Doppler vibrometer recorded the syrinx wall motion at 18 axial locations and three polar positions. Results indicated that the peak-to-peak amplitude of the SAS flow waveform in vivo was approximately tenfold that of the syrinx and in phase (SAS approximately 5.2 +/- 0.6 ml/s, syrinx approximately 0.5 +/- 0.3 ml/s). The in vitro flow waveform approximated the in vivo peak-to-peak magnitude (SAS approximately 4.6 +/- 0.2 ml/s, syrinx approximately 0.4 +/- 0.3 ml/s). Peak-to-peak in vitro pressure variation in both the SAS and syrinx was approximately 6 mm Hg. Syrinx pressure waveform lead the SAS pressure waveform by approximately 40 ms. Syrinx pressure was found to be less than the SAS for approximately 200 ms during the 860-ms flow cycle. Unsteady pulse wave velocity in the syrinx was computed to be a maximum of approximately 25 m/s. LDV measurements indicated that spinal cord wall motion was nonaxisymmetric with a maximum displacement of approximately 140 microm, which is below the resolution limit of MRI. Agreement between in vivo and in vitro MR measurements demonstrates that the hydrodynamics in the fluid filled coaxial elastic tube system are similar to those present in a single patient with syringomyelia. The presented in vitro study of spinal cord wall motion, and complex unsteady pressure and flow environment within the syrinx and SAS, provides insight into the complex biomechanical forces present in syringomyelia.     

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. This research was partially supported by the National Science Foundation; it was done in part by D. K. Whirlow in partial fulfillment of the requirements for the degree of Doctor of Philosophy at the Carnegie Institute of Technology.