DINAMICA DE FLUIDOS

El capitulo trata sobre el estudio de los fluidos, tomando en cuenta las condiciones de inicio de movimiento. Para
este estudio se analizaran lo fluidos como ideales y reales; es decir, en los ideales no se toman en cuenta conceptos como los
de viscosidad, considerando así solo fuerzas como las de corte y fricció n, en el caso de los reales se tomaran todos los
aspectos a modo de asemejar del mayor modo posible el aná lisis a la realidad.
El agua en movimiento presenta problemas de solució n má s difícil e incierta que los problemas tratados en la
está tica y cinemá tica de fluidos a causa de la presencia de resistencias debidas a la fricció n y otras causas perturbantes cuya
acció n muchas veces es difícil o imposible expresar matemá ticamente. La mayoría de los problemas concernientes a la
diná mica de fluidos; es decir, al flujo de fluidos en conductos y tubos implican la predicci ó n de las condiciones en una
secció n de un sistema, cuando se conocen las condiciones de alguna otra secció n (Fig4-1), la figura muestra una parte de un
sistema de distribució n de fluido con el flujo corriendo de la secció n 1, en el fondo, a la secció n 2, en la parte superior. En
cualquier secció n de dicho sistema, generalmente nos preocupa la presió n del fluido, la velocidad y la elevació n de la
secció n. Recuerde que la elevació n es el término usado para definir la distancia vertical desde algú n nivel de referencia a un
punto de interés, y se representa con la letra z. Cuando se tratan conductos y tubos, se mide la elevació n a la línea central de
la secció n de interés.

IV.1 OBJETIVOS DEL CAPITULO

El capitulo a continuació n abarca todo lo que implica la diná mica de fluido por lo que al finalizar el mismo al
alumno podrá :
i. Definir rapidez de flujo de volumen, de peso, de masa y sus respectivas unidades.
ii. Definir la diná mica de fluidos en fluidos incompresibles ideales y las diferentes ecuaciones.
iii. Escribir la ecuació n de Bernoulli, y sus aplicaciones.

148
 iv. Escribir la ecuació n de la energía.
v. Definir flujo por orificios y sus ecuaciones fundamentales.
vi. Escribir la ecuació n de cantidad de movimiento.
vii. Definir la diná mica de fluidos reales.
viii. Describir el efecto de la viscosidad y definir lo que es el nú mero de Reynolds.
ix. Escribir la ecuació n de Navier-stokes.
x. Definir flujo uniforme y turbulento, así como la teoría de la capa limite.
xi. Definir distribució n de velocidades y pérdidas de energía.

IV.2 DINAMICA DE FLUIDOS IDEALES INCOMPRESBLES

Estos flujos cumplen el llamado teorema de Bernoulli, enunciado por el matemá tico y científico suizo Daniel
Bernoulli. El teorema afirma que la energía mecá nica total de un flujo incompresible y no viscoso (sin rozamiento) es
constante a lo largo de una línea de corriente. Las líneas de corriente son líneas de flujo imaginarias que siempre son
paralelas a la direcció n del flujo en cada punto, y en el caso de flujo uniforme coinciden con la trayectoria de las part ículas
individuales de fluido. El teorema de Bernoulli implica una relació n entre los efectos de la presió n, la velocidad y la
gravedad, e indica que la velocidad aumenta cuando la presió n disminuye. Para el autor John Muller: "Este principio es
importante para la medida de flujos, y también puede emplearse para predecir la fuerza de sustentació n de un ala en vuelo.
En el caso de la diná mica de fluidos, el autor R.L Streeter, menciona que: "las ú nicas fuerzas de superficie son las
provocadas por la presió n, que sumadas a las demá s fuerzas, o de gravedad, son las responsables del movimiento del
fluido". Bajo estas condicione Newton represento su segunda ley, aplicada a un elemento fluido, o ecuació n de cantidad de
movimiento, la que se conoce como ecuació n de Euler.

IV.2.1 ECUACION DE EULER

Euler amplió y perfeccionó la geometría plana y de los só lidos, introdujo el método analítico a la trigonometría y a
él se debe el tratamiento moderno de las funciones log (x) y exp (x). Cre ó una consistente teoría de logaritmos de nú meros
negativos e imaginarios y descubrió que log (x) tiene un nú mero infinito de valores.
En el campo de la hidrá ulica podemos desarrollar las ecuaciones de movimiento de fluidos ideales (Fig.4-2), en la
figura se observan las fuerzas de presió n sobre un volumen de control, la presió n que siente el elemento en sus lados debido
al fluido que no tomamos en cuenta, para tomarla como fuerza debe multiplicarse por el á rea en contacto. Estas fuerzas
estará n en cada lado del mismo. Estas fuerzas de presione desbalanceadas en cada eje son:

149
  P   Py   P 
Fx    x x yz ; Fy    y xz ; Fz    z z yx
 x   y   x 

Dividiendo todas estas fuerza por unidad de masa: m  xyz

1 P 1 P 1 P
Fx   ; Fy   ; Fz  
 x  y  z

Pe  xyz 4.2.1

Tomando esta fuerza por unidad de masa tendremos:

Pc  xyz 
   g
m xyz 

1 P 1 P 1 P
Fx   ; Fy   ; Fz  
 x  y  z
La segunda ley de Newton dice: la sumatoria de las fuerzas que actú an sobre un cuerpo deben ser igual a su masa

por su aceleració n:  F  m  a , en el aná lisis previo pusimos todas las fuerzas por unidad de masa, por esa razó n las
fuerzas que hallamos precisamente será n iguales a la aceleració n en ese sentido:
V x V x v V x 1 P
 Vx  x V y  Vz   4.2.2
t x y z  x
V y V y v y V y 1 P
 Vx  Vy  Vz   4.2.3
t x y z  y
V z V z v V z 1 P
 Vx  z V y  Vz   4.2.4
t x y z  z

150
 Estas ú ltimas son las ecuaciones de Euler del movimiento de fluidos ideales. Estas ecuaciones se interpretan de la
siguiente manera: En los ejes X e Y las aceleraciones locales o convectivas se traducen en un cambio de presiones en esa
misma direcció n, es decir que una aceleració n en este sentido provoca presió n. En el sentido Z toda aceleració n provoca un
cambio en la presió n ayudada por el peso del cuerpo, por esto se introduce un cambio en la altura piezométrica P /   z

IV.3 RAPIDEZ DE FLUJO FLUIDO

La cantidad de flujo que fluye en un sistema por unidad de tiempo, se puede expresar mediante tres términos:
 Q La rapidez de flujo de volumen, es el volumen del flujo de fluido que pasa por una
secció n por unidad de tiempo.
 W La rapidez de flujo de peso, es la masa de fluido que fluye por una secció n por unidad
de tiempo.
 M La rapidez de flujo de masa, es la masa de fluido que fluye por una secció n por unidad
de tiempo.

El má s importante de estos tres es la rapidez de flujo de volumen, con la cual se establece la siguiente relació n:

Q  Av 4.3.1
[m3/s]
donde A es el á rea de la secció n y v es la velocidad promedio del flujo. Las unidades de Q se pueden derivar de la siguiente
manera:

Q  Av  m 2 * m / s  m 3 / s RAPIDEZ DE
FLUJO DE
VOLUMEN
el peso w, esta relacionada con Q y γ, donde las unidades derivan de la siguiente manera:
w  Q  N / m 3 * m3 / s  N / s

RAPIDEZ DE
FLUJO DE PESO
la masa m, esta relacionada con Q y la densidad ρ, y las unidades se derivan de la siguiente manera:
m  Q  kg / m 3 * m 3 / s  kg / s

RAPIDEZ DE
IV.4 TEOREMA DE BERNOULLI-(Ecuació n)
FLUJO DE MASA

Daniel Bernoulli, en 1738, demostró un teorema general, referente al movimiento de los fluidos, que es
probablemente el má s importante de la hidrá ulica. Toda la Hidrodiná mica reposa prá cticamente en dicho teorema; y gran
nú mero de problemas hidrá ulicos pueden ser resueltos con la ú nica ayuda de dicho teorema. Al enunciarlo y demostrarlo,
será mas simple despreciar, al principio, todas las resistencias de fricci ó n e investigar después de la influencia de estas
resistencias en los resultados.
De acuerdo con el teorema de Bernoulli, la carga absoluta en cualquier secció n transversal es igual a la carga
absoluta en una secció n de aguas abajo má s las perdidas de cargas intermedias. La carga de energía en todas las secciones
es constante para corriente uniforme y variable para corriente no uniforme.

151
 Estrictamente hablando, el teorema de Bernoulli solo puede aplicarse a la secció n transversal de una corriente
cuando la direcció n transversal cuando la direcció n de ésta sea horizontal. Su deducció n se basa en la hipó tesis de que todas
las velocidades son iguales y que cada punto de la secció n tiene el mismo contenido de energía.
Para efectos de demostració n, vamos a considerar un flujo ideal, irrotacional, permanente y unidimensional; en el
que, analizaremos una línea de corriente, que coincida con el eje de un filete o tubo de corriente, en la que consideramos un
elemento diferencial de fluido (Fig4-3).

En la vena liquida considerada, ésta tiene un elemento diferencial de á rea ´dA´ y un elemento longitudinal
diferencial ´ds´, por lo que el peso de ella será γdAds. La fuerza que produce la presió n en la base inferior del filete será
PdA y en el incremento ds la fuerza que produce la presió n en la base superior del filete será :

 p 
p ds dA
 s 
Entonces se puede hacer la sumatoria de fuerzas a través del eje ‘s’ considerado, y aplicamos la segunda ley de
Newton:

 p 
F  ma  pdA   p  ds dA  dA cos  dAds(a )
 s 
 p 
pdA  pdA   p  ds dA  dA cos  dAds(a )  0
 s 
p dz
   cos   (a )  0  cos 
s ds

152
 la aceleracio n se puede expresar como :
dV ds ds ds
a a  a V
dt ds dt dt
de lo que se obtiene finalmente :
 V2 
d 
 2 
a
ds
reemplazan do esta ecuacion de la aceleracio n y el coseno equivalente
tendremos :
 V2 
d  
p dz  2 
  0
s ds ds
dividiendo entre el peso especifico del fluido y ordenando
 p V2 
d    z 
 γ 2 g  0
 ds
p V2
  z  ctte Ec. Bernoulli 4.4.1
 2g

La ecuació n obtenida es la ecuació n de Bernoulli; la constante de integració n (conocida como constante de

Bernoulli) generalmente varía de una línea de corriente a otra, pero permanece constante a lo largo de una línea de corriente
en flujo permanente, sin fricció n e incompresible. Estas cuatro suposiciones con necesarias y se deben tener presentes al
aplicar este ecuació n. Cada término tiene dimensiones de (L/T2) o unidades de metros-newtons por kilogramo.
m  N m  kg  m / s 2 m 2
  2
kg kg s

Debido a que 1N=1Kg*m/s2, por consiguiente la Ec-1 se interpreta como energía por unidad de masa.
Puede interpretarse como energía por unidad de peso, metros-newtons por newton. Esta forma es particularmente
conveniente para desarrollar problemas de líquidos con una superficie libre.
Cada uno de los términos de la ecuació n de Bernoulli puede interpretarse como una forma de energía disponible.

Esta ecuació n también se conoce como la ecuació n de conservació n de energía mecá nica. Es particularmente
importante notar que esta ecuació n de energía se dedujo de la ecuació n de momentum.
Las perdidas de energía debida a la fricció n y a la transferencia de calor solamente pueden incorporarse a la
ecuació n diferencial de energía completa.
Al aplicar la ecuació n a dos puntos sobre una línea de corriente se obtiene,

Ec. Bernoulli p1 v12 p v2

z1    z2  2  2 4.4.2
desarrollada  2g  2g

Esta ecuació n (4.4.2) muestra que lo importante es la diferencia en energía potencial, energía de flujo y energía
cinética. Por consiguiente, z1-z2 es independiente del nivel de referencia particular, al igual que la diferencia en la elevació n
de los dos puntos. Similarmente, p1/γ y –p2/γ es la diferencia en las cabezas de presió n, expresada en unidades de longitud
del fluido fluyendo, y no se altera por la presió n de referencia particular seleccionada. Debido a que los términos de
velocidad son no lineales, su nivel de referencia es fijo.

153
 Por tanto los términos de la ecuació n de Bernoulli se conocen, a menudo, como cabezas refiriéndose a una altura
por encima de un nivel de referencia. El termino Z se le llama cabeza de elevació n, y el termino v 2 / 2 g se le conoce
como cabeza de velocidad. La suma de las tres se la conoce como cabeza total. Debido a que cada t érmino representa una
altura, un diagrama como se presenta en la Fig4-4, la cual es la utilidad para visualizar la relaci ó n entre los tres tipos de
energía. Se observa que, debido a la suposició n de que no se pierde o se agrega energía, la cabeza total permanece a un
nivel constante. Entonces, la altura relativa de cada término de cabeza varía segú n lo establecido por la ecuació n de
Bernoulli.

Como se muestra en la Fig4-3, la cabeza de velocidad en la secció n 2 será menor que en la secció n 1. Esto se puede
mostrar mediante la ecuació n de continuidad,

A1v1  A2 v 2
4.4.3
v 2  v1 ( A1 / A2 )

Puesto que A1<A2, y v 2  v1 , como la velocidad está al cuadrado en el término correspondiente a la cabeza de
2 2
velocidad, v2 / 2 g es mucho menor v1 / 2 g . Así mismo cuando se expande el tamañ o de la secció n como pasa en la
Fig4-4, la cabeza de presió n aumenta debido a que disminuye la cabeza de velocidad. Sin embargo, el cambio real también
se ve afectado.
En resumen, la ecuació n de Bernoulli explica el cambio en las cabezas de elevació n, de presió n y de velocidad
entre dos puntos en un sistema de flujo de fluido. Se supone que no existen pérdidas o ganancias de energía entre los dos

154
puntos, de modo que la cabeza total permanece constante. Es esencial que la presió n en los dos puntos de referencia se
expresen, ambas como presiones absolutas o como presiones manométricas. Es decir, deben tener las dos la misma presió n
de referencia.

IV.4.1 RESTRICCIONES A LA ECUACION DE BERNOULLI

Aunque la ecuació n es aplicable a un gran nú mero de problemas, es necesario mencionar ciertas limitaciones con
el fin de aplicar correctamente la ecuació n.
 Está será valida solo para fluidos incompresibles, dado que el peso especifico es el mismo para las
dos secciones.
 No pueden haber dispositivos mecá nicos entre las dos secciones tomadas en cuenta para el aná lisis,
estos pueden agregar o disminuir energía.
 No puede haber transferencia de calor hacia adentro o fuera del fluido.
No debe haber pérdidas de energía debido a la fricció n.
En la realidad ninguno de los sistemas satisfacen todas las restricciones. Sin embargo existen muchos sistemas
para los cuales solamente se tendrá un error pequeñ o y despreciable al aplicar la ecuació n.

PROCEDIMIENTO PARA LA APLICACIÓ N DE LA ECUACION

DE BERNOULLI EN FLUJO PERMANENTE

 Determine qué elementos son conocidos y qué se va a encontrar.

 Decidir cuales dos secciones del sistema se utilizaran cuando se escriba la ecuaci ó n. Se escoge una
secció n de la cual se conocen muchos datos.
La segunda es, por lo general, la secció n en la cual se debe calcular algo.
 Escriba la ecuació n para las dos secciones elegidas en el sistema. Es importante que la ecuació n que la
ecuació n se escriba en la direcció n del flujo. Es decir, el flujo debe ir de la secció n de la parte izquierda
de la ecuació n a la parte derecha.
 Simplifique la ecuació n, si es posible, mediante la cancelació n de los términos cuyo valor es cero o de
aquellos que son los mismos en ambos lados de la ecuació n.
 Resuelva la ecuación algebraicamente para el término deseado.
 Sustituya las cantidades conocidas y calcule el resultado, tome la precaución de asegurar el uso de
unidades consistentes a lo largo del cálculo.

EJEMPLO DE APLICACIÓN La figura 4-4 muestra agua fluyendo a 10 ºC de la sección 1 a la sección 2. En

la sección 1, que tiene 25mm de diámetro, la presión manométrica es de 345 Kpa y la velocidad de flujo es
de 3.0 m/s. La sección 2, que tiene 50mm de diámetro, esta a 2.0 m sobre la sección 1. Suponiendo que no
hay pérdidas de energía en el sistema, calcule la presión P 2.
Solución.

155
D1  25mm
D2  50mm
v1  3.0m / s
p 2  345 Kpa
Z 2  Z 1  2 .0 m
Debe encontrarse la presión P2. Es decir debemos calcular la presión en la sección 2, que es diferente de la
sección 1, debido al cambio de elevación y de área de flujo entre las dos secciones.
Se conoce P1, v1 y z1. Tenemos la ecuación de Bernoulli

P1 v2 P v2
 z1  1  2  z 2  2
 2g  2g
Despejando P2 de la ecuación tenemos:

P2 P1 v2 v2
  z1  1  z 2  2
  2g 2g
P v2 v2 
P2    1  z1  1  z 2  2 
 2g 2g 

Sin embargo, es conveniente agrupar las cabezas de elevación y las cabezas de velocidad. También, como

 ( P1 /  )  P1 , la solución final para P2 será:

 v 2  v22 
P2  P1    z1  z 2  1 
 2 g 

Sabemos que g=9.81Kn/m3, además sabemos que esta fluyendo agua a 10ºC en el sistema, por tanto γ
=9.81KN/m3. Nos queda por determinar v 2 , para esto recurrimos ala ecuación de continuidad vista en el
capitulo anterior.

A1  πD12 / 4  π  25mm  / 4  491mm 2

2
A1v1  A2 v2
A2  ππ 22 / 4  π  50mm  / 4  1963mm 2
2

v2  v1  A1  v1  3.0m/s
 A2 

v2  0.75m/s
Ahora sustituyendo en la ecuación despejada en un principio tendremos:

P2  345Kpa 
9.81Kn 
  2.0m 
 3.0m / s    0.75m / s 
2 2


m 3 
 2 9.81m / s 2  


156
Se sabe que z1  z 2  2.0m , debido a que z2 es mayor a z1
P2  345Kpa  15.4 KN / m 2  345Kpa  15.4 Kpa

P2  329.6 Kpa
la presión encontrada es manométrica debido a que fue calculada con respecto de P 1, que también es una
presión manométrica. En siguientes ejercicios supondremos que estas presiones son manométricas a menos
que se diga otra cosa.

EJEMPLO DE APLICACIÓN Por la tubería mostrada en la figura 4-5, circula agua, siendo la relación entre
el diámetro en el punto 1 y el diámetro en el punto 2 igual a 2.
En 1 la presión es de 0.5Kg/cm y la elevación 100 m. En 2 la presión es 3.38kg/cm 2 y la elevación 70 m.
2

Calcular la velocidad en dichos puntos despreciando las pérdidas por rozamiento.

Solución.

Por continuidad se tiene :

v1  A1  v 2  A2
de donde :
4d 22
v1  v 2 ( )
4d12
simplificando :
2
d2  1 
v1  v 2 ( 22 ) 2  v1  v 2  
d1  2
luego :
1
v1  (1)
2v 2
por Bernuolli :
v12 p v2 p
 1  z1  2  2  z 2
2g w 2g w
trasponiendo : p1
 5m de agua
v2 v1 p p w
 1  2  z1  z 2 (2) p2
2g 2g w w  33.8m de agua
Segun los datos del problema : w
p1  0.5kg / cm 2
p 2  3.38kg / cm 2
z1  100m
z 2  70m

157
 reemplazando estos datos y (1) en (2)
v 22 v2
 2  5  33.8  100  70  1.2
2 g 8g
4v 22  v 22  8 g( 1.2 )
este valor en (1) :
3v 22  9.6 g
9.8  9.6
v2   31.36  5.60 m / seg v2  5.60 m / seg
3
v1  2.80 m / seg
APLICACIÓN DEL TEOREMA DE BERNOULLI EN FLUJO
NO PERMANENTE

Se requiere a veces determinar el tiempo necesario para vaciar un reservorio, o el tiempo empleado para
que el nivel del liquido baje una cierta altura; es por ello que es necesario determinar una expresión matemática que
nos que nos determine el mencionado tiempo. Para efectos didácticos se elige un orificio con carga variable.
Como ejemplo de aplicación de este caso podemos mencionar la aplicación del teorema en el gasto de un
orificio con carga variable.
Sea h1 la carga del orificio en el momento inicial, como se muestra en la figura 4-6, y h 2 la carga, al final
de cierto tiempo “t”. Supondremos como constante la sección transversal horizontal del reservorio.

En cualquier instante determinado, la velocidad del chorro será:

V  Cv 2 gh
donde h es la carga variable en el instante considerado. El gasto que seria descargado por el orificio será:
Q  C d a 2 gh
En un tiempo infinitesimal “dt” se descargará un volumen dVol, por lo que el caudal que desciende en el
reservorio será:
dVol
Q
dt

Pero en ese mismo instante, la carga habrá disminuido una cantidad infinitesimal dh; entonces este caudal
puede ser expresado por la siguiente ecuación:

158
 Adh
Q
dt
igualando las dos anteriores
Adh
  C d a 2 gh
dt
A
 dt   C a 2 gh  h 2 dh
1

si en esta ultima ecuación, integramos en el tiempo de 0 a t e integramos h de h 1 a h2, tendremos entonces:

h1  h2
Tiempo de descarga T  2A 4.4.4
Cd a 2g [s]

EJEMPLO DE APLICACIÓN En los primeros cinco minutos un recipiente que tiene un orificio en el
fondo desciende 1 metro de su nivel inicial y en los siguientes 6 minutos desciende otro metro. Determinar
cual es el nivel inicial del recipiente.

Solución. Se considerará como punto de referencia la altura intermedia del recipiente entre los dos tiempos
dados por el problema

300 
Cd a
2A
2g
 H  1 1
2
H
1
2

360 
Cd a
2A
2g
H 1
2
  H  1
1
2


Dividiendo las dos ecuaciones se obtiene:

11H1/2-5(H-1)1/2= 6(H+1)1/2
despejando esta ecuación obtendremos la altura del nivel inicial:
H  2.82m

IV.5 FLUJO EN ORIFICIOS-TEOREMA DE TORRICELLI

El termino orificio, según se usa en hidráulica, se aplica a cualquier abertura, con perímetro cerrado, practicada en
una pared o tabique que permite el derrame del agua contenida en un recipiente. Los orificios entran en el diseño de muchas
estructuras hidráulicas; y se usan frecuentemente en la medición de caudales de las corrientes fluidas.
Los mas usados son lo circulares, El agua que fluye por orificios de conforme va abandonando el orificio, el chorro
va contrayéndose gradualmente, hasta formar un chorro cuya área transversal es algo menor que el área transversal del
orificio. Esto se debe a la convergencia de las trayectorias seguidas por las diferentes partículas, conforme se acercan al
orificio. Supóngase un deposito de liquido (Fig4-7), este tiene en la parte inferior un orificio por el que sale el liquido; el
área del orificio es pequeño y el de el deposito es suficientemente grande, y siendo el flujo permanente, de manera que el
gasto que sale por el orificio es igual al gasto que entrara en el deposito, por lo que se tendrá una altura del liquido “h”.

159
 Suponiendo el flujo ideal, podemos partir de la ecuación de Bernoulli entre el punto “0” en la superficie del deposito
y el punto de salida del orificio.
P0 V0 P V
  z0  s  s  z s
 2g  2g
Vemos que P0=Ps; V0 =es aproximadamente cero y zs =0, tendremos:
V s2
h
2g
Vs  2g  z0  z s 
Vs  2 gh [m/s] 4.5.1
esta última ecuación (4.5.1) se conoce como Teorema de Torricelli.

Donde la velocidad es igual a la que adquiriría una partícula de fluido al caer libremente desde una altura “h”
(velocidad ideal) y esta es independiente del peso especifico del fluido; con alcohol y mercurio la velocidad seria la misma.

En un orificio situado en un plano horizontal (Fig4-8), todas sus partes estarán bajo las mismas cargas, y la
velocidad de todas las partículas será igual, al llegar a la sección contraída. en este caso el chorro de agua se elevaría hasta
una altura igual a la carga que produce. Por supuesto, que la resistencia del aire impide que esto suceda en la realidad, como
también lo impiden la fricción entre el agua y el orificio, y la fricción de las partículas de agua. Los experimentos han
enseñado que, para cargas bajas (< 2 y 2.5 metros), la discrepancia es muy pequeña, aumentando conforme aumenta la
carga.
Para desarrollar este caso podemos partir de Bernoulli. Primero obtenemos una expresión para la velocidad del
chorro en el punto 2.

160
 P1 v2 P v2
 z1  1  2  z 2  2
 2g  2g
v2  2 gh
ahora se escribe la ecuación de Bernoulli entre los puntos 2 y 3 en el nivel de la superficie libre del fluido:
P2 v2 P v2
 z 2  2  3  z3  3
 2g  2g

Pero P2  P3  0 . Entonces para v3 , tendremos:

v3  v 22  2 g ( z 2  z 3 )
v 22  2 gh
( z 2  z 3 )  h
Tenemos:
v3  2 gh  2 g (  h)  0 [m/s] 4.5.2
este resultado verifica que la corriente alcanza justamente la altura de la superficie libre de fluido en el tanque.
Para alcanzar una altura mayor (fuentes), se puede desarrollar una mayor presión por encima del fluido en el
recipiente, o se puede utilizar una bomba para obtener una mayor presión.
EJEMPLO DE APLICACIÓN Para el tanque de la figura, calcular la velocidad de flujo que sale de la
boquilla para una profundidad de fluido de 3.0 m.
Solución. Esta es una aplicación directa del teorema de Torricelli:
vs  2 gh
vs  2(9.81m / s 2 )(3.0m)
vs  7.67m / s

EJEMPLO DE APLICACIÓN Para el sistema ilustrado en la figura 4-9, calcular la presión de aire
requerida por encima del agua para hacer que el chorro suba 40.0 desde la boquilla. La profundidad es de
6.0 pies.

161
 Solución. Usamos la ecuación de Bernoulli para obtener una expresión para la velocidad de flujo en la
boquilla como función de la presión de aire:
P1 v2 P v2
 z1  1  2  z2  2
 2g  2g
vemos que v1  0 y que P2 =0, tendremos:
v2  2 g  ( P1 /  )  ( z1  z2 )
v2  2 g  ( P1 /  )  h
Esta ecuación es parecida a la del teorema de Torricelli. Por analogía, el sistema presurizado

hará que el chorro se eleve a una altura   P1 /    h . Entonces, en este problema, si deseamos una

altura de 40 pies y h = 6 pies

P1 /   40  6  34 pies
P1   (34 pies)
P1  (62.4lb / pies3 )(34 pies) /(144 pu lg 2 )

P1  14.73lb / pu lg 2 relativa

IV.5.1 TRAYECTORIA DE LA VENA LIQUIDA

En la Figura 4-10 se ilustra una vena que se derrama por un orificio vertical bajo una carga h. La abscisa y la
ordenada de un punto m, situado en la trayectoria del chorro son respectivamente x e y. Si v es la velocidad en la vena
contracta, al final del tiempo t, entonces
x  vt
por la ley de la caída de los cuerpos,

y  1 gt 2
2
igualando a t y reemplazando
2v 2
x2  y 4.5.3
g

Esta es la ecuación de una parábola con su vértice en el orificio, como

v  Cv 2 gh Ec. para calcular la velocidad real de la

partícula fluida.
la ecuación (4.5.3) podrá escribirse en la forma:

x 2  4C v2 hy Ec. de una 4.5.4

parábola

162
IV.5.2 COEFICIENTE DE VELOCIDAD

Este se obtiene generalmente haciendo una serie de mediciones de trayectoria del chorro. Si una partícula sale del
orificio con una velocidad V r, y después de “t” segundos se halla la posición en un punto determinado de la trayectoria de la
vena liquida (4.5.4), de donde tomando las ecuaciones del punto anterior,
vr  C v 2 gh
despejando obtendremos:
g
x
2y 4.5.5
Cv 
2 gh

o Vr ; velocidad real de la partícula fluida

o Cv ; coeficiente de velocidad

IV.5.3 COEFICIENTE DE CONTRACION

Como se comento anteriormente el chorro de liquido al salir del orificio se contrae de manera que el área efectiva
del chorro es ligeramente menor que el área del orificio, entonces podemos hacer una relación de áreas, esta relación se
denomina coeficiente de contracción (Cc), este coeficiente varía ligeramente con el tamaño del orificio y con la carga de
líquido, tienen un valor promedio de 0.62 aproximadamente, para orificios Standard. Su uso principal se halla en la
determinación del área transversal de la vena liquida, en la sección contraída (Fig4-11).

Este coeficiente será:

a
Cc   a  AC c 4.5.6
A

o a ; área de la vena liquida

o A ; área del orificio
o Cc ; coeficiente de contracción

Por relación entre las formulas vistas anteriormente, podemos relacionar directamente los tres coeficientes vistos:
C d  Cc Cv
se muestra a continuación una tabla donde se observan las variaciones de estos coeficientes con el uso de esta ecuación

163
 Cv 1 0,99 0,98 0,97 0,96 0,95

Cc 0,586 0,6 0,615 0,31 0,47 0,664

Cd 0,586 0,594 0,603 0,12 0,621 0,631

Fuente: Introducción a la Hidráulica-

Texto auxiliar de la materia

IV.5.4 COEFICIENTE DE GASTO O DESCARGA

La cantidad o volumen “Vol” que fluye del orificio, por unidad de tiempo “t”, puede expresarse por el producto del
área efectiva “a” de la sección contraída, por la velocidad efectiva “V r”, adquirida por el agua al pasar por dicha sección;
tenemos:

Qr  aVr  ( AC c )(C v 2 gh )  C c C v A 2 gh 4.5.7

Donde CcCv = Cd , por lo que la ecuación (4.5.7) puede escribirse como:

Qr  C d A 2 gh [m3/s] 4.5.8

Qr
Cd  4.5.9
Qt
o Qr ; caudal real
o Qt = A(2gh)1/2 ; caudal teórico

este caudal teórico tendría lugar si no hubiera fricción, ni contracción de la vena liquida.

IV.5.5 PERDIDA DE CARGA EN UN ORIFICIO

A causa del rozamiento y de la viscosidad, la velocidad del agua que se derrama por un orificio es menor que
2 gh , o sea:

v  Cv 2 gh
la carga total que produce el derrame es, por lo tanto,

1 v2 [m]
h 2 4.5.10
Cv 2 g

la carga perdida es igual a la cara total menos la carga de velocidad, o sea, si h o es la carga perdida,

1 v2 v2  1  v2
ho      1 4.5.11
Cv2 2 g 2 g  Cv2  2g [m]

como v  C v vt ; donde vt  velocidad teorica ,

2
ho  (1  C v )h [m] 4.5.12

164
 Si, para un orificio de arista viva, Cv =0.98, se obtiene la siguiente relación:

v2
ho  0.041  0.040h [m] 4.5.13
2g

IV.5.6 ORIFICIO SUJETO A CARGAS BAJAS

Se tiene que, Q  C c A 2 gh , bajo la hipótesis de que la carga que produce el derrame es la que obra en el
centro del orificio. Cuando en un orificio vertical éste es pequeño en relación con la altura del mismo, hay una diferencia

apreciable entre el derrame teórico verdadero y el dado por la ecuación Q  C c A 2 gh .

En la figura 4-12 se ilustra un orificio rectangular de ancho L y altura M. Las cargas respectivas sobre los bordes
superior e inferior del orificio son h 1 y h2. Despreciando la velocidad de llegada o acceso, el derrame teórico a través de una
faja elemental de área Ldy, bajo una carga y, es
dQt  L 2 gy dy 4.5.14
integrada entre los límites h1 y h2,

Qt  2 3 L 2 g (h23 / 2  h13 / 2 ) [m3/s] 4.5.15

Cuando h1 es cero tendremos:

Qt  2
3 L 2 g ( h23 / 2 ) [m3/s] 4.5.16

esta es la formula teórica, sin corrección Por velocidad de acceso, para el derrame de un vertedero. La ecuación (4.5.15) es
el derrame teórico para orificios rectangulares. Para orificios circulares se puede deducir una formula semejante pero mas
complicada.

IV.5.7 ORIFICIOS SUMERGIDOS O AHOGADOS

Si un orificio descarga completamente dentro del agua, recibe el nombre de “sumergido”. Estos orificios son de uso
frecuente en las obras de ingeniería como en las esclusas, canales de escape, compuertas de marea, y muchas otras
construcciones (Fig4-13).

165
 Supóngase dos recipientes como se muestra en la figura 4-12, expuestos a la presión atmosférica. Se tendrá entonces
para los puntos “o” y “s”

Vs2
0  10.33  H 1   10.33  H 2
2g
de donde:
Vs  2g (H1  H 2 )

Vs  2 gH [m/s] 4.5.17

Se demuestra entonces fácilmente que el valor teórico de la velocidad es 2 gH , donde H es la diferencia entre
los niveles de agua. El gasto será, como antes:

Q  C d A 2 gh [m3/s] 4.5.18

IV.5.8 ORIFICIO EN LAS TUBERIAS

En la figura 4-14 se ilustra un orificio circular de arista viva dentro de un tubo. La velocidad en la vena contracta
“2”, se obtiene por la formula,

p p v2 
v2  2 g  1  2  1  h1 
 w w 2g 

reemplazando el término de pérdida de energía, h1, por Cv , esta ecuación se transforma en:

166
 p p v2 
v 2  C v 2 g  1  2  1 
 w w 2g 
2
En los orificios de los tubos no puede omitirse el factor de corrección por velocidad de acceso ( v1 / 2 g ).
Nuevamente, haciendo A2  C c A y C  C c  C v , la expresión para el derrame o gasto se transforma en,

p p v2 
Q  CA 2 g  1  2  1  4.5.19
 w w 2g 

En esta forma la ecuación (4.5.19) requiere una solución por aproximaciones sucesivas. Sin embargo, sustituyendo
v 2
1
por su valor obtenido de la ecuación de continuidad  Q / A1  , donde A , es el área del tubo, se deduce la siguiente
1

expresión:

p p 
2g 1  2 
Q  CA w w
2 4.5.20
2 A 
1  C  
 A1 
El valor de C varia no solamente con los factores que afectan a los orificios de los depósitos, sino también con la
relación del diámetro del orificio, d, al del tubo, d 1, y con el lugar donde se efectúan las tomas de presión. Los experimentos
indican que la vena contracta esta aproximadamente a la mitad del diámetro del tubo, aguas abajo del orificio. Cuando el
orificio esta en el extremo del tubo, p2 es igual a cero en la formula (4.5.20) y solamente se requiere un injerto o toma de
presión. En la figura 4-15 se muestran dos curvas que ilustran la relación entre C y el número de Reynolds. Algunas partes
de la curva aparecen de trazos porque los valores de R son inciertos, debido a los escasos datos de que se dispone en esta
región. En los anexos del presente capitulo se presenta información adicional relativa a los coeficientes de gato o derrame
para las tomas de presión en la vena contracta y para los injertos de brida.

Fuente: Manual de Hidráulica-H.W. King,

E.F. Brater

167
IV.5.9 GASTO EN ORIFICIOS DE PARED GRUESA

Cuando el contorno del orificio en estudio no tiene aristas afiladas, se dice que el orificio es de pared gruesa. En
este tipo de orificios podemos notar que después de la sección contraída, el chorro aun tiene espacio dentro del tubo para
expandirse y llenar la totalidad de la sección (Fig4-16). Ver Anexos

Entre la vena contraída y el final del tubo ocurre un rápido descenso en la velocidad, acompañado de una
turbulencia y una fuerte perdida de energía, en este caso la velocidad será:

V  Cv 2 gh 4.5.21
Donde:

o Cv =0.82 para e/D=3

o Cc =1
o Cd =0.82

Por lo que el gasto es aproximadamente, un tercio mayor que en un orificio de pared delgada. Este se debe a que el
la sección contraída se forma un vació parcial con presión ligeramente menor a la atmosférica, lo que aumenta el valor
efectivo de la carga H.

IV.6 ECUACION GENERAL DE LA ENERGIA

Cuando se desarrollo la ecuación de Bernoulli, y se mencionaron algunas restricciones, por lo que se menciono un
pequeño error despreciable en el resultado obtenido. Sin embargo para un sistema donde existen, algunas perdidas y
adiciones de energía entre las secciones de interés, esto sucede generalmente cuando intervienen dispositivos mecánicos
tales como bombas, válvulas, conectores, etc. Para este tipo de sistemas la ecuación de Bernoulli ya no es valida. Con
respecto al efecto que tienen los dispositivo mecánicos sobre el flujo, pueden clasificarse como dispositivos que añaden
energía o succionan energía (Fig4-17)

168
 Los motores de fluido, turbinas, accionadores giratorios y lineales son ejemplos de dispositivos que toman energía
de un fluido y la transfieren en forma de trabajo, ocasionando la rotación de un eje o el movimiento lineal de un pistón.
Muchos motores de fluido tienen la misma configuración básica que las bombas; la principal diferencia entre una
bomba y un motor de fluido es que cuando se pone en funcionamiento un motor, el fluido pone elementos giratorios en
funcionamiento. Las bombas funcionan a la inversa (Fig4-18). En algunos diseños, como el de los engranajes la bomba
podría actuar como motor.

Un fluido ofrece una resistencia de fricción al fluido. Parte de la energía del sistema se convierte en energía térmica
(calor), el cual se disipa a través de las paredes del conducto en el que el fluido se desplaza. La magnitud de la pérdida de
energía depende de las propiedades del fluido, la velocidad de flujo, el tamaño del conducto, la rugosidad de la pared del
conducto y la longitud del tubo.
Los elementos que controlan la dirección o la rapidez de flujo de un fluido en un sistema, típicamente establecen
turbulencias locales en el fluido, ocasionando que la energía se disipe en forma de calor (Fig4-19), estas perdidas de energía
se presentan siempre que haya una restricción, un cambio de velocidad de flujo o un cambio en su dirección.

En un sistema grande, las pérdidas debidas a la presencia de válvulas y conectores es pequeña en comparación con
las pérdidas por fricción en los conductos, estas son denominadas perdidas menores.

169
IV.6.1 NOMENCLATURA DE PÉRDIDAS Y ADICIONES DE ENERGIA

Se explican las perdidas y las adiciones de energía en un sistema en términos de energía por unidad de peso o de
fluido que fluye en el sistema. A esto también se le conoce como “cabeza”. Como un símbolo para el término cabeza
usaremos la letra h, cuando se hable de pérdidas y adiciones de energía. Específicamente, utilizaremos los términos
siguientes:
hA = energía añadida o agregada al fluido mediante un dispositivo mecánico
como puede ser una bomba.
hR = energía removida o retirada del fluido mediante un dispositivo
mecánico como podría ser un motor de fluido.
hL = perdidas de energía por parte del sistema, debidas a fricción en los
conductos, o perdidas menores debida a la presencia de válvulas y
conectores.

No tomaremos en cuenta los efectos de la transferencia de calor hacia o fuera del fluido, ya que son despreciables en
los tipos de problemas que estamos tratando. La magnitud de las perdidas de energía producidas por muchos tipos de
válvulas y de conectores es directamente proporcional a la velocidad del fluido. Lo anterior se puede expresar de manera
matemática como:
h L  K (v 2 / 2 g ) [m] 4.6.1
El termino K es el coeficiente de resistencia, que por lo general se le encuentra experimentalmente.

IV.6.2 ECUACION DE LA ENERGIA

170
 La ecuación general de la energía, es una expansión de la ecuación de Bernoulli, que hace posible resolver
problemas en los que se presentan pérdidas y adiciones de energía. La interpretación lógica de la ecuación de energía se
puede ver en la figura 4-20, que representa un sistema de flujo.
Los términos E1' y E 2' denotan la energía que posee el fluido por unidad de peso en las secciones 1 y 2,
respectivamente. También se muestra las adiciones, remociones y pérdidas de energía, h A, hR y hL. Para tal sistema, la
expresión del principio de conservación de energía es:
E1'  hA  hR  hL  E 2'
la energía que posee el fluido por unidad de peso es:
P v2
E'  z 4.6.2
 2g

la ecuación 4.6.2 queda entonces:

P1 v2 P v2
 z1  1  hA  hR  hL  2  z 2  2 4.6.3
 2g  2g

Esta es la forma de la ecuación de energía que utilizaremos con más frecuencia en el libro. Al igual que la ecuación
de Bernoulli, cada término de la ecuación representa una cantidad de energía por unidad de peso de fluido que fluye en el
sistema. Las unidades SI típicas son N. m/N o metros. Las unidades en el Sistema Británico de Unidades son lb-pie/lb o
pie.
Es de suma importancia que la ecuación general de energía esté escrita en la dirección de flujo, es decir, desde el
punto de referencia, en la parte izquierda de la ecuación, al punto correspondiente, en el lado derecho. Los signos
algebraicos juegan un papel crítico, debido a que el lado izquierdo de la ecuación establece que un elemento de fluido que
tenga una cierta cantidad de energía por unidad de peso en la sección 1, pueda tener una adición de energía (+h A), una
remoción de energía (-hR) o una pérdida de energía (-h L), antes de que alcance la sección 2. En tal punto contiene una
cantidad diferente de energía por unidad de peso según lo indican los términos de la parte derecha de la ecuación.
Por ejemplo, en la figura 4-20, los puntos de referencia son 1 y 2, y en cada uno de estos se indican las cabezas de
presión, de velocidad y de elevación. Después de que el fluido abandona el punto 1, entra a la bomba, donde se le agrega
energía. Un movilizador principal, que podría ser un motor eléctrico, hace funcionar la bomba y su movilizador transfiere
la energía al fluido (+hA). Entonces el fluido fluye por un sistema de conductos compuesto por una válvula, codos y
tramos de conducto en los que la energía se disipa (-h L). Ante de alcanzar el punto 2, el fluido fluye a través de un motor
de fluido que retira algo de la energía para hacer funcionar un dispositivo externo (-h R). La ecuación general de la energía
toma en cuenta todas esas energías. En caso de no tener alguno de los dispositivos, estos serán cero, y pueden sacarse de
la ecuación.

171
EJEMPLO DE APLICACIÓN De un recipiente grande fluye agua con una rapidez de 1.2 pies 3/s a través
de un sistema de conductos (Fig4-21). Calcular la cantidad total de energía perdida en el sistema debido a
la presencia de la válvula, los codos, la entrada del tubo y la fricción del fluido.

Utilizando un planteamiento similar al usado con la ecuación de Bernoulli, elija las dos secciones de
interés (1 y 2).
Las secciones en las que conocemos la mayor información sobre presión, velocidad y elevación son la
superficie del recipiente y la corriente libre del fluido a la salida del conducto, determinamos cuales
términos de la ecuación de la energía son cero:
P1 v12 P2 v 22
 z1   h A  hR  hL   z2 
 2g  2g
o P1 = 0 ; superficie del recipiente expuesta a la atmósfera
o P2 = 0 ; corriente libre de fluido expuesta a la atmósfera
o V1= 0 ; (aproximado) el área superficial del recipiente es grande
o hA = hR = 0 ; no hay dispositivos mecánicos
Entonces nuestra ecuación se resume a:
v 22
z1  hL  z 2 
2g
Puesto que estamos buscando la perdida de energía total del sistema hallamos h L
v 22
hL  ( z1  z 2 ) 
2g
z1  z 2  25 pies
v 2  Q / A2

Puesto que Q = 1.2pies3/s y el área de un chorro de 3 pulgadas de diámetro es de 0.0491 pies 2,

Q 1.2 pies 3 1
v2   *  24.4 pies / s
A2 s 0.0491 pies 2

172
 v 22 (24.4) 2
  9.25 pies
2 g 2(32.2)
Entonces la cantidad total de energía perdida por el sistema será:

hL  ( z1  z 2 )  v22 / 2 g  25 pies  9.25 pies

hL  15.75 pies

IV.7 FLUJO LAMINAR Y FLUJO TURBULENTO

Cuando se analiza un fluido en una corriente de flujo, es importante ser capaces de determinar el carácter de este
flujo. En algunas condiciones, el fluido parecerá que fluye en capas, de una manera uniforme y regular. Además se puede
observar este fenómeno cuando se abre un grifo de agua lentamente, hasta que el chorro es uniforme y estable. Este tipo de
flujo se lo conoce como flujo laminar.

Dirección del flujo por láminas

Experimento para determinar flujo laminar

En la figura 4-24, se muestra una manera de ver al flujo laminar, en un conducto circular. Anillos concéntricos de
fluido se trasladan siguiendo una trayectoria recta y uniforme. Hay poca mezcla del fluido a través de los “limites” de capa
a capa, conforme el flujo se desplaza por el conducto.

173
 En realidad, el fluido esta conformado por un numero infinito de capas. Otra forma de ver al flujo laminar es en un
tubo de vidrio donde fluye agua, donde se inyecta una partícula o una corriente de otro fluido (Fig4-25), esta se desplazara
en una línea recta no se mezclara con el volumen del fluido siempre y cuando este flujo siga siendo laminar. El flujo
laminar en un canal abierto (Fig4-26), se le llama flujo tranquilo. Un flujo tranquilo sobre un muro aparece como una
lamina uniforme de fluido, lo cal se utiliza a menudo en fuentes de ornato.

Flujo laminar en canal abierto

Experimento para determinar flujo turbulento

El flujo laminar esta gobernado por la ley que relaciona la tensión cortante con la velocidad de deformación angular,
es decir, la tensión cortante es igual al producto de la viscosidad del fluido por el gradiente de las velocidades o bien
  dv / dy . La viscosidad es la magnitud física predominante y su acción amortigua cualquier tendencia a la tubería.
El flujo turbulento es un movimiento no uniforme de las partículas, de forma desordenada en todas las direcciones.
Es imposible conocer la trayectoria de una partícula individualmente (Fig4-27). La tensión cortante en el flujo turbulento
puede expresarse como:
dv
     4.7.1
dy
Donde η(eta) = un factor que depende de la densidad del fluido y de las características del movimiento. El primer
término entre paréntesis (μ) representa los efectos debidos a la turbulencia.
Mediante los resultados obtenidos experimentalmente puede obtenerse la solución de las tensiones cortantes en el
caso de flujos turbulentos. Prandtl sugirió:
2
 dv 
  l 2   4.7.2
 dy 
Para expresar las tensiones cortantes en flujos turbulentos. Esta formula tiene el inconveniente de que la longitud de
mezcla (l) es función de y. Cuanto mayor es y, distancia de la pared a la tubería, mayor es el valor de l. Posteriormente Von
Karman sugirió:

174
  y
   0 1    k 2 2
 dv / dy  4
 d v / dy 2  2
4.7.3
 r0 
Aunque k no es una constante, este número adimensional se mantiene aproximadamente igual a 0.40.

IV.8 NÚMERO DE REYNOLDS

El comportamiento de un fluido, particularmente con respecto a las perdidas de energía, depende bastante de si el
flujo es laminar o turbulento. Por esta razón deseamos tener medios para predecir el tipo de flujo sin necesidad de
observarlo. En efecto, la observación directa es imposible para fluidos que se encuentran en conductos opacos. Se puede
mostrar experimentalmente y analíticamente que el carácter del flujo en un conducto redondo depende de cuatro variables.
Este número viene de la relación entre las fuerzas inerciales (FI) y las fuerzas de viscosidad (FV):

d L
L4

 L 
 L3
d
FI  Ma  L  2   2  V
3
;  dv  d T 2 d L2
d
FV  A    A   L 
T  T T
 dy  L T

L3
FI V T VL VL VL
  Re   4.8.1
FV L2   

T
donde:
L ; longitud característica del campo de flujo, en caso de tuberías el diámetro, en caso de placas paralelas, la
separación entre ellas.
V ; velocidad media del campo de flujo.
ρ ; densidad del fluido.
μ ; viscosidad del fluido.
υ ; viscosidad cinemática del fluido.

IV.8.1 NUMEROS DE REYNOLDS CRITICOS

Para aplicaciones practicas en flujos de conductos, tenemos que si el numero de Reynolds para el flujo es menor que
2000, el flujo será laminar. Tenemos también que si el numero de Reynolds es mayor que 4000, se puede suponer que el
flujo es turbulento. En el intervalo de números de Reynolds comprendido entre 2000 y 4000, es imposible predecir que tipo
de flujo existe; por consiguiente, este intervalo se lo conoce como región critica. Las aplicaciones típicas involucran flujos
que e encuentran bien colocados en el intervalo de los flujos laminares o en el intervalo de los flujos turbulentos, de modo
que la existencia de esta región de incertidumbre no ocasiona gran dificultad.

Mediante una cuidadosa minimización de las perturbaciones externas, es posible mantener un flujo laminar para
números de Reynolds hasta de 50000. Sin embargo, cuando N R es mayor que aproximadamente 4000, una perturbación
menor ocasionaría un cambio súbito en el flujo de laminar a turbulento. Es por esta razón, y por que estamos tratando con
aplicaciones prácticas supondremos:
Si NR < 2000, el flujo es laminar
Si NR > 4000, el flujo es turbulento

EJEMPLO DE APLICACIÓN Determinar si el flujo es laminar o turbulento, si fluye glicerina a 25 ºC en

un conducto cuyo diámetro inferior es de 150 mm. La velocidad promedio de flujo es de 3.6 m/s.
Solución. Se debe evaluar Reynolds,

175
 v  3.6m / s
d d  0.15m
NR 
   1258kg / m3 ( anexos)
  9.60  101 Pa  s (anexos)
Entonces tendremos:
(3.6)(0.15)(1258)
NR   708
9.60  101
debido a que NR = 708, menor que 2000, el flujo será laminar. Obsérvese que todos los términos fueron
convertidos a unidades SI antes de evaluar NR.
IV.9 DINAMIC A DE FLUIDOS REALES

Al incluir el análisis del movimiento de los fluidos, se debe tomar en cuenta varios factores que interviene en el
movimiento de estos fluidos, uno de estos factores es la viscosidad, está influye en la velocidad directamente, de tal manera
que los flujos son rotacionales, por lo tanto las deformaciones mencionadas anteriormente juegan un papel importante en el
estudio de estos.
Los fluidos reales se mueven generalmente bajo dos tipos de flujo, estos son el laminar o viscoso, donde el flujo es
dominado por las acciones de la viscosidad, es decir , esta se mueve en capa paralelas, y el otro tipo de flujo turbulento (se
vieron anteriormente en el punto IV.9 del presente), además de las fuerzas debidas a la viscosidad, actúan otras originadas
por el intercambio aleatorio y permanente de cantidad de movimiento dentro del campo de flujo, produciendo así un
movimiento inestable.

IV.9.1 EFECTO DE LA VISCOSIDAD

Si se considera la deformación de dos fluidos newtonianos diferentes, por ejemplo, glicerina y agua, se encontrará
que se deforman con diferente rapidez para una misma fuerza cortante. La glicerina ofrece mucha mayor resistencia a la
deformación que el agua; se dice entonces que es mucho más viscosa.
En la mecánica de fluidos se emplea muy frecuentemente el cociente de la viscosidad absoluta, “μ”, entre la
densidad, “ρ”. Este cociente recibe el nombre de viscosidad cinemática y se representa mediante el símbolo v. Como la
densidad tiene dimensiones [M/Lt], las dimensiones que resultan para v son [L2/t]. En el sistema métrico absoluto de
unidades, la unidad para v recibe el nombre de stoke = cm 2/s.
La viscosidad es una manifestación del movimiento molecular dentro del fluido. Las moléculas de regiones con alta
velocidad global chocan con las moléculas que se mueven con una velocidad global menor, y viceversa. Estos choques
permiten transportar cantidad de movimiento de una región de fluido a otra. Ya que los movimientos moleculares aleatorios
se ven afectados por la temperatura del medio, la viscosidad resulta ser una función de la temperatura.

176
 Cuando el fluido se mueve, se desarrolla en él una tensión de corte, cuya magnitud depende de la viscosidad del
fluido. La tensión de corte, denota con “τ” (tao), puede definirse como la fuerza requerida para deslizar una capa de área
unitaria de una sustancia sobre otra capo de la misma sustancia. Así pues, “τ” es una fuerza dividida entre un área y puede
medirse en unidades de newtons por metro cuadrado (Cap I). En un fluido como el agua, aceite, alcohol o cualquier otro
liquido común, encontramos que la magnitud de la tensión de corte es directamente proporcional al cambio de velocidad
entre diferentes posiciones del fluido.

Se muestra en la figura 4-28 el concepto de cambio de velocidad en un fluido mediante la exhibición de una capa
delgada del fluido situada entre dos superficies, una de las cuales esta estacionaria, mientras que la otra se esta moviendo.
Una condición fundamental que se presenta cuando tenemos un fluido real en contacto con una superficie frontera,
es que el fluido tiene la misma velocidad que la frontera, el fluido que esta en contacto con la superficie inferior tiene
velocidad cero y el que esta en contacto con la superficie superior tiene velocidad v. Si la distancia entre las dos placas es
pequeña, entonces la rapidez de cambio de velocidad con respecto de la posición “y” es lineal. El gradiente de velocidad es

una medida del cambio de velocidad y se define como v

y .
Conocida como rapidez de corte, el hecho de que la tensión de corte del fluido es directamente proporcional al
gradiente de velocidad puede establecerse, matemáticamente como:

    v  4.9.1
 y 
En la que la constante de proporcionalidad “μ”, se conoce como viscosidad dinámica del fluido. La acción de
revolver hace que se cree un gradiente de viscosidad en el fluido. Se requiere una mayor fuerza para revolver un aceite frió,
que tiene una viscosidad mayor, que la requerida para revolver agua, cuya viscosidad es menor. Esto es una indicación de la
mayor tensión de corte en el aceite frío.

IV.9.2 ECUACION DE NAVIER-STOKES

Navier y Stokes terminaron la deducción relacionando el campo de esfuerzos con la deformación del campo
resultante del campo de velocidad variable en el espacio y tiempo. Aquí se invoca la ley de viscosidad de Stokes, una
generalización de la ley de viscosidad de Newton. Si se supone que el fluido es incompresible, entonces se pueden mantener
las siguientes relaciones:
  
 xx  2  ;  yy  2  ;  zz  2
x y z

177
  u v   w u 
 xy   yx  2 z      ;  xz   zx  2  y     
 y x   x z 
 v w 
 yz   zy  2  x     
 z y 

Por consiguiente, el núcleo de la relación de esfuerzo cortante es la dependencia lineal del esfuerzo cortante con
respecto a la tasa de deformación, en donde el coeficiente de proporcionalidad es el coeficiente de viscosidad de Newton.

Dv  v 
    (v  )v    gh  p   2 v 4.9.2
Dt  t 

en forma de componentes:

D  u u u u  h p   2u  2u  2u 
    u v  w    g     2  2  2 
Dt  t x y z  x x  x y z 

Dv  v v v v  h p   2v  2v  2v 
    u v  w    g     2  2  2 
Dt  t x y z  y y  x y z 

Dw  w w w w  h p  2w 2w 2w 

    u v w    g     2   2 
Dt  t x y z  z z  x y 2 z 

Si todo el movimiento del fluido se detiene (v=0 ; a=0 ) y “z” se selecciona verticalmente hacia arriba en la línea de
acción de la gravedad, la ecuación hidrostática emerge como un caso especial del caso general de un movimiento completo
de fluidos. Ciertamente, el estudio de las ecuaciones de Navier-Stokes se debe concentrar en conjuntos de soluciones
especializadas tales como esta para condiciones de flujo especificas por la siguiente razón. Debido a la no linealidad de
estas ecuaciones, existe una desconcertante variedad de posibles resultados, tanto, que estas ecuaciones nunca han sido
resueltas en forma completa de manera analítica general.

IV.9.3 FLUJO ENTRE PLACAS PARALELAS

La figura 4-29 muestra dos placas paralelas separadas una distancia B e inclinadas un ángulo “α”. En el sentido
normal a la figura, las placas tienen una longitud infinita, entre ellas se mueve un fluido de densidad y viscosidad constante
y el gasto también constante. Las líneas de corriente serán rectas y paralelas a los bordes y entre sí.

178
 Dentro este flujo aislamos un volumen de control diferencial bidimensional de dimensiones ds por dn. El equilibrio
según las líneas de corriente:
 P    
Pdn   P  ds dn  dsdn sen  d    dn ds  0
 s   n 
P  z
 dsdn  sendsdn  dnds  0 ; sen  
s n s
P 
 ( P  z )dsdn  dnds  0
s n
si no existe aceleración normal entonces el valor P  z es constante, si no hay movimiento en el sentido normal no habrá
variación longitudinal del esfuerzo de corte, por lo que la ecuación será:
dh d P
  ; h  z  Altura piezométrica
ds dn 
Integrando la ecuación:
d dh
 dn dn    ds
dn

dh
la expresión no varia con n, entonces se tendrá:
ds
dh
  nC 4.9.3
ds

El esfuerzo cortante varia linealmente con “n”, puesto que el flujo entre las dos placas es simétrico, no queda otra
que las magnitudes de los esfuerzos cortantes en los contornos sean iguales, dado que en el punto medio n=B/2, el esfuerzo
cortante es nulo:
dh B
 0 C
ds 2
dh B
C  
ds 2

Sustituyendo tendremos:
dh  B 
     n 4.9.4
ds  2 

de la ecuación diferencial de la viscosidad y sustituyendo en la ecuación tendremos:

dv
 
dy
dv dh  B 
     n
dn ds  2 
 dh  B 
dv     n dn
 ds  2 
 dh  B 
 dv       n dn
ds  2 
 dh
v
2  ds
 Bn  n 2   C1 4.9.5

179
la constante de integración puede obtenerse de la condición que en los contornos la velocidad de estos y la del fluido es la
misma, es decir, si son estos estacionarios, v= nula para n =0 ó B, por lo que C 1 también será cero, entonces:

 dh
v
2  ds
 Bn  n 2  4.9.6

Esto implica que la ecuación será parabólica.

Como la distribución de velocidades es parabólica, la velocidad máxima se dará en la mitad de la separación “B/2”,
por lo que la ecuación máxima será:

 dh 2
v max  
8 ds
B  4.9.7

en el caso de una distribución parabólica la velocidad media será 2/3 de la máxima:

12 V 12 V
2 2
 dh 2
v
2
v max  
12 ds
B  dh

B 2   dh   ds
1 B
3 ds 2
1

Considerando que la velocidad media es un valor constante:

12 V
 h  h1  h2  4.9.8
B 2
esta ecuación permite calcular la variación de la altura piezométrica a lo largo del conducto. Si el flujo es uniforme
h  la perdida de energía .

IV.9.4 FLUJO EN TUBERIAS CIRCULARES DE SECCION CONSTANTE

Tomando un volumen de control concéntrico (Fig4-30) con el tubo de forma cilíndrica, la ecuación de equilibrio en
el sentido del flujo será:
 P  2
P r 2   P  s  r   r 2 s sen   2 rs  0
 s 

180
 El esfuerzo cortante a una distancia r del centro es constante, por cuanto el flujo es axialmente simétrico, por tanto

sen   dz , por lo que la ecuación se convierte en:

ds
r d r dh
  ( P  z )   
2 ds 2 ds
igualando la ecuación:
dv r dh
  
dn 2 ds
 dh
dv    dn 4.9.9
2  ds

integrando a lo largo de “n”, donde dh es constante y dr  dn .

ds

IV.9.5 TEORIA DE LA CAPA LIMITE

La teoría de la capa limite fue introducida por Prandtl. Esta teoría establece que, para un fluido en movimiento, toda
las perdidas por fricción tienen lugar en una delgada capa adyacente al contorno del solidó (llamada capa limite), y que el
flujo exterior a dicha capa puede considerarse como carente de viscosidad.
La distribución de velocidades en la zona próxima al contorno es influenciada por la tensión cortante en el contorno.
En general, la capa limite es muy delgada en la parte de aguas arriba del contorno y va aumentando su espesor hacia aguas
abajo por la acción continuada de las tensiones cortantes.
Para números de Reynolds bajos, toda capa limite es gobernada por la acción de las fuerzas viscosas y en su interior
el flujo es laminar. Para valores intermedios del numero de Reynolds la caspa limite es laminar cerca de la superficie del
contorno y turbulenta en las zonas algo mas alejadas. Para valores del números Reynolds muy elevados la capa limite es
totalmente turbulenta.

IV.10 TURBULENCIA

El fenómeno turbulento es ocasionado por la inestabilidad del flujo laminar, creando pequeños remolinos que se
mueven de manera aleatoria a lo largo y ancho del campo de flujo. Esta situación ocasiona un cambio constante de la
magnitud y dirección del vector velocidad en cualquier punto. La turbulencia es un intercambio continuo y aleatorio de
masa entre las diferentes zonas del campo de flujo que propicia la mezcla. Esto implica que materia de mayor energía
cinética que pasa por el centro de la tubería pase a las zonas laterales y viceversa ocasionado una mayor uniformidad de las
velocidades promedio en sentido del movimiento general (Fig4-31)

181
 El flujo turbulento representa un incremento sustancial en la perdida de energía. En resumen, la turbulencia se
caracteriza por su condición aleatoria en el tiempo y en el espacio, un rápido proceso de mezcla, la fluctuación
tridimensional de las velocidades y la alta disipación de energía, y por eso un fenómeno controlado por las características
del flujo como por las del fluido. La turbulencia se presenta para números de Reynolds elevados y es un movimiento
macroscópico de pequeños remolinos.
Para la determinación de esfuerzos cortantes en flujo turbulento se parte de la ecuación,
dV
 
dh
pero esta deja de tener validez, por lo que debe definirse como un promedio, pues tiene características aleatorias,
dV
  4.10.1
dh
o η ; viscosidad del remolino
o V ; velocidad promedio

La naturaleza de esta viscosidad de remolino “η”, presenta toda la dificultad del análisis de flujo turbulento, pues
este será función no solo del fluido, sino también de las características del flujo.
Para situaciones intermedias donde la viscosidad y la turbulencia tiene influencia, el esfuerzo cortante se puede
expresar como:
dV
  (   ) 4.10.2
dh
o μ ; viscosidad dinámica ya conocida

IV.11 DISTRIBUCION DE VELOCIDADES

La distribución de velocidades en una sección recta seguirá una ley de variación parabólica en el flujo laminar. La
velocidad máxima tiene lugar en el eje de la tubería y es igual al doble de la velocidad media.
Para el flujo turbulento resulta una distribución de velocidades mas uniforme. A partir de los datos experimentales
de Nikuradse y otros investigadores, se dan a continuación las ecuaciones de los perfiles de velocidades en función de
la velocidad en el eje de la tubería “ v c ” o en función de la velocidad de corte “ v corte ”

Una formula experimental es:

v  vc ( y / ro ) n 4.11.1
donde: n=1/7; para tuberías lisas hasta RE=100000
n=1/8; para tuberías lisas RE= 100000 a 400000

(a) Para tuberías lisas,

v  v corte (5.5  5.75 log yvcorte / v) 4.11.2

(b) Para tuberías lisas (5000 < R E < 3,000.000) y para tuberías rugosas en la zona de exclusiva
influencia de la rugosidad,

(v c  v)  2.5 v 0 /  ln y / r0 4.11.3

182
en función de la velocidad media, Vennard ha sugerido que
V puede escribirse como:
vc

1
(v c  v )  4.11.4
1  4.07 f /8

(c) Para tuberías rugosas,

v  v c (8.5  5.75 log y /  4.11.5
Donde:
 ; rugosidad absoluta de la pared de la tubería.
(d) Para contornos rugosos o lisos,

v V
 2 log y / r0  1.32 4.11.6
V f
donde: f ; Coeficiente de rozamiento de Darcy.

183