Fluidos Newtonianos

CAPTULO 7
FLUIDOS VISCOSOS NEWTONIANOS
7.1 CONCEPTOS GENERALES
En los slidos, tales como los metales, se ha observado que su deformacin es proporcional
a las solicitaciones aplicadas; sin embargo, medios como el agua y el aire presentan
comportamientos muy diferentes, ya que stos no son capaces de soportar ni siquiera los
esfuerzos de corte, producto de su propio peso. Por ejemplo, al aplicarse una solicitacin a
corte entre dos placas (figura 7.1), el fluido continuar su deformacin a corte mientras la
solicitacin permanezca. Queda claro, entonces, que cualquier fluido ser incapaz de
soportar solicitaciones de corte sin deformarse de manera permanente. La velocidad de
desplazamiento ser proporcional a la solicitacin aplicada y al eliminarse la carga, la
deformacin permanecer.
FIGURA 7.1 SOLICITACIN A CORTE ENTRE DOS PLACAS ENTRE LAS

QUE SE ENCUENTRA UN FLUIDO
En presencia de la gravedad un fluido como el agua, tomar la forma del recipiente que la
contiene, resultando imposible que mantenga su forma al retirar las paredes del recipiente.
Lo anterior significa que no soporta esfuerzos cortantes generados por su propio peso,
siendo entonces que en presencia de estos cortantes el fluido se deformar de manera
permanente y continua, limitada esta deformacin solo por la tensin superficial existente
INTRODUCCIN A LA MECNICA DEL MEDIO CONTINUO
entre el fluido y la superficie sobre la cual se extiende. Con base en las condiciones de
movimiento del fluido, se define a ste como un medio idealizado que durante su movimiento
como cuerpo rgido (considerando el propio estado de reposo) es incapaz de soportar
cualquier tipo de solicitaciones a corte. Asimismo, se tiene que para algunos casos la
densidad del fluido es aproximadamente constante. Esta situacin aplica, por ejemplo, para
el agua, en la cual en condiciones de carga muy variadas se considera que su densidad no
se altera (por lo tanto se describe como incompresible), por otra parte, el aire, como todos
los gases, se analiza sobre la premisa de que al variar la presin su densidad tambin se ve
afectada. Sin embargo, la descripcin de compresibilidad o invariabilidad de la densidad en
un fluido depende de las condiciones del estudio; por ejemplo, el aire a bajo nmero de
Match se le analiza como si se tratara de un fluido incompresible; por lo contrario, cuando se
estudia la propagacin de ondas elsticas en el agua se describe a sta como un fluido
compresible.
En los fluidos se observa que la resistencia al flujo depende de la velocidad y, por

consecuencia, de su velocidad de deformacin, esto de manera anloga a lo que sucede en
los slidos con relacin a su deformacin. Anlisis ms detallados revelan que existen fluidos
en los que la relacin de las cargas aplicadas con la velocidad de deformacin es lineal; tal
como pasa en los slidos de Hooke con la deformacin. Por otra parte, fluidos como la miel o
la propia sangre no presentan relaciones lineales. Es entonces que se pueden clasificar a los
fluidos como:
1. Fluidos newtonianos. Son aquellos en los que la relacin esfuerzo de corte-

velocidad de deformacin es lineal = c . A esta relacin de proporcionalidad se le
denomina como viscosidad, razn por la cual este tipo de fluidos se describen como
linealmente viscosos.
2. Fluidos no newtonianos. En este caso la relacin es no lineal, presentndose

fenmenos de almacenamiento de energa a la vez de los disipativos caractersticos
de los fluidos, a este tipo de medios se les denomina como fluidos viscosos no
lineales ( = c ; n 1)
n
Desde el punto de vista de la variacin de su densidad se describen como:
344
CAPTULO7. FLUIDOS VISCOSOS NEWTONIANOS
a) Fluidos compresibles
D
0
Dt
b) Fluidos incompresibles. Se puede considerar idealmente que la densidad del fluido

bajo estudio permanece constante
D
=0 v = 0
Dt
u = 0
Como ya ha sido mencionado un fluido es un medio idealizado, el cual en cualquier punto,

durante movimiento de cuerpo rgido o en reposo, no es capaz de soportar esfuerzos de
corte, por lo que con cualquier base que se analice el sistema, el estado de esfuerzos
siempre se presentar como
11 0 0
ij = 0 22 0

0 0 33

Esto debido a que el fluido en reposo o en movimiento de cuerpo rgido no presenta ningn
esfuerzo de corte; por otro lado, partiendo de la misma lgica se tiene que para un elemento
diferencial cualquiera en el seno del fluido al cortar ste con cualquier plano, solamente se
presentarn esfuerzos normales, lo que se expresa entonces como:
Tn = n
Para cualquier n , el esfuerzo en cualquier punto y para cualquier plano es normal al plano.
Considerando que el punto (elemento diferencial de fluido) es cortado por dos planos
cualesquiera cuyas normales son n1 , n2 , entonces se cumplir que
Tn1 = 1n1
Tn2 = 2 n2
345
n1Tn1 = 1 n1 n2
Por otra parte,
n1Tn2 = 2 n1 n2
dado que T = T T
n1Tn2 n2Tn1 = 0 = (2 1 )n1 n2
Existen entonces dos posibilidades:
a) n1 n2 = 0
b) 2 1 = 0
Ya que los planos no necesariamente son perpendiculares n1 n2 0
En otras palabras, en todos los planos que pasan a travs del punto no slo no existen
esfuerzos de corte, adems todos los esfuerzos normales son iguales, a estos se les
denomina como esfuerzos hidrostticos y representan una componente esfrica. Como los
esfuerzos en el seno del fluido deben ser compresivos, entonces
H = p
p = presin atmosfrica
T = pI
En notacin ndice se expresa
ij = p ij (7.1)
donde el escalar p representa la magnitud de los esfuerzos normales compresivos y por

consecuencia se define como presin hidrosttica.
346
7.2 FLUIDOS COMPRESIBLES E INCOMPRESIBLES
Como ya fue explicado, en diversos fluidos, como el agua, al variar la presin hidrosttica su
densidad se modifica en magnitudes tan pequeas que se definen a stos como
incompresibles, por lo cual se cumple que
D
=0
Dt
razn por la que la ecuacin de conservacin de la masa
D v
+ i =0
Dt xi
se reduce a
v = 0
O en notacin ndice
vi
=0
xi
Un fluido, aun cuando se considere como incompresible, no necesariamente deber

presentar uniformidad espacial en la densidad, de ser as se considerar como homogneo.
Resulta evidente que la suposicin de incompresibilidad simplificar el anlisis y por
consecuencia la solucin de los problemas, simplemente es necesario en cada caso evaluar
la conveniencia de considerar al fluido como compresible o incompresible, de tal forma que
se obtengan soluciones simples y con un alto grado de aproximacin.
7.3 ECUACIONES DE LA HIDROSTTICA
Un fluido, como todo medio continuo, deber cumplir con las ecuaciones generales, entre
stas la ecuacin de Cauchy, la cual considerando condiciones de equilibrio queda como
ij
+ Bi = 0
x j
347
Por otra parte, se ha definido que el estado de esfuerzos para un fluido en reposo se
representa por
ij = p ij
por lo que sustituyendo sta en la Ecuacin de Cauchy, se tiene
p
ij + Bi = 0
x j
p
+ Bi = 0
xi
p
= Bi (7.2)
xi
Lo que en notacin general se expresa como
p = B
EJERCICIO 7.1 Determine la variacin de la presin de un objeto que se encuentra sumergido

en un lquido de densidad , si se conoce que la presin sobre la superficie del fluido se
describe como p0 .
Con base en la ecuacin 7.2, se puede determinar la variacin de la presin en funcin de la

profundidad a la que se encuentra inmerso el medio al interior del fluido. Considrese que la
nica fuerza de cuerpo es producto del campo gravitacional y que el eje x3 corresponde con
la vertical, mientras que el plano horizontal est dado por x1 x2 ; la aceleracin producto de
las fuerzas de cuerpo ser
Bi = 0e1 + 0e2 + ge3 .
El fluido como cualquier otro medio deber cumplir con las ecuaciones de Cauchy en
equilibrio:
ij
+ Bi = 0
x j
348
ij = p ij
p
= Bi
xi
p = B
Considerando la aceleracin gravitatoria
B1 = 0, B2 = 0, B3 = g
p0 0
p x2
=0
x1
x1 x3
p
=0
x2
p
= g
x3
Por lo tanto, se tiene que p ( x1 , x2 ) = cte para cualquier x1 , x2 ; por otra parte:
p = gx3 + c
dado que en x3 = 0 p = p0 , se tiene que
p = gh + p0
EJERCICIO 7.2 Cuando se presentan diferencias de altitud menores, se puede considerar que
la atmsfera se encuentra a temperatura constante. Con base en lo anterior, determine las
ecuaciones que describan la variacin de presin y densidad de la atmsfera.
p
Con base en la ecuacin de la hidrosttica + Bi = 0 , se tiene que, para el eje x1
xi
p
=0 p = cte
x1
349
Para el eje x2
x3
p
=0 p = cte
x2
Para el eje x3
p x1 x2
= g
x3
Considerando la ecuacin de estado para un gas ideal, y adems que la temperatura es

constante
pV = mR
p = R , con = cte
p
=
R
Sustituyendo resulta entonces que

p p
= g
x3 R
Resolviendo la ecuacin diferencial se tiene
dp gdx3
=
p R
g
Ln p = x3 + c
R
Para x3 = 0 p = p0 (presin de referencia)
g
Ln p0 = (0) + c
R
g
p = p0 exp x3
R
350
De forma anloga se puede proceder para determinar la razn de variacin de la densidad

con la altura
p = R
R
= g
x3
d g
= dx3
R
g
Ln = x3 + c
R
Para x3 = 0 = 0 (densidad inicial o de referencia)
g
= 0 exp x3
R
7.4 MOVIMIENTO DE CUERPO RGIDO DEL FLUIDO
En los casos anteriores se consider que el fluido se encontraba en reposo o en condiciones

de equilibrio, por lo que se defini que su aceleracin era igual a cero, ahora se analizar
considerando que se trata de un movimiento de cuerpo rgido, por lo cual v = 0 , es decir,
las deformaciones y velocidades de deformacin son cero, por lo que la ecuacin de Cauchy
para este caso se expresa como
p
ij + Bi = ai
x j
p v
+ Bi = i
xi t
p + B = a (7.3)
351
EJERCICIO 7.3 Un recipiente con un fluido incompresible en su interior se mueve

verticalmente hacia arriba con una aceleracin constante a . Determine la presin en un
punto que se encuentra a una profundidad H de la superficie. Considere que en la superficie
H = 0 la presin est dada por p0 . La densidad se expresa como 0 .
El sistema coordenado se define considerando que la direccin positiva del eje x3 es hacia
abajo, mientras que los ejes x1 , x2 corresponden al plano horizontal. La aceleracin producto
de las fuerzas de cuerpo queda (slo se considera el campo gravitacional).
FIGURA 7.2 RECIPIENTE MOVINDOSE EN DIRECCIN VERTICAL
Bi = 0e1 + 0e2 + 0e3
Movimiento de cuerpo rgido
vi
=0
x j
Eje x3
p
+ g = a3
x3
La solucin de la ecuacin diferencial queda
p = ( g + a3 ) x3 + p0
352
EJERCICIO 7.4 Un vehculo arrastra una pipa de seccin cilndrica (figura 7.3), la cual tiene
una divisin central. El tanque tiene una longitud de 15 m por 2 m de dimetro. El fluido
dentro del tanque ocupa un 50 % del volumen de ste. Al ponerse la luz del semforo en
verde el vehculo debe acelerar con una magnitud constante (aceleracin en direccin
horizontal). Considerando movimiento de cuerpo rgido, determine el ngulo de la superficie
libre del tanque ( ) con relacin a la horizontal; asimismo, desarrolle la ecuacin que define
la presin para cualquier punto del tanque, tal que p = p ( x1 , x2 ) . Calcule la altura mxima
que alcanza el fluido al chocar con la pared vertical si a = g / 4 .
t = t0 t > t0
FIGURA 7.3 MOVIMIENTO DEL AGUA AL PRESENTARSE UNA ACELERACIN
La superficie del agua es normal a la resultante de la fuerza, si solo existe la aceleracin

producto del campo gravitacional, la superficie permanecer horizontal; sin embargo, al
acelerar el vehculo el agua, por efecto de su inercia, tiende a desplazarse hacia atrs dando
lugar a una superficie inclinada con ngulo con respecto a la horizontal. Esta superficie
ser perpendicular a la componente de la aceleracin resultante, por lo que
a 1
= tan 1 = tan 1 = 14

g
4
353
A partir de lo anterior y considerando las dimensiones del tanque ( , ) , donde representa

la longitud de la cmara y su dimetro, se calcula fcilmente la altura que alcanza el
lquido durante la aceleracin del vehculo. Adems, a partir de la ecuacin de movimiento
de cuerpo rgido de un fluido se puede determinar la ecuacin que defina la presin para
cualquier coordenada.
Tomando como referencia la superficie del fluido
p ( x1 , x2 ) = p0 + gx2 + gx1 tan
7.5
hmx = h + h = h + tan = 1 + tan14 = 1.9375 m
2 2
7.5
hmn = h h = h tan = 1 tan14 = 0.065 m
2 2

pmx = p0 + g (hmn + tan )
2
EJERCICIO 7.5 Demuestre que para un flujo unidireccional (figura 7.4) que corre en un plano
inclinado con relacin a la vertical, la cabeza piezomtrica h es constante en cualquier punto
dentro del flujo (se considera que el cambio de alturas z es mucho mayor que la dimensin
en direccin del eje x3 ).
p
h= +z
g
- densidad
g - aceleracin gravitacional
h - cabeza piezomtrica
354
FIGURA 7.4 FLUJO UNIAXIAL EN DIRECCIN DEL EJE x1
v1 0 v2 = v3 = 0
p
+ B2 = 0 (7.4)
x2
p
+ B3 = 0
x3
Como z apunta hacia arriba (figura 7.4), las fuerzas de cuerpo por unidad de masa se
expresan como
B = 0er + 0e gez
Entonces la proyeccin del campo gravitacional en direccin del eje x2 est dada por
B2 = B e2 = g ( ez e2 )
El vector de posicin de cualquier elemento diferencial del fluido queda
r = ( x1e1 + x2 e2 + x3 e3 )
355
Entonces la proyeccin del vector de posicin con respecto del eje vertical es
z = ez i r = ( ( ez ie1 ) x1 + ( ez ie2 ) x2 + ( ez ie3 ) x3 )
Es evidente que
( ez ie3 ) = 0
Por otra parte la variacin de

ez i r = ez i e2
x2
Se puede escribir entonces que
z gz
B2 = g ( ez ie2 ) = g =
x2 x2 x2
Sustituyendo en la ecuacin 7.4
p
( gz ) =0
x 2 x2

( gz + p ) = 0
x2
Entonces, para todos los puntos de un mismo plano el cual es perpendicular a la direccin de
flujo se tiene
p
z+ =0
x2 g
Por lo tanto
p
z+ = cte
g
356
7.5 FLUIDO NEWTONIANO
Como ya ha sido mencionado, se define como fluido newtoniano al medio continuo que se
caracteriza porque la relacin de los esfuerzos de corte a la velocidad de deformacin es
lineal, por consecuencia el estado de esfuerzos se podr describir como
ij = p ij + ij (7.5)
donde el tensor ij depende de la velocidad de deformacin y representa la componente

viscosa del estado de esfuerzos; en un fluido newtoniano, en cualquier punto asociado al
medio continuo, el esfuerzo ij depende en forma lineal de las componentes del tensor de
rapidez de deformacin ( Dij )
ij = ( Dij )
1 vi v j
Dij = +
2 x j xi

No existe ninguna razn experimental por la cual se pueda considerar que las propiedades
del fluido se modifican con la posicin, as como que stas dependen de la direccin. De lo
antes mencionado, se describe a ste como un fluido homogneo e isotrpico. Por analoga
con un slido elstico, homogneo, lineal e isotrpico, se define un fluido viscoso,
homogneo, lineal e isotrpico (fluido newtoniano); por lo tanto
ij = Cijkl kl
Realizando las mismas consideraciones que fueron efectuadas para el slido elstico
homogneo e isotrpico, se tiene que:
El tensor de esfuerzos y el de rapidez de deformacin son simtricos
ij' = 'ji ; kl = lk
357
El tensor de constantes viscosas es simtrico
Cijkl = C jikl ; Cijkl = Cijlk
Cijkl = Cklij
Por otra parte, el tensor de constantes viscosas Cijkl es isotrpico, esto es no se modifica
bajo cualquier base, razn por la que el sistema se puede representar en la forma
ij = kk ij + 2 ij
lo cual representa que solo existen dos constantes viscosas linealmente independientes:
f t
, =
l2
Como a travs de la Segunda ley de Newton se encuentran relacionadas la fuerza, masa,

longitud y tiempo, se tiene que
f t ml t m
, = = =
2 2 l t
l 2
t l
donde representa la razn de proporcionalidad entre el esfuerzo de corte y la rapidez con

la que decrece el ngulo entre dos lneas materiales mutuamente perpendiculares. A esta
constante se le denomina como primer coeficiente de viscosidad o simplemente viscosidad.
Por su parte, no tiene un significado fsico, ambas en el Sistema Internacional tienen
unidades de pascal-segundo, unidades que para la mayora de las aplicaciones prcticas
g
resultan muy elevadas. Es por lo anterior que se define el Poise , en sistema ingls
cm s
lb f -s slug
las unidades son o .
pie 2 pie-s
En la mecnica de fluidos es muy comn el empleo de la viscosidad cinemtica
358
entidad que proviene del cociente de la viscosidad , con respecto a la densidad . La

l2
viscosidad cinemtica tiene unidades de . En sistema mtrico la unidad de la
t
cm 2 4 m
2
viscosidad cinemtica recibe el nombre de Stoke = 10 .
s s
Retomando la expresin a travs de la cual se define el estado de esfuerzos, el tensor total

de esfuerzos queda
ij = p ij + Dkk ij + 2 Dij
11 = p + Dkk + 2 D11
22 = p + Dkk + 2 D22
33 = p + Dkk + 2 D33
ii = 3 p + ( 3 + 2 ) Dkk
2
H = p + + Dkk
3
Donde
1 D
Dkk = ( dV )
dV Dt
representa la rapidez de cambio de volumen y est dada por
vi
Dkk = v = = 11 + 22 + 33
xi
o Dkk = D11 + D22 + D33
Si el fluido es incompresible Dkk = 0 , el esfuerzo hidrosttico estar dado por la presin, sin
embargo, en el caso de que Dkk 0 , la presin p representar solo parte de la presin
hidrosttica, siendo necesario definir el coeficiente de viscosidad volumtrica
2
K =+
3
359
2
Por otra parte, si el fluido es compresible y H = p , solamente si + = 0 , lo cual se
3
define como Condicin de Stokes.
Fluido newtoniano incompresible
Para el caso de un fluido newtoniano incompresible, el estado de esfuerzos se expresa como
ij = p ij + 2 Dij
Por lo que
kk = 3 p + 2 Dkk
Dkk = 0 kk = 3 p
H = p
Es por tanto que en un fluido viscoso incompresible la presin hidrosttica no depende de

ninguna cantidad cinemtica y es indeterminada con relacin al comportamiento mecnico
de ste; por otra parte, para un fluido viscoso incompresible se podr superponer cualquier
presin al fluido sin que esto afecte su comportamiento mecnico. Por consecuencia, la
presin resulta indeterminada desde el punto de vista de las ecuaciones constitutivas que
caracterizan a un fluido viscoso incompresible. Retomando la ecuacin constitutiva
ij = p ij + Dkk ij + 2 Dij y considerando que
1 vi v j
Dij = +
2 x j xi

se tiene entonces que

vi v j
ij = p ij + +
x j xi

360
v1
11 = p +
x1
v2
22 = p +
x2
v3
33 = p + 2
x3
v1 v2
12 = +
x2 x1
v2 v3
23 = +
3 x2
x
v3 v1
31 = +
x1 x3
Ecuaciones de Navier-Stokes para fluidos incompresibles
Sustituyendo el estado de esfuerzos definido para un fluido newtoniano incompresible en la

ecuacin de Cauchy, se tiene
ij
+ Bi = ai
x j
ij Dvi
+ Bi =
x j Dt
ij v v
+ Bi = i + v j i
x j t x j

v v j v v
p ij + i + + Bi = i + v j i
x j x j xi t x j

361
p 2 vi 2v j v v
+ + + Bi = i + v j i
xi x j x j xi x j t x j

p 2 vi v j v v
+ + + Bi = i + v j i
xi x j x j xi x j t x j

Para un fluido incompresible
v j p 2vi v v
=0 + + Bi = i + v j i
x j xi x j x j t x j

Lo cual en notacin general queda
v
p + v + B = + vv (7.6)
t
A la ecuacin anterior se le conoce como ecuacin de Navier-Stokes, la cual describe el

movimiento de fluidos newtonianos incompresibles. Para este sistema vectorial de
ecuaciones, las incgnitas estn representadas por el campo de velocidades y la presin
(v1 , v2 , v3 , p ) , por esta razn se requiere una cuarta ecuacin que es la de la continuidad
(Conservacin de masa)
v = 0
En coordenadas rectangulares las ecuaciones de Navier-Stokes quedan en la forma
p 2 2 2 v v v v
+ 2 + 2 + 2 v1 + B1 = 1 + 1 v1 + 1 v2 + 1 v3
x1 t x1 x2 x3
x1 x2 x3
p 2 2 2 v v v v
+ 2 + 2 + 2 v2 + B2 = 2 + 2 v1 + 2 v2 + 2 v3
x2 x t x1 x2 x3
1 x2 x3
p 2 2 2 v v v v
+ 2 + 2 + 2 v3 + B3 = 3 + 3 v1 + 3 v2 + 3 v3
x3 x t x1 x2 x3
1 x2 x3
362
EJERCICIO 7.7 Si todas las partculas de un fluido tienen sus vectores de velocidad paralelos
a una direccin fija, el flujo se definir como uniaxial.
Demuestre que para un flujo uniaxial de un fluido newtoniano incomprensible el esfuerzo

normal compresivo en cualquier punto de cualquier plano paralelo y perpendicular a la
direccin de flujo es p .
Considere que la direccin de flujo corresponde con el eje x3 , por lo que las velocidades en
direccin de x1 , x2 sern igual a cero
v3 0; v1 = v2 = 0
De la ecuacin de la continuidad
v1 v2 v3
v = 0 + + =0
x1 x2 x3
v3
=0 v3 = c
x3
v3 = v ( x1 , x2 ; t )
D11 = D22 = D33 = 0
ij = p ij + 2 Dij
11 = 22 = 33 = p
Ecuaciones de Navier-Stokes en coordenadas cilndricas ( r , , z )
Dado que en este caso las incgnitas son vr , v , vz , p se requerir emplear las ecuaciones
de Navier-Stokes y la Ecuacin de la continuidad, las cuales quedan
D
+ ( v ) = 0
Dt
363
Como el fluido es incompresible
D
=0 v = 0
Dt
Por lo que
vr 1 v v
v = + + vr + z
r r z
1 ( rvr ) 1 v vz
v = + + z
r r r
El gradiente del campo de velocidades queda
vr 1 vr vr
v
r r z
v 1 v v
v = + vr
r r z
v 1 vz vz
z
r r z
Por su parte, la ecuacin de Navier-Stokes en notacin general se expresa
v
p + ( v ) + B = + v ( v )
t
La divergencia de un tensor de segundo orden ( T )r , , z en coordenadas cilndricas es
T 1 Tr Trr T Trz
( divT )r = rr + r + +
r r z
T 1 T Tr + T r T z
( divT ) = r + r + +
r r z
T 1 Tz Tzz Tzr
( divT ) z = zr + + +
r r z r
donde, por analoga, se puede considerar que T = v
364
Sustituyendo lo antes expuesto en la ecuacin de Navier-Stokes se tiene
p 2v 1 2 vr 2 vr 1 vr 2 v vr
+ 2r + 2 + 2 + 2 2 + Br =
r r 2
r r r
r z r
v v v v v v 2
= r + vr r + r + vz r
t r r z r

1 p 2v 1 2v 2v 1 v 2 vr v
+ 2 + 2 + 2 + + + B =
r r 2 r r r 2 r 2
r z
v v v v v vv
= + vr + + vz + r
t r r z r
p 2vz 1 2vz 2vz 1 vz

+ 2 + 2 + + + Bz =
z r r 2 z 2 r r

v v v v v
= z + vr z + z + vz z
t r r z
Ecuaciones de Navier-Stokes en coordenadas esfricas ( r , , )
Dado que en este caso las incgnitas son vr , v , v , p se requerir usar las ecuaciones de
Navier Stokes y la Ecuacin de la continuidad, las cuales quedan:
D
+ v = 0
Dt
Como el fluido es incompresible
D
=0 v = 0
Dt
Sea v = v ( r , , ) una funcin vectorial, entonces
vr 1 v 1 v vr v cot
v = div v = + + vr + + +
r r r sen r r
365
1 r 2vr 1 v sen 1 v
v = + +
r 2 r r sen r sen
vr 1 vr 1 vr
v v sen
r r r sen

v 1 v 1 v
v = + vr v cos
r r r sen

v 1 v 1 v vr v cot
r + +
r r sen r r
Por su parte, dado que la Ecuacin de Navier-Stokes en notacin general se expresa
v
p + ( v ) + B = + v ( v )
t
La divergencia de un tensor de segundo orden ( T )r , , en coordenadas esfricas est

dada por
( div T )r = 2
(
1 r Trr
2
) +1
( Tr sen ) +
1 Tr T + T

r r r sen r sen r

( div T ) = 3
(
1 r T r
3
) + 1
( T sen ) +
1 T Tr T r T cot
+
r r r sen r sen r

( div T ) = 3
(
1 r T r
3
) + 1
(
T sen +
1 T Tr T r + T cot
) +
r r r sen r sen r

Sustituyendo lo antes expuesto en la ecuacin de Navier-Stokes se tiene
366
p
+
r
1 2 1 vr 1 2vr 2 (v sen ) 2 v
2 ( r vr
) + sen + +
r r r r 2 sen r 2 sen 2 2 r 2 sen r 2 sen
v vr v vr v vr v2 + v2
Br = r + vr + +
t r r r sen r

1 p
+
r
1 2 v 1 1 1 2v 2 vr 2 cot v
r + (
2 sen v sen ) r 2 sen 2 2 + r 2 r 2 sen +
+
r r r r
2

v v v v v v vr v v2 cot
B = + vr + + +
t r r r sen r r

1 p
+
r sen
1 v 1 1 2 v vr 2 cot v
2 r 2
+
r r r r sen
2
v sen ( 1
) +
2
r 2 sen 2 2 r 2 sen + r 2 sen +
+

v v v v v v v vr v v cot
B = + vr + + + +
t r r r sen r r
7.6 LNEAS DE TRAYECTORIA Y LNEAS DE CORRIENTE
Para la mecnica de fluidos es necesario disponer de una representacin visual de un

campo de flujo, tal representacin se puede realizar a travs de la descripcin de las
trayectorias de las partculas del fluido.
Una lnea de trayectoria est constituida por la curva que se traza a travs del movimiento
de una partcula, esto es, si se pretendiera determinarla fsicamente habra que tomar una
imagen del flujo (empleando por ejemplo tinta o colorante) durante un intervalo de tiempo,
367
obteniendo as la descripcin de la trayectoria. Matemticamente, la lnea de trayectoria de

una partcula la cual estaba en X para un tiempo t0 puede ser obtenida a partir del campo
de velocidad v ( x, t ) . Si x = x ( t ) es la lnea de trayectoria, entonces
dx
= v ( x, t )
dt
x ( t0 ) = X
EJERCICIO 7.8 Para el campo de velocidades descrito por
2 x1
v1 = ; v2 = 0; v3 = cx3
( 2t + 1)
Determine la lnea de trayectoria de una partcula x = x ( X i , t0 ) , esto para el tiempo de

referencia t0 .
A partir de la definicin de velocidad y considerando la relacin entre coordenadas

lagrangianas y eulerianas, se tiene que
dx1 2 x1 x1dx1 t 2dt

dt
= v1 =
( 2t + 1)
X1 x1 t0 ( 2t + 1)
=
x1 ( 2t + 1) ( 2t + 1)
Ln = Ln x1 = X 1
X1 ( 2t0 + 1) ( 2t0 + 1)
dx2
= v2 = 0 x2 = X
dt
dx3 x3 dx3 t
dt
= v3 = cx3 X 3 x3
= cdt
t0
x3 = X 3e ( 0 )
x3 c t t
Ln = c ( t t0 )
X3
368
Las lneas de corriente son lneas dibujadas en el campo de flujo, de tal manera que en un
instante dado se encuentran siempre tangentes a la direccin de flujo en cada punto del
campo que describe ste. En otras palabras, una lnea de corriente en un tiempo t es una
curva cuya tangente en cualquier punto tiene la direccin instantnea del vector velocidad de
la partcula en ese momento. Experimentalmente, estas lneas son obtenidas adicionando
partculas reflejantes en el fluido al efectuar una toma fotogrfica sobre la superficie del fluido
con una apertura corta del obturador. Cada elemento reflejante produce una lnea corta
sobre el flujo, la cual se aproxima a la tangente de la lnea de corriente. Matemticamente,
las lneas de corriente pueden ser obtenidas a partir del campo de velocidades v ( x, t ) .
Considere que x = x ( s ) representa a la ecuacin paramtrica para la lnea de corriente al

tiempo t , la cual pasa a travs de un punto dado x0 , entonces una s puede ser escogida tal
que
dx
= v ( x, t ) x ( 0 ) = x0
ds
El perfil de las lneas de corriente puede cambiar de un instante a otro si la velocidad del flujo
es funcin del tiempo (flujo no establecido). Dado que las lneas de corriente son tangentes
al vector velocidad de cada punto del flujo, el fluido nunca puede cruzar una lnea de
corriente.
EJERCICIO 7.9 Para el campo de velocidades descrito para el ejercicio 7.7, determine las
lneas de corriente que pasan en el punto p :
dxi 2 x1
p ( p1 , p2 , p3 ) = v ( x, t ) = e + 0e + cx e
ds ( 2t + 1) 1 2 3 3
2s
x1 dx1 s 2ds ( t +1)
p1 = 2
x1 = p1e
x1 0 ( 2t + 1)
v2 = 0 x2 = p2
dx3 x3 dx3 s
v3 = cx3 = = = cds
ds p3 x 0
3
x3 = p3ecs
369
Las lneas de trazo se definen al observar un punto en el espacio correspondiente al flujo,

identificando todas las partculas que pasan a travs de ste y que al unirlas dan lugar a
dicha lnea. Una lnea de trazo a travs de un punto fijo x0 es la lnea a un tiempo t, formada
por todas las partculas que han pasado a travs de x0 para < t . Considere que X = X (x,
t) denota la funcin inversa de x = x (X, t), entonces la partcula que estaba en x0 para un
tiempo tiene las coordenadas materiales dadas por X = X (x0,), por lo que la misma
partcula est en x = x ( X ( x0, ), t ) al tiempo t. Es entonces que la lnea de trazo a un
tiempo t est dada por x = x ( X ( x0, ), t ) para un tiempo fijo t y variable .
EJERCICIO 7.10 Para el campo de velocidades descrito en los ejemplos anteriores
2 x1
v1 = ; v2 = 0; v3 = cx3
( 2t + 1)
Determine la lnea de trazo formada por las partculas que pasan a travs de la posicin
espacial p ( p1 , p2 , p3 ) .
Las ecuaciones de trayectoria fueron determinadas en el ejercicio 7.7, quedando:
( 2t + 1)
x1 = X 1
( 2t0 + 1)
x2 = X 2
x3 = X 3 e ( 0 )
c t t
las cuales tienen las funciones inversas
( 2t0 + 1)
X1 = x1
( 2t + 1)
X 2 = x2
X 3 = x3e ( 0 )
c t t
370
Entonces, la partcula que pasa a travs de la coordenada p ( p1 , p2 , p3 ) a un tiempo est

dada por
( 2t0 + 1)
X1 = p1
( 2 + 1)
X 2 = p2
X 3 = p3e ( 0 )
c t
Sustituyendo en las ecuaciones de trayectoria, se obtiene la ecuacin paramtrica para las

lneas de trazo
( 2t + 1)
x1 = p1
( 2 + 1)
x2 = p2
x3 = p3 e ( )
c t
7.7 FLUJO ESTABLECIDO Y FLUJO TRANSITORIO
El flujo denominado como establecido es aquel en el cual para cualquier localizacin fsica
las condiciones no cambian en el tiempo. En un flujo establecido las lneas de corriente, de
trayectoria y de trazo no cambian con el tiempo. Un flujo transitorio por su parte se
caracteriza en que ya sea la velocidad, aceleracin o temperatura cambian con el tiempo.

Por ejemplo, sea una variable cualquiera en un flujo transitorio 0.
t
7.8 FLUJO LAMINAR Y FLUJO TURBULENTO
Los flujos viscosos se pueden clasificar en laminares o turbulentos, de acuerdo con sus
condiciones de movimiento. El flujo laminar se caracteriza en que el movimiento se da en
371
capas, de tal forma que no existe mezcla entre las lneas de corriente. Se trata de un flujo
muy ordenado en el cual las partculas que fluyen se desplazan formando finas capas. Este
tipo de flujo se relaciona fundamentalmente con condiciones de baja velocidad.
Experimentalmente en este caso, si se inyecta tinta en el seno de la corriente se tendr que
no existe mezcla, mantenindose en el tiempo, separado el marcador del resto de la
corriente. En un flujo turbulento se presentan movimientos tridimensionales aleatorios de las
partculas del fluido. Si se inyecta tinta a un flujo turbulento, se tendr que sta rpidamente
se dispersa en todo el campo del flujo.
Para el caso de agua fluyendo en un conducto de seccin circular, Reynolds observ que un
parmetro adimensional (conocido ahora como nmero de Reynolds) N R = v d 1 [v ,
velocidad promedio; d , dimetro del tubo; , densidad y viscosidad del fluido] describa
las condiciones del flujo, de tal forma que para valores menores del orden de 2100, el flujo se
presentaba en capas, y donde un filamento de colorante introducido en el fluido permaneca
intacto. En dicho caso, cualquier alteracin en el flujo es rpidamente eliminada. Si NR se
incrementa, el flujo se vuelve sensible a la presencia de pequeas perturbaciones, de tal
manera que se llega al punto en que el filamento de colorante introducido al flujo se rompe,
mezclndose con la totalidad del flujo. A este fenmeno se le denomina como flujo
turbulento. Las condiciones en las cuales se presenta un flujo turbulento son difciles de
precisar con exactitud, dependiendo entonces de la alteracin que se genere.
Dada la complejidad de su anlisis, en primera instancia se proceder a estudiar la dinmica

de los fluidos considerando flujos laminares, por consecuencia, los resultados obtenidos
estarn limitados a dichas condiciones y definidos por su nmero de Reynolds caracterstico.
7.9 FLUJO DE COUETTE
En este apartado se va a analizar un flujo uniaxial entre dos placas planas semi-infinitas
(vase la figura 7.5). El flujo se considera uniaxial y se presenta en direccin de x1 cuando la
placa superior se mueve a una velocidad v . Dado que el movimiento del fluido es uniaxial,
entonces v1 = f ( x2 ) , v2 = v3 = 0 . Se trata de un fluido Newtoniano incompresible por lo que
la ecuacin de la continuidad queda v = 0 ; y a partir de las ecuaciones de la continuidad y
de Navier-Stokes, se tiene que:
372
2 v1
=0 v1 = ax2 + b
x22
Resulta por dems evidente que los valores de las constantes a y b estarn determinados
por la posicin en que se coloque el sistema de referencia; si el origen se coloca en la placa
inferior, tal como se muestra en la figura, se tiene que v1 = 0 para x2 = 0 y v1 = v para
x2 = d , es entonces que
vx2
v1 ( x2 ) =
d
vx2
v ( x2 ) =
d
FIGURA 7.5 FLUJO UNIAXIAL ENTRE PLACAS PLANAS SEMI-INFINITAS

INDUCIDO POR VELOCIDAD ( FLUJO DE COUETTE)
7.10 FLUJO UNIAXIAL PRODUCIDO POR PRESIN (FLUJO DE POISEUILLE)
El flujo denominado como de Poseuille se describe como un flujo uniaxial inducido por
presin, el cual se produce entre dos placas planas semi-infinitas (figura 7.6), donde la
presin p = p ( x1 ) (figura 7.7).
De acuerdo con lo observado, se tiene que:
v1 = f ( x2 )
v2 = v3 = 0
v1 (0) = vmx
v1 (a ) = v1 ( a) = 0
373
FIGURA 7.6 FLUJO UNIAXIAL ENTRE PLACAS SEMI-INFINITAS, INDUCIDO POR PRESIN
Ecuacin de Navier-Stokes considerando flujo establecido:
p 2 v1
+ 2 =0
x1 x2
p
=0
x2
p
=0
x3

FIGURA 7.7 PERFIL DE PRESIN
No se consideran fuerzas de cuerpo B1 = B2 = B3 = 0
p
=
x1
2 v1
=
x22
374
v1 2 2
= x2 + c1 0= a + 0 + 2c2 c2 = a
x2 2
2
v1 = x2 + c1 x2 + c2 0 = 0 + 2c1a + 0 c1 = 0
2
v1
Para la condicin de frontera en x2 = 0, v1 = vmx =0
x2
0 = 0 + 2c1a + 0 c1 = 0
y en x2 = a, v1 = 0
2 2
0= a + 0 + 2c2 c2 = a
2
2
v1 (a) = 0 = a + c1a + c2
2
2
v1 (a) = 0 = a c1a + c2
2
De lo anterior se tiene que
2 2
v1 ( x2 ) = x2 a
2 2
2 2 2 2
v1 ( x2 ) =
2
(x2 a ) v1 ( x2 ) =
2
(
a x2 )
1 2 2 p
v1 ( x2 ) =
2
(
a x2 )
x1
p
De lo cual se concluye que la velocidad del flujo v1 ( x2 ) es proporcional a y se presenta
x1
en sentido contrario a la direccin de crecimiento del gradiente; por otra parte, la velocidad
375
mxima del flujo se presenta en el plano medio entre ambas placas semi-infinitas y est
a 2 p
dada por; vmx = , por lo que el gasto volumtrico V se puede describir como
2 x1
a 2a3
V = wv1 ( x2 )dx2 = w
a 3
donde w representa el ancho de las placas, y la velocidad promedio v se describe a travs de
1 a a2
a 0
v= v ( x ) dx =
1 2 2
3
EJERCICIO 7.11 Para el mismo flujo uniaxial descrito en la figura 7.4, determine la presin en
el punto A .
pA p
+ z A = B + zB
g g
FIGURA 7.8 CAMPO DE VELOCIDADES DEL FLUJO EN UN CANAL INCLINADO
De la figura queda
p A = patm + g ( z B z A )
p A = patm + gh cos
donde patm representa a la presin atmosfrica.
376
Por otra parte, no existe restriccin al movimiento del fluido en la superficie superior, por lo
que su velocidad ser mxima. Cmo son las velocidades para el fondo del canal inclinado
y en la superficie libre del flujo?
El fondo del canal no presenta movimiento por lo que la velocidad en este plano es nula; por
Cmo ser el perfil de velocidades?
Dado que las paredes del canal estn inmviles, lo que genera el movimiento del fluido es la
diferencia de presin debida a la inclinacin del canal, razn por la que se trata de un flujo de
Poiseuille (inducido por presin) por lo que el perfil de velocidades es parablico, tal como
se observa en la figura 7.8.
7.11 FLUJO INDUCIDO POR PRESIN A TRAVS DE UN CONDUCTO DE

SECCIN CIRCULAR (TUBO)
Este flujo tambin denominado como de Hagen-Poiseuille, se describe como un flujo

uniaxial, axisimtrico inducido por presin. Resulta evidente que el anlisis se facilitar si
ste se realiza en coordenadas cilndricas. De lo antes expuesto se tiene que:
vr = 0, v = 0, vz = vz (r )
vz
v(0) = vmx =0
r
v( R) = 0
FIGURA 7.9 PERFIL DE VELOCIDADES PARA UN FLUJO DE HAGEN-POISEUILLE.

LA PRESIN ES FUNCIN LINEAL DE z
377
Ecuacin de la continuidad, considerando que se trata de un fluido incompresible
D
+ v = 0 v = 0
Dt
Resulta entonces que la divergencia del campo de velocidades se expresa
vr 1 v v
v = + + vr + z
r r z
De lo antes expuesto se concluye que
vz
vr = v = 0; =0
z
En virtud de que el gradiente del campo de velocidades descrito en coordenadas cilndricas
est dado por
vr 1 vr vr
v
r r z
v 1 v v
v = + vr
r r z
v 1 vz vz
z
r r z
Considerando la ecuacin de Navier Stokes correspondiente al eje z
p 2v 1 2 vz 2 vz 1 vz v v v v v
+ 2z + 2
r
+ 2 + + Bz = z + vr z + z + vz z
z r 2
z r r t r r z
sta se simplifica al ser un flujo establecido ya que
vz v v v v
+ vr z + z + vz z = 0
t r r z
Adems de que vz = v ( r ) .
Por consecuencia la ecuacin de Navier-Stokes para el eje z se reduce a
p 2 vz 1 vz
+ 2 + = 0
z r r r

378
La ecuacin se puede simplificar considerando
2vz 1 vz 1 vz 1 r vz 2vz
+ = r = +
r r r r r
r
r 2 r r r r 2
Por lo que la ecuacin se presenta como
p vz
+ =0
z r r r
r
El gradiente de presin se describe
p
=
z
Entonces, despejando
vz r
=
r r
r

e integrando
vz r2
r = + c1
r 2
v z r c1
= +
r 2 r
r2
vz = + c1 Ln r + c2
4
De las condiciones de frontera, tenemos
v z
Para r = 0 = 0, dado que v z = vmx
r
r 2 2 c1
0= + 0 = r 2 2 c1
2 r 2 r
c1 = 0
r2
vz ( r ) = + c2
4
379
vz = 0 para r =R
R2
c2 =
4
de lo cual se concluye que

vz ( r ) =
4
(
R2 r 2 =
1
4
)
R2 r 2 (
p
z
)
La mxima velocidad se presenta en el centro del tubo r = 0 , por lo que est dada por
d 2 p
vz mx =
16 z
donde d es el dimetro del tubo.
La velocidad promedio queda
4 d 2 p vz mx
v=
d2
vz dA = =
32 z 2
A
Por lo que el gasto volumtrico se describe por V = vA
d2 d 4 p
V= v=
4 128 z
EJERCICIO 7.11 Anlisis del flujo inducido por velocidad (Coutte) para dos capas de fluidos
newtonianos incompresibles, cuyo movimiento corresponde a un flujo laminar en rgimen
permanente.
Se tienen dos capas, con espesores d1 , d 2 (figura 7.10) que corresponden a fluidos
newtonianos incompresibles con viscosidad 1 , 2 , los cuales se encuentran entre dos
placas paralelas y semi-infinitas. Considerando que las placas se desplazan entre s,

determine el perfil de velocidades correspondiente. Para el anlisis desprecie el efecto de las
fuerzas de cuerpo.
380
Para ambas capas
v1 = f ( x2 )
v2 = v3 = 0
p = 0
Considere que el origen del sistema coordenado se establece en la interfase entre ambas
capas.
FIGURA 7.10 MOVIMIENTO UNIAXIAL INDUCIDO POR VELOCIDAD EN DOS

CAPAS DE FLUIDOS NEWTONIANOS INCOMPRESIBLES
Para la capa 1, las ecuaciones de Navier-Stokes se reducen a
2v11 v11
1 =0 = cte
x22 x2
v11 = ax2 + b
para x2 = d1 v11 = 0
0 = ad1 + b b = ad1
v11 ( x2 ) = ax2 + ad1
Realizando el mismo anlisis para la capa 2
2v12
2 =0 v12 ( x2 ) = cx2 + f
x22
381
Para x2 = d 2 v12 = v p (velocidad de la placa)
Analizando en la interfase ( x2 = 0) las velocidades debern ser iguales (no existe

deslizamiento):
v12 (0) = v11 (0) 0
Para la capa 2 en contacto con la pared mvil
v12 (d2 ) = v p = cd2 + f f = v p cd2
( )
v12 ( x2 ) = v p cd2 + cx2 = c( x2 d2 ) + v p
v12 (0) = v11 (0) c (0 d 2 ) + v p = a ( 0 + d1 )
v p = cd2 + ad1
En virtud de existir equilibrio de fuerzas en la interfase, los esfuerzos de corte definidos en el

seno de cada fluido (evaluados estos en la interfase) sern iguales:
12
1
= 12
2
1 v1 v12
12
1
= 2 1D12
1
= 2 1 1 +
x1
2 x2
D - velocidad de deformacin
12
2
= 22 D12
2
v11
12
1
= 1 ; v11 = ax2 + ad1
x2
v11 12
1
= a; 12
1
= 1a a=
x2 1
382
v12 d1
12
2
= 2 ; v12 = cx2 +
x2 a
v12 12
2
= c; c=
x2 2
En la interfase los esfuerzos de corte son iguales, por lo que:
12
1
= 12
2
= a 1 = c2
1
c= a
2
Para la capa 1
v11 = a( x2 + d1 )
Para la capa 2
1
v12 = a( x2 d 2 ) + v
2
Como ya se mencion, en la interfase las velocidades son iguales ya que se considera que
no hay deslizamiento entre los fluidos
1
v11 = v12 = a( x2 + d1 ) = a( x2 d 2 ) + v
2
1 vp
ad1 = ad + v a d1 + 1 d 2 = v p a=
2 2 p 2 1
d1 + d
2 2
v p 2
v11 ( x2 ) = ( x2 + d1 )
2d1 + 1d 2
1 v p 2
y v12 ( x2 ) = (x d ) + vp
2 2 d1 + 1d 2 2 1
v p 1
v12 = ( x2 d1 ) + v p
2 d1 + 1d 2
383
7.12 FLUJO INDUCIDO POR VELOCIDAD ENTRE DOS CILINDROS

CON LONGITUD INFINITA
Al flujo laminar bidimensional, en estado estable, de un fluido newtoniano incompresible

entre cilindros concntricos de longitud infinita se le denomina como flujo de Couette para
una base en coordenadas cilndricas. Para el anlisis de este caso considere dos cilindros
concntricos, de radios r1 , r2 y altura infinita h , tal como se muestra en la figura 7.11; stos
giran a velocidades angulares constantes ( 1 , 2) . De un primer anlisis resulta evidente
que el flujo se puede describir como uniaxial e inducido por velocidad, de tal forma que si el
origen se define en el centro del arreglo se tendr que
vr = vz = 0; v = v ( r )
FIGURA 7.11 MOVIMIENTO DE UN FLUIDO NEWTONIANO INCOM-

PRESIBLE ENTRE DOS CILINDROS COAXIALES
Al tratarse de un fluido incompresible, la ecuacin de la continuidad se expresa como
vr 1 v v
v = + + vr + z = 0
r r z
De lo que se tiene que
v
= 0,

384
situacin que resulta evidente al tratarse de un flujo establecido; por otra parte, de las
ecuaciones de Navier-Stokes se tiene para el eje
1 p 2v 1 2 v 2 v 1 v 2 vr v
+ 2 + 2 + + + 2 + B =
r r 2
2 r r 2
r
r z r
v v v v v vv
+ vr + + vz + r
t r r z r
El sistema se simplifica en virtud de que se realizan las siguientes consideraciones:
a. Flujo establecido
b. Se desprecia el efecto de las fuerzas de cuerpo
c. El flujo es inducido por velocidad
d. El gradiente de presin es igual a cero
Es por consecuencia que se reduce a
2 v 1 v v
+ =0
r 2 r r r 2
La solucin de la ecuacin diferencial est dada por
B
v ( r ) = Ar +
r
Por lo que
v B
= A 2
r r
2 v 2B
=
r 2
r3
Sustituyendo en la ecuacin diferencial se comprueba la validez de la solucin planteada
2B 1 B 1 B
+ A 2 2 Ar + = 0
r r
3
r r r
2B A B A B
3
+ =0
r r r3 r r3
385
Demostrada la solucin se tiene que las velocidades en cada cilindro concntrico estn
dadas por
B
r = r1 , v ( r1 ) = 1r1 v ( r1 ) = 1r1 = Ar1 +
r1
B
r = r2 , v ( r2 ) = 2 r2 v ( r2 ) = 2 r2 = Ar2 +
r2
Por lo que, considerando las condiciones de frontera, se llega a que las constantes de
integracin estn dadas por
w2 r22 w1r12 r12 r22 ( w1 w2 )

A= y B=
r22 r12 r22 r12
Por otra parte, a travs de lo expuesto se pueden calcular fcilmente los esfuerzos de corte
en la pared de cualquiera de los cilindros y, por consecuencia, es posible calcular el
momento transmitido para hacer girar stos. Se concluye, por ltimo, que dicho momento es
funcin de la viscosidad del fluido y de las caractersticas geomtricas y dimensionales del
arreglo, por lo que el sistema en anlisis puede ser empleado fcilmente para determinar la
viscosidad del fluido.
La velocidad de deformacin en coordenadas cilndricas se expresa como
vr 1 1 vr v v 1 vr vz
+ +
r 2 r r r 2 z r
1 v 1 v 1 vz
= D r + vr +
2 z r
Dr z
r
vz
Dzr Dz
z
Como para el caso en estudio
vr = 0, vz = 0, v = v ( r )
386
Para el caso analizado, el tensor de velocidad de deformacin se reduce a
1 v v
0 r r 0
2
1 v v
Dr z = 0 0
2 r r
0 0 0

Dado que el campo de velocidades se expresa como
B
v ( r ) = Ar +
r
de lo que se tiene
v B
= A 2
r r
Por consecuencia,
1 B 1 B 2B B
Dr = A 2 Ar + = 2 = 2
2 r r r 2r r
B
Dr =
r2
Por otra parte, los esfuerzos viscosos se expresan
r z = ( v ) I + 2 D
Entonces, el tensor de esfuerzos queda
v v
0 0
r r
v v
'r z = 0 0
r r
0 0 0

387
Por lo que
B
' r = 'r = 2
r2
Considerando, por ejemplo, el cilindro exterior, el esfuerzo de corte que existe entre ste y el
fluido viscoso est dado por
B
r = 2
r22
La fuerza aplicada estar dada por la integral de dicho esfuerzo a travs del rea de
contacto, de tal forma que
2 z 2 B
f = r dA = 0 rd dz
A 0 r22
B
f = 4 z
r2
Por consecuencia, el par aplicado al sistema para producir el movimiento ser
B
T = fd = 4 r2 z
r2
T = 4 Bz
donde ste queda en funcin de los parmetros geomtricos a travs de los cuales se define
la constante B , por consecuencia se puede determinar la viscosidad a travs de evaluar el
par aplicado
T
=
4 Bz
7.13 FLUJO ROTACIONAL E IRROTACIONAL
Durante su movimiento los fluidos pueden describir un flujo rotacional, donde, la rapidez de
rotacin de una partcula de un fluido se define como la velocidad angular promedio de
dos segmentos de lnea, los cuales son mutuamente perpendiculares.
388
Dado que el gradiente de velocidad v es un tensor de segundo orden, entonces ste se

puede descomponer en su parte simtrica, la que describe al tensor de rapidez de
deformacin D , y su componente no simtrica, la cual representa la rapidez de rotacin
v =
1
2
( 1
) (
v + ( v ) + v ( v )
T
2
T
)
1
2
(
v + ( v ) = D
T
)
1
2
(
v ( v ) =
T
)
A partir de lo anterior se define el vector de vorticidad que, en notacin ndice y con auxilio
del permutador, se representa de la forma siguiente:
i = 2 ijk jk
i = 2 (23e1 + 31e2 + 12 e3 )
v3 v2 v1 v3 v2 v1
i = e1 + e2 + e3
x2 x3 x3 x1 x1 x2
Por tanto, el vector de vorticidad est dado por el rotacional del campo de velocidades
(v) . Para un sistema coordenado rectangular (cartesiano) se tiene que el rotacional del
campo de velocidades se describe como
e1 e2 e3

i = Rot vi = v =
x1 x2 x3

v1 v2 v3
v v v1 v3 v2 v1
= 3 2 e1 + e2 + e3
x2 x3 x3 x1 x1 x2
389
Por su parte, a partir de que el rotacional est dado por la componente antisimtrica del
gradiente de velocidades se tiene que en coordenadas cilndricas se expresa como
1 v v v v v 1 v v
r , , z = z er + r z e + r + ez
r z z r r r r
Mientras que, en coordenadas esfricas
1 v v cot 1 v 1 vr v v
r , , z = + er + e

r r r sen r sen r r
v v 1 v
+ + r e
r r r
Flujo irrotacional
Se caracteriza porque v = 0 , esto es, las partculas que constituyen el flujo se mueven
sin presentar rotacin, es por consecuencia que la condicin antes expresada se deber
cumplir, respectivamente, en cualquiera que sea la base, cilndricas y esfricas. Por tanto, se
tiene que
v3 v2 v1 v3 v2 v1
= = =0

2x x
3 3x x
1 1x x2
1 vz v vr vz v 1 vr v
z = z r = r r + r =0
r
1 v v cot 1 v 1 vr v v v v 1 vr
r + = = + =0
r r sen r sen r r r r r
Es entonces que el vector de vorticidad asociado a un flujo irrotacional deber ser nulo
i = 0
390
Los flujos irrotacionales se caracterizan en virtud de que pueden ser descritos con auxilio de
( )
una funcin escalar = ( x1, x2 , x3 ) , tal que
v =

vi = ei
xi

v1 = e1; v2 = e2 ; v3 = e3
x1 x2 x3
El vector de vorticidad (para coordenadas rectangulares) queda

i =
1

e +
e +
2 x x x x 3
e
x2 x3 x3 x2 x3 x1 x1 x3 1 2 2 1
Por lo tanto, siempre que el flujo se pueda representar por una funcin escalar tal que

vi = ei , entonces ser un flujo irrotacional, donde la funcin ( x1 , x2 , x3 , t ) define un
xi
flujo irrotacional.
Considerando la ecuacin de la continuidad para un fluido incompresible e irrotacional, sta

se expresa como:

v = 0 =0
xi xi
2 = 0
2 2 2
+ + =0
x12 x22 x32
391
Estado de esfuerzos para un flujo irrotacional de un fluido

incompresible de densidad homognea
El estado de esfuerzos en un fluido Newtoniano se expresa como
= pI +
ij = p ij + ij
ij = kk ij + 2 ij
donde representa a los esfuerzos viscosos, si el fluido es incompresible kk = 0 , por lo

que el estado esfuerzos para un fluido newtoniano incompresible queda:
ij = p ij + 2 ij
1 v v j
ij = p ij + 2 i
+ ij
2 x j xi

Si se considera que el flujo es irrotacional, entonces es posible describir el campo de

velocidades a travs de una funcin escalar v = , tal como se explic anteriormente,
entonces sustituyendo ahora en la ecuacin de Navier-Stokes
v
p + ( v ) + B = + v ( v )
t
se tiene que
p 2 vi v
+ + Bi = + vj i
xi x j x j xi t x j
p 2 v v
+ + Bi = i + v j i
xi xi x 2j x j
t
392
Dado que el fluido es incompresible, la ecuacin de la continuidad se expresa
2
= 0; o 2 = 0
xi2
De lo que se desprende que el trmino que representa el efecto de la viscosidad del fluido es
igual a cero, por lo que la ecuacin se reduce a
p v v
+ Bi = i + v j i
xi
t x j
En notacin general queda
v
p + B = + v ( v )
t
Dicha ecuacin corresponde con la ecuacin de Euler, la cual describe el comportamiento de

un fluido viscoso e incompresible.
7.14 FUNCIONES DISIPATIVAS EN FLUIDOS NEWTONIANOS
DW
La rapidez de trabajo (potencia) = P desarrollada por las fuerzas (de cuerpo y de
Dt
superficie) aplicadas sobre el MC se ha demostrado que se puede descomponer en:
DW D
= ( EC ) + PS dV
Dt Dt
vi
PS = ij
x j
donde el trmino PS representa la potencia desarrollada para cambiar la forma y

dimensiones de la partcula de volumen diferencial (elemento diferencial de volumen), y
D
( EC ) representa la rapidez de cambio de la energa cintica asociada al elemento.
Dt
393
En el caso de un fluido incompresible, el trmino PS se puede expresar
ij = p ij + ij'
vi v vi
ij = p i + ij'
x j xi x j
Dado que se trata de un fluido incompresible
vi
=0
xi
vi v
ij = ij' i
x j x j
ij' = Dii ij + 2 Dij
Antes de continuar es conveniente recordar que para describir la velocidad de deformacin

se emplea en este texto
Dij ij
Como se trata de un fluido incompresible
Dii = 0
vi v
ij' = 2 Dij i
x j x j
Por otra parte desarrollando, se tiene que
vi v v v v v v
Dij = D11 1 + D12 1 + D13 1 + D21 2 + D22 2 + D23 2
x j x1 x2 x3 x1 x2 x3
v3 v v
+ D31 + D32 3 + D33 3
x1 x2 x3
1 vi v j
Dij = +
2 x j xi

394
Reagrupando y considerando que Dij = D ji (el tensor de rapidez de deformacin es

simtrico)
vi v 1 v v 1 v1 v3
Dij = D11 1 + D12 1 + 2 + D13 +
x j x1 2 x2 x1
2 x3 x1
1 v v v2 1 v v
+ D21 2 + 1 + D22 + D23 2 + 3
2 x1 x2 x2 2 x3 x2
1 v v 1 v v v3
+ D31 3 + 1 + D32 3 + 2 + D33
2 x1 x3 2 x2 x3 x3
vi
Dij = Dij Dij
x j
vi vi
ij = ij' = 2 Dij Dij
x j x j
Por lo que
(
PS = 2 D11
2
+ D22
2
+ D33
2
+ 2 D12
2
+ 2 D23
2
+ 2 D31
2
)
El trmino PS representa el trabajo realizado por unidad de volumen y por unidad de tiempo
(el cual siempre ser positivo), para cambiar la forma del elemento diferencial. Esta parte del
trabajo realizado se acumula con el tiempo en funcin de la variacin de la velocidad de
deformacin Dij , entonces, la funcin indi s (funcin de disipacin) representa la velocidad
con la que el trabajo se convierte en calor para un fluido newtoniano incompresible.
indis = 2 D11
2
(+ D22
2
+ D33
2
+ 2 D12
2
+ 2 D23
2
+ 2 D31
2
)
= 2 Dij Dij
395
Funcin disipativa para un fluido newtoniano compresible
Como
ij = p ij + ij + 2 Dij
= v
La rapidez de trabajo desarrollado por los esfuerzos queda
vi
ij = p + 2 + indis = p +
x j
donde = ( D11 + D22 + D33 ) 2 + indis
Lo anterior debido a que
vi v v v v
ij = p i ij + i i ij + 2 Dij i
x j x j xi x j x j
vi v v v
ij = p i + i i + 2 Dij Dij
x j xi xi xi
vi
ij = p + 2 + 2 Dij Dij
x j
donde discom es la funcin disipativa para fluido compresible, la cual se puede representar
tambin como
discom = ( D11
2
+ D22
2
+ D33
2
) + indis
indis = 2 ( D11
2
+ D22
2
+ D33
2
+ 2D12
2
+ 2D23
2
+ 2D31
2
)
2 2
discom = + ( D11 + D22 + D33 )2 + ( D11 D22 )2 + ( D33 D11 )2 + ( D22 D33 )2
3 3
+ 4 ( D12
2
+ D13
2
+ D23
2
)
396
Resulta evidente que la funcin de disipacin (rapidez con la que el trabajo se transforma en
calor) asociada a un fluido compresible ser siempre mayor que la correspondiente para un
fluido incompresible.
7.15 DIFUSIVIDAD TRMICA
En el caso de difusin de una especie en un solvente se ha demostrado experimentalmente

que existe una relacin entre la concentracin C del soluto en funcin de la posicin x , y del
tiempo, a dicha relacin se le conoce como Segunda Ley de Fick y se expresa como
C C
= D
t xi xi
A partir de considerar que el coeficiente de difusin D no depende de la posicin, esta

ecuacin se simplifica a
C 2C C
=D 2 o = D 2C
t xi t
De manera anloga se define una ecuacin diferencial en derivadas parciales que relaciona
la temperatura con la posicin y el tiempo, empleando para esto un parmetro definido como
difusividad trmica . Para su determinacin considrese la ecuacin de conservacin de
energa, la cual en notacin general se expresa
Du
T v + q qs = 0
Dt
donde representa la densidad, mientras que u es la energa interna; T o es el tensor
de esfuerzos, q es el calor por conduccin y qs o h representa la rapidez de generacin de
calor en el interior del medio continuo. Si solamente existe un flujo de calor por conduccin,
el cual sigue la Ley de Fourier q = , donde representa la temperatura, es la
397
constante de conductividad trmica y la rapidez de generacin de calor al interior del medio

continuo es igual a cero; a partir de lo anterior se tiene que
Du v 2
ij i =0
Dt x j xi xi
En notacin general se tiene que
Du
T v 2
Dt
Por otra parte, considerando la funcin de disipacin para un fluido incompresible indis y
que la energa interna es una funcin de la temperatura u = c , donde c es el calor
especfico, entonces
D
c indis 2 = 0
Dt
donde la funcin de disipacin representa la energa disipada como calor por efecto de las
fuerzas viscosas. Existen, por otra parte, una gran cantidad de condiciones en las cuales la
energa disipada es muy pequea cuando se compara con la magnitud del flujo de calor,
razn por la cual la ecuacin de simplifica a
D
c 2 = 0
Dt
Reordenando la ecuacin, sta se puede representar como
D 2
=0
Dt c

Sustituyendo = (difusividad trmica), por lo que
c
D
= 2
Dt
Ecuacin similar a la de Fick, cuya solucin, permite determinar la temperatura en funcin

de posicin y tiempo. La difusividad trmica es una nueva constante de proporcionalidad a la
constante de conductividad trmica e inversamente proporcional al producto de la densidad
por el calor especfico (energa interna por unidad de volumen).
398
7.16 FLUJO IRROTACIONAL DE UN FLUIDO NO VISCOSO DE DENSIDAD HOMOGNEA
Como ha sido mencionado anteriormente para el caso de un fluido no viscoso, el estado de

esfuerzos Tij se puede representar como Tij = p ij , de lo que la ecuacin de conservacin
de cantidad de movimiento se expresa como
p v v
+ Bi = i + v j i
xi
t x j
O, en notacin general
v
p + B = + vv
t
descripcin que es conocida como Ecuacin de Movimiento de Euler.
Ecuacin de Bernoulli
Cuando se estudia un fluido no viscoso, incompresible y cuya densidad es homognea, ser

siempre posible que se presente un flujo irrotacional. Considere que los campos son
conservativos, de tal forma que el vector de fuerzas de cuerpo se puede desarrollar a partir
de una funcin de potencial , tal que; B = . El caso ms comn es aquel en que el
medio continuo se encuentra exclusivamente bajo el efecto de la fuerza de gravedad,
considerando que el eje x3 coincide con la vertical se tiene que = gx3 , por lo que
Bi = 0e1 + 0e2 ge3 , por otra parte al tratarse un flujo irrotacional, v = ( v ) .
T
x3
x2
x1
FIGURA 7.12 PARA LA DEFINICIN DE LAS FUERZAS DE CUERPO

EL EJE x3 SE HACE COINCIDIR CON LA VERTICAL
399
Sustituyendo en la ecuacin de Euler se tiene
p v v j
= i + vj
xi xi t xi
Por lo que, reordenando
p v j v
vj = i
xi xi xi t
Como
v j 1 v j v j 1 2 1
vj = = vj = v jv j
xi 2 xi 2 xi 2 xi
Adems, para un flujo irrotacional el campo de velocidades se puede expresar a travs de

una funcin , tal que

vi =
xi
Es entonces que
1 2
p + + v =0
xi 2 t
Por lo tanto, se puede concluir que
1
p + + v2 =k
2 t
donde k = constante con la posicin
1
p + + v2 = f ( t )
2

Si el flujo es establecido (estacionario) = cte , por lo que
t
p v2
+ + = cte
2
400
Lo que entonces permite deducir la ecuacin de Bernoulli:
v2
p + gx3 + = C
2
Siendo C una constante.
Esta expresin es muy til en aquellos casos en que la viscosidad pueda ser despreciada o
tambin cuando, como se demostrar ms adelante, el flujo se idealice como irrotacional.
Ecuacin de Torricelli
Considere ahora un fluido incompresible y no viscoso. De acuerdo con lo expuesto, si el flujo

es estable se cumplir con la ecuacin de Bernoulli:
v2
p + + gx3 = C

2
Ahora bien, si el nivel del recipiente (figura 7.13) permanece constante y la presin sobre el
nivel del fluido es la misma que la de descarga, se cumplir que
v12 v22
p1 + + gh1 = p2 + + gh2
2 2
donde la condicin 1 se refiere al plano de entrada (superficie del lquido en el tanque),

mientras que la 2 se refiere a la descarga; adems de lo anterior, se considera que el nivel
del fluido en el tanque es constante por lo que v1 = 0 , por consecuencia
p1 = p2 = patm presin atmosfrica
v1 = 0 h1 h2 = h diferencia de nivel
v2 = 2 g h
401
La cual es conocida como Ecuacin de Torricelli.
FIGURA 7.13 UN TANQUE CUYO NIVEL DE FLUIDO, CON RESPECTO A LA

POSICIN DE DESCARGA, ES h PRESENTAR UNA VELOCIDAD DE
SALIDA, LA CUAL EST DADA POR v = 2 gh
Flujos irrotacionales como solucin a la ecuacin de Navier-Stokes
Considrese en primera instancia la ecuacin de Navier Stokes para un fluido incompresible
v
p + ( v ) + B = + vv
t
En notacin ndice, sta queda
p 2 vi v v
+ + Bi = i + v j i
xi x j x j t x j

El campo de velocidades para un flujo irrotacional se puede describir a partir de una funcin
escalar , tal que

= vi
xi
entonces, sustituyendo
2vi 2 2
= =
x j x j x j x j xi xi x j x j

402
Pero se ha demostrado que la ecuacin de la continuidad (fluido incompresible) se expresa

en la forma
2
= =0
2

x j x j

Por lo tanto, para un flujo irrotacional de un fluido Newtoniano incompresible se tiene
p 2 v v
+ + Bi = i + v j i
xi xi x j x j

t
x j
Por lo que el trmino que representa los esfuerzos viscosos se hace cero, entonces la
ecuacin constitutiva se reduce a
p v v
+ Bi = i + v j i
xi t x j

la cual corresponde con la ecuacin de Euler (para un fluido no viscoso).
Por lo tanto, todos los resultados desarrollados para flujos no viscosos corresponden
tambin al caso de flujos irrotacionales. Sin embargo, en todo problema real existirn
fronteras fsicas en las cuales la velocidad del fluido ser de cero (o la velocidad de la
frontera) en virtud de que el fluido se adhiere a sta. Es por consecuencia que la condicin
v = no se podr cumplir en las condiciones de frontera.
7.17 ECUACIN DE TRANSPORTE DE VORTICIDAD PARA UN FLUIDO VISCOSO

INCOMPRESIBLE DE DENSIDAD HOMOGNEA
La imposibilidad de que exista una funcin , la cual se cumple para las paredes (frontera)
que confinan el movimiento del fluido, da lugar a la existencia de vorticidad confinada a una
capa adyacente a la frontera.
403
De nueva cuenta, retomando la ecuacin de Navier-Stokes
1 p 2vi v v
+ + B = i +vj i
xi x j x j
t x j
Asimismo, sustituyendo de tal forma que = / (viscosidad cinemtica) y que

B = , = gx3
p 2vi v v
+ + = i +vj i
xi x j x j t x j

Empleando el operador mni (tomando rotacional de ambos lados de la ecuacin
xn
anterior):
vi vi
mni t =
t mni = m
xn xn t
vi v j vi v
mni vj = mni + vj mni i
xn x j xn x j x j xn
v j vi
= mni + vj m
xn x j x j
2 p p
mni + = 0 + = 0
xn xi
2 vi 2 vi 2 m
mni = mni =
xn
x j x j x j x j xn x j x j
La ecuacin de Navier Stokes toma entonces la siguiente forma
m v j vi 2 m
+ v j m + mni =
t x j xn x j x j x j
404
Por otra parte, se puede demostrar que el tercer trmino del lado izquierdo de la ecuacin es
igual a cero, por lo que
v j vi
mni =0
xn x j
D m vm 2 m
= n +
Dt xn x j x j
Lo cual en forma invariante se expresa
D
= ( v ) +2
Dt
En una frontera slida el fluido se adhiere, lo que da lugar a que las velocidades en esta
interfase estn definidas por la superficie. Los vrtices son generados en la superficie
difundindose al flujo. En algunos casos los vrtices quedan confinados a una capa delgada
en la vecindad de la frontera, de tal forma que fuera de esta capa el flujo es irrotacional.
7.18 EL CONCEPTO DE CAPA LMITE
De lo que ha sido discutido con antelacin se ha comprobado que las funciones que
describen el comportamiento en un fluido viscoso y no viscoso son iguales, sin embargo,
debido a la presencia de esfuerzos cortantes en el seno del fluido viscoso, la condicin a ser
satisfecha en las superficies rgidas de fronteras S , en contacto con el fluido viscoso son
diferentes que el caso no viscoso. Para el caso del flujo del fluido viscoso, en la superficie de
frontera S , la velocidad estar dada por vS que representa la velocidad a la que se mueve la
superficie. Si sta se encuentra en reposo es evidente que vS = 0 . Sin embargo, las
condiciones impuestas al fluido implican que la componente normal de la velocidad de ste
sea la misma que la de la superficie slida (en el punto de contacto), lo cual representa en s
una restriccin a la componente tangencial. Esto en consecuencia representa que el fluido
en contacto con la superficie slida se deba mover en conjuncin con dicha superficie, lo
cual representa que el fluido est adherido a la superficie y por consecuencia no
405
puede deslizarse sobre sta. Esta condicin fue primero propuesta por Stokes, y es
conocida como condicin de no deslizamiento. Con la intencin de satisfacer esta condicin
de frontera, Prandtl en 1905 propuso la hiptesis que en una zona muy cercana, adyacente a
la superficie de la frontera, la velocidad relativa del fluido se incrementa muy rpido desde
cero (en la frontera slida) hasta la del flujo en la zona exterior de dicha zona. Esta delgada
capa es denominada como capa lmite, al interior de la cual los efectos de la viscosidad son
predominantes. Fuera de sta, las condiciones se pueden considerar como de un fluido no
viscoso. Por consecuencia, los fenmenos disipativos se presentarn en dicha capa.
Ecuacin de transporte de vorticidad para fluidos viscosos

incompresibles de densidad constante (homognea)
Se asume que las fuerzas de cuerpo son derivables a partir de una funcin de potencial ( )

Bi =
xi
condicin que aplicada a las ecuaciones de Navier-Stokes se expresa
vi v p 2 vi
+ vj i = + +
t x j xi x j x j

como ya se mencion, el trmino = representa la denominada viscosidad cinemtica.

Flujo irrotacional como solucin de las ecuaciones de Navier-Stokes
Si bien en las ecuaciones de Navier-Stokes al considerar la descripcin de un flujo

irrotacional para el cual vi = y dado que a partir de la ecuacin de la continuidad se
xi
2 2 2
debe cumplir con 2 = 0; + +
x12 x22 x32
406
Esto se traduce en que las ecuaciones de Navier-Stokes se expresen en la forma
1 v
p + ( v ) + B = + vv
t
En el caso de un fluido viscoso en flujo irrotacional, la ecuacin, como se ha demostrado con

antelacin, se transforma en la ecuacin de Euler, la cual corresponde a un flujo no viscoso.
1 v
p + B = + vv
t
Los resultados indican que un flujo irrotacional no es factible (dinmicamente posible) para
una situacin en donde se presentan fronteras slidas. Un fluido viscoso se adhiere a las
fronteras de tal forma que las componentes normal y tangencial de la velocidad del fluido
correspondern a la frontera. Esto representa que las componentes de la velocidad estn
predefinidas en la frontera. Por ejemplo, si y = 0 representa a la frontera slida la cual se
encuentra en reposo, entonces en sta las componentes tangenciales v x = v z = 0 y la
componente normal vy = 0 . Para un flujo irrotacional las condiciones preestablecidas
(funcin de flujo) en la frontera son = constante para y = 0 lo mismo que v x = v z = 0
=k

=0
y
v2 = 0 y = 0
Pero es conocido que en general no existe solucin para la ecuacin de Laplace = 0

2

cuando = k y n = = 0 en las fronteras del sistema. En consecuencia, a menos que
n
las condiciones en las fronteras del slido tiendan a ser consistentes con las condiciones de
irrotacionalidad, se presentar la formacin de vrtices en las fronteras, los cuales tendern
a propagarse en el seno del fluido de acuerdo con ciertas restricciones. En condiciones
adecuadas, la vorticidad generada por las fronteras slidas es confinada a una capa delgada
de fluido en la vecindad de la frontera, de tal forma que el exterior de la capa de flujo es
irrotacional. Dicha capa a la cual se limita la presencia de vrtices se denomina como capa
lmite.
407
Demostracin de la imposibilidad de cumplimiento

de la ecuacin de Laplace
2 = 0
Considere un problema de conduccin de calor en estado estable. Un fluido se encuentra en

reposo entre dos placas de dimensiones infinitas. La placa inferior se encuentra a una
temperatura constante l y la superior, a n .
FIGURA 7.14 LA TEMPERATURA VARA EN FORMA LINEAL ENTRE LAS DOS PLACAS
La distribucin de temperaturas en estado estable, de acuerdo con la ecuacin de Laplace
2 2 2
es + + =0
x12 x22 x32
que en el caso en estudio se reduce a
2
=0 =c = c1x2 + c2
x22 x2
Empleando las condiciones de frontera = l para x2 = 0 ; = n para x2 = d , entonces

las constantes de integracin quedan
l = c1 (0) + c2 l = c2
n l
n = c1d + l c1 =
d
x2 ( n l )
= + l
d
408
Como se puede observar, cuando los valores de son predefinidos en las placas, los
d
valores de en ellas estn completamente determinados.
dx2
d n l
=
dx2 d
Esto permite ilustrar que en un problema de conduccin de calor en estado estable

(gobernado por la ecuacin de Laplace) en general, no es posible prescribir los valores de
y de las normales a las derivadas de sta en los mismos puntos de la frontera, a menos que
estos resulten consistentes uno con otro.
En una frontera slida el fluido se adhiere, lo que da lugar a que las velocidades en esta
interfase estn definidas por la superficie. Los vrtices son generados en la superficie
difundindose al flujo. En algunos casos los vrtices quedan confinados a una capa delgada
en la vecindad de la frontera, de tal forma que fuera de esta capa el flujo es irrotacional
(figura 7.15). Por ejemplo, en las alas de un avin la capa lmite se extiende en un espesor
no mayor a un centmetro de la superficie del slido, esto es, las velocidades varan
rpidamente desde la correspondiente al avin (en la superficie del ala) hasta la del medio
(velocidad del viento) quedando limitados a esta zona los elevados valores del nmero de
Reynolds, reducindose rpidamente stos conforme el flujo se aleja del ala.
La viscosidad es la responsable de la generacin de vrtices en la regin adyacente al

slido, su efecto depender de la velocidad del flujo v0 . A elevados valores de la velocidad
v0 , la influencia del slido se confina a sus zonas adyacentes, mientras que a bajas
velocidades su efecto se extiende en el fluido en todas direcciones.
Por consecuencia, a elevadas velocidades, el efecto de formacin de vrtices es confinado a

una pelcula delgada alrededor del obstculo a la cual se denomina como capa lmite. A las
afueras de esta capa el flujo es irrotacional. Este concepto permite, al plantear la solucin de
un problema, dividir el flujo en una regin externa irrotacional y una capa lmite viscosa. Esto
simplifica la complejidad de aplicar las ecuaciones de Navier.
409
FIGURA 7.15 CONCEPTO DE CAPA LMITE. LA INFLUENCIA DE LA VISCOSIDAD DEPENDE DE

LA VELOCIDAD DEL FLUJO. A ELEVADAS VELOCIDADES LA VISCOSIDAD (SUS
EFECTOS) SON CONFINADOS A UNA CAPA DELGADA (CAPA LMITE)
7.19 FLUIDO NEWTONIANO COMPRESIBLE
Como ya ha sido mencionado con antelacin, aquellos flujos en los que las variaciones de
densidad son insignificantes se pueden describir como incompresibles; situacin que
favorece la solucin del problema al reducir el nmero de variables. Cuando las variaciones
de densidad en el flujo son importantes y su efecto no se puede despreciar es necesario
definir al fluido como compresible. No se puede generalizar de entrada y relacionar con el
estado de la materia (lquido o gas), de tal forma que se considere a los lquidos siempre
como incompresibles y a los gases como compresibles. Dicha generalizacin, si bien
corresponde con la mayora de los casos prcticos, presenta sus limitaciones ya que los
gases se pueden describir como incompresibles cuando el flujo se caracteriza por
velocidades muy por debajo de la del sonido en el flujo. Definiendo al nmero de Mach M
como la relacin existente entre la velocidad del fluido v y la del sonido vs de tal forma que
v
M = ,
vs
es entonces que se ha demostrado que los cambios de densidad son del orden del 2% para
M < 0.3 , esto representa que para el aire a temperatura ambiente se puede considerar
como incompresible a velocidades menores de 100 m/s. Por otra parte, existen una infinidad
de aplicaciones de ingeniera para las cuales los efectos de la compresibilidad de gases y
lquidos son fundamentales para la correcta descripcin de los fenmenos.
410
En un fluido compresible para ser consistente con el estado de esfuerzos para el reposo y
movimiento la presin p no depender explcitamente de algn trmino cinemtica, por lo
que p = p ( , ) .
Por ejemplo, para los gases ideales se considera que
p = R
donde R es la constante universal del gas ideal ( R = 8.31 J mol K )
ij = p( , ) + ij + 2 Dij
11 = p + Dkk ij + 2 D11
22 = p + Dkk ij + 2 D22
33 = p + Dkk ij + 2 D33
ii = 3 p + (3 + 2 ) Dkk
En el caso de que el fluido sea compresible

ii
p
3
La presin no representa entonces a los esfuerzos compresivos totales.
Por otra parte, se define la compresibilidad volumtrica como
2
k =+ compresibilidad volumtrica
3
Cuando se trata de gases monoatmicos
2
+ =0
3
Por consecuencia, la ecuacin de esfuerzos para un fluido compresible se expresa
2
ij = p ij Dkk ij + 2 Dij + kDkk ij
3
411
Sustituyendo en la ecuacin de movimiento, queda entonces
ij Dv
+ Bi =
x j Dt
p D 1 v v j v v
ij + kk ij + 2 i + + Bi = i + v j i
x j x j x j 2 x j xi t x j

p D 2v v j v v
+ kk + i
+ + Bi = i + v j i
xi xi x j x j xi x j

t
x j
Reordenando
p 2 vi v v
+k Dkk + Dkk + + Bi = i + v j i
xi xi 3 xi x j x j t x j

vi v1 v2 v3
Dkk = = + + = v = div v
xi x1 x2 x3
En forma general, la ecuacin constitutiva queda
Dv
p + ( v ) + ( v ) + k ( v ) + B =
3 Dt
Por consecuencia, se tienen 6 incgnitas v1 , v2 , v3 , p, , , por lo que a las anteriores 4

ecuaciones (ecuacin de estado escalar y ecuacin de movimiento vectorial) habr que
adicionar la Ecuacin de Continuidad
D v
+ i =0
Dt xi
donde la sexta ecuacin es la de la Energa, sin considerar el calor que se genere al interior
del MC
Du v 2
= ij i +
Dt x j x j x j
412
De la ecuacin anterior surgen nuevas incgnitas, ya que es necesario determinar el estado

de esfuerzos ij y la energa interna u :
u = Cv
donde Cv representa el calor especfico a volumen constante; en general, se cumplir que
u = u ( , )
Se concluye entonces que el anlisis para un fluido compresible ser mucho ms complejo
que para el caso de que ste se declare como incompresible; las incgnitas estn dadas por:
i) Campo de velocidades i
ii) Densidad
iii) Presin p ; p = p ( , )
iv) Temperatura
v) Energa interna u
7.20 ONDAS ACSTICAS
La propagacin del sonido puede ser aproximada, considerando la propagacin de disturbios

infinitesimales en un fluido compresible no viscoso. Para un fluido no viscoso negando el
efecto de fuerzas de cuerpo, la ecuacin de movimiento queda
Dv
p =
Dt
p v v
= i + vj i
xi
t x j
De otra forma
1 p vi v
= + vj i
xi t x j
413
Considerando que se parte de condiciones de reposo
= 0 , vi = 0, p = p0
donde y p0 son la densidad y presin en un tiempo inicial.
Se supone que a partir del reposo se producen los disturbios
= 0 + '( xi , t ), vi = v '( xi , t ), p = p0 + p '( xi , t )
FIGURA 7.16 LA PROPAGACIN DE ONDAS ACSTICAS GENERA PEQUEAS

VARIACIONES DE PRESIN Y DENSIDAD
Sustituyendo en la ecuacin de movimiento
vi v 1 p
+ j i =
t x j ( 0 + ) xi
Como los disturbios generados son muy pequeos (infinitesimales)
vi
vj 0 0 + 0
x j
vi 1 p
= (7.6)
t 0 xi
De manera anloga, si se considera la ecuacin de conservacin de masa se tendr
D
+ ( v) = 0
Dt
414
v
+ vi + i =0
t xi xi
vi
0 1 + + + vi =0
xi t xi
Al ser un disturbio infinitesimal

vi 0, 0
xi 0
vi
0 = (7.7)
xi t
Derivando 7.6 con respecto a xi
2vi 1 2 p
= (7.8)
xi t 0 xi xi
Derivando 7.7 con respecto a t
2vi 1 2
=
txi 0 t 2 (7.9)
Igualando de 7.8 y 7.9, se tiene que
2 p 2
=0
xi xi t 2
Se puede considerar sin incurrir en un error mayor, que el flujo es a temperatura constante,
de tal forma que la presin depende explcitamente solo de la densidad. Considerando que
p1V1 p2V2
=
1 2
p = R
p = p( )
415
Expandiendo la ecuacin en series de Taylor y considerando que se trata de un fluido

barotrpico p = p ( ) , entonces la presin se puede expresar como
p 1 n p
p = p0 + ( ) + ...... + ( 0 )n
0 n ! n
0
0
Si despreciamos el efecto de trminos de orden superior, se puede definir
p = p0 + p ( x, t )
p
p = p0 + ( 0 )
0
p = C02
p
C02 =
0
= 0
p = p p0
Entonces, para un flujo barotrpico (la densidad slo depende de la presin) se tiene que,
partiendo de la ecuacin diferencial en derivadas parciales (ecuacin 7.10):
2 p 2 p
=
xi2 t 2
2 2
C02 =
xi2 t 2
Esta ecuacin es equivalente a la de la onda elstica (slido elstico, homogneo, lineal e

isotrpico) que ha sido definida en el captulo 6; razn por la cual se concluye que las
variaciones de presin y de la densidad se propagan en el fluido con una velocidad
dp
C0 =
d 0
416
Dicha velocidad C0 es definida como la del sonido en condiciones de estancamiento, donde
dp
la velocidad local del sonido C est dada por C =
d
Cp
Para un flujo isoentrpico se tiene p = k , donde = , siendo estos los calores
Cv
especficos del fluido a presin y volumen constante. En este caso la velocidad del sonido se
P
expresa C = .

Con la finalidad de comprobar la relacin antes desarrollada con la establecida para SEHLI,
considerando desplazamientos infinitesimales (ecuacin de Navier), se tiene
2 ui jj 2 ui
0 = ( + ) + + 0 Bi
t 2 xi x j x j
Como el fenmeno es lineal, se pueden estudiar por separado las ondas longitudinales y
transversales. Si se considera una onda longitudinal,
2
u1 = a sen ( x1 v t ) ; u2 = 0 ; u3 = 0
Considerando un fenmeno lineal, la ecuacin diferencial que gobierna queda
2u1 2u1
= v2 ; u1 = u ( x1, t )
t 2 x12
+ 2
v2 = y u2 = u3 = 0
0
p
por analoga v2 =
0
1/2
p
Por lo tanto, la velocidad del sonido est dada por C0 =
0
417
EJERCICIOS RESUELTOS
1. Un elemento cilndrico de radio r , longitud l y densidad , est sujeto al fondo de un

recipiente de acuerdo con la siguiente figura. El fluido en el que se encuentra inmerso el
elemento cilndrico presenta una densidad l . Con base en lo antes expuesto, determine
el nivel de esfuerzos en el cable C cuyo dimetro es d y su lmite elstico es 0 .
FIGURA 7.17 RECIPIENTE CON DENSIDAD INMERSO EN UN FLUIDO
CON DENSIDAD l Y SUJETO AL FONDO
SOLUCIN
Tomando en cuenta el principio de Arqumedes:
d = 2r
f = ( l ) Vg
d2
V = l
4
d2
f = l ( l ) g
4
FIGURA 7.18 SI LA DENSIDAD DEL CUERPO ES MAYOR
QUE LA DEL LQUIDO LA CARGA SOBRE EL
f d l ( l ) g
2
CABLE SER DE COMPRESIN
= =
A d2
2
d
= ( l ) gl
d
418
2. Un recipiente con un fluido incompresible se mueve verticalmente hacia abajo con una
aceleracin constante a . Determine la presin en un punto que se encuentra a una
profundidad h de la superficie. Considere que en la superficie ( h = 0 ) la presin est
dada por p0 y que la densidad se expresa como .
SOLUCIN
Al moverse el recipiente hacia arriba, la inercia del fluido genera un efecto adicional.
FIGURA 7.19 EL RECIPIENTE QUE CONTIENE AL FLUIDO DE DENSIDAD SE

MUEVE HACIA ABAJO CON UNA ACELERACIN a
para un fluido se tiene: p + B = a
p
Eje x1 =0 p = cte. con respecto a x1
x1
p
Eje x2 =0 p = cte. con respecto a x2
x2
B3 = g a3 = a
p
Eje x3 + B3 = a3
x3
p = ( a3 g ) x3 + p0
419
3. Un vehculo arrastra un contenedor rectangular separado por una divisin central. El

contenedor tiene una longitud l de 6 m por 2 m de ancho e igual de altura, el fluido en
reposo dentro del contenedor alcanza una altura de 1 m. Al ponerse la luz del semforo
en rojo, el vehculo debe desacelerar con una magnitud constante ( a -aceleracin en
direccin horizontal). Considerando movimiento de cuerpo rgido, determine el ngulo
de la superficie libre del tanque; asimismo, determine la ecuacin que define la presin
para cualquier punto del tanque. Calcule la altura mxima que alcanza el fluido al chocar
con la pared vertical ( a = g /3).
FIGURA 7.20 AL ACELERAR EL VEHCULO POR EFECTOS DE LA INERCIA EL LQUIDO TIENDE A

MOVERSE HACIA ATRS (FIGURA SUPERIOR); POR OTRA PARTE AL PRESENTAR UNA
ACELERACIN NEGATIVA (FIGURA INFERIOR) EL FLUIDO SE DESPLAZA HACIA ADELANTE
Para el caso en que el vehculo acelera a partir del reposo, el lquido por efecto de la inercia
se inclina hacia atrs, en caso de acelerar la inclinacin ser en sentido opuesto, situacin
que corresponde al caso en estudio.
a1 = a, a2 = a3 = 0, B1 = B3 = 0, B2 = g
SOLUCIN
FIGURA 7.21 EL NGULO DE INCLINACIN DE LA SUPERFICIE

ES FUNCIN DE LA ACELERACIN
420
g
= angtan 3 = 18.43
g

B = 0e1 + ge2 + 0e3
a=g e1 + 0e2 + 0e FIGURA 7.22 LOS EJES SE COLOCAN EN EL VERTICE

3
DE LA SUPERFICIE DEL FLUIDO
p + B = a
p
=0 p = cte. con respecto a x3
x3
p
+ g = 0 p ( x2 ) = gx2 + f ( x1 ) + p0
x2
p ( x1, x2 ) = gx2 + ax1 + p0
x2 x2
= tan x1 =
x1 tan
El fluido es incompresible por lo que el volumen ocupado por ste antes y despus de
su desplazamiento son iguales, asimismo debe recordarse que el tanque tiene una
divisin en el centro
h1 + h2 l l
b = h0 b
2 2 2
l l 6
h2 = h1 + tan h2 = h0 + tan h2 = 1 + tan18.43 = 1.5m
2 4 4
l
V = cte ; V0 = h0 b
2
b = cte ; l = cte
h1 + h2 l
Vf = b h1 = h0 h = 0.5 m
2 2
421
4. El tensor ij describe el estado de esfuerzos para un punto X1 de un fluido

incompresible. Para ij determine el tensor ij que depende exclusivamente de la
velocidad de deformacin, asimismo determine la presin hidrosttica asociada.
45 6 10
ij = 6 15 15 103 [Pa]
10 15 0

Si la densidad del fluido es de 1000 Kg/m3 , indique la profundidad (de ser esto posible) a
la que se encuentra el punto en estudio.
SOLUCIN
ij = p ij + ij
ij = p ij + kk ij + 2 ij
Fluido incompresible
ij = p ij + 2 ij
11 = p ij + 211
22 = p ij + 2 22
33 = p ij + 2 33
ii = 3 p + 2 ii
Como ii = 0 ii = 3 p
ii
H = = p = 20 103 [ Pa] p = 20 103 [ Pa]
3
kg
H = gh ; = 1000 3
m
422
m
g = 9.8
s2
h 2m
ij = ij + p ij
20 0 0

p =H = 0 20 0 103 [ Pa ]
0 0 20

25 6 10
ij = 6
5 15 103 [ Pa ]
10 15 20

5. La compuerta AB de la figura es rectangular de 1 m de ancho y 4 m de largo, sta gira

sobre el punto A . Si la masa de la compuerta es de 400 kg, encuentre las reacciones en
A y B . El fluido en el tanque es agua.
FIGURA 7.23 COMPUERTA SUMERGIDA EN UN TANQUE CON AGUA
SOLUCIN
Para analizar el problema se puede considerar a la compuerta como una viga simplemente
apoyada en A, B ; la carga se incrementa linealmente de acuerdo con la profundidad
f A = pA = ghA A , mientras que en B , la carga ser f B = f A + g hA . Entre estos dos
puntos la carga vara de acuerdo con f = f A + g x sen . Donde x se mide en la direccin
423
de la compuerta a partir de A ( x = 0 ) y hasta B ( x = 4 m ) . Si se considera el peso de la

compuerta, ste deber colocarse en el centro de cargas, tal como se observa en la figura 7.24:
FIGURA 7.24 DISPOSICIN DE LA COMPUERTA
p A = ghA ; hA = 3m
pB = ghB = g ( hA + h )
h = l sen B = 4sen 70 = 3.75 m
hB = hA + h = 3 + 3.75 = 6.75 m
FIGURA 7.25 DESCOMPOSICIN DE CARGAS
424
Dado que se trata de un sistema que se encuentra en equilibrio, entonces

f = 0; M = 0 , lo momentos en los apoyos es nulo, por lo tanto,
g h
f = 0 mg cos + ( ghA ) wl + wl = RA + RB
2
l l g h 2
M A = 0 0 = RBl ( mg cos ) ( ghA ) wl wl l
2 2 2 3
4 RB (2681.43 + 235200 + 196000) = 0
RB = 108470 N
RA = 192440 RB
RA = 83970.7 N
6. Desarrolle (en coordenadas rectangulares) la ecuacin de conservacin de cantidad de

movimiento para un fluido newtoniano incompresible. Indquela en forma general
SOLUCIN
Ecuacin de Cauchy
Dv
div + B =
Dt
v
T + B = + vv
t
T = pI + T
T - esfuerzos viscosos
T = I + 2 D
1
2
( T
T = pI + ( v ) I + 2 v + ( v )

)
425
Sustituyendo en la ecuacin de Cauchy:

(
1
2
)
T

v
pI + ( v ) I + 2 v + ( v ) + B = + vv
t
v
p + ( v ) + ( v ) + ( v ) + B = + vv
t
Considerando incompresibilidad
v
p + ( v ) + B = + vv
t
7. El tensor de esfuerzos asociado a un punto q de un fluido newtoniano incompresible est

dada por:
60 3 4

ijq = 3 100 8 [ kPa ]
4 20
8
x1
x2 h
x3
FIGURA 7.26 UN OBJETO SE ENCUENTRA SUMERGIDO A UNA PROFUNDIDAD h
Si la densidad del fluido es de 1000 kg/m3 y Bi = 0e1 + 0e2 ge3 donde g es la aceleracin
de la gravedad, determine la profundidad a la que se encuentra inmerso el cuerpo en el seno
del fluido. Si la viscosidad es de 9.82x10-4 Pa-s, determine el tensor de rapidez de
deformacin y sus valores principales.
426
SOLUCIN
60 3 4
ij = 3 100 8

[ kPa ]
4 20
8
Fluido incompresible
ij = p ij + 2 ij
11 = p + 2 ij ; 22 = p + 2 ij ; 33 = p + 2 ij
kk = 3 p + 2 ( 11 + 22 + 33 ) ; kk = 0
kk = 3 p = 3 H
kk
H = = 60 kPa
3
p = gh h = p
g
60 kPa
h= 6m
1000 9.8
60 0 0

ijH = 0 60 0 [ kPa ]
0 0 60

ij = p ij + ij
ij = kk ij + 2 ij ij = 2 ij
ij
ij =
2
N
= 9.82 104 -s
m2
427
0 3 4
ij = 3 40 8

[ kPa ]
4 8 40

0 1527.5 2036.6
ij = 1527.5 20366 4073

2036.6 4073 20366

20570 0 0

ijP = 0 307 0 s 1
0 0 20877

8. Sea ( xi , t ) una funcin escalar la cual define el campo de velocidades del medio
continuo en el intervalo analizado, como

vi = ; v =
xi
a) Qu caractersticas deber cumplir la funcin para que el flujo en el intervalo sea

irrotacional?
b) Si la funcin ( xi , t ) est asociada a un medio incompresible Cmo quedar
expresada la ecuacin de conservacin de masa, en el intervalo en estudio?
SOLUCIN

= ( xi , t ) tal que vi =
xi
a) El rotacional se expresa como:
1 v v j
wij = i
2 x j xi
sustituyendo
1 1 2 2
wij = = +
2 x j xi xi x j 2 x j xi xi x j

428
wij = 0 si ( xi , t ) tal que

vi =
xi
b) Medio incompresible
v = 0 Coordenadas rectangulares
vi
=0
xi

=0
xi xi
2
=0 2 = 0
xi2
( )
9. Analice si la funcin escalar = A x12 x22 + 2 x32 t describe un flujo irrotacional, si
t s, x m, A = 1 s 2 . Asimismo, determine el campo de esfuerzos asociado
considerando que la viscosidad est dada por . Verifique si el campo de esfuerzos
satisface las ecuaciones de Navier-Stokes. Las fuerzas de cuerpo se pueden despreciar.
SOLUCIN
(
= A x12 x22 + 2 x32 t )
Flujo irrotacional
w=
1
2
(
v vT = 0 )
v1 = 2 x1 At ;
v2 = 2 x2 At ;
v3 = 4 x3 At
2 = 0 At (2 + 2 4) = 0
Por lo tanto, se cumple la condicin de incompresibilidad.
429
2 0 0

v = 0 2 0 A t
0 0 4

v = ( v )
T
w=0
Por consecuencia, se cumple la condicin de irrotacional.
El estado de esfuerzos est dado por
ij = p ij + ij
donde ij = esfuerzos viscosos
Por analoga con SEHLI se tiene que la ecuacin constitutiva para un fluido est dada por
vk 1 v v j
ij = ij + 2 i +
xk 2 x j xi

vk
= 2 = 0
xk
Por lo tanto, un fluido irrotacional e incompresible la ecuacin se reduce a
vi v j
ij = +
x j xi

vi v j
Por lo que ij = p ij + +
x j xi
p + 4 At 0 0

ij = 0 p + 4 At 0

0 0 p 8 At
430
10. Considere que entre las placas A y B existe un fluido newtoniano incompresible. La placa
A se desplaza a una velocidad de 5 m/s, mientras que la placa B permanece en reposo.
Si la distancia entre ambas placas es de 1 m, determine la velocidad del fluido a 0.2 m de
la placa A. Considere que las placas estn horizontales y que sus dimensiones son muy
grandes.
A
v = 5 m/s
B
x3
h = 1 m , x = 0.2 m
FIGURA 7.27 FLUJO INDUCIDO POR VELOCIDAD (COUETTE)
SOLUCIN
De las ecuaciones de NavierStokes
p vi Dv
+ i = i
xi x j x j Dt
Considerando que es un flujo inducido por velocidad y establecido
p Dvi
=0 y =0
xi Dt
De la ecuacin de conservacin de masa
v1 v2 v3
v = + + =0
x1 x2 x3
v1
Como v2 = v3 = 0 =0
x1
v1 es constante en x1 y v1 = v1 (x3)
431
En la ecuacin de NavierStokes para el eje x1
2v1
+ B1 = 0
x32
Despreciando el efecto de las fuerzas de cuerpo
2v1
=0 v1 = x3 +
x32
De las condiciones de frontera en x3 = 0 , v1 = vA = 5 =5
En x3 = h = 1 , v1 = 0 =5
v1 = 5x3 + 5
Para x3 = 0.2 v = 4 m/s
11. Para el flujo que se describe en la figura, determine el perfil de velocidades en la zona del
conducto que va de AA .
FIGURA 7.28 FLUJO INDUCIDO POR PRESIN
Considere que se trata de un fluido newtoniano incompresible con viscosidad . El canal

es rectangular con un ancho igual a 2e y altura d , donde d << e y por consecuencia se
puede despreciar el efecto en las paredes laterales. El flujo en el canal es uniaxial de tal
forma que v2 = v3 = 0 . Determine el gasto volumtrico.
432
SOLUCIN
Se trata de un flujo de Poiseville, es decir, un flujo uniaxial inducido por presin.
FIGURA 7.29 EL FLUJO VA EN SENTIDO OPUESTO AL GRADIENTE DE PRESIN. EN

LAS PAREDES DEL CONDUCTO LA VELOCIDAD ES NULA (LA DE LAS
PAREDES) Y EN EL CENTRO ES MXIMA
v1 = v ( x2 )
v2 = v3 = 0
p 2 2 2 v v v v
+ 2 + 2 + 2 v1 + B1 = 1 + v1 1 + v2 1 + v3 1
x1 x1 x2 x3 t x1 x2 x3
No considere el efecto de fuerzas de cuerpo, y adems se considera flujo establecido

v1
=0
t
Condiciones de frontera
d
para x1 = v1 = 0
2
d
para x1 = v1 = 0
2
para x1 = 0 v1 = vmx
433
Ecuacin de Conservacin de masa
v1 v2 v3
+ + =0
x1 x2 x3
v1
=0 El flujo no se acelera a travs del conducto
x1
De las ecuaciones de Navier-Stokes
p 2v 2 v1 p 1
+ 21 = 0 =
x1 x2 x22 x1
v1 p x2
= + c1
x2 x1
p x22
v1 = + c1 x2 + c2
x1 2
En un mximo, la primera derivada es cero
v1 p x2
=0= + c1 c1 = 0
x2 x2 = 0 x1
d
Cuando x2 = v1 = 0
2
( )
2
d
p 2
0= + c2
x1 2
p d 2
c2 =
x1 8
1 d p
2
v1 ( x2 ) = x22
2 2 x1
El perfil de velocidades es parablico para un flujo inducido por presin, y ste va en

sentido contrario al del gradiente de la presin.
434
12. Sea ( xi , t ) una funcin escalar a travs de la cual se pretende definir el campo de
velocidades en un medio continuo, de acuerdo con

vi =
xi
a) Qu caractersticas deber cumplir la funcin para que el flujo en el intervalo sea

irrotacional?
b) Si la funcin ( xi , t ) est asociada a un medio incompresible Cmo quedar
expresada la ecuacin de conservacin de masa en el intervalo en estudio?
SOLUCIN
a) La funcin deber ser continua y continuamente derivable en el intervalo.

vi = o v =
xi
( ) = 0
b) La ecuacin de conservacin de masa para el intervalo de estudio queda expresada
v = 0 2 = 0
13. Un continuo con una ecuacin constitutiva de la forma
ij = kk ij + 2 ij
presenta un flujo irrotacional e incompresible, el cual se describe a travs de una funcin

escalar , que permite describir el campo de velocidades a travs de v = .
a) Con base en lo anterior, determine la funcin disipativa WV = ij ij o WV =
b) Determine la ecuacin de Cauchy para el material en cuestin considerando la

descripcin de su campo de velocidades y de su ecuacin constitutiva.
435
SOLUCIN
ij = kk ij + 2 ij

v = vi = e1 e2 e3
x1 x2 x3
WV = ij ij
Para un flujo irrotacional e incompresible
2 = 0
Y la ecuacin constitutiva
ij = 2 ij
WV = ij ij = 2 ij ij
(
WV = ij ij = 2 11
2
+ 22
2
+ 33
2 2
(
+ 2 12 + 23
2
+ 31
2
))
Ecuacin de Conservacin de masa
2 2 2
+ + =0
x1x1 x2 x2 x3 x3
2 = 0
14. Dado el siguiente campo de velocidades de un fluido viscoso newtoniano incompresible
v1 = kx1 v2 = kx2 v3 = 0
a) Demuestre si el campo de velocidades es irrotacional.

b) Determine el tensor de esfuerzos.
c) Determine el campo de aceleraciones.
d) Demuestre si el campo de velocidades satisface la ecuacin de Navier-Stokes, y esto
permite determinar la funcin de presin (distribucin de presiones). Para lo anterior
desprecie las fuerzas de cuerpo y considere que p = p0 en el origen.
e) Determine la rapidez de disipacin de energa mecnica en calor.
436
f) Si x2 = 0 representa una frontera fsica, qu condicin no es satisfecha por el

campo de velocidades?
SOLUCIN
a)
Irrotacional

vi =
xi
e1 e2 e3
v v v v v v
v = = 3 2 e1 + 1 3 e2 + 2 1 e3
x1 x2 x3 x2 x3 x3 x1 x1 x2
v1 v2 v3
v3 v2 v v v2 v1
= 0; 1 = 3 = 0; = =0
x2 x3 x3 x1 x1 x2
v = 0 Campo irrotacional
ij = p ij + ij
ij = kk + 2 ij
v1 v2 v3
kk = + + = k k +0=0
x1 x2 x3
k 0 0
ij = 0 k 0
0 0 0
k 0 0
ij = 0 k 0 2
0 0 0
437
b)
p + 2 k 0 0
ij = 0 ( p + 2 k ) 0
0 0 0
Dv1 vi vi
= + vj
Dt t x j
Dv1 v1 v1 v v
= + v1 + 1 v2 + 1 v3
Dt t x1 x2 x3
Dv1
= k 2 x1
Dt
Dv2 v2 v2 v v
= + v1 + 2 v2 + 2 v3
Dt t x1 x2 x3
Dv2
= k 2 x2
Dt
c)
Dvi
ai = = k 2 ( x1e1 + x2 e2 + 0e3 )
Dt
Ecuacin de Navier-Stokes
v
p + ( v ) + B = + ( v ) v
t
p 2 vi v v
+ + Bi = + i vj
xi t x j
x j x j
Despreciando fuerzas de cuerpo y considerando incompresibilidad
Eje x1
p 2v v v
+ 21 + B1 = 1 + 1 v1
x1 x
1 t x1
438
p
= k 2 x1
x1
k 2 x12
p = + f ( x2 , x3 ) + cte1
2
Eje x2
p 2v v v
+ 22 + B2 = 2 + 2 v2
x2 x
2 t x2
p
= k 2 x2
x2
k 2 x22
p = + f ( x1, x3 ) + cte2
2
Eje x3
p 2v
+ 23 + B3 = 0
x3 x
3
p ( x3 ) = cte3
k 2 2
p =
2
( )
x1 + x22 + p0 Incompresible
d)
vi vi
W = ij
x j
(
= p ij + 2 ij
x j
)
v vi 2
W = ij i + 2 ij ij = 0 Irrotacional =0
xi xi xi2
W = 2 ij ij
W = 2 11
2
(+ 22
2
+ 33
2
+ 2 12
2
(
+ 23
2
+ 31
2
))
439
e) Al sustituir ij
( )
W = 2 k 2 + k 2 = 4 k 2
f) No se cumple con la ecuacin de conservacin de masa (Ecuacin de la continuidad)
, x2 = k , para cualquier x2 en x2 = 0, v2 0
15. Demuestre si la ecuacin de movimiento para un fluido newtoniano compresible se puede

expresar como:
v
p + k ( v ) + ( v ) + (v ) + B = + vv
3 t
a) Cuntas incgnitas se presentan en este sistema?

b) Cules son stas?
c) A qu otras ecuaciones se puede acudir para resolver el sistema?
d) Exprese la ecuacin en notacin ndice.
SOLUCIN
Fluido newtoniano compresible
ij = p ij + kk ij + 2 ij
Ecuacin de Cauchy
Dv ij Dvi
+ B = + Bi =
Dt x j Dt
440
v 1 v v Dv
p ij + j ij + 2 i + j + Bi = i
x j
x j
2 x
xi
j
Dt

ii = ( 3 + 2 ) kk
2 v j
H = + kk p = p + k
3 x j
p v j 2 vi v v
+ ( + ) + + Bi = i + i v j
xi xi x j x j x j t x j

p v j v j 2 vi Dv
+k + + = i
xi xi x j 3 xi x j

x j x j

Dt
v
p + k ( v ) + ( v ) + ( v ) = + ( v ) v
3 t
Nmero de incgnitas: 5 y son v1 , v2 , v3 , y p
Para resolver las ecuaciones se debe recurrir a las ecuaciones de conservacin de masa y
conservacin de energa.
441
EJERCICIOS PROPUESTOS
Fluidos newtonianos. Hidrosttica. Movimiento de cuerpo rgido
1. Un fluido en estado de reposo presenta una ecuacin de estado de la forma
p = k
donde y k son constantes, y es la densidad.
La nica fuerza de cuerpo al que est sometido el sistema es la gravedad. Considerando

una orientacin de ejes como se muestra en la figura y que la presin en la referencia
est dada por p0 , determine:
x1
a) Variacin de la presin en los ejes x1 , x2 . x2
b) Variacin de la presin en funcin de x3 .
x3
2. La compuerta AB de la figura es rectangular de 40 cm de ancho y 3 m de largo, la cual

rota sobre el punto A . Si el peso de la compuerta es de 400 kg, encuentre las reacciones
en A y B . El fluido en el tanque es agua.
3. La compuerta de la siguiente figura tiene 6 m de largo por 4 m de ancho. Si el peso de la

compuerta es de 1000 N, determine el nivel del agua h al cual la compuerta empieza a
caer.
442
4. Un sistema hidrulico soporta las masas M1 y M 2 en cada uno de los vasos

comunicantes, los cuales contienen fluidos de densidad y viscosidad 1 , 2 , 1 , 2 .
Determine la masa M 2 que garantiza el equilibrio del sistema, esto a partir de
1 , 2 , M1 , h . Considere que los fluidos son newtonianos incompresibles y no existe
mezcla.
M1
M2
1
h 2
5. Qu caracteriza a un fluido newtoniano? D algunos ejemplos de fluidos que se pueden

modelar como newtonianos.
6. Cmo se describe el comportamiento de un fluido no newtoniano? D algunos ejemplos

de fluidos que se pueden modelar como no newtonianos.
7. Desarrolle la ecuacin de conservacin de cantidad de movimiento para un fluido

newtoniano incompresible. Indquela en forma general y en notacin ndice. Desarrllela
en coordenadas rectangulares, cilndricas y esfricas.
8. Los lquidos en los vasos comunicantes de la siguiente figura estn en equilibrio.

Determine h2 como una funcin de 1 , 2 , 3 , h1 , h3 .
443
Los lquidos no se pueden mezclar.
3 h3
h1 2
h2
2
9. Un recipiente con agua se mueve verticalmente hacia arriba con una aceleracin a .
Determine la presin a la que se encuentra un punto en la profundidad h .
10. En aplicaciones de astrofsica, una atmsfera que tiene una relacin de la forma
n
P
=
P0 0
donde P0 y 0 son la presin y densidad de referencia, se denomina como atmsfera

politrpica. Para este tipo de atmsferas determine la distribucin de presin y densidad.
11. Defina el concepto de capa lmite.
12. La ecuacin de estado de un fluido barotrpico presenta la forma p = , donde y k

k
son constantes; siendo el flujo isoentrpico. Verifique si para este caso la ecuacin de
Bernoulli en estado estable est dada por
kp 1
+ + v 2 = cte
(k + 1) 2
Ln 1
y si el flujo es isotrmico la ecuacin de Bernoulli queda + p + v 2 = cte
2
13. Verifique que el gasto volumtrico para un proceso de extrusin de un polmero, el cual
se desarrolla mediante un tornillo extrusor est dado por
D h3 sen 2
Q = 0.5 2 D 2 N h sen cos p
12 l
444
Para lo cual, considere que el fluido se puede modelar como newtoniano y que las
condiciones son isotrmicas.
h = espesor de la pelcula
p = presin D = dimetro del can o barril
= viscosidad N = velocidad de rotacin del husillo
vA= velocidad del cojinete (tangencial) h = profundidad del canal del husillo
vB= velocidad de la flecha (tangencial) = ngulo entre la hlice y la direccin
perpendicular al husillo
P = presin de descarga del husillo
l = longitud del husillo
= viscosidad
Husillo
h
w
b
Z
D
X
Canal de flujo Barril

p
b
Modelo simplificado del canal
Y , Y
Vx
Vx
X X
Z Z
Flujo de arrastre (QA) Flujo de presin (QP)
Para lo anterior, se define que el gasto volumtrico neto Q se puede expresar como la suma
de un flujo de arrastre por velocidad QA (Couette) menos un flujo de presin QP , este ltimo
generado por el incremento de presin que se produce hacia la zona de descarga. El flujo de
arrastre se desplaza hacia adelante y ocurre por el movimiento de la superficie del husillo en
contacto con el fluido, mientras la otra permanece fija.
445
Para facilitar el anlisis al calcular el flujo de arrastre, suponga que el barril gira y el husillo
permanece inmvil; adems, considere que el flujo est dado por Q = vAc , donde v es la
velocidad promedio y Ac la seccin transversal del canal.
w - velocidad angular del husillo

p - paso del husillo
Considere que Q = QA QP
14. Demuestre que la ecuacin de movimiento para un fluido newtoniano compresible se

puede expresar como
v
p + k (div v ) + (div v ) + div(v ) + B = + vv
3 t
a) Cuntas incgnitas se presentan en este sistema?

b) Cules son stas?
c) A qu otras ecuaciones se puede acudir para resolver el sistema?
d) Exprese la ecuacin en notacin ndice.
15. Para un fluido newtoniano incompresible con flujo irrotacional, deduzca la ecuacin de
Torricelli a partir de las ecuaciones de Navier-Stokes. Considere que slo existe el campo
gravitacional y que el flujo es establecido.
16. Deduzca la ecuacin que describa el comportamiento de un fluido newtoniano

incompresible, presntela en forma invariante (notacin general) y en notacin ndice
para coordenadas rectangulares, en este ltimo caso exprese las ecuaciones a las que
se da lugar.
a) Cules son las incgnitas y qu otras ecuaciones se emplearn para poder resolver
las incgnitas?
b) Exprese las ecuaciones complementarias.
c) Qu pasa en el caso de que el fluido newtoniano sea compresible?
d) Cules sern las incgnitas adicionales, qu otras ecuaciones son empleadas para
la solucin del problema?
446
17. Considere que entre las placas A y B existe un fluido newtoniano incompresible. La
placa A se desplaza a una velocidad de 1 m/s, mientras que la placa B permanece en
reposo. Si la distancia entre ambas placas es de 1 m, determine la velocidad de placa del
fluido a 0.3 m de la placa A . Considere que las placas estn horizontales y que sus
dimensiones son muy grandes.
2
18. Al trmino + en un fluido newtoniano se le denomina viscosidad volumtrica. En
3
un fluido newtoniano, en una coordenada ( x1 , x2 , x3 ) a un tiempo t el estado de
esfuerzos est dado por ij . Si el fluido presenta una viscosidad volumtrica igual a cero,
determine la profundidad a la que se encuentra inmerso, si las fuerzas de cuerpo estn
dadas por Bi . La densidad es de 1000 kg/m3.
9.8 4 6
x3
ij = 4 39.2 2

kPa
g
6 9.8
2
x2
x1
Asimismo, determine el tensor de esfuerzos viscosos. Considere que:
Bi = 0e1 + 0e2 ge3 g = 9.8 m/s 2
19. El tensor ij describe el estado de esfuerzos para un punto X i de un fluido
incompresible. Para ij determine el tensor ij que depende exclusivamente de la

velocidad de deformacin; asimismo, determine la presin hidrosttica asociada.
447
25 8 5
ij = 8 15 12
5 12 5

20. Para un flujo irrotacional de un fluido no viscoso, homogneo e isotrpico, deduzca la

ecuacin de Bernoulli a partir de las ecuaciones de Navier-Stokes.
21. Determine el tensor de esfuerzos asociado a los siguientes campos de velocidad,

considerando que se trata de un fluido viscoso.
a) v1 = 0, v2 = 0, v3 = x2
b) v1 = 0, v2 = x23 , v3 = 2 x2 x3
c) v1 = v1 ( x1 , x2 ) , v2 = v2 ( x1 , x2 ) , v3 = 0
22. Para un flujo axisimtrico estable inducido por presin a travs de una tubera de
dimetro d , compruebe que el gasto msico est dado por
d4
M =
128
considerando que se trata de un fluido newtoniano incompresible de densidad , donde

representa el gradiente de presin
p
=
z
donde p es la presin y z es la direccin de flujo.
23. Demuestre que el campo de velocidades dado por
x2 x2 2 x1 x2
v1 = 1 4 2 , v2 = , v3 = 0
R R4

donde representa una constante diferente de cero y R = x1 + x2 0 obedece a la

2 2 2
ecuacin de movimiento de Euler. Asimismo, determine la distribucin de presiones

asociada al campo de velocidades.
448
24. Para el siguiente campo de velocidades en coordenadas cilndricas:
vr = v(r , ) , v = 0 , vz = 0
f ( )
A partir de la ecuacin de la continuidad verifique que vr =
r
2 f ( f )2
En ausencia de fuerzas de cuerpo, compruebe que +4f + =0 y
2
f
p = 2 + c , donde c representa una constante.
r2
25. Explique usted, con base en conceptos fundamentales de la MMC, por qu un chorro de
agua lanzado verticalmente hacia abajo tiende a adelgazarse (reduce aparentemente su
seccin) conforme se desplaza hacia el suelo. Qu pasa ahora cuando el chorro se
lanza verticalmente hacia arriba?
26. Determine la ecuacin constitutiva para un fluido newtoniano para el cual su condicin de
Stokes es cero.
27. Demuestre que para el campo de velocidades v1 = v ( x2 , x3 ) , v2 = v3 = 0 , las ecuaciones

de Navier Stokes, despreciando el efecto de las fuerzas de cuerpo, se reducen a
2v1 2v1 1 p
+ = = ; = cte .
x22 x32 x1
28. Determine el gasto (inducido por presin) de un fluido Newtoniano incompresible a travs
de un tubo de seccin transversal elptica: El tubo tiene un semieje menor y un
semieje mayor .
29. Dado el campo de velocidades para un fluido Newtoniano incompresible

v1 = k ( x x ) ; v2 = 2kx1 x2 ; v3 = 0
2
1
2
2
Desarrolle la ecuacin de Navier Stokes.

Determine el tensor de deformaciones y de esfuerzos asociado.
449
30. Si va usted a analizar el flujo de la sangre en venas y arterias, indique cuales

consideraciones deber realizar. Para dicho anlisis, es factible emplear las ecuaciones
de Navier-Stokes?
31. Explique el principio de sustentacin de un objeto ms pesado que el aire. la velocidad

de despegue de un avin, por cules variables estar determinada?
32. Para el diseo estructural de un submarino, qu parmetros debern ser considerados?
450

Fluidos Newtonianos

Cargado por

Copyright:

Formatos disponibles

Fluidos Newtonianos

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Fluidos Newtonianos

Cargado por

Copyright:

Formatos disponibles

CAPTULO 7

FLUIDOS VISCOSOS NEWTONIANOS

7.1 CONCEPTOS GENERALES

FIGURA 7.1 SOLICITACIN A CORTE ENTRE DOS PLACAS ENTRE LAS

En los fluidos se observa que la resistencia al flujo depende de la velocidad y, por

1. Fluidos newtonianos. Son aquellos en los que la relacin esfuerzo de corte-

2. Fluidos no newtonianos. En este caso la relacin es no lineal, presentndose

Desde el punto de vista de la variacin de su densidad se describen como:

b) Fluidos incompresibles. Se puede considerar idealmente que la densidad del fluido

Como ya ha sido mencionado un fluido es un medio idealizado, el cual en cualquier punto,

cualesquiera cuyas normales son n1 , n2 , entonces se cumplir que

Por otra parte,

n1Tn2 n2Tn1 = 0 = (2 1 )n1 n2

Existen entonces dos posibilidades:

Ya que los planos no necesariamente son perpendiculares n1 n2 0

En notacin ndice se expresa

donde el escalar p representa la magnitud de los esfuerzos normales compresivos y por

7.2 FLUIDOS COMPRESIBLES E INCOMPRESIBLES

razn por la que la ecuacin de conservacin de la masa

Un fluido, aun cuando se considere como incompresible, no necesariamente deber

7.3 ECUACIONES DE LA HIDROSTTICA

por lo que sustituyendo sta en la Ecuacin de Cauchy, se tiene

Lo que en notacin general se expresa como

EJERCICIO 7.1 Determine la variacin de la presin de un objeto que se encuentra sumergido

Con base en la ecuacin 7.2, se puede determinar la variacin de la presin en funcin de la

Bi = 0e1 + 0e2 + ge3 .

Considerando la aceleracin gravitatoria

dado que en x3 = 0 p = p0 , se tiene que

Considerando la ecuacin de estado para un gas ideal, y adems que la temperatura es

Sustituyendo resulta entonces que

Resolviendo la ecuacin diferencial se tiene

Para x3 = 0 p = p0 (presin de referencia)

De forma anloga se puede proceder para determinar la razn de variacin de la densidad

Para x3 = 0 = 0 (densidad inicial o de referencia)

7.4 MOVIMIENTO DE CUERPO RGIDO DEL FLUIDO

En los casos anteriores se consider que el fluido se encontraba en reposo o en condiciones

EJERCICIO 7.3 Un recipiente con un fluido incompresible en su interior se mueve

FIGURA 7.2 RECIPIENTE MOVINDOSE EN DIRECCIN VERTICAL

Bi = 0e1 + 0e2 + 0e3

Movimiento de cuerpo rgido

La solucin de la ecuacin diferencial queda

FIGURA 7.3 MOVIMIENTO DEL AGUA AL PRESENTARSE UNA ACELERACIN

La superficie del agua es normal a la resultante de la fuerza, si solo existe la aceleracin

A partir de lo anterior y considerando las dimensiones del tanque ( , ) , donde representa

Tomando como referencia la superficie del fluido

p ( x1 , x2 ) = p0 + gx2 + gx1 tan

FIGURA 7.4 FLUJO UNIAXIAL EN DIRECCIN DEL EJE x1

El vector de posicin de cualquier elemento diferencial del fluido queda

z = ez i r = ( ( ez ie1 ) x1 + ( ez ie2 ) x2 + ( ez ie3 ) x3 )

Por otra parte la variacin de

Se puede escribir entonces que

Sustituyendo en la ecuacin 7.4

7.5 FLUIDO NEWTONIANO

donde el tensor ij depende de la velocidad de deformacin y representa la componente

El tensor de esfuerzos y el de rapidez de deformacin son simtricos

El tensor de constantes viscosas es simtrico

Cijkl = C jikl ; Cijkl = Cijlk

Como a travs de la Segunda ley de Newton se encuentran relacionadas la fuerza, masa,