Ou 4
Ou 4
Ecuaciones básicas
del flujo de fluidos
(ρu)x (ρu)x+∆x
∆z
∆y
(x, y, z)
∆x
FIGURA 4.1
Región del volumen x y z fijo en el espacio a través del cual el fluido se está moviendo.
[(ρu) x ] [
− ( ρu) x + ∆x ∆y ∆z + ( ρυ ) y − ( ρυ ) y + ∆y ∆x ∆z ]
∂ρ
[ ]
+ ( ρw ) z − ( ρw ) z + ∆z ∆x ∆y = ∆x ∆y ∆z
∂t
(4.1)
∂ρ ∂( ρu) ∂( ρυ ) ∂( ρw )
= − + + = −(∇ ⋅ ρV) (4.3)
∂t ∂x ∂y ∂z
∂ρ ∂ρ ∂ρ ∂ρ ∂u ∂ υ ∂w
+u +υ +w = −ρ + + (4.4)
∂t ∂x ∂y ∂z ∂x ∂y ∂z
Dρ ∂u ∂υ ∂w
= −ρ + + = − ρ (∇ ⋅ V ) (4.5)
Dt ∂x ∂y ∂z
donde Dρ/Dt es la derivada sustancial o la derivada que sigue del movimiento. Ésta es la
velocidad del cambio de la densidad que podría notar un observador si se mueve hacia
la salida a la velocidad del fluido. Observe que ∂ρ/∂t en la ecuación (4.3) es la velocidad
de cambio observada desde un punto fijo. En un estado estacionario ∂ρ/∂t = 0 como se
indica en las ecuaciones (4.3) y (4.4).
∂u ∂υ ∂w
∇⋅V = + + =0 (4.6)
∂x ∂y ∂z
Cuando un fluido pasa por un ducto o una tubería es conveniente expresar la ecuación
de continuidad (y otras ecuaciones) en coordenadas cilíndricas. La relación entre las
coordenadas cilíndricas y las rectangulares se muestra en la figura 4.2a; las que se dan entre
las variables x, y y z en el sistema rectangular, y las del sistema cilíndrico r, θ y z son
sen
(4.7)
a) b)
FIGURA 4.2
a) Coordenadas cilíndricas; b) coordenadas esféricas.
sen sen
sen
sen sen
(4.8)
donde uØ es la velocidad en la dirección axial.
∫
m� = ρ u dS
s
(4.10)
–
La velocidad V también es igual a la velocidad de flujo volumétrico total, dividida entre
el área de la sección transversal del conducto, y de hecho, generalmente se calcula de esta
forma. Se considera como el flujo de volumen en m3/m2 · s o ft3/ft2 · s. Por lo tanto,
q
V = (4.12)
S
donde q es la velocidad de flujo volumétrico.
Para el caso especial importante donde el flujo se lleva a cabo a través de canales de
sección transversal circular
m� = 14 π Da2 ρ a Va = 14 π Db2 ρ b Vb
de la cual
2
ρ a Va Db
=
ρ b Vb Da
(4.14)
†
Si las paredes del conducto son permeables, como en los tubos de membrana estudiados en el capítulo
26, estas ecuaciones no son aplicables.
Velocidad de masa
A partir de la ecuación (4.9) se escribe
m�
Vρ = ≡G (4.13)
S
Esta ecuación define la velocidad de masa G, que se calcula dividiendo la velocidad de
flujo de masa, entre el área de la sección transversal del canal. En la práctica, la velocidad
de masa se expresa en kilogramos por metro cuadrado por segundo, libras por pie cua-
drado por segundo, o libras por pie cuadrado por hora. La ventaja de utilizar G consiste
.
en que es independiente de la temperatura y la presión cuando el flujo es estacionario (m
constante) y la sección transversal no varía (S constante). Este hecho es especialmente
–
útil, cuando se consideran fluidos compresibles, en los cuales tanto V como ρ varían
con la temperatura y la presión. Ciertas relaciones aparecen después en este libro en las
–
que interviene el producto de V y ρ, de tal manera que la velocidad de masa representa
el efecto neto de ambas variables. La velocidad de masa G se describe como la densi-
dad de corriente o flujo de masa, donde el flujo se define generalmente como cualquier
magnitud que pasa a través de una unidad de área en la unidad de tiempo. La velocidad
–
promedio V , tal como se muestra en la ecuación (4.10), se expresa como el flujo de
volumen del fluido.
EJEMPLO 4.1 A través del sistema de tuberías que se muestra en la figura 4.3 fluye el
petróleo crudo, que tiene una gravedad específica = 0.887 a 60 °F/60 °F. La tubería A es de
2 in. (50 mm) Norma 40, la tubería B es de 3 in. (75 mm) Norma 40, y cada una de las tuberías
C tiene 1–12 in. (38 mm) Norma 40. A través de cada una de las tuberías C fluye una cantidad
igual de líquido. El flujo a través de la tubería A es de 30 gal/min (6.65 m3/h). Calcule: a) la
velocidad de flujo de masa en cada tubería, b) la velocidad lineal promedio en cada tubería,
y c) la velocidad de masa en cada tubería.
Solución Las dimensiones y las áreas de las secciones transversales de la tubería estándar
se dan en el apéndice 3. Las áreas de las secciones transversales que se requieren son, para
la tubería de 2 in., 0.0233 ft2; para la tubería de 3 in., 0.0513 ft2; y para la tubería de 1–12 in.,
0.01414 ft2.
a) La densidad del fluido es
Puesto que hay 7.48 gal en 1 ft3 (apéndice 1), la velocidad de flujo volumétrico total es
30 × 60
q= = 240.7 ft 3 /h
7.48
La velocidad de flujo de masa a través de cada tubería C es la mitad del total o 13 300/2 =
6 650 lb/h (0.8379 kg/s).
B 1
12 in.
A
2 in. C
3 in. 1
12 in. C
FIGURA 4.3
Sistema de tuberías para el ejemplo 4.1.
240.7
VA = = 2.87 ft/s
3 600 × 0.0233
a través de la tubería B es
240.7
VB = = 1.30 ft/s
3 600 × 0.0513
240.7
VC = = 2.36 ft/s
2 × 3 600 × 0.01414
GA =
13 300
0.0233
(
= 571 000 lb/ft 2 ⋅ h 744 kg/m 2 ⋅ s )
a través de la tubería B es
GB =
13 300
0.0513
(
= 259 000 lb/ft 2 ⋅ h 351 kg/m 2 ⋅ s )
y a través de cada una de las tuberías C es
GC =
13 300
2 × 0.01414
(
= 470 000 lb/ft 2 ⋅ h 637 kg/m 2 ⋅ s )
EJEMPLO 4.2 El aire a 20 °C y a 2 atm de presión absoluta entra a un calentador de vapor con
aletas a través de un tubo de 50 mm a una velocidad promedio de 15 m/s. Sale del calentador
a través de un tubo de 65 mm a 90 °C y a 1.6 atm de presión absoluta. ¿Cuál es la velocidad
promedio del aire a la salida?
Las densidades a la entrada y a la salida son entonces ρa = Mpa /(RTa ) y ρb = Mpb /(RTb). Por
lo tanto
ρa paTb
=
ρb pbTa
V ρ D2 V p T D2
VbV=ρ aD2a 2a V= p aT aD2b 2a
Vb = a aρ aD = a ap bT aD
b b b a b
ρb Db2 pb Ta Db2
15 × 2 × 02.052 × 363.16
= × 2 × 0.05 ×2363.16
15
= 1.6 × 0.065 × 293.16
1.6 × 0.0652 × 293.16
= 13.74 m/s
= 13.74 m/s
(τzx )z+∆z
z
(τyx )y+∆y
(τyx )y
y
FIGURA 4.4
Elemento del volumen x y z
con flechas que indican la direc- (x, y, z)
ción en la cual el componente x (τzx )z
del momento se transporta a través
de las superficies. x
En la figura 4.4 el fluido está moviéndose a través de las seis caras del elemento de
volumen en cualquier dirección arbitraria. Como la velocidad es un vector, la ecuación
(4.16) tiene componentes en cada una de las direcciones coordenadas x, y, z. Primero
considere sólo el componente x de cada término en la ecuación (4.16); los componentes
y y z se tratan en forma análoga.
Considere las velocidades de flujo del componente de momento x dentro y fuera del
elemento de volumen de la figura 4.4. El momento entra y sale del elemento de volumen
en parte por convección del flujo del fluido global, y en parte por la acción viscosa como
resultado de los gradientes de velocidad. La velocidad de entrada por convección del
componente de momento x en la cara x es (ρuu)x y z, y la velocidad de salida en x
x es (ρuu)x+∆x y z. La entrada en y a través de la cara x z tiene una velocidad (ρυu)y
x z. Se escriben expresiones similares para las otras tres caras. Por consiguiente, existe
flujo convectivo del momento x a través de las seis caras del elemento. El flujo convectivo
neto dentro del elemento de volumen es
[ ] [
∆y ∆z ( ρuu) x − ( ρuu) x + ∆x + ∆x ∆z ( ρυu) y − ( ρυu) y + ∆y ]
[
+ ∆x ∆y ( ρwu) z − ( ρwu) z + ∆z ] (4.17)
[ ]
∆y ∆z (τ xx ) x − (τ xx ) x + ∆x + ∆x ∆z τ yx
( ) y − (τ yx ) y + ∆y
[
+ ∆x ∆z (τ zx ) z − (τ zx ) z + ∆z ] (4.18)
Entonces τxx es el esfuerzo normal sobre la cara x, y τyx es la tensión tangencial o esfuerzo
cortante dirigida hacia x, sobre la superficie de la cara y, resultante de las fuerzas vis-
cosas. Se desarrollan expresiones similares para el flujo de momento en las direcciones
y y z. Los esfuerzos cortantes resultan de la deformación del elemento de volumen; la
tensión normal está relacionada principalmente con el cambio en ∂u/∂x y la dilatación
del elemento.2a, 5
En muchos casos, las fuerzas importantes que actúan sobre el sistema se originan
de la presión del fluido p y la fuerza gravitacional por unidad de masa g. La resultante de
las fuerzas en la dirección x es
∆y ∆z( px − px + ∆x ) + ρg x ∆x ∆y ∆z (4.19)
∂ ∂ ∂ ∂
ρu = − ρuu + ρυu + ρwu
∂t ∂x ∂y ∂z
∂ ∂ ∂ ∂p
− τ xx + τ yx + τ zx − + ρg x (4.20)
∂x ∂y ∂z ∂x
Esta ecuación se reordena con la ayuda de la ecuación de continuidad [ecuación
(4.3)] para dar
Du ∂p ∂τ xx ∂τ yx ∂τ zx
ρ =− − + + + ρg x (4.21)
Dt ∂x ∂x ∂y ∂z
Se derivan ecuaciones similares para los componentes y y z. Considerando los tres com-
ponentes vectoriales se obtiene
DV
ρ = −∇p − [∇ ⋅ τ ] + ρ g (4.22)
Dt
Las tensiones (o esfuerzos) en cualquier punto dependen de los gradientes de veloci-
dad y de las propiedades reológicas del fluido. Por ejemplo, para los fluidos newtonianos,
los componentes en la dirección x del tensor de tensión son2a
∂u 2
τ xx = −2 µ + µ − κ (∇ ⋅ V ) (4.23)
∂x 3
∂u ∂υ
τ xy = τ yx = − µ + (4.24)
∂y ∂x
∂u ∂υ ∂w
∇⋅V = + + (4.25)
∂x ∂y ∂z
donde κ es la viscosidad global. Hay mucha incertidumbre acerca del valor de κ. Para los
gases monoatómicos, es cero y probablemente de menor importancia para los líquidos
y gases densos.2a Las ecuaciones para las tensiones en las direcciones y y z (y todas las
ecuaciones en los sistemas coordenados cilíndricos y esféricos) están dadas por Bird,
Stewart y Lightfoot.2c
Las ecuaciones generales del movimiento para un fluido newtoniano cuando varían
la densidad y la viscosidad se ejemplifican en la siguiente ecuación para la dirección x,
obtenida por sustitución de las ecuaciones (4.23) y (4.24) en la ecuación (4.21)
Du ∂p ∂ ∂u 2
ρ =− + 2µ − µ − κ ( ∇ ⋅ V )
Dt ∂x ∂x ∂x 3
∂ ∂u ∂υ
+ µ +
∂y ∂y ∂x
∂ ∂w ∂u
+ µ + + ρg x (4.26)
∂z ∂x ∂z
Dυ ∂p ∂ ∂υ ∂u
ρ =− + µ +
Dt ∂y ∂x ∂x ∂y
∂ ∂υ 2
+ 2µ − µ − κ ( ∇ ⋅ V )
∂y ∂y 3
∂ ∂w ∂υ (4.27)
+ µ + + ρg y
∂z ∂y ∂z
Dw ∂p ∂ ∂w ∂u ∂ ∂w ∂υ
ρ =− + µ + + µ +
Dt ∂z ∂x ∂x ∂z ∂y ∂y ∂z
∂ ∂w 2
+ 2µ − µ − κ (∇ ⋅ V ) + ρgz (4.28)
∂z ∂z 3
∂u ∂u ∂u ∂u
ρ + u +υ +w
∂t ∂x ∂y ∂z
∂ 2 u ∂ 2 u ∂ 2 u ∂p
=µ 2 + 2 + 2 − + ρg x (4.29)
∂x ∂y ∂z ∂x
∂υ ∂υ ∂υ ∂υ
ρ +u +υ +w
∂t ∂x ∂y ∂z
∂ 2υ ∂ 2υ ∂ 2υ ∂p
=µ 2 + 2 + 2− + ρg y (4.30)
∂x ∂y ∂z ∂y
∂w ∂w ∂w ∂w
ρ +u +υ +w
∂t ∂x ∂y ∂z
∂ 2 w ∂ 2 w ∂ 2 w ∂p
=µ 2 + 2 + 2− + ρgz (4.31)
∂x ∂y ∂z ∂z
∂u r ∂u r u θ ∂u r u2 ∂u r
ρ + ur + − θ +u
∂t ∂r r ∂θ r ∂z
∂ 1 ∂(r u r ) 1 ∂2ur 2 ∂u θ ∂2ur ∂ρ
=µ + − + − + ρgr (4.33)
∂r r ∂r r 2 ∂θ 2 r 2 ∂θ ∂z 2 ∂r
∂u θ ∂u θ u θ ∂u θ ur u θ ∂u θ
ρ + ur + − +u
∂t ∂r r ∂θ r ∂z
∂ 1 ∂(ru θ ) 1 ∂2uθ 2 ∂u r ∂2uθ 1 ∂ρ
=µ + − 2 + − + ρgθ
∂r r ∂r r ∂θ
2 2 r ∂θ ∂z 2 r ∂θ
(4.34)
∂u ∂u u θ ∂u ∂u
ρ + ur + +u
∂t ∂r r ∂θ ∂z
1 ∂ 1 ∂u 1 ∂2u ∂2u ∂ρ
=µ + + 2 − + ρgz
r ∂r r ∂r r 2 ∂θ 2 ∂z ∂z (4.35)
∂u r ∂u r u θ ∂u r u φ ∂u r u 2θ + u 2φ
ρ + ur + + −
∂t ∂r r ∂θ r sen θ ∂φ r
2 2 ∂u θ 2 2 ∂u φ ∂p
= µ ∇ 2ur − ur − 2 − 2 u θ cot θ − 2 − + ρgr
r2 r ∂θ r r senθ ∂φ ∂r (4.36)
∂u θ ∂u θ u θ ∂u θ u φ ∂u θ ur u θ u 2φ cot θ
ρ + ur + + + −
∂t ∂r r ∂θ r senθ ∂φ r r
2 ∂u r uθ 2 cos θ ∂u φ 1 ∂p (4.37)
= µ ∇ 2uθ + − 2 2 − 2 2 − + ρgθ
r 2 ∂θ r sen θ r sen θ ∂φ r ∂θ
∂u φ ∂u φ u θ ∂u φ u φ ∂u φ u φ ur uθ uφ
ρ + ur + + + + cot θ
∂t ∂r r ∂θ r senθ ∂φ r r
uφ 2 ∂u r 2 cos θ ∂u θ 1 ∂p
= µ ∇ 2uφ − + 2 + 2 2 − + ρgφ
r 2 sen2 θ r senθ ∂φ r sen θ ∂φ r senθ ∂φ
(4.38)
1 ∂ ∂ 1 ∂ ∂ 1 ∂2
∇2 = r2 + sen θ +
r 2 ∂r ∂r r senθ ∂θ
2 ∂θ r sen2 θ
2 ∂φ 2
McCabe 04.indd 82 14/12/06 23:37:20
r ∂θ r sen θ r sen θ ∂φ r ∂θ
∂u φ ∂u φ u θ ∂u φ u φ ∂u φ u φ ur uθ uφ
ρ + ur + + + + cot θ
∂t ∂r r ∂θ r senθ ∂φ r r
uφ 2 ∂u r 2 cos θ ∂u θ 1 ∂p
= µ ∇ 2uφ − + 2 + 2 2 − + ρgφ
r 2 sen2 θ CAPÍTULO 4 Ecuaciones básicas del∂φ
r senθ ∂φ r sen θ ∂φ r senθ flujo de fluidos 83
donde
1 ∂ ∂ 1 ∂ ∂ 1 ∂2
∇2 = r2 + sen θ + (4.39)
r 2 ∂r ∂r r 2 senθ ∂θ ∂θ r 2 sen2 θ ∂φ 2
Ecuación de Euler. Para una densidad constante y viscosidad cero, como es el caso
del flujo potencial, se utiliza la ecuación del movimiento, conocida como la ecuación de
Euler, que es la siguiente:
DV
ρ = −∇p + ρ g (4.40)
Dt
EJEMPLO 4.3 Un fluido newtoniano está confinado entre dos placas verticales paralelas
anchas, separadas por una distancia B, como se muestra en la figura 4.5. La placa sobre la
izquierda es estacionaria; y sobre la derecha está moviéndose verticalmente en forma ascen-
dente a una velocidad constante υ0. Suponga que el flujo es laminar y obtenga la ecuación
para el perfil de la velocidad en el estado estacionario del fluido.
dυ x dp (4.42)
− + ρg = C1
dx µ dy
Integrándola de nuevo se obtiene
x 2 dp
υ− + ρg = C1 x + C2 (4.43)
2 µ dy
v0
FIGURA 4.5
x Flujo entre las placas verticales del ejemplo 4.3.
υ=−
1 dp
2 µ dy
( )
+ ρg Bx − x 2 + υ 0
x
B
(4.44)
Flujo de Couette
Cuando las placas en un sistema como el del ejemplo 4.3 son horizontales (o en cual-
quier situación donde la fuerza de gravedad sea despreciable), la velocidad del fluido
varía linealmente con la distancia de la placa estacionaria, y el gradiente de velocidad
es constante. La viscosidad entonces está relacionada con la tensión de corte (o esfuerzo
cortante) Fs /A por la ecuación siguiente:
Fs B
µ= (4.45)
Aυ 0
donde A es el área de cada placa. El fluido bajo estas condiciones se conoce como flujo
de Couette.
La capa de flujo con una superficie libre se estudia en la siguiente sección.
∑
� �
F = Mb − M a (4.46)
Ma Mb
a)
FIGURA 4.6 x dirección
Balance de momento. b)
Para encontrar β para cualquier situación de flujo dado, se debe conocer la variación de
u con la posición en la sección transversal.
Por lo tanto, la ecuación (4.46) se escribe de la manera siguiente:
∑ F = m� (β V b b − β a Va ) (4.51)
Al utilizar esta expresión, hay que tener cuidado en identificar e incluir en Σ F todos los
componentes de las fuerzas que actúan sobre el fluido en la dirección del componente de
velocidad en la ecuación. Quizá intervengan varias de estas fuerzas: 1) cambio de presión
en la dirección de flujo; 2) esfuerzo cortante en el límite entre la corriente de fluido y
el conducto, o bien (si el conducto en sí mismo se considera como parte del sistema),
las fuerzas externas que actúan sobre la pared sólida; 3) si la corriente está inclinada, el
componente apropiado de la fuerza de gravedad. Si se considera el flujo unidimensional
en la dirección x, una situación típica se representa por la ecuación siguiente:
∑F = p S a a − pb Sb + Fw − Fg (4.52)
τA
φ
r
Fg
δ
jo
Flu
FIGURA 4.7
Fuerzas que actúan sobre el elemento de líquido
en el flujo de capa.
extremos del volumen son iguales y en dirección opuesta. Por consiguiente desaparecen.
También, por suposición, el corte en la superficie superior del elemento es despreciable. Las
dos fuerzas restantes son entonces la fuerza de corte sobre la superficie inferior del volumen
de control y el componente de gravedad en la dirección del flujo. Entonces
Fg cos φ − τ A = 0 (4.53)
µ ∫
rdr rdr
0 0 δ δ
ρg cos φ
u = u = ρg(cos
δ −φ rδ ) − r )
2µ 2µ (
2 22 2 (4.56)
donde δ es el espesor total de la capa líquida. La ecuación (4.56) muestra que en el flujo
laminar sobre una placa, la distribución de la velocidad es parabólica.
Ahora considere un elemento diferencial del área de la sección transversal dS, donde
.
dS = b dr. La velocidad de flujo de la masa diferencial dm a través de este elemento es
igual a ρub dr. La velocidad de flujo de la masa total del fluido es entonces
δ
m� = ∫ ρub dr
0
(4.57)
m� δ 3 ρ 2 g cos φ
= =Γ (4.58)
b 3µ
.
donde Γ ≡ m/b y se denomina la carga del líquido. Las unidades de Γ son los kilogramos
por segundo por metro de ancho o libras por segundo por pie de anchura.
Reordenando la ecuación (4.58) se obtiene, para el espesor de la capa,
1/ 3
3µΓ
δ = 2 (4.59)
ρ g cos φ
El número de Reynolds para el flujo debajo de una placa plana está definido por la
ecuación
4rH V ρ m� ρ 4Γ
Re = = 4δ = (4.60)
µ ρL pδ µ µ
donde rH = radio hidráulico, definido por la ecuación (5.10). Para el flujo de un líquido
debajo de la entrada o salida de la tubería, en general el espesor de la capa es una fracción
muy pequeña del diámetro de la tubería y el número de Reynolds es el mismo que para
una placa plana, como se indica en la ecuación (4.60).
La ecuación (4.59) para el espesor de una película laminar de caída fue presentada
primero por Nusselt,4 quien usó el resultado para predecir los coeficientes de transfe-
rencia de calor para la condensación de vapores. Mediciones del espesor de la película
en una superficie vertical (cos φ = 1) muestran que la ecuación (4.59) es aproximada-
mente correcta para el Re ≈ 1 000, pero el espesor realmente varía alrededor de 0.45
de la potencia del número de Reynolds, y las capas son más delgadas que las previstas
a un valor inferior de Re y más gruesas a un Re superior a 1 000. Las desviaciones
se deben a ondulaciones en las películas, las cuales pueden aparecer incluso a números de
Reynolds muy bajos.
La transición del flujo laminar al turbulento no se detecta con tanta facilidad como
el flujo en la tubería, ya que la película es muy delgada y las ondulaciones hacen difícil
observar turbulencia en la película. Un número de Reynolds crítico de 2 100 ha sido
utilizado a menudo para el flujo de capa, pero las mediciones6 del espesor de la película
indican una transición al Re ≈ 1 200. Por encima de este punto, el espesor aumenta al-
rededor de 0.6 la potencia de la velocidad del flujo.
Fθ r = T (4.61)
V
ur
P uθ
r2
Q
Dirección
O r1 de la rotación
FIGURA 4.8
Momento angular del líquido en movimiento.
Suponga que la figura 4.8 representa parte del impulsor de una bomba centrífuga o
.
turbina a través de la cual circula un fluido a una velocidad de masa constante m . Entra
por el punto Q cerca del centro de rotación a una distancia radial r1 desde el punto O
y sale a una distancia radial r2. Sus velocidades tangenciales en estos puntos son uθ1 y
uθ2, respectivamente. La fuerza tangencial Fθ que actúa sobre el fluido en el punto P, es
proporcional a la variación de la velocidad del momento angular del fluido; por lo tanto,
a partir de la ecuación (4.61), el momento de torsión viene dado por la relación
Una ecuación que describe las interconversiones de la energía que ocurren en un fluido
en movimiento, puede derivarse para formar el producto escalar de la velocidad local V
con la ecuación del movimiento. Al aplicar este procedimiento a la ecuación (4.22) se
obtiene una ecuación general que establece que la velocidad de aumento en la energía
cinética por unidad de masa es igual a la velocidad neta de entrada de la energía cinética
por convección menos lo siguiente: 1) velocidad de trabajo hecha por la presión de los
alrededores; 2) velocidad de la conversión reversible para la energía interna; 3) velocidad
de trabajo efectuada por las fuerzas viscosas; 4) conversión irreversible para la energía
interna; 5) velocidad de trabajo hecho por la gravedad (este término puede ser positivo
o negativo).2b Como se estudiará más adelante, también es posible adicionar la energía
mecánica al fluido mediante una bomba o soplador.
Aquí las derivaciones son inicialmente restringidas al flujo de fluidos unidireccional
de densidad constante y viscosidad cero, utilizando la ecuación de Euler.
∂u ∂u ∂p
ρu + u = −u + ρugx
∂t ∂x ∂x
o ρ
(
∂ u 2 /2 )
+u
(
∂ u 2 /2 ) = −u ∂p + ρug (4.64)
x
∂t ∂x ∂x
Ésta es la ecuación de la energía mecánica para el flujo potencial unidireccional de fluidos
de densidad constante cuando la velocidad de flujo varía con el tiempo.
Considere ahora un elemento de volumen de un tubo de corriente dentro de una
corriente mayor de fluido, que circula con flujo estacionario, como se muestra en la figura
4.9. Suponga que la sección transversal del tubo aumenta continuamente en la dirección
del flujo, y que el eje del tubo es recto e inclinado hacia arriba formando un ángulo φ con
la vertical. Represente la presión, la velocidad del fluido y la elevación a la entrada por
pa, ua y Za, respectivamente, y sean las correspondientes magnitudes a la salida pb, ub y Zb .
Considere el eje x paralelo al eje del tubo.
pb
gx
c ción φ g
Dire x
pa e l flujo
d xb
Zb
Za
FIGURA 4.9
Flujo potencial a través de un tubo de corriente xa
inclinado.
u
(
d ρu 2 /2 ) + u dp + ρug cos φ = 0 (4.65)
dx dx
Entonces para el flujo estacionario es posible dividir entre la velocidad u. Al hacer esto y
también al dividir entre ρ y sustituir por cos φ, la ecuación (4.65) se convierte en
(
d u 2 /2 ) + 1 dp + g dZ = 0 (4.66)
dx ρ dx dx
La ecuación (4.66) es el punto formado de la ecuación de Bernoulli sin fricción.
Aunque derivadas para el caso especial de una sección transversal en expansión y con
flujo ascendente, esta ecuación es aplicable a secciones transversales de contracción
o constantes y flujo horizontal o descendente (el signo de la diferencial dZ corrige el
cambio en la dirección).
Cuando la sección transversal y la densidad son constantes, u no cambia con la
posición, el término d(u2/2)/dx es cero, y la ecuación (4.66) se vuelve idéntica a la ecua-
ción (2.2) para un fluido estacionario. Entonces, en el flujo potencial unidireccional a
una velocidad constante, la magnitud de la velocidad no afecta a la caída de presión en
el tubo; la caída de presión depende sólo de la velocidad del cambio de elevación. En
consecuencia, en un tubo recto horizontal, no hay caída de presión en el flujo potencial
de velocidad constante estacionario.
Integrando la ecuación (4.66) sobre el sistema mostrado en la figura 4.9 da
pa u2 p u2
+ gZa + a = b + gZb + b (4.67a)
ρ 2 ρ 2
En unidades fps,
pa gZa ua2 p gZ u2
+ + = b + b + b (4.67b)
ρ gc 2 gc ρ gc 2 gc
La ecuación (4.67) es conocida como la ecuación de Bernoulli sin fricción. Es una
forma particular de un balance de energía mecánico, pero debido a las condiciones
especiales que permiten dividir la ecuación (4.65) entre la velocidad u para formar la
ecuación (4.66), la ecuación (4.67) también puede derivarse de un balance de concha
de momento en el elemento de la figura 4.9. Sin embargo, esto no es cierto para las
formas más completas de la ecuación de Bernoulli que se estudiarán más adelante en
este capítulo.
Cada término en la ecuación (4.67) es escalar y tiene las dimensiones de la energía
por unidad de masa, lo que representa un efecto de la energía mecánica basado en una
unidad de masa del fluido en movimiento. Los términos gZ y u2/2 son la energía poten-
cial y la cinética, respectivamente, de una unidad de masa de fluido; y p/ρ representa el
trabajo mecánico realizado sobre el fluido por las fuerzas, externas a la corriente, que
lo empujan dentro del tubo o por el trabajo recuperado del fluido que sale del tubo. La
ecuación (4.67) muestra que en la ausencia de fricción, cuando se reduce la velocidad
u, la altura Z o la presión p, o ambas, se incrementan. Cuando la velocidad aumenta, lo
hace sólo a costa de Z o p. Si la altura se modifica, tiene que haber una compensación
mediante una variación de la presión o la velocidad.
La ecuación de Bernoulli tiene un amplio intervalo de validez que la que se desprende
de su deducción. Por otra parte, aunque en la deducción se ha hecho la suposición de
que el tubo de corriente es recto, el principio de la conservación de la energía permite
ampliar la ecuación al flujo potencial en tubos de corriente curvos. Si el tubo es curvo,
la dirección de la velocidad cambia y en la ecuación de Bernoulli se utiliza un valor
escalar de la velocidad, en vez del vector velocidad. En todas las situaciones reales hay
algunas pérdidas por fricción en el fluido y algunas variaciones de velocidad dentro de
una sección transversal del tubo, pero en algunos casos son suficientemente pequeñas
para ser ignoradas. En otras situaciones, al emplear factores de corrección, puede modifi-
carse la ecuación para su utilización en el flujo de capa límite, donde existen variaciones
de velocidad dentro de una sección transversal y se producen efectos de fricción. Estas
correcciones se estudian en las secciones siguientes.
Para aplicar la ecuación de Bernoulli a un problema específico, es esencial identificar
la línea de corriente o el tubo de corriente y elegir unos puntos definidos de corriente de
salida o entrada. Los puntos a y b se seleccionan con base en su conveniencia y, por lo
general, se toman en localizaciones donde se dispone de la mayor información acerca de
las presiones, velocidades y alturas.
EJEMPLO 4.4 Por el fondo de un gran tanque abierto se vacía salmuera, con gravedad
específica = 1.15 a 60 °F/60 °F a través de una tubería de 50 mm. La tubería de drenado
termina en un punto que está a 5 m por debajo de la superficie de la salmuera en el tanque.
Considerando una línea de corriente que inicia de la superficie de la salmuera en el tanque y
pasa a través del centro de la línea de drenado hasta el punto de descarga, y suponiendo que
la fricción a lo largo de la línea de corriente es despreciable, calcule la velocidad de flujo a lo
largo de la línea de corriente al punto de descarga de la tubería.
5g = 12 ub2
La velocidad total de flujo de masa está dada por las ecuaciones (4.8) y (4.13), y la ener-
gía cinética por unidad de masa de fluido en movimiento, la cual reemplaza a u2/2 en la
ecuación de Bernoulli, es
∫ u dS = ∫ u dS
1 3 1 3
�
Ek 2
S
2
S
= (4.69)
m�
∫ u dS
S
VS
�
αV 2 Ek
≡ � = S ∫
u 3 dS
2 m 2VS
∫ u dS
3
S
α= (4.70)
V 3S
Si se conoce α, puede utilizarse la velocidad promedio para calcular la energía cinética
–
empleando α V 2/2 en lugar de u2/2. Para calcular el valor de α a partir de la ecuación
(4.70), debe conocerse la velocidad local, como una función de la localización en la sec-
ción transversal, de tal manera que puedan evaluarse las integrales de dicha ecuación. El
mismo conocimiento de la distribución de la velocidad es necesario para calcular el valor
–
de V por medio de la ecuación (4.11). Como se indica en el capítulo 5, α es 2.0 para el
flujo laminar y cerca de 1.05 para el flujo altamente turbulento.
pa α V2 p α V2
+ gZa + a a = b + gZb + b b + h f (4.71)
ρ 2 ρ 2
Las unidades de hf y las de los términos restantes en la ecuación (4.71) son la energía por
unidad de masa. El término hf representa toda la fricción que se produce por unidad de
masa del fluido (y por consiguiente, toda la conversión de la energía mecánica en calor)
que tiene lugar en el fluido entre los puntos a y b. Ésta difiere de los términos restantes
en la ecuación (4.71) en dos aspectos:
El signo de hf, como se define en la ecuación (4.71), es siempre positivo. Por supuesto,
en el flujo potencial, es igual a cero.
La fricción se produce en las capas límite, debido a que el trabajo realizado por las
fuerzas de corte para mantener los gradientes de velocidad, tanto en el flujo laminar como
en el turbulento, se convierte finalmente en calor por la acción viscosa. La fricción que se
genera en las capas límite no separadas se llama fricción de superficie. Cuando las capas
límite se separan y forman estelas, se produce una disipación adicional de energía en la
estela, y a la fricción de este tipo se le llama fricción de forma, ya que es una función de
la posición y de la forma del sólido.
En una situación determinada, pueden actuar la fricción de superficie y la fricción
de forma a la vez en grado variable. En el caso que se representa en la figura 3.9a, la
fricción es completamente de superficie; en el de la figura 3.9b, la fricción es en gran
parte fricción de forma, debido a que la estela es grande y la fricción de superficie es
relativamente poco importante. La fricción total hf en la ecuación (4.71) incluye ambos
tipos de pérdidas por fricción.
EJEMPLO 4.5 Agua con una densidad de 998 kg/m3 (62.3 lb/ft3) entra en un accesorio de
tubería de 50 mm (1.969 in.), de manera horizontal como se muestra en la figura 4.10, con una
velocidad estacionaria de 1.0 m/s (3.28 ft/s) y una presión manométrica de 100 kN/m2 (14.48
lbf /in.2). El agua sale horizontalmente del accesorio, a la misma elevación, formando un ángulo
de 45° con respecto a la dirección de entrada. El diámetro a la salida es 20 mm (0.787 in.).
Suponiendo que la densidad del fluido es constante, que los factores de corrección de la energía
cinética y del momento, tanto a la entrada como en la salida son la unidad, y que la pérdida
por fricción en el accesorio es despreciable, calcule: a) la presión manométrica a la salida del
accesorio y b) las fuerzas en las direcciones x y y ejercidas por el accesorio sobre el fluido.
Solución
–
a) V a = 1.0 m/s. De la ecuación (4.14),
2 2
D 50
Vb = Va a = 1.0 = 6.25 m/s
Db 20
pa = 100 kN/m 2
V̄b
V̄b,y
pb
V̄b,x
V̄a θ
pa
FIGURA 4.10
Flujo a través de un accesorio de reducción, visto
desde arriba, ejemplo 4.5.
pa − pb V b2 − V a2
=
ρ 2
A partir de la cual
pb = pa −
(
ρ V b2 − V a2 ) = 100 − 998(6.252 − 1.02 )
2 1 000 × 2
= 100 − 18.99 = 81.01 kN/m 11.75 lb f /in.2
2
( )
b) Las fuerzas que actúan sobre el fluido se obtienen combinando las ecuaciones (4.51) y
(4.52). Puesto que Fg = 0 para el flujo horizontal, en la dirección x se obtiene
( )
m� βbVb, x − β aVa, x = pa Sa, x − pb Sb, x + Fw, x (4.72)
donde Sa,x y Sb,x son las áreas proyectadas de Sa y Sb sobre planos normales a la dirección
inicial de flujo. (Recuerde que la presión p es una magnitud escalar.) Ya que el flujo entra en
– –
la dirección x, Va, x = Va y
π
Sa, x = Sa = 0.050 2 = 0.001964 m 2
4
De la figura 4.10
Por lo tanto,
( )
Fw, y = m� βbVb, y − β aVa, y − pa Sa, y + pb Sb, y
= 1.96( 4.42 − 0) − 0 + 81.01 × 0.000222 × 1 000
(
= 8.66 + 17.98 = 26.64 N 5.99 lb f )
La energía mecánica distribuida al fluido es, por lo tanto, ηWp, donde η < 1. La ecuación
(4.71) corregida para el trabajo de bomba es
pa α V2 p α V2
+ gZa + a a + ηWp = b + gZb + b b + h f (4.74)
ρ 2 ρ 2
La ecuación (4.74) es una expresión final de la ecuación de Bernoulli para el tratamiento
de problemas sobre el flujo de fluidos no compresibles.
EJEMPLO 4.6 En la instalación que se muestra en la figura 4.11, se bombea una disolu-
ción de gravedad específica igual a 1.84 desde un tanque de almacenamiento a través de
una tubería de acero de 3 in. (75 mm) Norma 40. La eficiencia de la bomba es de 60%. La
velocidad en la línea de succión es de 3 ft/s (0.914 m/s). La bomba descarga a través de
una tubería de 2 in. (50 mm) Norma 40, hasta un tanque elevado. El extremo de la tubería
tubería de 2 in.
50 ft
tubería de 3 in.
a
FIGURA 4.11
Diagrama de flujo para el ejemplo 4.6.
g V2
Wpη = Zb + b + h f
gc 2 gc
De acuerdo con el apéndice 3, las áreas de las secciones transversales de las tuberías de 2 y 3
pulgadas son 0.0513 y 0.0233 ft2, respectivamente. La velocidad en la tubería de 2 in. es
3 × 0.0513
Vb = = 6.61 ft/s
0.0233
Entonces
g 6.612
0.60Wp = 50 + + 10 = 60.68
gc 64.34
y
60.68
Wp = = 101.1 ft ⋅ lb f /lb
0.60
La presión desarrollada por la bomba se encuentra al aplicar la ecuación (4.74) sobre la bomba.
El punto a está en la conexión de succión, y el punto b está en la descarga de la bomba. La
diferencia de altura entre los niveles de succión y descarga es despreciable, de tal manera que
Za = Zb, y la ecuación (4.74) se transforma en
pb − pa V a2 − V b2
= + Wpη
ρ 2 gc
32 − 6.612
pb − pa = 1.84 × 62.37 + 60.68
2 × 32 . 17
= 6 902 lb f /ft 2 o
6 902
144
(
= 47.9 lb f /in.2 330 kN/m 2 )
La potencia utilizada por la bomba es el producto de Wp por la velocidad de flujo de masa,
dividida entre el factor de conversión 1 hp = 550 ft · lbf /s. La velocidad de flujo de masa es
y la potencia es
�
mWp 17.66 × 101.1
P= = = 3.25 hp (2.42 kW )
550 550
SÍMBOLOS
A Área, m2 o ft2
B Distancia entre las láminas, ejemplo 4.3, m o ft
b Ancho de la capa límite, m o ft
C1, C2 Constantes en las ecuaciones (4.33) y (4.34)
D Diámetro del canal circular, m o ft; Da., en el punto a; Db, en el punto b
Ek Energía cinética del fluido, J o ft · lbf; Ek, flujo de energía cinética por unidad de
tiempo, J/s o ft · lbf/s
F Fuerza, N o lbf; Fg, componente de la fuerza de gravedad; Fw, fuerza neta de la
pared del canal sobre el fluido; Fw,x, componente de Fw en la dirección x; Fw,y,
componente en la dirección y; Fθ, fuerza tangencial o componente de fuerza
G Velocidad de masa, kg/m2 · s o lb/ft2 · s
g Vector de la fuerza de gravedad, N o lbf; gx, componente en la dirección x
g Aceleración de la gravedad, m/s2 o ft/s2
gc Factor de proporcionalidad de la ley de Newton, 32.174 ft · lb/lbf · s2
h Pérdida de fricción, J/kg o ft · lbf /lb; hf, pérdida de la fricción en el conducto
entre los puntos a y b; hfp, pérdida de fricción total en la bomba
L Longitud, m o ft .
M Peso molecular; también momento, kg · m/s o ft · lb/s; M,. velocidad de flujo de .
momento por unidad de tiempo, kg · m/s2 o ft · lb/s2; M a, en el punto a; M b,
en el punto b
.
m Masa, kg o lb; m, velocidad de flujo de masa, kg/s o lb/s
n Número de moles
P Potencia, kW o hp
p Presión, N/m2 o lbf /ft2; pa, en el punto a; pb, en el punto b
q Velocidad volumétrica de flujo, m3/s o ft3/s
Re Número de Reynolds, adimensional
r Distancia radial, m o ft; r1, en el punto 1; r2, en el punto 2; también, distancia
de la superficie líquida, m o ft
rH Radio hidráulico, m o ft
S Área de la sección transversal, m2 o ft2; Sa en el punto a; Sb, en el punto b; Sa,x,
Sb,x, proyecciones de Sa y Sb sobre los planos normales al eje x; Sa,y, Sb,y, proyec-
ciones sobre los planos normales al eje y
T Momento de torsión, N · m o ft · lbf; también temperatura, °C o °F; Ta, tempe-
ratura en el punto a; Tb, en el punto b
t Tiempo, s
Letras griegas
PROBLEMAS
4.1. Por una tubería de 75 mm circula un líquido con flujo estacionario. La velocidad local varía
con la distancia desde el eje de la tubería como se muestra en la tabla 4.1. Calcule: a) la velo-
–
cidad media V , b) el factor de corrección de la energía cinética α, y c) el factor de corrección
de momento β.
4.2. Para un fluido y una distancia dados entre las placas en un sistema como el que se muestra en
la figura 4.5, la velocidad υ0 debe ser suficientemente grande para contrarrestar los efectos de la
gravedad; de lo contrario, algo del fluido se desplazaría hacia abajo. En un sistema particular,
la distancia entre las placas es de 1 mm, y el fluido es un aceite con una densidad de 900 kg/m3
y una viscosidad de 50 mPa · s. La caída de presión dp/dy es despreciable comparada con el
término ρg. a) ¿Cuál es la velocidad ascendente mínima de la placa en movimiento, así como
de todos los fluidos de movimiento ascendente? b) Si υ0 se establece en su valor mínimo,
TABLA 4.1
Datos para el problema 4.1
¿cuál es la velocidad del fluido a la mitad de las placas? c) ¿Cuál es la velocidad de corte en
el fluido en una placa estacionaria, en una placa en movimiento y a la mitad de ellas? Utilice
las ecuaciones del ejemplo 4.3.
4.3. a) Un tanque con agua tiene 30 ft de diámetro, y la profundidad normal es de 25 ft. La sali-
da está al fondo de una tubería horizontal de 4 in. Si esta tubería es interrumpida cerca del
tanque, ¿cuál es la velocidad inicial de flujo del agua desde el tanque? (Ignore la pérdida
de fricción en el pedazo de tubería.) b) ¿Cuánto tiempo le tomará al tanque para vaciarse?
c) Calcule la velocidad promedio de flujo y compárela con la velocidad de flujo inicial.
4.4. Se bombea agua a 20 °C con una velocidad constante de 9 m3/h desde un gran depósito
situado en el suelo hasta la parte superior abierta de una torre experimental de absorción. El
punto de descarga está a 5 m por encima del piso, y las pérdidas por fricción en la tubería
de 50 mm desde el depósito hasta la torre son de 2.5 J/kg. ¿A qué altura ha de mantenerse
el nivel de agua en el depósito si la potencia que puede desarrollar la bomba es tan sólo de
0.1 kW?
4.5. En un codo de 90 ° y 100 mm de diámetro interno DI, situado en un plano horizontal, entra
agua con una velocidad de 6 m/s y una presión manométrica de 70 kN/m2. Despreciando la
fricción, ¿cuál será la magnitud y la dirección de la fuerza que ha de aplicarse al codo para
que permanezca inmóvil?
4.6. Un reactor cilíndrico vertical de 2.5 m de diámetro y 4 m de altura es enfriado por pulverización
del agua sobre la parte superior y permitiéndole al agua fluir debajo de la pared de salida. El
agua circula a una velocidad de 0.15 m3/min, y la temperatura promedio del agua es de 40 °C.
Estime el espesor de la capa de agua.
4.7. Una solución acuosa, con una densidad de 1 040 kg/m3 y una viscosidad de 1.25 mPa · s, se
aspersa sobre la parte superior de una columna empacada con anillos Raschig de 2.5 2.5
cm. La velocidad de flujo del líquido es de 4 500 kg/m2 · h. a) Estime el espesor promedio
de la película líquida que fluye hacia abajo sobre las piezas del empaque. Asuma un valor de
45º para el ángulo de inclinación y que el líquido se distribuye uniformemente sobre la capa
superior del empaque. Véase tabla 18.1 para obtener información del empaque; b) ¿qué frac-
ción del volumen de la columna es alcanzada por los líquidos?; c) si sólo se moja la mitad de
la superficie total del empaque, ¿cuál es el espesor promedio de la película líquida?
4.8. Un dispositivo de almacenaje por bombeo extrae agua de un río por la noche, cuando la de-
manda es baja, y la bombea hasta un recipiente en la cima de una colina que se encuentra a 500
ft (pies) por encima del nivel del río. El agua se distribuye de regreso pasando por turbinas
durante el día para cumplir con la demanda pico. a) Para dos tuberías de 30 in, cada una de
2 500 ft de longitud que acarrean 20 000 gal/min, ¿qué potencia de bombeo se necesita, si la
eficiencia de las bombas es de 85%? La pérdida por fricción se estima en 15 ft de agua;
b) ¿cuánto poder potencia generar las turbinas usando la misma velocidad de flujo total?, y c)
¿cuál es la eficiencia total de esta instalación como sistema de almacenaje de energía?
4.9. Se obtiene una aproximación del perfil de velocidades de una corriente utilizando un modelo
de tres secciones, en el que la velocidad es 1.6 m/s en 15% del área transversal, 3.2 m/s en
35% del área y 4.5 m/s en el resto del área. a) ¿En cuánto difiere el momentum total del fluido
del producto de la velocidad de flujo de masa y de la velocidad promedio?, b) ¿Puede usted
pensar una situación donde exista un perfil de velocidades en la corriente de flujo y que el
momentum del flujo sea exactamente igual al producto del flujo de masa por la velocidad
promedio?
REFERENCIAS
1. Bennett, C.O. y J.E. Myers. Momentum, Heat, and Mass Transfer. 3a. ed. Nueva York: McGraw-
Hill, 1982.
2. Bird. R.B., W.E. Stewart y E.N. Lightfoot. Transport Phenomena. Nueva York: Wiley, 1960:
a) p. 79. b) p. 81. c) pp. 84-91.
3. Geankoplis, C.J. Transport Processes and Unit Operations. 3a. ed. Englewood Cliffs, NJ:
Prentice-Hall, 1993.
4. Nusselt, W. VDIZ. 60: 541, 569 (1916).
5. Perry, R.H. y D.W. Green, eds. Perry’s Chemical Engineers’ Handbook. 7a. ed. Nueva York:
McGraw-Hill, 1997, p. 6-7.
6. Portalski, S. Chem. Eng. Sci. 18: 787 (1963).
7. Streeter, V.L. y E.B. Wylie. Fluid Mechanics. 8a. ed. Nueva York: McGraw-Hill, 1985.