Next: Axisymmetric Stokes Flow In Up: Incompressible Viscous Flow Previous: Axisymmetric Stokes Flow

# Axisymmetric Stokes Flow Around a Solid Sphere

Consider a solid sphere of radius that is moving under gravity at the constant vertical velocity through a stationary fluid of density and viscosity . Here, gravitational acceleration is assumed to take the form . Provided that the typical Reynolds number,

 (10.95)

is much less than unity, the flow around the sphere is an example of axisymmetric Stokes flow. Let us transform to a frame of reference in which the sphere is stationary, and centered at the origin. Adopting the standard spherical coordinates , , , the surface of the sphere corresponds to , and the surrounding fluid occupies the region . By symmetry, the flow field outside the sphere is axisymmetric (i.e., ), and has no toroidal component (i.e., ). The physical boundary conditions at the surface of the sphere are

 (10.96) (10.97)

that is, the normal and tangential fluid velocities are both zero at the surface. A long way from the sphere, we expect the fluid velocity to asymptote to . In other words,

 (10.98) (10.99)

Let us write

 (10.100)

where is the Stokes stream function. As we saw in the previous section, axisymmetric Stokes flow is characterized by

 (10.101)

Here, the differential operator is specified in Equation (10.89). The boundary conditions (10.96)-(10.99) reduce to

 (10.102) (10.103) (10.104)

Equation (10.104) suggests that can be written in the separable form

 (10.105)

In this case,

 (10.106) (10.107)

and Equations (10.101)-(10.104) reduce to

 (10.108) (10.109) (10.110)

Let us try a test solution to Equation (10.108) of the form . We find that

 (10.111)

which implies that , , , . Hence, the most general solution to Equation (10.108) is

 (10.112)

where , , , are arbitrary constants. However, the boundary condition (10.110) yields and , whereas the boundary condition (10.109) gives and . Thus, we conclude that

 (10.113)

and the stream function becomes

 (10.114)

(See Figure 10.4.) From Equation (10.88), the fluid vorticity is

 (10.115)

(See Figure 10.5.) Moreover, from Equation (10.82),

 (10.116)

Hence,

 (10.117) (10.118)

which implies that the effective pressure distribution within the fluid is

 (10.119)

where is an arbitrary constant. (See Figure 10.6.) However, , where . Thus, the actual pressure distribution is

 (10.120)

From Section 1.20, the radial and tangential components of the force per unit area exerted on the sphere by the fluid are

 (10.121) (10.122)

Now, . Moreover, because , it follows from Equation (1.168) that . Finally, Equation (10.88) yields . Hence,

 (10.123) (10.124)

Thus, the force density at the surface of the sphere is

 (10.125)

It follows that the net vertical force exerted on the sphere by the fluid is

 (10.126)

which reduces to

 (10.127)

By symmetry, the horizontal components of the net force both average to zero. We can recognize the second term on the right-hand side of the previous equation as the buoyancy force due to the weight of the fluid displaced by the sphere. (See Chapter 2.) Moreover, the first term can be interpreted as the viscous drag acting on the sphere. Note that this drag acts in the opposite direction to the relative motion of the sphere with respect to the fluid, and its magnitude is directly proportional to the relative velocity.

Vertical force balance requires that

 (10.128)

where is the sphere's mass. In other words, in a steady state, the weight of the sphere balances the vertical force exerted by the surrounding fluid. If the sphere is composed of material of mean density then . Hence, in the frame in which the fluid a large distance from the sphere is stationary, the steady vertical velocity with which the sphere moves through the fluid is

 (10.129)

where is the fluid's kinematic viscosity. Obviously, if the sphere is more dense than the fluid (i.e., if ) then it moves downward (i.e., ), and vice versa. Finally, the typical Reynolds number of the fluid flow in the vicinity of the sphere is

 (10.130)

For the case of a grain of sand falling through water at , we have and (Batchelor 2000). Hence, , where is measured in meters. Thus, expression (10.129), which is strictly speaking only valid when , but which turns out to be approximately valid for all Reynolds numbers less than unity, only holds for sand grains whose radii are less than about microns. Such grains fall through water at approximately . For the case of a droplet of water falling through air at and atmospheric pressure, we have and (Batchelor 2000). Hence, , where is measured in meters. Thus, expression (10.129) only holds for water droplets whose radii are less than about microns. Such droplets fall through air at approximately .

At large values of , Equations (10.106), (10.107), and (10.113) yield

 (10.131) (10.132)

It follows that

 (10.133)

and

 (10.134)

Hence,

 (10.135)

where is the Reynolds number of the flow in the immediate vicinity of the sphere. [See Equation (10.95).] Our analysis is based on the assumption that advective inertia is negligible with respect to viscosity. However, as is clear from the previous expression for the ratio of inertia to viscosity within the fluid, even if this ratio is much less than unity close to the sphere--in other words, if --it inevitably becomes much greater than unity far from the sphere: that is, for . In other words, inertia always dominates viscosity, and our Stokes flow solution therefore breaks down, at sufficiently large .

Next: Axisymmetric Stokes Flow In Up: Incompressible Viscous Flow Previous: Axisymmetric Stokes Flow
Richard Fitzpatrick 2016-03-31