## NON-NEWTONIAN FLUIDS

A Newtonian Fluid is one where there is a linear relationship between

*stress*and*strain-rate*: the ratio of stress to strain-rate is the viscosity of the fluid. A*Hookean solid*is one where there is a linear relationship between*stress*and*strain*: the ratio of stress to strain is the*modulus*of the solid. Many materials have intermediate properties between those of a Newtonian fluid and a Hookean solid. If the properties are predominantly solid-like, the materials are called*non-Hookean*and the materials are described as*viscoelastic*. If they are predominantly fluid-like, they are called*non-Newtonian*and the materials are described as*elasticoviscous*.
A non-Newtonian fluid has, therefore, a simultaneously elastic and viscous nature. In fact, all fluids are non-Newtonian on an appropriate time-scale, though for many common fluids such as air and water the time-scale is extremely short. When the time-scale of a flow t

_{f}is much less than the*relaxation time*t_{r}of an elasticoviscous material, elastic effects dominate. This typically happens when there are abrupt changes in flow geometry. When on the other hand t_{f}is much greater than t_{r}elastic effects relax sufficiently for viscous effects to dominate. This typically happens when there are not abrupt changes in flow geometry. The ratio of t_{f}to t_{r}is a dimensionless number of particular significance in the study of flow of non-Newtonian fluids: depending on the circumstances, this number is called the Deborah Number or the Weissenberg Number. If the Deborah or Weissenberg number is small, elastic effects can be neglected and the non-Newtonian fluid treated as a purely viscous material, albeit with a nonconstant viscosity.
The mass and linear momentum conservation equations which govern the flow of any fluid are respectively [see Richardson (1989)]:

(1)

where ρ denotes the density of the fluid, u velocity, p pressure, τ

*deviatoric stress*, g gravitational acceleration and t time. Note that it is assumed that τ is a symmetrical tensor, thus ensuring that angular momentum is conserved. Determination of the flow of any fluid requires solution of these conservation equations, subject to appropriate boundary and initial conditions, together with*constitutive equations*for the fluid. It is in the constitutive equations that the nature of a given fluid is manifested. For a Newtonian fluid τ is given by:
(3)

where e denotes the strain-rate given by:

(4)

and η denotes the viscosity of the fluid, which is independent of e (though it may, of course, depend on temperature T and perhaps on pressure p). A non-Newtonian fluid is one for which Eq. (3) does not hold or one for which Eq. (3) holds but only with η dependent on e.

A flow in a long duct such as a pipe is susceptible to rather simpler analysis than a more general flow because it is essentially a

*shear flow*(which may loosely be regarded as one in which the velocity varies in a direction perpendicular to the direction of flow) as opposed to an*extensional flow*(which may loosely be regarded as one in which the velocity varies in the direction of flow). In a shear flow, the relevant rheological or flow property is*viscosity*η and τ is given by Eq. (3). A fluid for which (3) holds but η depends on e is called a*generalized Newtonian fluid*.
The viscosity η is a scalar; the strain-rate e is a tensor. If η depends on e, then it can be shown [see Richardson (1989) and Bird, Armstrong and Hassager (1977)] that η must be a function of the three scalar invariants I

_{e}, II_{e}and IIIe of e which are given by:
(5)

(6)

(7)

where tr( ) and det( ) denote trace and determinant, respectively. For an incompressible fluid:

(8)

A shear flow is one for which:

(9)

Thus II

_{e}is the only nonzero invariant for a shear flow of an incompressible fluid. The shear-rate γ is defined by:
(10)

Hence the constitutive equation for the stress t of an incompressible generalized Newtonian fluid is:

(11)

For a Newtonian fluid, η is independent of I

_{e}, II_{e}and III_{e}and hence of γ.
Many different semi-empirical relationships for η(γ) exist [(see Bird, Armstrong & Hassager (1977), Skelland (1967) and Wilkinson (1960)]. One of the most widely used is the

*power-law fluid*:
(12)

where η

_{c}denotes consistency or unit shear-rate viscosity and n the power law exponent. A Newtonian fluid is one for which n = 0. Many non-Newtonian fluids are*shear-thinning*or*pseudoplastic*, that is the viscosity h) decreases with increasing shear-rate γ (see**Figure 1**). The magnitude of the shear stress always increases with increasing shear-rate, however, so the product η γ always increases with increasing γ. Thus −1 < n ≤ 0; typically, n is between −0.8 and 0. If n = 0, the viscosity is, of course, independent of shear-rate; if n = −0.8, an increase in g by a factor of 10 causes a decrease in η by a factor of about 6. For a*shear-thickening*or*dilatant*fluid, in contrast, η increases with increasing γ (see**Figure 1**).
A relationship which is applicable to fluids which have a yield stress is the

*Bingham plastic*or*Bingham fluid*:
(13)

(see

**Figure 1**), where η_{p}denotes the plastic viscosity and τ_{γ}the yield stress. If |τ| ≤ τ_{γ}, that is the stress in the fluid is less than its yield stress, the fluid does not deform and so γ = 0 but η = ∞ since |τ| ≠ 0 in general. In contrast if |τ| > τ_{γ}the fluid deforms and γ is linearly related to |τ|. If |τ| >> τ_{γ}then η - η_{p}. For a Newtonian fluid, τ_{γ}= 0.
A somewhat more complicated relationship, which combines features of the power-law fluid and Bingham plastic, so that it is often called a

*yield power-law fluid*, is the*Herschel-Bulkley fluid*:
(14)

which reduces to the Bingham plastic if n = 0 and to the power-law fluid if τ

_{y}= 0 and is identified with η_{c}.
A large number of other generalized Newtonian fluids exist [see Bird, Armstrong & Hassager (1977), Skelland (1967), Tanner (1985) and Wilkinson (I960)]. These can, for example, involve the introduction of a characteristic time or a limiting viscosity at high or low shear-rates and hence generally involve further parameters. Thus the Newtonian fluid involves one parameter, η; the power-law fluid involves two, η

_{c}and n; the Bingham plastic involves two, η_{p}and τ_{y}; the Herschel-Bulkley fluid involves three, η_{p}, τ_{y}and n.
In some fluids, the viscosity h varies with time t when the shear-rate g is held constant. Two different types of behavior can be distinguished [see Skelland (1967) and Wilkinson (1960)]. For a

*thixotropic fluid*, η decreases with increasing t; for a*rheopectic fluid*, η increases with increasing t (see**Figure 2**). The reason for the variation of η with t is associated with changes in the structure of the fluid. When straining of the fluid breaks down that structure, for example, by destroying local linkages within the fluid, η decreases with increasing t and the fluid is thixotropic. A similar argument can be advanced for the decrease in η with increasing γ for a pseudoplastic fluid. In contrast, when straining builds up the structure, for example by causing local alignment in the fluid, η increases with increasing t and the fluid is rheopectic. Again, a similar argument can be advanced for the increase in η with increasing γ for a dilatant fluid. There is, therefore, an analogy between thixotropy and pseudoplasticity on the one hand and between rheopexy and dilatancy on the other. Nevertheless, variation of viscosity with time is different from, and should not be confused with, variation with shear-rate.
Constitutive equations for non-Newtonian fluids which are not purely viscous are considerably more difficult to develop and, because of their complexity, are considerably more difficult to use. The constitutive equations in general define what is called the

*extra stress*τ_{E}: the deviatoric stress τ in**Equation (2)**is related to τ_{E}by:
(15)

where I denotes the unit tensor. There are five main principles which can be used to deduce equations for τ

_{E}. In a purely mechanical theory (which thus ignores thermal effects and so on) for a purely viscous fluid, the first three of these principles are:- translational invariance: τ
_{E}depends on u and not on u; - local action: τ
_{E}at a given spatial position x depends only on the flow in the neighborhood of x, - reference-frame indifference or rotational invariance: it can be shown [see Leigh (1968)] that this means that τ
_{E}depends on e given by (4) and not on vorticity w given by:

(16)

It can then be shown [(see Richardson (1989) and Leigh (1968)] that:

(17)

where A, B and C are constants depending on I

_{e}, II_{e}and III_{e}given by (5), (6) and (7), respectively. For an incompressible fluid, the term A I can be subsumed into pressure p. Generalized Newtonian fluids are thus fluids for which B = η and C = 0. A fluid for which B = η and C ≠ 0 is called a*Reiner-Rivlin fluid*. For elasticoviscous fluids, a fourth principle is added:- determinism: τ
_{E}for a given material point X at time t depends on the deformation of X at all times t* ≤ t and not t* > t. A typical constitutive equation might then be:

(18)

where m( ) denotes a memory function and D( ) a deformation measure. An elastic material has a perfect memory in the sense that, however it is deformed, it will always try to revert to its original shape. A viscous material, in contrast, has no memory and never tries to revert to its original shape. An elasticoviscous material has a fading memory; it remembers the recent past better than the distant past. A typical functional form for m( ) is:

(19)

where t

_{r}denotes a relaxation time. Constitutive equations can be written in integral form, such as Eq. (18), or in differential form. A widely used one in differential form is the*Oldroyd fluid*:
(20)

where t

_{r1}and t_{r2}denote relaxation times and:
(21)

where a, b and c are constants. Thus is a generalization of the substantial derivative D/Dt given by:

(22)

D/Dt allows for translation following a given material point X: allows for rotation and stretching as well as translation following X. A special example of the Oldroyd fluid is the

*Maxwell fluid*for which t_{r2}= 0 and t_{r1}= t_{r}. For very small deformations, it can be shown that:
(23)

and hence that, for a Maxwell fluid:

(24)

For a nearly steady flow, that is one for which :

(25)

which is the constitutive equation for a Newtonian fluid. In contrast for a rapidly changing flow, that is one for which :

(26)

where G denotes the

*shear modulus*. Thus Eq. (26) can be integrated to yield:
(27)

But the time-integral of the strain-rate e is just the strain so Eq. (22) is the constitutive equation for a Hookean solid. Thus the Maxwell fluid can be seen to possess both viscous and elastic characteristics.

The study of the flow of non-Newtonian fluids is called

*rheology*: measurements are made in*rheometers*[see Walters (1975)]. In addition to viscosity, or strictly*shear viscosity*, η, three other parameters are of interest. They are*extensional viscosity*η_{ext},*first normal stress difference coefficient*Ψ_{1}and*second normal stress difference coefficient*Ψ_{2}. This is because, in addition to**Shear Stress**which is normally important for Newtonian fluids,**Normal Stress**is also important for non-Newtonian fluids. For a simple shear flow (see**Figure 3)**:
(28)

where i

_{x}denotes a unit vector in the x-direction. The stress τ comprises nine components. For a Newtonian fluid, seven of these (τ_{xz}, τ_{zx}, τ_{yz}, τ_{zy}, τ_{xx}, τ_{yy}and τ_{ZZ}) vanish for a simple shear flow. For a non-Newtonian fluid in contrast, only four of these (τ_{xz}, τ_{zx}, τ_{yz}and τ_{zy}) vanish. Of those which do not vanish:
(29)

(30)

(31)

For a pseudoplastic fluid such as a molten polymer, η is positive and decreases as γ increases; Ψ

_{1}is positive and decreases as γ increases and Ψ_{2}is negative and increases as γ increases (see**Figure 4**).
For a simple extensional flow with extension-rate e (see

**Figure 5**):
(32)

where i

_{z}and i_{r}denote unit vectors in the z-direction and r-direction, respectively. By definition:
(33)

For a Newtonian fluid:

(34)

For a pseudoplastic fluid, η

_{ext}is greater than 3 η and increases as ε increases: η_{ext}can be several orders of magnitude larger than η at high ε.
Non-Newtonian fluids exhibit a number of effects not shown by Newtonian fluids. A Newtonian fluid issuing at high Reynolds number from an orifice of diameter D will tend to form a jet of diameter 0.87D, whereas a non-Newtonian fluid issuing in the same way, though normally at low Reynolds number because of the higher viscosity involved, will tend to form a jet of typical diameter 3D. This so-called

*extrudate swell*has important consequences, for example, in the formation of droplets.*Rod-climbing*, as the name suggests, means that a non-Newtonian fluid will, because of normal stress effects, tend to climb up a rotating rod, whereas a Newtonian fluid will tend to move down the rod because of centrifugal effects. This has important consequences, for example, in the stirring of reactors.*Drag-reduction*is the significant reduction in friction factor, typically by two orders of magnitude, for turbulent flow of certain non-Newtonian fluids such as polymer solutions in water compared with that of the corresponding Newtonian fluid such as water. Unfortunately, the reduction in friction factor is accompanied by a similar reduction in heat transfer coefficient, so there is no obvious application of this effect to flows in heal exchangers.