Capítulo 4

Fluidos no-Newtonianos

Este capítulo abarca a los fluidos Newtonianos generalizados (ley de potencia y viscoplás-
ticos) en flujos tangenciales y en tubos.

4.1. La ecuación de Weissenberg-Rabinowitch (W-R).

Esta ecuación se utiliza para obtener la viscosidad no-Newtoniana de datos que se obtienen
de un viscosímetro capilar. Se parte de la definición de flujo volumétrico:
Z R
Q = 2π vz r dr (4. 1)
0

Realizando un cambio de variable: u = vz ,

R
r2 
  Z R 2  
∂ vz r ∂ vz
du = dr, se obtiene: Q = 2 π vz  − 2π dr
∂r 2 0 0 2 ∂r
 
∂ vz
Definiendo la velocidad de corte: γ̇w = − , obtenemos:
∂r
Z R
Q = π γ̇w r2 dr (4. 2)
0

Por balance de momentum lineal:

   
∂p r ∆p r
τrz = = (4. 3)
∂z 2 L 2

y definiendo el esfuerzo en la pared:

 
∆p R
τw = (4. 4)
L 2

1 R
obtenemos: dr = dτrz = dτrz
∆p/2L τw
Z τw 3
R 2


luego, la ecuación 4. 2 puede expresarse como: Q = π γ̇w τrz dτrz
0 τw3

1
o similarmente:
Q τw3
Z τw
2


= γ̇w τrz dτrz (4. 5)
π R3 0

Derivando con respecto a τw :

τw3 dQ 3Q 2
τw = γ̇w τw2


3
+ 3
π R dτw πR

 
1 dQ
τw + 3Q = γ̇w , obtenemos:
π R3 dτw

   
4Q 3 1 τw dQ 4Q 3 1 d lnQ
γ̇w = + = + (4. 6)
π R3 4 4 Q dτw π R3 4 4 d lnτw

que es la ecuación W-R.

En el caso de un fluido Newtoniano, y por la ley de Hagen-Poiseuille, el esfuerzo en la

pared está dado por:

   
∆p R 4Q
τw = = µ (4. 7)
L 2 π R3
4Q
En este caso d ln τw = d ln Q, por lo que la ecuación 4. 6 se reduce a γ̇w = π R3
. Con
esta ecuación es posible obtener la viscosidad no Newtoniana del fluido:

 
τw
η = (4. 8)
γ̇w

4.2. Flujos no Newtonianos en tuberías.

a) Ley de potencias.
En este caso, la ley de Newton se generaliza a:
∂ vz n
 
τrz = K (4. 9)
∂r

Donde K (índice de consistencia) y n (coeficiente de la ley de la potencia) son constantes.

Por balance de momentum,
∂ vz n
   
∆p r
τrz = = (4. 10)
L 2 ∂r

2
  1/n Z
r /n dr + C1 , con condiciones de frontera: vz = 0 en r = R.
∆p 1
De donde vz = 2KL
Aplicando las condiciones de frontera obtenemos:
1/n 1
nR( n +1 )  r  1 +1 
 
∆p n
vz = 1− (4. 11)
2KL n+1 R

En el caso Newtoniano (K = µ y n = 1) se obtiene:

∆p R2
   r 2 
vz = 1− (4. 12)
4µL R

Lo que recupera el perfil parabólico. El gasto volumétrico correspondiente es:

1/n 1
nR( n +1 )  r  1 −1 
ZZ Z 2π Z R  Z R
∆p n
Q= v · n dS = vz r dr dθ = 2π 1− r dr
0 0 2KL n+1 0 R

Z R  r  1 −1  n+1
n
Ya que 1− r dr = R2 , entonces:
0 R 2(3n + 1)

1/n
n π R3

∆p R
Q = (4. 13)
2KL 3n + 1

Nuevamente, si K = µ y n = 1, obtenemos la ley de Poiseuille:

π R4 ∆p
Q = (4. 14)
8µ L

La ecuación 4. 13 puede ser expresada como:

4 Q 3n + 1 n
  
∆p R
τw = = K (4. 15)
L 2 π R3 4n

En donde:
 
4Q 3n+1
γ̇w = (4. 16)
π R3 4n

Comparando la ecuación 4. 16 con la ecuación 4. 6 se obtiene:

3 1 d ln Q 3n+1
+ = (4. 17)
4 4 d ln τw 4n

debido a que por la ecuación 4. 15 d ln τw = (n) d ln Q.

3
b. Fluidos de Bingham (viscoplásticos).
Estos fluidos poseen un esfuerzo característico, el esfuerzo de cedencia, el cual es necesario
superar para que el líquido fluya. Para valores menores a este esfuerzo, el material es un
sólido que no fluye. En un tubo, el esfuerzo de cedencia se encuentra a una distancia radial
determinada (r0 ), de acuerdo a la figura 4.1.

Figura 4.1. Vista del plano r −z de un tubo.

La zona sólida (r < r0 ) se encuentra en la
parte central y la zona líquida (r > r0 ) se
encuentra en la vecindad de la pared.

En este caso, la componente del tensor de esfuerzos τrz está dada por:
∂ vz
τrz = τ0 + µ0 (4. 18)
∂r

donde τ0 es el esfuerzo de cedencia. Por el balance lineal de momentum se tiene que:

 
∆pr
τrz = si r ≥ r0
L 2
 
∆p r0
τrz = τ0 = si r ≤ r0
L 2

Existen dos fases en el líquido: una fase sólida cuando τrz < τ0 y una fase fluida cuando
τrz > τ0 . Esto es:

∂ vz
τrz = τ0 + µ0 si τrz > τ0 (4. 19)
∂r

∂ vz
= 0 si τrz < τ0 (4. 20)
∂r

4
Perfil de velocidad en la fase líquida (1). Se obtiene integrando la ecuación 4. 18,
Z Z
∆p τ0
vz1 = r dr − dr +C1 . La constante C1 se evalúa con la condición:
2µ0 L µ0

vz (r=0) = 0,
resultando en:
∆p R2
   r 2  τ0 R h ri
vz1 = 1− − 1− para r ≥ r0 (4. 21)
4µ0 L R µ0 R

En la fase sólida (2), si r = r0 , entonces

∆p R2
   r 2   ∆p r  R h
0 0 r0 i
vz2 = 1− − 1−
4µ0 L R L 2 µ0 R

resultando en la siguiente simplificación:

∆p R2 h
 
r0 i2
vz2 = 1− para r ≤ r0 (4. 22)
4µ0 L R

El flujo volumétrico se evalúa por partes:

R
r2 
Z R Z R 2   Z R 
r ∂ vz ∂ vz
Q = 2π vz r dr = 2 π vz  − 2π dr = − π r2 d r
0 2 0 0 2 ∂r r0 ∂r
Z r0  
∂ vz
ya que r2 dr = 0. Por la ecuación 4. 19, tenemos que:
0 ∂r
Z R   
∆p r τ0 2 ∆p 4 4

τ0 3 3


Q = π − r dr = π R − r0 − R − r0
r0 µ0 L 2 µ0 8µ0 L 3µ0

∆p R4 r04 ∆p R4 τ0
      
4µ0 3 3


= π 1− 4 − R − r0
8µ0 L R 8µ0 L 3µ0 τw R3

∆p R4 τw R3
debido a que = , entonces:
8µ0 L 4 µ0

∆p R4 r04 4τ0 r03

    
4τ0 r0 τ0
Q = π 1− 4 − + 3
. Como = , entonces:
8µ0 L R 3τw 3τw R R τw

∆p R4 τ4 4 τ04 ∆p R4 1 τ04
    
4 τ0 4 τ0
Q = π 1 − 04 + 4
− = π 1+ −
8µ0 L τw 3 τw 3 τw 8µ0 L 3 τw4 3 τw
(4. 23)

La ecuación 4. 23 es conocida como la ecuación de Buckingham-Reiner.

5
4.3. Flujo torsional en fluidos no Newtonianos.
a) Ley de potencias.
Se hace referencia aquí al flujo similar al de un viscosímetro de Couette, en donde el fluido
está contenido entre dos cilindros concéntricos, y uno de los cilindros gira a una velocidad
angular determinada. Se obtendrá el torque en función de la velocidad angular.
La expresión de la ley de la potencia en este caso es la generalización al caso no Newtoniano
de la ley de potencia. El único componente del tensor de esfuerzos diferente de cero es:
d  vθ  n
 
τr θ = K r (4. 24)
dr r
junto con el balance de momentum:
1 d 2

− 2 r τr θ = 0 , (4. 25)
r dr
de donde se obtiene:

d  v n
θ
r2 r = C1
dr r

1 C1 /n
 1 Z  1/n
d  vθ  vθ 1 C1
luego: = 2
. Integrando: = dr + C2
dr r r r r r r2

con condiciones de frontera: vθ = 0 en r = κ R y vθ = R Ω0 en r = R

El resultado de determinar las constantes por medio de las condiciones de frontera es el
perfil de velocidades:
vθ 1 − ( κrR )2/n
= (4. 26)
Ω0 r 1 − κ 2/n
El torque resultante es:
Z L Z 2π
|T |θ = r τrz r dθ dz = r2 τrθ (2 π L) (4. 27)
0 0
2 n
κ R /n
 

d  vθ  n
  d 1− r 
donde r2 τrθ 2
= Kr r 2
= K r r  .
 
2/n
dr r  dr 1−κ 

 2/n !  
d κR 2/n2 (− 2 − 1)
Ahora: 1− = − (κR) − r n . Por lo tanto,
dr r n
"   #n n
2/n 
2 2 2 (κR) ( 2
− n − 1) 2 2
2/n
r τrθ = K r r r = Kκ R , lo que da como
n 1 − κ 2/n 1 − κ 2/n
resultado:
 n
2/n
n 2 2
|T |θ = 2 π L K Ω0 κ R (4. 28)
1 − κ 2/n

6
b) Fluidos de Bingham.
En este caso, la ecuación constitutiva de Bingham con el componente del tensor de
esfuerzos es:
∂  vθ 
τrθ = τ0 + µ0 r (4. 29)
∂r r
Este caso involucra, como uno anterior, al fluido comprendido entre dos cilindros concén-
tricos.
En la figura 4.2 el fluido presenta dos fases: una fase sólida en la región de menor esfuerzo
(que es la zona adyacente a la superficie en movimiento) y una fase líquida en la región de
mayor esfuerzo (zona adyacente a la superficie estática).

κR
R
r
0

Figura 4.2. Flujo torsional visto en el

plano (r − θ ). El cilindro interior está
estático y el exterior se mueve a una
velocidad angular Ω0 .

Las condiciones de frontera para la velocidad son:

vθ = 0 en r = κR y vθ = R Ω0 en r = R (4. 30)

Además, la condición de interfase:

vθ = r0 Ω0 en r = r0 (4. 31)

Substituyendo la ecuación 4. 27 en la 4. 29 obtenemos:

T τ0 d  vθ 
= + , integrando:
2 π L µ0 r3 r µ0 dr r
vθ T τ0
= − 2
− ln (r) + C1 Utilizando la condición 4. 31 obtenemos:
r 4 π L µ0 r µ0
T τ0
Ω0 = − 2
− ln (r0 ) + C1 Esto implica que:
4 π L µ0 r0 µ0
 
vθ T 1 1 τ0  r0 
= Ω0 + − 2 + ln
r 4 π L µ0 r02 r µ0 r

7
Por lo que el campo de velocidades entre los dos cilindros es el siguiente:

r02
   
vθ T τ0 r
= Ω0 + 1 − − ln κ R ≤ r ≤ r0 (4. 32)
r 4 π L µ0 r02 r2 µ0 r0

vθ
= Ω0 r0 ≤ r ≤ R (4. 33)
r
Si r0 → R, entonces la ecuación 4. 32 con la condición 4. 30 se modifica a:

R2
 
vθ T τ0  r 
= Ω0 + 1− 2 − ln (4. 34)
r 4 π L µ0 R2 r µ0 R

Utilizando la condición de velocidad cero en el cilindro interior se obtiene:

 
T 1 τ0
Ω0 = − 1 + ln κ (4. 35)
4 π L µ0 R2 κ 2 µ0

La ecuación 4. 35 es la ecuación de Reinier-Rivlin. Se utiliza para determinar el esfuerzo

de cedencia y la viscosidad a partir de datos viscosimétricos.

8
TAREA 6

1. Los siguientes datos fueron obtenidos en un tubo capilar de un material x.

a) En caso de que x fuera un fluido Newtoniano, calcular su viscosidad.

b) En caso de que x fuera un fluido de ley de potencia, calcular K y n.
c) Graficar la curva de viscosidad en función de la rapidez de corte (log − log)
para x.

4 Q/π R3 τw (lb f/ f t 2 )
11 1.46
24 1.97
42 2.49
75 3.30
125 4.20