Understanding Rheology: Ross Clark
Understanding Rheology: Ross Clark
Understanding Rheology: Ross Clark
Ross Clark
Distinguished Research Fellow
San Diego R&D
Page 1
Background
Page 3
Types of deformation
Page 4
Basic Deformations
Extension (Tension)
Bending
Simple Shear
Page 5
Laminar Shear Flow
1 meter
separation
1meter / second 1
Shear rate
1meter second
Since viscosity is defined as shear stress/shear rate, the final units (in the SI system) for
viscosity are Newtonsecond / m2. This can also be given as a Pascalsecond since a
Pascal is one Newton / m2.
In the more traditional physics units, the units of viscosity are dynesseconds / cm2. This
is defined as 1 Poise
Density is sometimes included for gravity driven capillary instruments. This unit is a
centiStoke.
Page 6
More terms
Page 7
Viscoelasticity measurement
Viscoelastic Response of a Perfectly Elastic Sample
Strain
Stress
0 10 0
Strain
Stress
Viscous materials, like water or thin oils, will
always have stress and strain shifted 90 from
each other. This is because the most
resistance to movement occurs when the rate
-1 . 2
of the movement is the greatest.
0 10 0
1. 2
Viscoelastic Response of a Xanthan Gum Sample
Strain
Stress
Most of the world is viscoelastic in nature and
so shares characteristics of elastic and viscous
materials. The phase shift ( ) will always be
between 0 and 90.
-1 . 2
0 10 0
Page 8
Where do G and G come from?
Viscoelastic Response of a Perfectly Elastic Sample Viscoelastic Response of a Perfectly Viscous Sample
1 . 2 1 . 2
Viscoelastic Response of a Xanthan Gum Sample
1 . 2
- 1 . 2 - 1 . 2
- 1 . 2
0 1 0 0 0 1 0 0
0 1 0 0
G (viscous)
G (viscous)
G (elastic) G (elastic) G (elastic)
Page 9
What is the use of viscoelasticity?
Page 10
Still more terms for steady shear
Suspension
Stabilization
Pouring
Viscosity
Mouthfeel
Dispensing
Pumping
Coating
Spraying
10 -6 10-5 10-4 10-3 10-2 10-1 1 10 10+2 10+3 10+4 10+5 10+6
Shear Rate
Page 12
Couette Rheometers, Design
Page 13
Couette Calculations
T
=
2* * Rb2 * Hb
. 2* * rpm Rc2
=2* *
60 Rc2 - Rb2
Cup Wall Bob Wall Cup Wall Bob Wall Cup Wall Bob Wall
Page 15
Couette, nonNewtonian Corrections
For all cases except where the gap between the cup & bob is very small, that is Rc/Rb > 0.95, we must
correct for the flow field in the gap. When a more pseudoplastic fluid is tested, the shear rate in the gap
tends to be the highest near the rotating member (the bob in the case of the Brookfield). This is because the
shear stress is at a maximum at this point and the fluid tends to flow faster under high shear stress values.
In any event, we need to correct for this flow profile in the gap. This is most commonly done by assuming
that the material will obey the power law, that is = K n.
Below, a step by step procedure is listed for correction of the shear rate for non-Newtonian fluids in Couette
viscometers: .
Step #1, calculate the shear stress, on the bob using the equation for Newtonian fluids.
Step #2, calculate the value of for each of the speeds used.
Step #3, make a loglog plot of and . Calculate the value of n (slope) and b (intercept).
Step #4, determine the correction factor, as: = 1+ [ln (s)/ n] where s = Rc2
ln (s) Rb2
Step #5, determine the K value as K=b n
. .
Step #6, determine the shear rate, as =
.
Step #7, determine the viscosity, as = /
Page 16
Couette Errors
In the table given below, the shear rate and errors associated with the Newtonian calculation of shear rate
for two different theoretical pseudoplastic fluids are given. In each case the speed of the viscometer is the
same (0.3 rpm) and the Rb/Rc calculation is given for the three different viscometer gaps.
It can be seen that the errors become very large for pseudoplastic fluids in wide gap instruments. This is
one reason to avoid the use of small spindles in the Brookfield small sample adapter or any other
viscometer.
Newtonian Pseudoplastic Error Pseudoplastic Error
Couette Gap n=1 n=0.658 (%) n=.261 (%)
Narrow 0.967 0.974 0.99 1.64% 1.063 9.14%
An example of some of the dimensions for various Brookfield attachments is given in the table below. As you
can see, the error associated with even the best conditions (#18 bob with the small sample adapter or the
UL) is still significant, especially when the degree of pseudoplasticity of most of our fluids is taken into
account.
Bob Radius Cup Radius
(mm) (mm) Rb/Rc
Small Sample Adapter, #18 bob 8.7175 9.5175 0.916
Page 17
Effect of Couette corrections
100
10 Nearly
2x higher
(s -1 )
for xanthan!
Shear rate,
0.25% Keltrol
Glycerin / Water
0.1
1 10 100
Spindle speed (rpm)
Page 18
Cone/Plate Rheometers, Design
Page 19
Cone/Plate Calculations
T
=
2/3* * Rc3
. 2* * rpm / 60
=2*
sine ( )
Top View
Torque . 2 (rpm / 60)
2/3 r3 sine
10 mm If we assume a speed of 12 rpm or 1.26 rad/s:
By the equation given above for shear rate we get 48.0 s-1.
Page 21
Capillary Rheometers, Design
Page 22
Capillary Calculations
p * Rc
=
2 * Lc
. 4*Q
=2*
* Rc3
Page 24
Newtonian
.
Equation: =K*
Newtonian
3 6
10 10
For this graph:
5 K = 150
10
Shear stress
4
10
rare. Low molecular
Viscosity
2
10
2 1
10 -1 0 1 2 3
10
10 10 10 10 10
Rate
Page 25
Bingham Plastic
.
Equation: = o + (K *
For this graph:
Bingham K = 120
3 6
10 10
o = 20
5
10
This modification of the
Newtonian model allows
for a yield stress. This
Shear stress
4
10
Viscosity
Page 26
Casson
.
Equation: = o + K*
For this graph:
Casson model
10
3
10
5 K = 40
o = 20
4
10
This is a variation of the
Bingham model. It is
Shear stress
Viscosity
frequently used to
2 3
10 10 extrapolate to a yield
stress from low shear
rate data. Has been
2 successfully used with
10
chocolate.
1 1
10 -1 0 1 2 3
10
10 10 10 10 10
Rate
Page 27
Power Law
. n
Equation: =K*
For this graph:
Power Law
3 4
10 10 K = 150
n = 0.5
2 3
10 10 The power law model
Shear stress
is the most frequently
Viscosity
Page 28
Ellis (Power Law with Yield Stress)
. n
Equation: = o + K * (
For this graph:
Ellis or power law + yield stress
10
3 4
10 K = 120
n = 0.5
o = 25
2 3
10 10
Shear stress
Viscosity
This modification of
the power law model
1 2
allows for a yield
10 10
stress like the
Bingham. It is difficult
to fit; 2 known
parameters, 3
0 1
10 -1 0 1 2 3
10 unknowns.
10 10 10 10 10
Rate
Page 29
Cross Model
Equation: = + o
.n For this graph:
1+ K * (
Cross model = 0.002
3 4
10 10
o = 500
3
10
10
2
K = 50
2
10
1
n = 0.8
10 1
10 This model fits
Shear stress
Viscosity
10
0 0
10
viscosity data rather
than shear stress
-1
-1
10 data. It allows for
10
-2
upper and lower
10
-2
Newtonian viscosity
10 -3
10 values. With 4
-3 -4
parameters, it
10 -7 -5 -3 -1 1 3 5
10 requires nonlinear
10 10 10 10 10 10 10
Rate methods.
Page 30
Time sweep test
1% xanthan in 0_05M N aC l AR ES time sweep
3
10
1% Xanthan in
0.05 M NaCl
50 mm parallel
plate, 0.5 mm
gap
)
23C
[dyn/cm ]
2
10 rad/s
G" (
50% strain 2
10
ARES instrument
)
[dyn/cm ]
2
G' (
1
10
5 0 .0 10 0 .0 1 5 0 .0 2 0 0 .0 2 5 0 .0 3 0 0 .0
0.0 350.0
tim e [s]
Page 31
What is learned from a time sweep?
Page 32
Strain sweep test
1% xanthan in 0_05M N aC l AR ES strain sweep
3
10
1% Xanthan in
0.05 M NaCl
50 mm parallel
plate, 0.5 mm
gap L i n e a r vi s co e l a sti c re g i o n
)
2
23C 10
[dyn/cm ]
2
10 rad/s
G" (
1 to 1000% strain
ARES instrument
)
[dyn/cm ]
2
G' (
1
10
0
10
1 2
0 10 10 3
10 10
Strain [% ]
Page 33
What is learned from the strain sweep?
Page 34
Frequency sweep test
1% xanthan in 0_05M N aC l AR ES freq sweep
3
10
1% Xanthan in
0.05 M NaCl
50 mm parallel
plate, 0.5 mm
gap
)
23C
[dyn/cm ]
2
50% strain 2
10
ARES instrument
)
[dyn/cm ]
2
G' (
1
10
-1 0 1
-2 10 10 10 2
10 10
Freq [rad/s]
Page 35
What is learned from a freq. sweep?
Page 36
Stress ramp test
0.216% Primacel (0.36% total) in STW control
1000.0
1% bacterial 9 0 0 .0
cellulose in tap
water
8 0 0 .0
50 mm parallel
plate, 0.5 mm
7 0 0 .0
gap
23C 6 0 0 .0
)
5 0 0 .0
over a 120
(t) (
second period 4 0 0 .0
SR-2000
Yield stress:(2 0.59 dynes/cm 2 @ 139.57% strain)
instrument 3 0 0 .0
2 0 0 .0
1 0 0 .0
0.0
5 .0 1 0 .0 1 5 .0 2 0 .0 2 5 .0
0.0 30.0
2
(t) [dyn/cm ]
Page 37
Why do a stress ramp test?
Page 38
Temperature sweep test
Figure 1. E ffect o f A P V G aulin H o m o genizatio n o n G ellan G um S et T em perature
3
10
C o n tro l Ts~ 97C
0.5% gellan gum T re a tm e n t # 1 T s 8 2 C
in 4 mM CaCl 2 T re a tm e n t # 2 T s 7 4 C
50 mm parallel
plate, 0.5 mm
gap
23C
Dynamic test at 10
2
)
10 rad/s and 5%
[dyn/cm ]
2
strain
G' (
ARES instrument
with Peltier
heating / cooling
S e t te m p e ra tu re d e fi n e d a s G ' va l u e
10
1 o f 1 0 d yn e s/cm 2
4.0
6 5 .0 7 5 .0 8 5 .0
55.0 95.0
< ------ C o o l i n g te st ru n fro m h i g h te m p to l o w < -------
Tem p [C]
Page 39
Why do a temperature sweep test?
Page 40
Creep test
1% xanthan in 0_05M N aC l creep test
3 2
10 10
1% Xanthan in
0.05 M NaCl
50 mm parallel
plate, 0.5 mm
gap
23C
stress(t) (
)
50 dynes / cm 2
[dyn/cm ]
stress
2
[% ]
Strain(t) (
10
2
)
1 1
10 10
1 0 0 .0 2 0 0 .0 3 0 0 .0 4 0 0 .0 5 0 0 .0 6 0 0 .0 7 0 0 .0 8 0 0 .0
0.0 900.0
tim e [s]
Page 41
What is learned from a creep test?
Page 42
Spring and dashpot models
Used to provide a physical model of a materials properties
Spring-represents a Dashpot-represents a
perfectly elastic element perfectly viscous element
2 3
10 10
1% Xanthan in
0.05 M NaCl
50 mm parallel
plate, 0.5 mm
gap
23C
)
s tre s s (t) (
50% strain
[d yn /c m
ARES instrument 1 2
S tra in (t) (
10 10
[% ]
2
]
)
0 1
10 10
0 .0 5 0 .0 1 0 0 .0 1 5 0 .0 2 0 0 .0 2 5 0 .0
tim e [s ]
W e d n e s d a y , S e p te m b e r 0 6 , 2 0 0 0 R o s s C la r
Page 44
Startup (stress overshoot) test
8.0
7.0
6.0
5.0
)
[Pa]
stress(t) (
4.0
3.0
2.0
1.0
0.0
0.0 100.0 200.0 300.0 400.0 500.0 600.0 700.0 800.0 900.0 1000.0
Strain [%]
Page 45
What is learned from a stress overshoot?
Page 46
Creep frequency transform
G ', G " from 1% xanthan in 0_05M N aC l creep test
3
10
1% Xanthan in
0.05 M NaCl
T h e se d a ta a re fro m a fre q u e n cy sw e e p cu rve
50 mm parallel
plate, 0.5 mm
gap
)
23C
[dyn/cm ]
2
Frequency sweep
G" (
data collected
with auto stress 10
2
adjust. Strain
)
[dyn/cm ]
controlled from
2
20 to 50%
G' (
Creep data
collected at 50
dynes / cm 2 stress
SR-2000
instrument
These data are from a transform ation of
the creep cu rve data.
1
10
-2 -1 0 1
-3 10 10 10 10 2
10 10
Freq [rad/s]
Page 47
Why use a creep freq. transform?
Page 48
Relaxation frequency transform
3
10
1% Xanthan in
0.05 M NaCl
50 mm parallel
plate, 0.5 mm
gap
)
23C
[dyn/cm ]
2
Frequency sweep
G" (
data collected
with auto stress 10 2
adjust. Strain
)
controlled from
[dyn/cm ]
2
20 to 50%
G' (
Creep data
collected at 50
dynes / cm 2 stress
Solid lines from relaxation data Dashed lines from frequency sweep
SR-2000
instrument
1
10
-3 10 -2 10 -1 10 0 10 1 10 2 3
10 10
Freq [rad/s]
Page 49
Why use a relaxation freq transform?
Page 50
Time Temperature Superposition
Page 51
TTS data
[Pg13]:TTS Session--TTS O verlay C urve
6
10
0.5% gellan gum
in 4 mM CaCl 2
25 mm parallel 5
plate, 0.5 mm 10
gap
-35 to +55C
)
[dyn/cm ]
2
Frequency sweep 10
4
data collected at
G" (
10% strain
ARES instrument
)
[dyn/cm ]
2
3
10
G' (
2
10
1
10
0 1 2
-1 10 10 10 3
10 10
Freq [rad/s]
Page 52
TTS data after shifting
[Pg18]:TTS Session--TTS O verlay C urve
7
10
minimization. 10
Cubic spline
interploation
)
[dyn/cm ]
2
5
10
G" ( )
[dyn/cm ]
2
4
10
G' (
3
10
2
10
-5 -4 -3 -2 -1 0 1
-6 10 10 10 10 10 10 10 2
10 10
Freq [rad/s]
For more information concerning minimization methods developed by Brent and Powell, see W. H.
Press, et. al., Numerical Recipes in C, Cambridge University Press, 1992, ISBN 0 521 43108
Page 53
TTS master curve
7
10 90.0
80.0
106
70.0
)
60.0
[dyn/cm ]
PhaseAngle (
2
105
G" (
50.0
[]
)
[dyn/cm ]
40.0
2
4
10
G' (
)
30.0
20.0
103
10.0
2
10 0.0
-6 10-5 10-4 10-3 10-2 10-1 100 101 2
10 10
Freq [rad/s]
Page 54
Practical tips 1
Page 55
Practical tips 2
Page 56
Practical tips 3
Page 58
Most overlooked tests
Page 59
Most overused tests
Page 60
In conclusion
Page 61