The oseen flow of finite clusters of spheres and the wake effect on the drag factor

In this paper, the wake effect on drag factor in the axisymmetric Oseen flow of the finite clusters of equally spaced spheres with same size is studied. Putting the Oseen lets on the centres of all the spheres, the series solution of the problem is obtained. By truncating the infinite series and applying the collocation method to solve a set of the linear algebraic equations, the approximate solution of the Oseen flow of finite clusters of spheres and the drag factor for each sphere are presented. The effect of the sphere number and spacing on the drag factor of each sphere under different Reynolds numbers are calculated and the wake effect as well as the shielding effect and the end effect are revealed. The influence of various parameters on the effects is considered and compared with the corresponding results of the Stokes flow. The convergence of the method is also studied numerically in this paper.     

球栅的Oseen流动及阻力的尾流效应

本文通过大小相同且等距放置的球栅的轴对称Oseen流动研究了阻力的尾流效应。在每个球的球心放置Oseen流子得到了问题的级数解。截断无穷级数并采用配置法解线性代数方程组求出了球栅Oseen流动的近似解及每个球所遭受的阻力。 在不同的球的个数,不同的球的间距以及不同的雷诺数下计算了各个圆球的阻力系数,发现除球栅的遮蔽效应和端缘效应外还存在着尾流效应。研究了上述参数对这些效应的影响并与Stokes流动的结果进行了比较。文章还对方法的收敛性进行了数值研究。     

THE CREEPING MOTION IN THE ENTRY REGION OF A SEMI-INFINITE CIRCULAR CYLINDRICAL TUBE

In 1969,Lew and Fung[1]considered the inlet flow into a se-mi-infinite circular cylinder at low Reynolds number.Dagan etal.[2]in1982 obtained a series solution for the creeping motionthrough a pore of finite length directly.The numerical resultsobtained in[1]also describe the entrance flow in a tube of afinite length as the Fourier integrals in the general solutions arereplaced by Fourier series.In the present paper,the Fourier in-tegralss are evaluated numerically and the velocity,pressure dis-tribution and the stream function in the entry region of a semi-infinite circular cylindrical tube is close to the factor1.3 sug-gested by Lew and Fung[1].The collocation technique applied inthe present paper is shown to converge rapidly and it should beuseful in other similar problems.     

Wall effects for motion of spheres in power-law fluids

The steady motion of spheres representing particles inside tubes filled with different fluids has been investigated using both a finite-element and a finite-volume method. The rheology of the fluids has been modelled by the power-law able to describe the shear-thinning (pseudoplastic) behaviour of a series of polymer solutions. New results have been obtained for a series of tube/sphere diameter ratios in order to investigate the wall effects on the drag exerted by the fluid on the sphere. The results agree well with previous simulations for an unbounded medium (infinite diameter ratio). Experimental investigations have also been carried out and simulated, and the results compare favourably with the experiments. The present simulations revealed the convergence of the drag coefficient to a constant value independent of tube-to-sphere diameter ratio when the power-law index approaches zero.     

Slow Motion of a Porous Sphere Translating Along the Axis of a Circular Cylindrical Pore Subject to a Stress Jump Condition

The coupled flow problem of an incompressible axisymmetrical quasisteady motion of a porous sphere translating in a viscous fluid along the axis of a circular cylindrical pore is discussed using a combined analytical–numerical technique. At the fluid–porous interface, the stress jump boundary condition for the tangential stress along with continuity of normal stress and velocity components are employed. The flow through the porous particle is governed by the Brinkman model and the flow in the outside porous region is governed by Stokes equations. A general solution for the field equations in the clear region is constructed from the superposition of the fundamental solutions in both cylindrical and spherical coordinate systems. The boundary conditions are satisfied first at the cylindrical pore wall by the Fourier transforms and then on the surface of the porous particle by a collocation method. The collocation solutions for the normalized hydrodynamic drag force exerted by the clear fluid on the porous particle is calculated with good convergence for various values of the ratio of radii of the porous sphere and pore, the stress jump coefficient, and a coefficient that is proportional to the permeability. The shape effect of the cylindrical pore on the axial translation of the porous sphere is compared with that of the particle in a spherical cavity; it found that the porous particle in a circular cylindrical pore in general attains a lower hydrodynamic drag than in a spherical envelope.     

The creeping motion of multiple spherical liquid drops

This paper deals with the drag factor of the multiple spherical liquid drops in the creeping motion by means of the Sampson singularities and collocation technique. The drag factors of the drops are calculated under distinct conditions: different number of liquid drops in the chain and different sphere spacing. From the results the influence of the viscosity ratio on the shielding effect and end effect are revealed. The convergence of the method is also studied in this paper.In this paper the collocation technique developed by Gluckman et al. in treating the rigid sphere case is applied to deal with the creeping motion of multiple spherical liquid drops which has improtant applications in bioengineering and chemical engineering. Writing the general solutions in inner and outer regions of the spheres and satisfying the kinematic and dynamic matching conditions at the collocation points on the interfaces, a set of linear algebraic equations is obtained to determine the unknown coefficients in the solutions. By means of any matrix inversion technique the approximate solutions are presented. In the first section of this paper the mathematic formulation of the problem is given and then in the second section the numerical results are introduced and analysed.     

A general formula for the drag on a sphere placed in a creeping unsteady micropolar fluid flow

In the present work, we investigate the creeping unsteady motion of an infinite micropolar fluid flow past a fixed sphere. The technique of Laplace transform is used. The drag formula is obtained in the physical domain analytically by using the complex inversion formula of the Laplace transform. The well known formula of Basset for the drag on a sphere placed in an unsteady viscous fluid flow and that of Ramkissoon and Majumdar for steady motion in the case of micropolar fluids are recovered as special cases. The obtained formula is employed to calculate the drag force for some micropolar fluid flows. Numerical results are obtained and represented graphically.     

Creeping motion of a slip spherical particle in a circular cylindrical pore

A combined analytical–numerical study for the creeping flow caused by a spherical fluid or solid particle with a slip-flow surface translating in a viscous fluid along the centerline of a circular cylindrical pore is presented. To solve the axisymmetric Stokes equations for the fluid velocity field, a general solution is constructed from the superposition of the fundamental solutions in both cylindrical and spherical coordinate systems. The boundary conditions are enforced first at the pore wall by the Fourier transforms and then on the particle surface by a collocation technique. Numerical results for the hydrodynamic drag force acting on the particle are obtained with good convergence for various values of the relative viscosity or slip coefficient of the particle, the slip parameter of the pore wall, and the ratio of radii of the particle and pore. For the motion of a fluid sphere along the axis of a cylindrical pore, our drag results are in good agreement with the available solutions in the literature. As expected, the boundary-corrected drag force for all cases is a monotonic increasing function of the ratio of particle-to-pore radii, and approaches infinity in the limit. Except for the case that the cylindrical pore is hardly slip and the value of the ratio of particle-to-pore radii is close to unity, the drag force exerted on the particle increases monotonically with an increase in its relative viscosity or with a decrease in its slip coefficient for a constant ratio of radii. In a comparison for the pore shape effect on the axial translation of a slip sphere, it is found that the particle in a circular cylindrical pore in general acquires a lower hydrodynamic drag than in a spherical cavity, but this trend can be reversed for the case of highly slippery particles and pore walls.     

Slow motion of a slip spherical particle perpendicular to two plane walls

The problem of the quasisteady motion of a spherical fluid or solid particle with a slip-flow surface in a viscous fluid perpendicular to two parallel plane walls at an arbitrary position between them is investigated theoretically in the limit of small Reynolds number. To solve the axisymmetric Stokes equation for the fluid velocity field, a general solution is constructed from the superposition of the fundamental solutions in both circular cylindrical and spherical coordinate systems. The boundary conditions are enforced first at the plane walls by the Hankel transform and then on the particle surface by a collocation technique. Numerical results for the hydrodynamic drag force exerted on the particle are obtained with good convergence for various values of the relative viscosity or slip coefficient of the particle and of the relative separation distances between the particle and the confining walls. For the motions of a spherical particle normal to a single plane wall and of a no-slip sphere perpendicular to two plane walls, our drag results are in good agreement with the available solutions in the literature for all relative particle-to-wall spacings. The boundary-corrected drag force acting on the particle in general increases with an increase in its relative viscosity or with a decrease in its slip coefficient for a given geometry, but there are exceptions. For a specified wall-to-wall spacing, the drag force is minimal when the particle is situated midway between the two plane walls and increases monotonically when it approaches either of the walls. The boundary effect on the particle motion normal to two plane walls is found to be significant and much stronger than that parallel to them.     

Acceleration of a sphere behind planar shock waves

The paper presents the results of an investigation on the motion of a spherical particle in a shock tube flow. A shock tube facility was used for studying the acceleration of a sphere by an incident shock wave. Using different optical methods and performing experiments in two different shock tubes, the trajectory and velocity of a spherical particle were measured. Based upon these results and simple one-dimensional calculations, the drag coefficient of a sphere and shading effect of sphere interaction with a shock tube flow were studied.     

Drag of a rotating sphere

We consider the problem of steady incompressible viscous fluid flow about a rotating sphere, with the flow specified on a sphere of finite radius, which reduces to the solution of the complete Navier-Stokes equations.The dimensionless stream functions and circulai velocity are sought in the form of series in powers of the Reynolds numbers, which converge for small values of this number. Recurrence formulas are derived for determining the coefficients of these series. The pressure, rotational resistance torque, and drag are determined. It is established that the rotating sphere has higher drag than a stationary sphere. The leading term of the series in powers of the Reynolds number for the drag and resistive torque is calculated.     

Viscous flow of a suspension of liquid drops in a cylindrical tube

Viscous flow in a circular cylindrical tube containing an infinite line of viscous liquid drops equally spaced along the tube axis is considered under the assumption that a surface tension, sufficiently large, holds the drops in a nearly spherical shape. Three cases are considered: (1) axial translation of the drops, (2) flow of the external fluid past a line of stationary drops, and (3) flow of external fluid and liquid drops under an imposed pressure gradient. Both fluids are taken to be Newtonian and incompressible, and the linearized equations of creeping flow are used.The results show that both drag and pressure drop per sphere increase as the spacing increases at fixed radius and also increase as the radius of the drop increases. The presence of the internal motion reduces the drag and pressure gradients in all cases compared to rigid spheres, particularly for drops approaching the size of the tube.     

On the slow translation of a solid submerged in a fluid with a surfactant surface film—II

In Shail & Gooden (1982) the problem of a solid particle translating in a semi-infinite fluid, whose surface is contaminated with a surfactant film, was examined in the quasi-steady Stokes flow régime. Various linearised models governing the variation of film concentration were considered, but the analysis was approximate in that the fluid motion generated was represented by that due to a Stokeslet situated at the centre of the particle. In this paper we remove the latter restriction and treat two specific solids, namely a rigid flat circular disk and a sphere, which move axisymmetrically perpendicular to the fluid surface. This surface is assumed to remain plane throughout the motion. The velocity field in the translating-disk problem is represented in terms of harmonic functions, and the resulting mixed boundary-value problems are reduced, for each of the film behaviours examined, to the solution of sets of simultaneous Fredholm integral equations of the second kind. These equations are solved both iteratively and numerically, and the drag on the disk is computed. For the sphere a stream-function formulation in bispherical coordinates is used. Application of the boundary conditions at the sphere and film results in infinite sets of simultaneous linear equations for the coefficients in the eigenfunction expansion of the stream function. These equations are solved by the method of truncation, and the drag on the sphere is determined.     

Slow motions of a circular cylinder experiencing slip near a plane wall

A theoretical study is presented for the two-dimensional creeping flow caused by a long circular cylindrical particle translating and rotating in a viscous fluid near a large plane wall parallel to its axis. The fluid is allowed to slip at the surface of the particle. The Stokes equations for the fluid velocity field are solved in the quasi-steady limit using cylindrical bipolar coordinates. Semi-analytical solutions for the drag force and torque acting on the particle by the fluid are obtained for various values of the slip coefficient associated with the particle surface and of the relative separation distance between the particle and the wall. The results indicate that the translation and rotation of the confined cylinder are not coupled with each other. For the motion of a no-slip cylinder near a plane wall, our hydrodynamic drag force and torque results reduce to the closed-form solutions available in the literature. The boundary-corrected drag force and torque acting on the particle decrease with an increase in the slip coefficient for an otherwise specified condition. The plane wall exerts the greatest drag on the particle when its migration occurs normal to it, and the least in the case of motion parallel to it. The enhancement in the hydrodynamic drag force and torque on a translating and rotating particle caused by a nearby plane wall is much more significant for a cylinder than for a sphere.     

A boundary element simulation of the unsteady motion of a sphere in a cylindrical tube containing a viscoelastic fluid

A boundary element method is used to simulate the unsteady motion of a sphere falling under gravity along the centreline of a cylindrical tube containing a viscoelastic fluid. The fluid is modelled by the upper-convected Maxwell constitutive equation. Results show that the viscoelasticity of the liquid leads to a damped oscillation in sphere velocity about its terminal value. The maximum sphere velocity, which occurs in the first overshoot, is approximately proportional to the square root of the Weissenberg number when the ratio of the sphere radius to the tube radius is sufficiently small. Particular attention is also paid to the wall effects. It is shown that a closer wall reduces the oscillatory amplitude of the sphere velocity but increases its frequency. The results suggest that the falling-ball technique, which is now widely used for viscosity measurement, might also be used for the determination of a relaxation time for a viscoelastic fluid.     

Stokes flow of a rotating sphere around the axis of a circular orifice or a circular disk

This paper presents an infinite series solution to the creeping flow equations for the axisymmetric motion of a sphere of arbitrary size rotating in a quiescent fluid around the axis of a circular orifice or a circular disk whose diameters are either larger or smaller than that of the sphere. Numerical tests of the convergence are passed and the comparison with the exact solution and other computational results shows an agreement to five significant figures for the torque coefficients in both cases. The torque coefficients are obtained for the sphere located up to a position tangent to the wall plane containing either the orifice or the disk. It is concluded that the torque coefficients of the sphere and the disk are monotonically increasing with the decrease of the distance from the disk or the orifice plane in both cases.     

Stokes flow past finite coaxial clusters of spheres in a circular cylinder

Solutions are presented for the Stokes flow past finite axial assemblages of up to 9 spheres in an infinitely long cylindrical tube for a wide range of sphere spacings and sphere to cylinder diameter ratios. General solutions are constructed from the fundamental solutions to the governing equation in both the cylindrical and spherical coordinate systems. No-slip boundary conditions are enforced on the tube surface by constructing the Fourier transform of the general disturbance created by the spheres, as detected on the cylinder wall. The boundary conditions are then applied on the sphere surfaces by a previously developed series truncation technique.The calculated drag forces and zero-drag velocities demonstrate the interparticle interaction effects, the sphere-wall interactions, and the effects of wall damping on the inter-particle shielding phenomenon.     

The settling of a sphere along the axis of a long square duct at low Reynolds' number

A theoretical treatment is presented for the determination of the drag upon a sphere settling along the axis of a long square duct under the condition that the creeping motion equations are applicable. In order to obtain the second reflection velocity field, it was necessary to develop a new general solution in cartesian coordinates to the creeping motion equation, applicable within the domain of a long square duct. Using the second reflection velocity field solution an Faxen's law, a third reflection (first correction to the Stokes' value) drag correction is obtained. The results show that the drag correction for a square container is quite close to (but smaller than) the drag correction produced by a cylinder whose diameter is the same as the duct width.     

Axisymmetric creeping motion of a slip spherical particle in a nonconcentric spherical cavity

A semianalytical study of the creeping flow caused by a spherical fluid or solid particle with a slip surface translating in a viscous fluid within a spherical cavity along the line connecting their centers is presented in the quasisteady limit of small Reynolds number. In order to solve the Stokes equations for the flow field, a general solution is constructed from the superposition of the fundamental solutions in the two spherical coordinate systems based on both the particle and cavity. The boundary conditions on the particle surface and cavity wall are satisfied by a collocation technique. Numerical results for the hydrodynamic drag force exerted on the particle are obtained with good convergence for various values of the ratio of particle-to-cavity radii, the relative distance between the centers of the particle and cavity, the relative viscosity or slip coefficient of the particle, and the slip coefficient of the cavity wall. In the limits of the motions of a spherical particle in a concentric cavity and near a cavity wall with a small curvature, our drag results are in good agreement with the available solutions in the literature. As expected, the boundary-corrected drag force exerted on the particle for all cases is a monotonic increasing function of the ratio of particle-to-cavity radii, and becomes infinite in the touching limit. For a specified ratio of particle-to-cavity radii, the drag force is minimal when the particle is situated at the cavity center and increases monotonically with its relative distance from the cavity center to infinity in the limit as it is located extremely away from the cavity center. The drag force acting on the particle, in general, increases with an increase in its relative viscosity or with a decrease in its slip coefficient for a given configuration, but surprisingly, there are exceptions when the ratio of particle-to-cavity radii is large.     

Accelerating motion of a sphere in a maxwell fluid

The translatory accelerating motion of a sphere due to an arbitrarily applied force in an unlimited Maxwell fluid is considered. The exact solutions for the velocity of the sphere for three particular types of accelerating motion are presented. The first is for a falling sphere; the second is for the decelerating motion of a sphere after the force which maintains the sphere with a constant velocity is removed; the third is for the motion of the sphere subjected to an impulsive force. The exact solutions are expressed in terms of real, regular, definite integrals which can be evaluated by numerical technique. Also presented are the asymptotic solutions for the velocity of the sphere in all three cases which are valid for small values of time.