Peristaltic motion of a non-Newtonian fluid

Summary To understand theoretically the flow properties of physiological fluids, we have considered as a model the peristaltic motion of a power law fluid in a tube, with a sinusoidal wave of small amplitude travelling down its wall. The solution for the stream function is obtained as a power series in terms of the amplitude of the wave. The stream function and the velocity components are evaluated by solving numerically two point boundary value problems with a singular point at the origin. The influence of the applied pressure gradient along with non-Newtonian parameters on the streamlines and velocity profiles are discussed in detail.     

NUMERICAL SIMULATION OF STEADY FLOW IN A COMPLIANT TUBE OR CHANNEL WITH TAPERED WALL THICKNESS

Flow through compliant tubes with linear taper in wall thickness is numerically simulated by finite element analysis. Two models are examined: a compliant channel and an axisymmetric tube. For verification of the numerical method, flow through a compliant stenotic vessel is simulated and compared to existing experimental data. Steady two-dimensional flow in a collapsible channel with initial tension is also simulated and the results are compared with numerical solutions from the literature. Computational results for an axisymmetric tube show that as cross-sectional area falls with a reduction in downstream pressure, flow rate increases and reaches a maximum when the speed index (mean velocity divided by wave speed) is near unity at the point of minimum cross-sectional area, indicative of wave-speed flow limitation or “choking” (flow speed equals wave speed) in previous one-dimensional studies. For further reductions in downstream pressure, the flow rate decreases. Cross-sectional narrowing is significant but localized. For the particular wall and fluid properties used in these simulations, the area throat is located near the downstream end when the ratio of downstream-to-upstream wall thickness is 2; as wall taper is increased to 3, the constriction moves to the upstream end of the tube. In the planar two-dimensional channel, area reduction and flow limitation are also observed when outlet pressure is decreased. In contrast to the axisymmetric case, however, the elastic wall in the two-dimensional channel forms a smooth concave surface with the area throat located near the mid-point of the elastic wall. Though flow rate reaches a maximum and then falls, the flow does not appear to be choked.     

Peristaltic flow of a Johnson-Segalman fluid through a deformable tube

To understand theoretically the flow properties of physiological fluids we have considered as a model the peristaltic motion of a Johnson–Segalman fluid in a tube with a sinusoidal wave traveling down its wall. The perturbation solution for the stream function is obtained for large wavelength and small Weissenberg number. The expressions for the axial velocity, pressure gradient, and pressure rise per wavelength are also constructed. The general solution of the governing nonlinear partial differential equation is given using a transformation method. The numerical solution is also obtained and is compared with the perturbation solution. Numerical results are demonstrated for various values of the physical parameters of interest.      

Nonlinear and viscous effects on wave propagation in an elastic axisymmetric vessel

In this paper, a power series and Fourier series approach is used to solve the governing equations of motion in an elastic axisymmetric vessel with the assumption that the fluid is incompressible and Newtonian in a laminar flow. We obtain solutions for the wave speed and attenuation coefficient, analytically where possible, and show how these differ under a number of different conditions. Viscosity is found to reduce the wave speed from that predicted by linear wave theory and the nonlinear terms to increase the wave speed in comparison to the linear solution. For vessels with a wall stiffness in the arterial range, the reduction in the wave speed due to the viscous terms is approximately 10% and the increase due to the nonlinear terms is approximately 5%. This difference between the linear and nonlinear wave speeds was found to be largely constant irrespective of the number of terms considered in the power series for the velocity profile. The linear wave speed was found to vary weakly with stiffness, whilst the nonlinear wave speed was found to vary significantly with the stiffness, especially at low values of stiffness. The 10% variation in the wave speed due to the viscous terms was found to be constant with wall stiffness whilst the 5% variation due to the nonlinear terms was found to vary with wall stiffness. The importance of the number of terms considered in the power series is discussed showing that only a relatively small number is required in the viscous case to obtain accurate results.     

DYNAMICS OF AN IN-LINE TUBE ARRAY SUBJECTED TO STEAM–WATER CROSS-FLOW. PART III: FLUIDELASTIC INSTABILITY TESTS AND COMPARISON WITH THEORY

This paper presents the results of comprehensive flow-induced vibration tests conducted on an in-line array in steam–water two-phase flow. The responses of three essentially isolated flexible cylinders at different depths within the array were simultaneously measured. The main test parameters were, ambient pressure (and saturation temperature), in the range 0·5–5·8 MPa, void fraction, 0·70–0·96, and phase flow velocity. Tests reported here were conducted simultaneously with the damping tests reported in Part I of this study. At the highest pressures (3·0 and 5·8 MPa), strong instabilities, in homogeneous flow akin to single-phase flow occurred. The test tube located in the central region of the array was the most susceptible to instability. This was attributed partly to reduced two-phase damping deep in the array, while differences in local fluid forces at different locations in the array are not ruled out. The flow at 0·5 MPa was a nonhomogeneous intermittent slug-type flow. Strong turbulence excitation obscured clear fluidelastic instability; intermittent instability was, however, ascertained. Stability boundary calculations were done using unsteady fluid forces presented in Part II of this series of papers. Results for the case of P=5·8 MPa show good agreement with the measured instability boundary.     

Wave propagation in an elastic tube: a numerical study

In this paper, a fluid–wall interaction model, called the elastic tube model, is introduced to investigate wave propagation in an elastic tube and the effects of different parameters. The unsteady flow was assumed to be laminar, Newtonian and incompressible, and the vessel wall to be linear-elastic, isotropic and incompressible. A fluid–wall interaction scheme is constructed using a finite element method. The results demonstrate that the elastic tube plays an important role in wave propagation. It is shown that there is a time delay between the velocity waveforms at two different locations and that the peak velocity increases while the low velocity decreases in the elastic tube model, contrary to the rigid tube model where velocity waveforms overlap each other. Compared with the elastic tube model, the increase of the wall thickness makes wave propagation faster and the time delay cannot be observed clearly, however, the velocity amplitude is reduced slightly due to the decrease of the internal radius. The fluid–wall interaction model simulates wave propagation successfully and can be extended to study other mechanical properties considering complicated geometrical and material factors.     

DYNAMICS OF AN IN-LINE TUBE ARRAY SUBJECTED TO STEAM–WATER CROSS-FLOW. PART II: UNSTEADY FLUID FORCES

In this paper, experimental results of unsteady fluid-force measurements are reported. Important deviations of the measured fluid forces from their single-phase flow counterparts were uncovered. Most importantly, the resulting force coefficients are not simple functions of the reduced flow velocity U/fD, as is the case for single-phase flow. Test results at 0·5 MPa challenge the basic assumption of the existence of a time-invariant linear transfer function between tube displacement and the resulting fluid forces. Time–frequency analysis using Wigner–Ville transforms shows that the phase difference between tube displacement and the fluid force (an indicator of stabilizing or destabilizing fluid effects) undergoes significant variation under what may be considered steady flow conditions. This variation may explain the previously reported phenomenon of intermittent fluidelastic instability in two-phase flows.     

Numerical studies of slow viscous rotating flow past a sphere. III

The flow of steady incompressible viscous fluid rotating about the z-axis with angular velocity ω and moving with velocity u past a sphere of radius a which is kept fixed at the origin is investigated by means of a numerical method for small values of the Reynolds number Re_ω. The Navier–Stokes equations governing the axisymmetric flow can be written as three coupled non-linear partial differential equations for the streamfunction, vorticity and rotational velocity component. Central differences are applied to the partial differential equations for solution by the Peaceman–Rachford ADI method, and the resulting algebraic equations are solved by the ‘method of sweeps’. The results obtained by solving the non-linear partial differential equations are compared with the results obtained by linearizing the equations for very small values of Re_ω. Streamlines are plotted for Ψ = 0·05, 0·2, 0·5 for both linear and non-linear cases. The magnitude of the vorticity vector near the body, i.e. at z = 0·2, is plotted for Re_ω = 0·05, 0·24, 0·5. The correction to the Stokes drag as a result of rotation of the fluid is calculated.     

Electrohydrodynamic peristaltic flow of a viscoelastic Oldroyd fluid with a mild stenosis: Application of an endoscope

The effect of a vertical alternating current, electric field, and heat transfer on a peristaltic flow of a dielectric viscoelastic Oldroyd fluid is studied. This analysis involves uniform and nonuniform annuli having a mild stenosis. The analytical solutions of equations of motion are based on the perturbation technique. This technique depends on two parameters: amplitude ratio and small wave number. Numerical calculations are performed to obtain the effects of several parameters, such as the electrical Rayleigh number, temperature gradient, Reynolds number, wave number, maximum height of stenosis, and Weissenberg numbers, on the distributions of velocity, temperature, electric potential, and wall shear stress. It is found that the above-mentioned distributions in the case of a convergent tapered tube are larger than those in the case of a non-tapered one as well as a diverging tapered tube.     

Oscillatory flow of second grade fluid in cylindrical tube

The unsteady oscillatory flow of an incompressible second grade fluid in a cylindrical tube with large wall suction is studied analytically. Flow in the tube is due to uniform suction at the permeable walls, and the oscillations in the velocity field are due to small amplitude time harmonic pressure waves. The physical quantities of interest are the velocity field, the amplitude of oscillation, and the penetration depth of the oscillatory wave. The analytical solution of the governing boundary value problem is obtained, and the effects of second grade fluid parameters are analyzed and discussed.     

FLOW VISUALIZATION AROUND A CIRCULAR CYLINDER NEAR TO A PLANE WALL

Flow visualization, particle image velocimetry and hot-film anemometry have been employed to study the fluid flow around a circular cylinder near to a plane wall for Reynolds numbers, based on cylinder diameter, between 1200 and 4960. The effect of changing the gap between the cylinder and the wall, G, from G=0 (cylinder touching the wall) to G/D=2, was investigated. It is shown that the flow may be characterized by four distinct regions. (a) For very small gaps, G/D≤0·125, the gap flow is suppressed or extremely weak, and separation of the boundary layer occurs both upstream and downstream of the cylinder. Although there is no regular vortex shedding, there is a periodicity associated with the outer shear-layer. (b) In the “small gap ratio” region, 0·125<G/D<0·5, the flow is very similar to that for very small gaps, except that there is now a pronounced pairing between the inner shear-layer shed from the cylinder and the wall boundary layer. (c) Intermediate gap ratios, 0·5<G/D<0·75, are characterized by the onset of vortex shedding from the cylinder. (d) For the fourth region, characterized by the largest gap ratios considered, G/D>1·0, there is no separation of the wall boundary layer, either upstream or downstream of the cylinder.     

Exact and numerical simulation of peristaltic flow of a non‐Newtonian fluid with inclined magnetic field in an endoscope

In the present article, we have studied the effects of inclined magnetic field on the peristaltic flow of Jeffrey fluid through the gap between two coaxial inclined tubes. The inner tube is rigid, whereas the outer tube has sinusoidal wave traveling down its wall. The governing equations are simplified using long wave length and low Reynolds number approximations. Exact and numerical solutions have been derived for velocity profile. The expressions for pressure rise and friction force are calculated using numerical integration. Graphical results and trapping phenomenon is presented at the end of the article to see the physical behavior of different parameters. Copyright © 2010 John Wiley & Sons, Ltd.     

The effect of the slip condition on Stokes and Couette flows due to an oscillating wall: exact solutions

Stokes and Couette flows produced by an oscillatory motion of a wall are analyzed under conditions where the no-slip assumption between the wall and the fluid is no longer valid. The motion of the wall is assumed to have a generic sinusoidal behavior. The exact solutions include both steady periodic and transient velocity profiles. It is found that slip conditions between the wall and the fluid produces lower amplitudes of oscillations in the flow near the oscillating wall than when no-slip assumption is utilized. Further, the relative velocity between the fluid layer at the wall and the speed of the wall is found to overshoot at a specific oscillating slip parameter or vibrational Reynolds number at certain times. In addition, it is found that wall slip reduces the transient velocity for Stokes flow while minimum transient effects for Couette flow is achieved only for large and small values of the wall slip coefficient and the gap thickness, respectively. The time needed to reach to steady periodic Stokes flow due to sine oscillations is greater than that for cosine oscillations with both wall slip and no-slip conditions.     

Permanent Waves in Slow Free-Surface Flow of a Herschel–Bulkley fluid

Unsteady flow of a viscoplastic fluid on an inclined plane is examined. The fluid is described by the three-parameter Herschel–Bulkley constitutive equation. The set of equations governing the flow is presented, recovering earlier results for a Bingham fluid and steady uniform motion. A permanent wave solution is then derived, and the relation between wave speed and flow depth is discussed. It is shown that more types of gravity currents are possible than in a Newtonian fluid; these include some cases of flows propagating up a slope. The speed of permanent waves is derived and the possible surface profiles are illustrated as functions of the flow behavior index.     

Stokes flow through a permeable tube

Pressure-driven Stokes flow through a circular tube with a permeable wall is considered as a model of blood flow through a capillary vessel. Fluid penetrates the tube wall over a test section according to Starling law relating the normal fluid velocity to the transmural pressure defined as the difference between the wall and the uniform ambient pressure. The problem is formulated using the integral representation for Stokes flow, and the solution is computed with high accuracy using a boundary-element method for specified values of the wall permeability and percentage of fluid escaping through the walls. The results illustrate the structure of the flow and validate the predictions of a model based on the assumption of locally unidirectional flow for sufficiently small permeability.     

FINITE ELEMENT SURFACE MODEL FOR FLOW AROUND VERTICAL WALL ABUTMENTS

A two-dimensional finite element surface model is developed to determine velocities, depths, and turning angles around vertical wall abutments. The model solves the Reynolds-averaged turbulent flow equations along a horizontal plane passing through the average water surface. This approach is an improvement over the depth-averaged flow models where dispersion terms reflecting vertical effects are neglected. In the model, vertical gradient effects are accounted for through the use of power law for the vertical distribution of the longitudinal velocity; a similar treatment is applied to lateral turbulent shear stresses. The model is capable of computing the dynamic pressure distribution, which in turn is converted to water elevation values. The model, being two dimensional, is computationally efficient and practical to use. The numerical model was successfully verified using experimental data from vertical wall abutments and groins with protrusion ratios (ratio of protrusion length perpendicular to direction of flow to total channel width) of 0·1, 0·2 and 0·3. The results show the occurrence of a high intensity velocity zone close to the upstream abutment nose similar to those observed experimentally. The effects of roughness, depth, and energy slope on the intensity of flow field is investigated and an analytical expression is developed. Numerical experiments indicate that grain roughness affects flow field around the abutment nose by controlling the magnitude of the lateral velocity component and by controlling the lateral extent of the affected zone. Velocity amplification at the abutment nose is found to be mainly related to the protrusion ratio and to the friction factor, and can be up to 1·75 times the approach velocities for protrusion ratios of 0·3. For a protrusion ratio of 0·3, for a typical range of roughness values the increase in nose velocities due to friction factor alone was found to be up to 20 percent.     

Effect of the rotation,magnetic field and initial stress on peristaltic motion of micropolar fluid

In this work, the effect of magnetic field, rotation and initial stress on peristaltic motion of micropolar fluid in a circular cylindrical flexible tube with viscoelastic or elastic wall properties has been considered. Runge–Kutta technique are used. Runge–Kutta method is developed to solve the governing equations of motion resulting from a perturbation technique for small values of amplitude ratio. The time mean axial velocity profiles are presented for the case of free pumping and analyzed to observe the influence of wall properties, magnetic field, rotation and initial stress for various values of micropolar fluid parameters. In the case of viscoelastic wall, the effect of viscous damping on mean flow reversal at the boundary is seen. The numerical results of the time mean velocity profile are discussed in detail for homogeneous fluid under the effect of wall properties, magnetic field, initial stress and rotation for different cases by figures. The results indicate that the effect of wall properties, rotation, initial stress and magnetic field are very pronounced. Numerical results are given and illustrated graphically.     

非对称槽道中涡旋波的特性研究

利用PIV流场显示技术，对振荡流体在非对称槽道中涡旋波的产生、发展和消失的规律进 行了实验研究和分析，测得了涡旋波流场的速度矢量图，阐明了涡旋波流场周期性变化的特 点. 结合涡动力学方程，深入分析并揭示了做周期性运动的流体能在槽道中产生波的特性这 一规律，从中发现：流体周期变化的非定常性和不对称的槽道结构是形成涡旋波流动的主要 因素. 本文对涡旋波流场中各个旋涡的速度分布和涡量进行了测量和计算，分析了涡旋波 强化传质的机理，并研究了Re数对涡旋波流动的影响     

Experimental study on critical heat flux in bilaterally heated narrow annuli

Critical heat flux (CHF) experiments using deionized water as working fluid have been conducted in a range of pressure from 0.6 to 4.2 MPa, mass flow velocity from 60 to 130 kg/m² s and wall heat flux from 10 to 90 kW/m² for vertical narrow annuli with annular gap sizes of 0.95 and 1.5 mm. We found that the CHF, occurring only on the inside tube, or on the outside tube or on both tubes of the annular channel, depends on the heat flux ratio between surfaces of the outside and inside tubes. The CHF, occurring on the surface of the inside tube, reaches the maximum value under the pressure of 2.3 MPa while it occurring on the surface of the outside tube keeps increasing with the increase of the pressure. The CHF, occurring on the inside or outside tubes, increases with the increase of the mass flow velocity and the annular gap size; and decreases with the increase of critical quality and the other tube wall heat flux. Empirical correlations, which agree quite well with the experimental data, have been developed to predict the CHF occurring on surfaces of the inside or outside tubes of the narrow annular channel on the conditions of low pressure and low flow.     

Propagation of pressure waves in two-phase flow

We deal with a pressure wave of finite amplitude propagating in a gas and liquid medium or in the fluid in an elastic tube. We study the effects of pipe elasticity on the propagation velocity of the pressure wave. Pressure waves of finite amplitude progressing in the two-phase flow are treated considering the void fraction change due to pressure rise. The propagation velocity of the two-phase shock wave is also investigated, and the behavior of the reflection of the pressure wave at the rigid wall is analyzed and compared to that in a pure gas or liquid. The results are compared to experimental data of a pressure wave propagating in the two-phase flow in a vertical shock tube.