Transport Phenomena by George Hirasaki - HTML preview

PLEASE NOTE: This is an HTML preview only and some elements such as links or page numbers may be incorrect.
Download the book in PDF, ePub, Kindle for a complete version.

Chapter 5Stress in Fluids

Topics Covered in this Chapter
  • Cauchy's stress principle and the conservation of momentum

  • The stress tensor

  • The symmetry of the stress tensor

  • Hydrostatic pressure

  • Principal axes of stress and the notion of isotropy

  • The Stokesian fluid

  • Constitutive equations of the Stokesian fluid

  • The Newtonian fluid

  • Interpretation of the constants λ and μ

Reading assignment
Chapter 1 in BSL
Chapter 5 in Aris

The only material property of the fluid we have so far discussed is the density. In the last chapter we introduced the rate of deformation or rate of strain tensor. The distinguishing characteristic between fluids and solids is that fluids can undergo unlimited deformation and yet maintain its integrity. The relation between the rate of deformation tensor and stress tensor is the mechanical constitutive equation of the material. An ideal fluid has a stress tensor that is independent of the rate of deformation, i.e., it has an isotropic component, which is identified as the pressure and has zero viscosity.

Cauchy's stress principle and the conservation of momentum

The forces acting on an element of a continuous medium may be of two kinds. External or body forces, such as gravitation or electromagnetic forces, can be regarded as reaching into the medium and acting throughout the volume. If the external force can be describes as the gradient of a scalar, the force is said to be conservative. Internal or contact forces are to be regarded as acting on an element of volume through its bounding surfaces. If an element of volume has an external-bounding surface, the forces there may be specified, e.g., when a constant pressure is applied over a free surface. If the element is internal, the resultant force is that exerted by the material outside the surface upon that inside. Let n be the unit outward normal at a point of the surface S and t(n) the force per unit area exerted there by the material outside S. Then Cauchy's principle asserts that t(n) is a function of the position x, the time t, and the orientation n of the surface element. Thus the total internal force exerted on the volume V through the bounding surface S is

(5.1)
_autogen-svg2png-0013.png

If f is the external force per unit mass (e.g. if 03 is vertical, gravitation will exert a force ge(3) per unit mass or ρge(3) per unit volume, the total external force will be

(5.2)
_autogen-svg2png-0017.png

The principle of the conservation of linear momentum asserts that the sum of these two forces equals the rate of change of linear momentum of the volume, i.e.,

(5.3)
_autogen-svg2png-0018.png

This is just a generalization of Newton's law of motion, which states that the rate of change of momentum of a particle is equal to the sum of forces acting on it. It has been extended to a volume that contains a number of particles.

From the form of these integral relations we can deduce an important relation. Suppose V is a volume of a given shape with characteristic dimension d. Then the volume of V will be proportional to d3 and the area of S to d2, with the proportionality constants depending only on the shape. Now let V shrink to a point but preserve its shape, then the volume integrals in the last equation will decrease as d3 but the surface integral will decrease as d2. It follows that

(5.4)
_autogen-svg2png-0028.png

or, the stresses are locally in equilibrium.

The stress tensor

To elucidate the nature of the stress system at a point P we consider a small tetrahedron with three of its faces parallel to the coordinate planes through P and the fourth with normal n (see Fig. 5.1 of Aris). If dA is the area of the slanted face, the areas of the faces perpendicular to the coordinate axis Pi is

(5.5) d Ai = nid A .

The outward normals to these faces are e(i) and we may denote the stress vector over these faces by t(i). (t(i) denotes the stress vector when +e(i) is the outward normal.) Then applying the principle of local equilibrium to the stress forces when the tetrahedron is very small we have

(5.6)
_autogen-svg2png-0038.png

Now let Tji denote the ith component of t(j) and t(n)i the ith component of t(n) so that this equation can be written

(5.7) t ( n ) i = Tjinj .

However, t(n) is a vector and n is a unit vector quite independent of the Tji so that by the quotient rule the Tji are components of a second order tensor T. In dyadic notation we might write

(5.8) t ( n ) = Tn .

This tells us that the system of stresses in a fluid is not so complicated as to demand a whole table of functions t(n)(x,n) at any given instant, but that it depends rather simply on n through the nine quantities Tji(x). Moreover, because these are components of a tensor, any equation we derive with them will be true under any rotation of the coordinate axis.

Inserting the tensor expression for the stress into the momentum balance and using the equation of continuity and Green's theorem we have

(5.9)
_autogen-svg2png-0055.png

Since all the integrals are now volume integrals, they can be combined as a single integrand.

(5.10)
_autogen-svg2png-0056.png

However, since V is an arbitrary volume this equation is satisfied only if the integrand vanishes identically.

(5.11)
_autogen-svg2png-0058.png

This is Cauchy's equation of motion and a is the acceleration. It holds for any continuum no matter how the stress tensor T is connected with the rate of strain.

The symmetry of the stress tensor

A polar fluid is one that is capable of transmitting stress couples and being subject to body torques, as in magnetic fluids. In case of a polar fluid we must introduce a body torque per unit mass in addition to the body force and a couple stress in addition to the normal stress t(n). The stress for polar fluids is discussed by Aris.

A fluid is nonpolar if the torques within it arise only as the moments of direct forces. For the nonpolar fluid we can make the assumption either that angular momentum is conserved or that the stress tensor is symmetric. We will make the first assumption and deduce the symmetry.

Return now to the integral linear momentum balance with the internal force expressed as a surface integral. If we assume that all torques arise from macroscopic forces, then not only linear momentum but also the angular momentum x×(ρv) are expressible in terms of f and t(n).

(5.12)
_autogen-svg2png-0065.png

The surface integral has as its ith component

(5.13)
_autogen-svg2png-0067.png

by Green's theorem. However, since xj,p=δjp, this last integrand is

(5.14)
_autogen-svg2png-0069.png

where T is the vector _autogen-svg2png-0071.png.

Since _autogen-svg2png-0072.png, applying the transport theorem to the angular momentum we have

(5.15)
_autogen-svg2png-0073.png

Substituting back into the equation for the angular momentum and rearranging gives

(5.16)
_autogen-svg2png-0074.png

However, the left-hand side vanishes for an arbitrary volume and so

(5.17) T× ≡ 0 .

The components of T× are _autogen-svg2png-0077.png, _autogen-svg2png-0078.png, and _autogen-svg2png-0079.png and the vanishing of these implies

(5.18) Tij = Tji

so that T is symmetric for nonpolar fluids.

Hydrostatic pressure

If the stress system is such that an element of area always experiences a stress normal to itself and this stress is independent of orientation, the stress is called hydrostatic. All fluids at rest exhibit this stress behavior. It implies that nT is always proportional to n and that the constant of proportionality is independent of n. Let us write this constant p, then

(5.19)
_autogen-svg2png-0086.png

However, this equation means that any vector is a characteristic vector or eigenvector of T. This implies that the hydrostatic stress tensor is spherical or isotropic. Thus

(5.20)
_autogen-svg2png-0088.png

for the state of hydrostatic stress.

For a compressible fluid at rest, p may be identified with the classical thermodynamic pressure. On the assumption that there is local thermodynamic equilibrium even when the fluid is in motion this concept of stress may be retained. For an incompressible fluid the thermodynamic, or more correctly thermostatic, pressure cannot be defined except as the limit of pressure in a sequence of compressible fluids. We shall see later that it has to be taken as an independent dynamical variable.

The stress tensor for a fluid may always be written

(5.21)
_autogen-svg2png-0090.png

and Pij is called the viscous stress tensor. The viscous stress tensor of a fluid vanishes under hydrostatic conditions.

If the external or body force is conservative (i.e., gradient of a scalar) the hydrostatic pressure is determined up to an arbitrary constant from the potential of the body force.

(5.22)
_autogen-svg2png-0092.png

Buoyancy (Deen, 1998)

A consequence of the fact that pressure increases with depth in a static fluid is that the pressure exerts a net upward force on any submerged object. To calculate this force, consider an object of arbitrary shape submerged in a constant-density fluid as shown in the figure. The net force on this object, Fp, is given by

(5.23)
_autogen-svg2png-0094.png
Figure 5.1

where the minus sign reflects the fact that positive pressure are compressive (i.e., pressure acts in the n direction). The pressure force is evaluated most easily by considering the situation in the figure where the solid has been replaced by an identical volume of fluid. Because the pressure in the fluid depends on depth only, this replacement of the solid does not affect the pressure distribution in the surrounding fluid. In particular, the pressure distribution on the fluid control surface in (b) must be identical to that on the solid surface in (a). Situation (b) has the advantage that we can apply the divergence theorem to the integral in the last equation, because p is a continuous function of position within the volume V; this is not necessarily true for (a), because we have said nothing about the meaning of p within a solid. Applying the divergence theorem and using the hydrostatic pressure field, we find that

(5.24)
_autogen-svg2png-0098.png

where ρf is the (constant) density of the fluid. Thus, the upward force on the object due to pressure equals the weight of an equivalent volume of liquid; this is Archimedes' law. An important implication of this equation is that Fp is independent of the absolute pressure (provided that the density is independent of pressure).

The buoyancy force is the net force on the object due to the same static pressure variation and gravity. Evaluating the gravitational force on the solid body by setting ρ=ρf, the buoyancy force is found to be

(5.25)
_autogen-svg2png-0102.png

Principal axes of stress and the notion of isotropy

The diagonal terms T11, T22, T33 of the stress tensor are sometimes called the direct stresses and the terms T12, T21, T31, T13, T23, T32 the shear stresses. When there are no external or stress couples, the stress tensor is symmetric and we can invoke the known properties of symmetric tensors. In particular, there are three principal directions and referred to coordinates parallel to these, the shear stresses vanish. The direct stresses with these coordinates are called the principal stresses and the axes the principal axes of stress.

An isotropic fluid is such that simple direct stress acting in it does not produce a shearing deformation. This is an entirely reasonable view to take for isotropy means that there is no internal sense of direction within the fluid. Another way of expressing the absence of any internally preferred direction is to say that the functional relation between stress and rate of deformation must be independent of the orientation of the coordinate system. We shall show in the next section that this implies that the principal axes of stress and rate of deformation coincide.

The Stokesian fluid

The constitutive equation of a non-elast