Masako Sugihara-Seki

The motion of an ellipsoid in tube flow at low Reynolds numbers

The motion of a rigid ellipsoidal particle freely suspended in a Poiseuille flow of an incompressible Newtonian fluid through a narrow tube is studied numerically in the zero-Reynolds-number limit. It is assumed that the effect of inertia forces on the motion of the particle and the fluid can be neglected and that no forces or torques act on the particle. The Stokes equation is solved by a finite element method for various positions and orientations of the particle to yield the instantaneous velocity of the particle as well as the flow field around it, and the particle trajectories are determined for different initial configurations. A prolate spheroid is found to either tumble or oscillate in rotation, depending on the particle tube size ratio, the axis ratio of the particle, and the initial conditions. A large oblate spheroid may approach asymptotically a steady, stable configuration, at which it is located close to the tube centreline, with its major axis slightly tilted from the undisturbed flow direction. The motion of non-axisymmetric ellipsoids is also illustrated and discussed with emphasis on the effect of the particle shape and size.

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.

The motion of an elliptical cylinder in channel flow at low Reynolds numbers

The motion of an elliptical cylindrical particle immersed in an incompressible Newtonian fluid in a narrow channel is examined numerically in the zero-Reynolds-number limit. It is assumed that no external forces or torques act on the elliptical cylinder, and the effects of inertia forces on the motion of the fluid and the particle are neglected. The Stokes equations are solved by a finite-element method for various positions and orientations of the cylinder, yielding the instantaneous velocities of the particle that satisfy the conditions of zero force and zero torque on the particle. Using the computed longitudinal, lateral, and angular velocities of the particle, the evolution of the particles position and orientation is determined for various initial configurations. An elliptical cylinder is found to either tumble or oscillate in rotation, depending on the particle-channel size ratio, the axis ratio of the elliptical cylinder, and the initial conditions. In the first case, the particle rotates continuously in one direction that is well approximated by Jefferys solution for an elliptical cylinder in unbounded shear flow with a so-called equivalent axis ratio; in the second case, the particle changes its direction of rotation during part of each period. In both cases, the particle translates with a periodically varying longitudinal velocity, accompanied by a considerable side drift due to the walls. The oscillatory motion is more likely to occur when the particle-channel size ratio or axis ratio is increased. The tumbling motion is inhibited for elliptic cylinders whose size ratios are larger than threshold values that depend on the axis ratio.

A mechanism for erythrocyte-mediated elevation of apparent viscosity by leukocytes in vivo without adhesion to the endothelium.

In spite of the relatively small number of leukocytes in the circulation, they have a significant influence on the perfusion of such organs as skeletal muscle or kidney. However, the underlying mechanisms are incompletely understood. In the current study a combined in vivo and computational approach is presented in which the interaction of individual freely flowing leukocytes with erythrocytes and its effect on apparent blood viscosity are explored. The skeletal muscle microcirculation was perfused with different cell suspensions with and without leukocytes or erythrocytes. We examined a three-dimensional numerical model of low Reynolds number flow in a capillary with a train of erythrocytes (small spheres) in off-axis positions and single larger leukocytes in axisymmetric positions. The results indicate that in order to match the slower axial velocity of leukocytes in capillaries, erythrocytes need to position themselves into an off-axis position in the capillary. In such off-axis positions at constant mean capillary velocity, erythrocyte axial velocity matches on average the axial velocity of the leukocytes, but the apparent viscosity is elevated, in agreement with the whole organ perfusion observations. Thus, leukocytes influence the whole organ resistance in skeletal muscle to a significant degree only in the presence of erythrocytes.

Numerical study of asymmetric flows of red blood cells in capillaries

Experimental studies have shown that red blood cells in capillaries may flow in single-file or multifile arrangements. To model multifile rheological behavior, the asymmetric flows of rigid circular cylinders in a two-dimensional channel are studied by numerical analysis. The rigid circular cylinders are arranged off-center in a channel in a row or two rows with equal spacings. The motion of the suspending fluid is analyzed by the finite element method applied to the Stokes equation, and the motions of the particles are simultaneously determined under the zero force and zero moment conditions appropriate to neutorally buoyant particles. The velocity difference between the particles and the bulk flow is significantly affected by the arrangement of the particles. The particle velocity is reduced as the particles are moved away from the centerline of the channel. At a constant concentration of the particles, the relative apparent viscosity of an off-center arrangement is considerably higher than that of a single-file flow of the particles located on the centerline of the channel. The present results suggest that changes of the radial distribution of red cells flowing through narrow vessels may lead to alterations of the Fahraeus and Fahraeus-Lindqvist effects.

Flow across microvessel walls through the endothelial surface glycocalyx and the interendothelial cleft

A mathematical model is presented for steady fluid flow across microvessel walls through a serial pathway consisting of the endothelial surface glycocalyx and the intercellular cleft between adjacent endothelial cells, with junction strands and their discontinuous gaps. The three-dimensional flow through the pathway from the vessel lumen to the tissue space has been computed numerically based on a Brinkman equation with appropriate values of the Darcy permeability. The predicted values of the hydraulic conductivity Lp, defined as the ratio of the flow rate per unit surface area of the vessel wall to the pressure drop across it, are close to experimental measurements for rat mesentery microvessels. If the values of the Darcy permeability for the surface glycocalyx are determined based on the regular arrangements of fibres with 6nm radius and 8nm spacing proposed recently from the detailed structural measurements, then the present study suggests that the surface glycocalyx could be much less resistant to flow compared to previous estimates by the one-dimensional flow analyses, and the intercellular cleft could be a major determinant of the hydraulic conductivity of the microvessel wall.

Asymmetric flows of spherical particles in a cylindrical tube

To study the rheological behavior of blood cells in various flow patterns through narrow vessels, we analyzed numerically the motion of blood cells arranged in one row or two rows in tube flow, at low Reynolds numbers. The particles are assumed to be identical rigid spheres placed periodically along the vessel axis at off-axis positions with equal spacings. The flow field of the suspending fluid in a circular cylindrical tube is analyzed by a finite element method applied to the Stokes equations, and the motion of each particle is simultaneously determined by a force-free and torque-free condition. In both cases of single- and two-file arrangements of the particles, their longitudinal and angular velocities are largely affected by the radial position and the axial spacing between neighboring particles. The apparent viscosity of the asymmetric flows in higher than that of the symmetric flow where particles are located on the tube centerline, and this is more pronounced when particles are placed farther from the tube centerline and when the axial distance between neighboring particles is reduced.

Transport of spheres suspended in the fluid flowing between hexagonally arranged cylinders

The motion of a spherical particle suspended in an incompressible Newtonian fluid flowing longitudinally between hexagonally arranged circular cylinders has been numerically analysed by a finite-element method in the Stokes flow regime. The results are applied to study the diffusive and convective transport of spherical solutes across the vascular endothelial surface glycocalyx, based on the quasi-periodic ultrastructural model. The obtained values of diffusive permeability and reflection coefficient of the solutes show a reasonable agreement with experimental observations, and conform to the hypothesis that the endothelial surface glycocalyx forms the primary size selective structure to solutes in microvascular permeability.

Force Acting on Spheres Adhered to a Vessel Wall

To evaluate the force and torque acting on leukocytes attached to the vessel wall, we numerically study the flow field around the leukocytes by using rigid spherical particles adhered to the wall of a circular cylindrical tube as a model of adherent leukocytes. The adherent particles are assumed to be placed regularly in the flow direction with equal spacings, in one row or two rows. The flow field of the suspending fluid is analyzed by a finite element method applied to the Stokes equations, and the drag force and torque acting on each particle, as well as the apparent viscosity, are evaluated as a function of the particle to tube diameter ratio and the particle arrangements. For two-row arrangements of adhered particles where neighboring particles are placed alternately on opposite sides of the vessel, the drag and the torque exerted on each particle are higher than those for single-row arrangements, for constant particle to tube diameter ratio and axial spacing between neighboring particles. This is enhanced for larger particles and smaller axial spacings. The apparent viscosity of the flow through vessels with adhered particles is found to be significantly higher than that without adhered particles or when the particles are freely floating through the vessels.

Effect of irregularities of vessel cross-section on vascular resistance

Irregularities of the vessel luminal geometry affect resistance to blood flow in microvessels. In the present study, the effect of the protrusion of endothelial cell nuclei into the vessel lumen on vascular resistance is numerically evaluated. It is assumed that nuclei of endothelial cells protrude regularly into the vessel lumen, which has an otherwise circular cross-section. The flow of an incompressible Newtonian fluid in these vessels is numerically analyzed by a finite element method applied to the Stokes equations, and the relationship of the flow rate and the pressure drop between upstream and downstream is used to calculate vascular resistance. It is found that the vascular resistance is constantly elevated compared to that for vessels without protrusions. The vascular resistance increases as the area of the vessel cross-section decreases, its shape is more distorted from the circular, or its variation along the vessel axis is more significant.