Non Newtonian
Non Newtonian
Non Newtonian
Non-Newtonian fluid
3-1. Introduction:
The study of the deformation of flowing fluids is called rheology; the
rheological behavior of various fluids is sketched in Figure 3-1. Newtonian
fluids, defined as fluids for which the shear stress is linearly proportional
to the shear strain rate. Newtonian fluids (stress proportional to strain
rate). Many common fluids, such as air and other gases, water, kerosene,
gasoline, and other oil-based liquids, are Newtonian fluids. Fluids for which
the shear stress is not linearly related to the shear strain rate are called
non-Newtonian fluids. Examples include slurries and colloidal
suspensions, polymer solutions, blood, paste, and cake batter. Some non-
Newtonian fluids exhibit a "memory"-the shear stress depends not only on
local strain rate, but also on its history. A fluid that returns (either fully or
partially) to its original shape after the applied stress is released is called
viscoelastic.
n=1
n<1
𝝉𝑩 n=1
n>1
Figure 3-4.
𝝆𝒖𝒅𝒊
𝑹𝒆 = (3.1)
𝝁
Another definition is based, not on the true shear rate at the wall, but on
the flow characteristic. This quantity, which may be called the apparent
viscosity for pipe flow, is given by:
For laminar flow, 𝝁𝒂𝒑 has the property that it is the viscosity of a Newtonian
fluid having the same flow characteristic as the non-Newtonian fluid when
subjected to the same value volumetric flow rate for the same pressure
gradient in the same pipe. This suggests that 𝝁𝒂𝒑 might be a useful quantity
for correlating flow rate-pressure gradient data for non-Newtonian flow in
pipes. This is found to be the case and it is on 𝝁𝒂𝒑 that a generalized
Reynolds number 𝑹𝒆 is based:
𝝆𝒖𝒅𝒊
𝑹𝒆′ = (3.4)
𝝁𝒂𝒑
The pipe flow apparent viscosity, defined by equation 3.3, is given by:
When this equation for 𝝁𝒂𝒑 is substituted into equation 3.4, the generalized
Reynolds number takes the from:
′ ′
′
𝝆𝒖𝟐−𝒏 𝒅𝒏 𝒊
𝑹𝒆 = ′ −𝟏 (3.7)
𝟖𝒏 𝑲′
Use of this generalized Reynolds number was suggested by Metzner and
Reed (1955). For Newtonian behavior, 𝑲′ = 𝝁 and 𝒏′ = 1 so that the
generalized Reynolds number reduces to the normal Reynolds number.
For laminar flow of a non-Newtonian fluid, the wall shear stress can be
expressed in terms of 𝑲′ and 𝒏′ as:
′
𝟖𝒖 𝒏
𝝉𝒘 = 𝑲′ (3.9)
𝒅𝒊
On substituting for 𝝉𝒘 equation 3.8, the Fanning friction factor for laminar
non-Newtonian flow becomes:
A stability analysis made by Ryan and Johnson (1959) suggests that the
transition from laminar to turbulent flow for inelastic non-Newtonian fluids
occurs at a critical value of the generalized Reynolds number that depends
on the value 𝒏′ . The results of this analysis are shown in Figure 3-5. This
relationship has been for shear thinning and for Bingham plastic fluids and
number with 𝒏′ .
has been found to be accurate. Over the range of shear thinning behavior
encountered in practice, 0.2 ≤ 𝒏′ < 1.0, the critical value of 𝑹𝒆′ is in the
range 2100 ≤ 𝑹𝒆′ ≤ 2400.
Experimental results for the Fanning friction factor for turbulent flow of
shear thinning fluids in smooth pipes have correlated by Dodge and
Metzner (1959) as a generalized form of the von K𝐚rm𝐚n equation:
Figure 3-6 Friction factor chart for purely viscous non-Newtonian fluids.
Having determined the value of the friction factor 𝒇 for specified flow
rate and hence 𝑹𝒆′ , the pressure gradient can be calculated in the normal
way using equation 2.12.
𝟐 𝟐𝒇𝑳𝝆𝒖𝟐
𝑳 𝝆𝒖
∆𝑷𝒇 = 4𝒇 = (3.12)
𝒅𝒊 𝟐 𝒅𝒊
For the conditions at the pipe wall, denoted by the subscript 𝒘, the
equation of power law fluid can be written as:
𝒏
𝐝𝒖
𝝉𝒘 = 𝑲 (− 𝜸𝒘 ) = 𝑲 ( − )𝒏𝒘 (3.14)
𝐝𝒚
The minus sign has been placed inside the parentheses recognizing the
fact that the shear rate 𝜸 (equal to 𝐝𝒖/𝐝𝒚) is negative. 𝜸𝒘 is the true shear
rate at wall and is related to the flow characteristic (𝟖𝒖/𝒅𝒊 ) by the
Rabinowitsch- Mooney equation:
𝟖𝒖 𝟑𝒏′ + 𝟏
− 𝜸𝒘 = (3.15)
𝒅𝒊 𝟒𝒏′
𝒏′ = 𝒏 (3.17)
𝟑𝒏+ 𝟏
𝑲′ = 𝑲 (3.18)
𝟒𝒏
Nomenclature
𝑲 consistency coefficient, Pa 𝐬 𝒏
𝑳 length of pipe, m
𝝁 dynamic viscosity, Pa s
𝝉 shear stress, Pa
𝝉𝑩 yield stress, Pa
Solution:
𝝆𝒖𝒅𝒊
𝑹𝒆′ =
𝝁𝒂𝒑
𝟖𝒖 𝟖 (𝟐.𝟎 𝐦/𝐬)
= = 210 s-1
𝒅𝒊 𝟎.𝟎𝟕𝟔𝟐 𝐦
and
′
𝟖𝒖 𝒏 −𝟏
= (210 s-1 ) (0.3-1.0) = 0.0237 s0.7
𝒅𝒊
Hence
′
(𝟗𝟔𝟏 𝐤𝐠/𝐦𝟑 ) 𝟐.𝟎 𝐦.𝐬 𝟎.𝟎𝟕𝟔𝟐 𝐦
𝑹𝒆 = = 4178
(𝟎.𝟎𝟑𝟓𝟏 𝐏𝐚 𝐬)
From Figure 3-6, the Fanning friction factor 𝒇 has a value 0.0047.
Example 3-2
Calculate the frictional pressure gradient ∆𝑷𝒇 /𝑳 for a time independent non-
Newtonian fluid in steady state flow in a cylindrical tube if:
𝟖𝒖 𝟖 (𝟏.𝟎 𝐦/𝐬)
= = 100 s-1
𝒅𝒊 𝟎.𝟎𝟖 𝐦
and
′
𝟖𝒖 𝒏 −𝟏
= (100 s-1) (0.5-1.0) = 0.1 s0.5
𝒅𝒊
Hence
and
′
(𝟏𝟎𝟎𝟎 𝐤𝐠/𝐦𝟑 ) 𝟏.𝟎 𝐦.𝐬 𝟎.𝟎𝟖 𝐦
𝑹𝒆 = = 400
(𝟎.𝟐 𝐏𝐚 𝐬)
𝟏𝟔
From Figure 3-6, at 𝑹𝒆′ = 400, the Fanning friction factor 𝒇 = :
𝑹𝒆′
𝟏𝟔
then 𝒇= = 0.04.
𝟒𝟎𝟎
Therefore the frictional pressure gradient ∆𝑷𝒇 /𝑳 is given by:
𝟐𝒇𝝆𝒖𝟐
∆𝑷𝒇 /𝑳 =
𝒅𝒊
𝒅𝒖 𝒏
𝝉𝒘 = 𝑲 (− )
𝒅𝒚
𝒅𝒖
𝝉𝒘 = 10 (− )0.2
𝒅𝒚
When
𝒅𝒖
= 10 s-1 , 𝝉𝒘 = 10 * (𝟏𝟎)𝟎.𝟐 = 15.85 N/𝐦𝟐
𝒅𝒚
𝒅𝒖
= 50 s-1 , 𝝉𝒘 = 10 * (𝟓𝟎)𝟎.𝟐 = 21.87 N/𝐦𝟐
𝒅𝒚
𝒅𝒖
𝝉𝒘 = 𝝉𝑩 + 𝝁𝒂 (− )
𝒅𝒚
15.85 = 𝝉𝑩 + 𝝁𝒂 (10),
21.8 7 = 𝝉𝑩 + 𝝁𝒂 (50),
𝒅𝒖
Thus, Bingham plastic equation is: 𝝉𝒘 = 14.35 + 0.15 (− )
𝒅𝒚
Solution:
𝒅𝒖 𝒏
𝝉𝒘 = 𝑲 (− )
𝒅𝒚
𝒅𝒖 𝒏−𝟏 𝒅𝒖
𝝉𝒘 = 𝑲 − ∗ − ………. (1)
𝒅𝒚 𝒅𝒚
𝝉𝒘
Since 𝝁𝒂 = 𝒅𝒖 ………. (2)
−
𝒅𝒚
𝒅𝒖 𝒏−𝟏
𝝁𝒂 = 𝑲 − ………. (3)
𝒅𝒚
𝟎.𝟑𝟑−𝟏
𝝁𝒂 = 𝑲 𝟏𝟎𝟎𝟎
−𝟎.𝟔𝟕
𝟎. 𝟏 = 𝑲 𝟏𝟎𝟎𝟎
𝑲 = 10 Pa sn
𝒅𝒖 𝟎.𝟑𝟑
Rheological equation is: 𝝉𝒘 = 10 (− )
𝒅𝒚