Folleto de Mecánica de Fluidos PROFESIONALIZACIÓN
Folleto de Mecánica de Fluidos PROFESIONALIZACIÓN
Folleto de Mecánica de Fluidos PROFESIONALIZACIÓN
INGENIERÍA ELECTROMECÁNICA
4
Donde:
liq, liq. – Densidad y peso específico a la temperatura dada respectivamente.
agua 40°C, agua 40°C – Densidad y peso específico del agua a 40°C respectivamente.
Compresibilidad
Propiedad del líquido de cambiar su volumen bajo la acción de la presión. La compresibilidad
también se caracteriza por el coeficiente de compresión volumétrica p, que no es más que el
inverso del módulo de elasticidad volumétrico, o sea:
1 ∆v (I.6)
β=− ∗
V ∆p
La magnitud inversa del módulo de compresión volumétrica β es el módulo volumétrico de
elasticidad K
El signo menos en la fórmula significa que al incremento positivo de la presión p le corresponde el
negativo es decir la reducción del volumen (V). Al despejar el volumen final de la ecuación I.6 se
tiene:
V = Vo (1 − βp ∆P) (I.7)
Dilatación Térmica
Se caracteriza por el coeficiente de expansión volumétrica βt que representa el cambio relativo del
volumen al variar la temperatura T
1 ∆v (I.8)
βt = ∗
Vo ∆t
Viscosidad
La viscosidad () de un fluido es una medida de su resistencia a la deformación cortante o angular.
Las fuerzas de fricción en un fluido en movimiento son el resultado de la cohesión y del intercambio
de la cantidad de movimiento entre moléculas. La viscosidad de los fluidos varía con la temperatura.
Al aumentar la temperatura la viscosidad de todo líquido disminuye, mientras que la viscosidad de
todos los gases aumenta. Esto se debe a que la fuerza de cohesión entre moléculas disminuye al
aumentar la temperatura. Esta fuerza es predominante en los líquidos por lo cual la viscosidad de
estos disminuye. Mientras que en lo gases es predominante el intercambio molecular entre capas
adyacentes de fluidos a diferentes velocidades. El incremento de los choques moleculares provoca
un aumento de la viscosidad.
5
Examínese el caso clásico de dos placas paralelas (Figura I. 1) suficientemente grande para que las
condiciones de contorno sean despreciables. Ellas están separadas una distancia Y, estando el
espacio entre ellas lleno de fluido. Se supone que la superficie inferior es estacionaria, mientras que
la superficie superior se mueve en dirección paralela, a una velocidad U provocada por la aplicación
de la fuerza F que se corresponde con el área de la placa móvil.
En los contornos las partículas de fluido se adhieren a las paredes, por lo que su velocidad es cero
con respecto a la pared. Esta condición, llamada condición de no deslizamiento, ocurre en todo
fluido viscoso. Como resultado en la figura I. 1 la velocidad del fluido en contacto con la placa inferior
tiene que ser cero mientras que la velocidad del fluido en contacto con la superficie superior tiene
que ser U. La forma de la variación de la velocidad con la distancia entre las dos superficies se
denomina perfil de velocidades. Si la separación entre las placas Y la velocidad relativa entre ellas U
no son demasiados grandes, y si el flujo neto de fluido a través del espacio entre las placas es nulo,
se produce un perfil de velocidades lineal.
Figura I.1. Perfil de velocidades del fluido al ser arrastrado por la placa plana superior.
Se ha demostrado en experimentos que una gran cantidad de fluidos en las condiciones anteriores
cumplen la siguiente relación.
F
AU
(I.9)
Y
Por semejanza de triángulos en la figura I.1 U/Y se puede reemplazar por el gradiente de
velocidades du/dy. Introduciendo una constante de proporcionalidad (), el esfuerzo cortante ()
entre dos capas finas de fluido cualesquiera se puede expresar como:
F U du
(I.10)
A Y dy
6
La ecuación 5 se denomina Ecuación de Viscosidad de Newton.
El coeficiente de proporcionalidad () se denomina Coeficiente de Viscosidad, viscosidad absoluta,
viscosidad dinámica o simplemente viscosidad del fluido.
(I.11)
du
dy
Las dimensiones de la viscosidad absoluta son:
F
2 F
L 2 t (I.12)
du L L
dy t
L
En el Sistema Internacional las unidades son: N*s/m2 (Pa/s) y lb*s/pies2 en el Sistema Inglés.
En muchos problemas relacionados con la viscosidad, esta aparece dividida por la densidad. Esta
relación se conoce como viscosidad cinemática () así denominada porque la fuerza no está
involucrada en las dimensiones, quedando únicamente la longitud y el tiempo. De esta forma
(I.13)
Generalmente se expresa en m2/s (SI) o pies2/s en el Sistema Ingles. Otra unidad muy usada en el
Sistema Métrico es el cm2/s, también denominada Stoke (St).
1.2 Problemas resueltos
Problema 1
Un barril de aceite pesa 1,5 kN, calcule el peso específico(𝛾), la densidad (𝜌) y la densidad relativa
(𝜌r )de este aceite. La capacidad del barril es de 159 litros o 0,159 m3. El peso propio del barril es de
110 N.
Datos
Peso del barril (W) = 1,5 kN
Volumen (V) = 159 l o 0,159 m3
Peso propio del barril (W) = 110N
Calcule
𝛾, 𝜌, 𝜌r
7
Solución.
Primero se determina el peso neto del aceite
Peso neto del aceite = Pesototal − PesoPropio del Barril
Peso neto del aceite = 1,5 x103 − 110
Peso neto del aceite = 1390 N
Haciendo uso de la definición de peso específico se obtiene:
Peso neto del aceite 1390
γaceite = =
V 0,159
N
𝛾aceite = 8742,1
m3
La densidad del aceite queda determinada por:
Masa (m)
ρaceite =
V
Peso neto 1390
g 9,81
ρaceite = =
V 0,159
kg
ρaceite = 891,14
m3
La densidad relativa se obtiene
ρaceite 891
ρr = =
ρagua 4C 1000
ρr = 0,891
Problema 2
En un recipiente de 30 litros de capacidad se encuentra aire bajo una presión de 150 at y a +50 0C.
Determinar el peso del aire en la botella teniendo en cuenta que el peso específico del aire a la
kgf
temperatura de +15 0C y a la presión atmosférica es γ = 1,225 m3
Datos
V= 30 l
P= 150 at
T= 50 0C
kgf
γ15 0C = 1,225
m3
T0= 15 0C
8
Solución
Considerando la Constante Universal de los Gases Ideales constante:
P0 V0 P1 V1
= =R
T0 T1
P0 P1
=
γ0 T0 γ1 T1
El peso específico despejando la relación anterior queda:
P1 T0
γ1 = γ0
P0 T1
150 15 + 273
γ1 = 1,225 ∗
1 50 + 273
kgf
γ1 = 163,8 3
m
El peso del aire en la botella
W = γ1 ∗ V1 = 163,8 ∗ 0,03
W = 4,91 kgf
Problema 3
Determinar cómo cambiaran las tensiones tangenciales en el líquido bombeado por una tubería, si el
kgf
keroseno se sustituye por la gasolina. El peso específico del keroseno es γkeroseno = 800 y el
m3
kgf
de la gasolina es γgasolina = 700 . El coeficiente cinemático de viscosidad del keroseno se
m3
9
𝜐 = μρ
μk = ρk ∗ υk = 800 ∗ 2,2 x10−6
μk = 0,00176 Pa ∗ s
μg = ρg ∗ υg = 700 ∗ 0,84 x10−6
μk = 0,00058 Pa ∗ s
La tensión tangencial por su parte queda expresada como:
dV
τ𝑘 = μ k ( )
dy k
La relación entre los valores de viscosidades dinámicas del keroseno y de la gasolina nos dará la
disminución de las tensiones tangenciales.
μk 0,00176
= = 3,03
μg 0,00058
Las tensiones tangenciales disminuirán 3 veces.
Problema 4
En un sistema hidráulico fue echado 250 l de líquido; 25 l quedaron en un depósito abierto.
Determinar los volúmenes mínimos y máximo del líquido en el depósito, si la temperatura puede
variar en los límites de +60 0C a -60 0C, y se rellenó a +20 0C. El coeficiente de expansión
1
volumétrica del líquido es β = 0,00072 grad
10
Solución
Partiendo de conocer el coeficiente de expansión volumétrica βt (ecuación I.8), se puede obtener
entonces la variación de volumen
1 ∆v
βt = ∗
Vo ∆t
∆v = βt ∗ Vo ∗ ∆t
La primera de las variaciones de volumen se determina para un ∆t = 40℃
1
∆v = 0,00072 ∗ 250 l ∗ 40 ℃ = 7,20 l
grad
∆v = βt ∗ Vo ∗ ∆t
para un ∆t = 80℃
1
∆v = 0,00072 ∗ 250 l ∗ 80 ℃ = 14,40 l
grad
Vmax = 25 l + 7,20 = 32,2 l
Vmin = 25 l − 14,40 = 10,6 l
Problema 5
Un recipiente de 10 l (Litros) de cabida está lleno de líquido. Calcular la cantidad de líquido que
escurrirá del recipiente si bajamos la presión de este desde 200 at hasta la atmosférica. El
2
coeficiente de compresión volumétrica es βp = 0,77. 10−4 cm ⁄kgf
Datos
Vo=10 L
P=200 at hasta la presión atmosférica
2
βp = 0,77. 10−4 cm ⁄kgf
Solución
A partir de la ecuación 6, el coeficiente de compresión volumétrica queda expresado como:
∆V
βp =
Vo ∆P
Despejando la variación de volumen para determinar cuánto de este fue escurrido se tiene:
∆V = βp Vo ∆P
2 kgf⁄
= 10 L ∗ 0,77. 10−4 cm ⁄kgf ∗ 200 cm2
11
= 0,154 L = 154 cm3
Problema 6
En el ensayo hidráulico una tubería de diámetro d=500 mm y longitud L=2 km se llenó de agua bajo
la presión de 75 at. Al final del ensayo la presión bajó hasta 70 at. Despreciando la deformación de
la tubería, determinar el volumen de agua derramada durante el ensayo a través de las uniones
estancadas de los tubos.
Datos
Diámetro de la Tubería d=500 mm
Longitud (L) L=2 km
P1= 75 at.
P2= 70 at.
1
K = β = Módulo volumétrico de elasticidad
p
kgf
K agua = 20 000
cm2
Partiendo del inverso del módulo de compresión volumétrica tenemos:
1
∆V =K
Vo ∆P
kgf 1
20 000 = ∆V
cm2
1,962 x109
medios serán mayores los esfuerzos tangenciales y en cuántas veces (siendo iguales las
demás condiciones).
Solución
Las tensiones tangenciales en el agua serán mayores 64,5 veces.
3. En la figura I.3 está representado el esquema del aparato para la calibración de los
manómetros. La presión del aceite en la cámara que se transmite a ambos manómetros se crea
atornillando el émbolo buzo de diámetro d=1 cm. Determinar cuántas revoluciones habrá que
tener el émbolo buzo para formar una presión de 250 at, si el paso del tornillo es t = 2 mm y el
volumen de la cámara que puede ser considerado invariable al cambiar la presión es igual a
cm2
300 cm3. El coeficiente de compresión volumétrica del aceite es βp = 0,47 x10−4 kgf
13
La potencia está ligada con la tensión tangencial 𝜏 en la superficie del vástago por la fórmulaN =
2𝛿𝑁
τLπdw, de donde μ = πL𝑤2 𝑑2
5. En una tubería llena de agua, por medio del émbolo se produjo un cambio de presión 𝑑𝑝 =
1000 kgf/m2 . Calcular el cambio de densidad provocado por dicho efecto.
Solución
La oscilación de la densidad del agua es 5,01 kgf*s2/m4, y la del aire, a 0,0086 kgf*s2/m4
6. La presión del aceite en los espacios libres entre las partes móviles y fijas de la máquina
hidráulica aumentó desde 10 at hasta 700 at, mientras que la temperatura del mismo y la
velocidad de la corriente quedaron invariables. ¿En cuántas veces cambiarán las tensiones
tangenciales?
Solución
Caso 1: 0,1785; para el caso 2: 1,455
14
Capítulo II. Estática de los Fluidos
2.1. Estática de los fluidos
La estática de los fluidos consiste en el estudio de los problemas en los cuales no existe movimiento
relativo entre los elementos de estos fluidos. Al no existir este, entre los elementos individuales, no
puede existir esfuerzo de corte, sea cual sea la viscosidad, por lo que no tiene efecto sobre los
problemas de estática.
La presión en un punto es la relación entre la fuerza normal y el área, por lo que cumple con una
serie de supuestos o propiedades como:
1. La Presión en un punto en el seno de un fluido en reposo es igual en todas las direcciones.
2. La Presión en todos los puntos situados en un mismo plano horizontal en el seno de un
fluido en reposo tiene el mismo valor.
3. Las leyes de la estática de los fluidos reales, no se diferencia de las leyes de la estática
obtenida para los fluidos ideales.
4. La Fuerza de presión en un fluido en reposo se dirige siempre hacia el interior del fluido.
5. La superficie libre de un líquido en reposo es siempre horizontal.
La presión tiene unidades de fuerza por unidad de área, las cuales pueden ser newton por metro
cuadrado, llamadas pascales (Pa), libras por pie cuadrado (psf) o libras por pulgada cuadrada (psi)
F (II.1)
P=
A
En el SI
N Kg. m Kg
P= 2
= 2 2= = Pa
m s .m m. s2
La presión puede expresarse con respecto a cualquier nivel de referencia arbitraria, los más usuales
son el cero absoluto y la presión atmosférica local. Cuando la presión se expresa como una
diferencia entre su valor y un vacío completo, se conoce como presión absoluta, Cuando se expresa
entonces como la diferencia entre su valor y la presión atmosférica local, se conoce como presión
manométrica.
15
Figura II.1. Unidades y escalas para la medida de la presión.
Según Vennard & Street (1992), la ecuación de la hidrostática queda:
P2 + ρgZ2 = (P0 + ρgZ0 = C (II.2)
Fuerza de presión del líquido sobre una pared plana
La fuerza total de presión del líquido sobre una pared plana es igual al producto de la superficie de la
pared por el valor de la presión hidrostática en el centro de gravedad de esta superficie.
Fuerza sobre superficies planas sumergidas
16
𝑑𝐹 = 𝛾𝑙𝑑𝐴 𝑠𝑒𝑛 𝛼
Y la fuerza total sobre el área A se encuentra por la integración de esta expresión sobre el área, lo
que queda
A
𝐹 = γ senα ∫ dA (II.3)
A
Aquí∫ dA es el momento del área A, con respecto a la línea O-O, el que también se da por el
producto del área A y la distancia perpendicular lc, desde O-O hasta el centroide (C) del área. Así
A
∫ l dA = lc A
Sustituyendo entonces en la ecuación anterior
F = γAlc senα
Sin embargo
hc = lc sen α
Lo que se reduce la ecuación anterior a
F = γhc A (2.4)
Indica que la magnitud de la fuerza resultante sobre (a un lado de) cualquier área plana sumergida,
se puede calcular multiplicando el área, A, por la presión en su centroide,𝛾ℎ𝑐 .
Fuerza sobre superficies curvas sumergidas
Las fuerzas sobre superficies curvas sumergidas se pueden determinar con facilidad por medio del
cálculo de sus componentes horizontales y verticales. En figura II. 3 se puede determinar la
magnitud de la fuerza resultante F con componentes FH y FV, según se muestra.
17
El análisis de cuerpo libre de fluido ABC, permite el cálculo de FH´ y FV´, las componentes de la
fuerza resultante ejercida por la superficie AB sobre el fluido, y de las respectivamente iguales y
opuestas FH y FV.
Del equilibrio estático del cuerpo libre.
∑ Fx = FHC − FH´ = 0
FH´ = FBC
FV´ = WABC + FAC
2.2 Problemas resueltos
Problema 1
Determinar la presión excesiva (Figura II. 4) del aire en un depósito de carga según la lectura del
manómetro de mercurio compuesto de dos tubos en U. Las marcas de los niveles se dan en metros
a partir del suelo.
18
Problema 2
Calcular la altura h (Figura II. 5) de la elevación del mercurio en un tubo de cristal, si la presión
absoluta del aire en el depósito es P1=0,15 at. La presión atmosférica es igual a 775 mm c.m. La
altura H=1,0 m.
19
Datos
T= 2 Kg
H=20 m
D=300 mm
kgf
γaire = 1,29
m3
kgf
γhid = 0,1 3
m
P + γaire (H) − γhidrogeno (H) = P0
Kgf Kgf
P − P0 = γaire (H) − γhidrogeno (H) = 1,29 ∗ 20 m − 0,1 ∗ 20 m
m3 m3
Kgf Kgf Kgf
P − P0 = 25,8 2
− 2 2 = 23,8 2
m m m
Problema 4
En la figura II. 7 se puede encontrar la presión en A, en Pascales ¿Cuál es la presión del aire en el
tubo?
20
ρliq1
ρrel1 =
ρagua 40 C
ρliq1 = ρrel ∗ ρagua 40 C de donde
Kg Kg
ρliq1 = 0,90 ∗ 1000 = 900
m3 m3
ρliq2
ρrel2 =
ρagua 40 C
ρliq2 = ρrel ∗ ρagua 40 C de donde
Kg Kg
ρliq2 = 2,94 ∗ 1000 3
= 2940 3
m m
PA + ρliq1 gh1 − ρliq2 gh2 − ρliq3 gh3 = Patm
Considerando Patm = 0 entonces
PA = −ρliq1 gh1 + ρliq2 gh2 + ρliq3 gh3
Kg m Kg m Kg m
PA = (−900 3
∗ 9,81 2 ∗ 0,6 m) + (2940 3 ∗ 9,81 2 ∗ 0,5 m) + (1000 3 ∗ 9,81 2
m s m s m s
∗ 0,3 m)
PA = 12066,3 Pa
Determinación de la presión de aire dentro del tubo.
PA + ρliq1 gh1 − ρliq2 gh2 = Paire
Paire = PA + ρliq1 gh1 − ρliq2 gh2
Paire = 12066,3 Pa + 5297,4 Pa − 14420,7 Pa
Paire =2943 Pa
Problema 5
Hallar la fuerza que soporta la presa (Figura 2. 8), si la profundidad del agua delante de la presa es
H=4m, el ángulo de inclinación del borde de presión de la presa es α=600. Llevar a cabo el cálculo
para un metro lineal de longitud ¿Cómo depende la fuerza incógnita del ángulo φ?
21
Figura 2.8. Esquema de análisis.
Datos
Profundidad H=4m
Ángulo de inclinación: φ=600
De la expresión para el cálculo de la fuerza sobre superficies planas tenemos:
Fp = (ρghc )A
Determinación del área para un metro lineal de longitud.
1∗H
A=
2
1 ∗ 4m
A= = 2m
2
Kg m
Fp = (1000 3 ∗ 9,81 2 ∗ 4m) ∗ 2m
m s
Fp = 78480 N = 8000 kgf
Problema 6
La figura 2.9 representa un aliviadero automático de presa AOB. El ángulo AOB es rígido; OA=150
cm; OB= 180 cm. La hoja OA tiene una masa de 3000 kg y la hoja OB tiene una masa de 3600 kg.
La dimensión normal al dibujo es 4 m. Desprecie el rozamiento en O y B. W es un contrapeso cuyo
centro de gravedad se encuentra a una distancia de 165 cm de O. El aliviadero está en equilibrio
cuando el nivel de agua se encuentra como en la figura.
22
Figura 2.9. Esquema de análisis.
Determinar:
a) Fuerza debido a la presión de agua sobre OA y OB
b) Centro de presión sobre OA
c) Valor del contrapeso W
Según el enunciado del problema se tiene
m𝑂𝐴 = 3000 kg
m𝑂𝐵 = 3600 kg
𝑂𝐴 = 1,50 m
𝑂𝐵 = 1,80 m
𝑏 = 4,0 m
𝑂𝐺𝑊 = 1,65 m
a) Fuerza sobre el elemento OA
OA
FpOA = ρghp A = ρg cos 300 ∗ AOA
2
FpOA = 1000 ∗ 9,81 ∗ 0,75 ∗ 0,866 ∗ 6
FpOA = 38 230 N
Fuerza sobre el elemento OB
OB
FpOB = ρghp A = ρg ( ∗ cos 600 + OA ∗ cos 300 ) ∗ 7,2
2
𝐹𝑝𝑂𝐵 = 1000 ∗ 9,81 ∗ (0,9 ∗ 0,5 + 1,5 ∗ 0,866) ∗ 7,2
𝐹𝑝𝑂𝐵 = 123 535 𝑁
23
b) Llamando 𝐴𝐶𝑂𝐴 a la distancia del centro de presiones sobre OA, medida a partir de A, se
tendrá que:
∫ y 2 dA b OA3 /3 2 2
ACOA = = 2
= OA = 1,5
∫ y dA b OA /2 3 3
ACOA = 1 m
Por tanto, la distancia OCOA pedida será:
OCOA = OA − ACOA = 1,5 − 1
OCOA = 0,5 m
Como en la pregunta b, llamando DCOB a la distancia del centro de presiones medida a partir de D:
∫ y 2 dA
DCOB =
∫ y dA
Aplicando el teorema de Steiner para el cálculo de estas integrales se tiene:
4
∫ y 2 dA = (4,3983 − 2,5983 )
3
4
∫ y dA = (4,3982 − 2,5982 )
3
2(4,3983 − 2,5983 )
DCOB = = 3,575 m
3(4,3982 − 2,5982 )
OCOB = DCOB − OD = 3,575 − 2,598
OCOB = 0,977 m
c) Cálculo del valor del contrapeso W
El equilibrio de la compuerta exige que la suma de los momentos de todas las fuerzas que sobre ella
actúan con relación al punto O sea igual a 0, es decir, tomando como positivos los momentos en el
sentido de las agujas del reloj.
m𝑂𝐴 ∗ 0,5 𝑂𝐴 ∗ cos 600 + F𝑃𝑂𝐴 ∗ 𝑂𝐶𝑂𝐴 + 𝑊𝑂𝐴 ∗ 0,5 ∗ 𝑂𝐵
∗ cos 300 − F𝑃𝑂𝐵 ∗ 𝑂𝐶𝑂𝐵 + 𝑊 ∗ 1,65 ∗ cos 300 = 0
𝑊
123 535 ∗ 0,977 − 3000 ∗ 9,81 ∗ 0,5 ∗ 1,5 ∗ 0,5 − 38 230 ∗ 0,5 − 3600 ∗ 9,81 ∗ 0,5 ∗ 1,80 ∗ 0,866
=
1,65 ∗ 0,866
𝑊 = 44 101,8 N
Problema 7
24
El tanque de fuel horizontal de la figura 2. 10, de sección transversal circular, d=2,6 m y longitud
kg
l=9,6 m, se halla totalmente lleno de fuel de densidad ρ = 900 . La presión en el exterior del
m3
tanque es la atmosférica. Calcular la fuerza total (módulo y dirección) que ejerce el fluido en la mitad
ABC del tanque.
Fpz
𝛼 = arc tg = 38,150
Fpx
2.3. Problemas propuestos
1. Calcular la presión Px en la figura 2. 11 si l = 760 mm, h = 500 mm, el líquido 1 es agua (1 =
1000 kg/m3) y el líquido 2 es mercurio (2 = 13 600 kg/m3).
25
Figura 2.11. Esquema de análisis.
Solución
𝐏𝑿 = 59 252,4 Pa
2. Si el barómetro de la figura 2. 12 está lleno de aceite de silicio de densidad relativa r=0,86,
calcular la altura h si la presión barométrica absoluta es de 101,3 kPa.
Solución:
p2=71,4 at.
26
Figura 2.13. Multiplicador hidráulico.
4. Un recipiente de gas (Figura 2.14), fabricado con chapa de acero (𝛿 = 7,85), de 9 mm de
espesor, tiene la forma de un cilindro invertido de 15 m de diámetro y 9 m de alto. Los refuerzos o
roblones, añaden un 2% al peso del recipiente h`=75 mm. La densidad del gas 0,58 kg/m3, y del
aire, que se supondrá constante, 1,28 kg/m3. En la conducción del gas a h=120 m se instala un
manómetro con agua que marca Δh.
Calcular:
a) Δh
b) El contrapeso C que mantiene en equilibrio el sistema.
27
Figura 2. 15. Depósito para combustible líquido.
Solución:
La fuerza de presión del líquido sobre las paredes frontales es igual a 𝜋𝛾𝑎𝑏 2 , la fuerza de presión
𝜋
que dilata la pared lateral es igual a (1 − 4 ) 𝛾𝑎𝑏𝑙
28
Capítulo III. Hidrodinámica y aplicación de la Ecuación de Bernoulli
Para un fluido con todas las partículas moviéndose a la misma velocidad, su energía cinética será
también:
𝟏
EC = 𝟐 m0 v 2 (3.2)
Donde:
m0 =∑ m.- Sumatoria de las masas del conjunto de partículas en movimiento.
Expresada esta Energía Cinética de forma específica (Energía / Unidad de peso) se tiene:
1 1 γ
EC mv2 ( )Vυ2 υ2
= 2 γV = 2 g
= 2g (3.3)
Peso γV
29
𝐿 2
EC υ2 ( 𝑡 )
= = 𝐿 =𝐿
Peso 2g
𝑡2
Donde:
= 1 .- Para una distribución uniforme de velocidades
= 1,02 – 1,15 .- Para flujo turbulento
= 2 .- Para Flujo laminar
Es insignificante el error cuando consideramos = 1 en la mayoría de los cálculos de ingeniería
para flujo turbulento.
Energía Potencial
La energía potencial de una partícula de fluido depende de su altura (z) por encima de un plano
arbitrario de referencia.
EP = 𝑊 ∗ Δ𝑍 (3.5)
Donde:
W .- Peso de la partícula de fluido
Z .- Distancia por encima del plano de referencia
La energía potencial por unidad de peso es:
EP 𝑊∗Δ𝑍
= = Δ𝑍 (3.6)
𝑃𝑒𝑠𝑜 𝑊
Dimensionalmente
EP 𝐹∗𝑙
= =l
𝑃𝑒𝑠𝑜 𝐹
En el Sistema Internacional la unidad es el metro (m).
Energía Interna
La energía interna es energía almacenada y está asociada con el estado molecular o interno de la
sustancia.
30
Se puede almacenar en muchas formas:
Térmica.
Nuclear.
Química.
Electrostática.
En el curso solo se estudiará la térmica, ésta depende de:
El movimiento molecular.
Las fuerzas de atracción intermolecular.
“La Energía Interna (Ei) es una función de la temperatura”.
Para un proceso a volumen constante la variación de la energía interna por unidad de masa es:
ΔEi = Cv ∗ ΔT (3.7)
Donde:
Cv.- Calor específico de la sustancia a volumen constante. Sus unidades en el Sistema Internacional
N m
son:
Kg K
La ecuación 3.8 expresa la Energía Interna por unidad de masa, dividiendo esta expresión entre la
aceleración gravitatoria se obtiene la Energía Interna por unidad de peso.
N∗m N∗m
ΔEi = kg∗g = =𝑚 (3.9)
N
La Ecuación general de la energía para el flujo estacionario de cualquier fluido (3.10) y la Ecuación
de Continuidad son ecuaciones claves en la solución de Problemas de Flujo de Fluido.
Para el flujo compresible es necesario utilizar una tercera ecuación que es la Ecuación de Estado de
los Gases (Ideales o reales).
La Ecuación General de la Energía para Flujo Estacionarios Incompresibles puede simplificarse en
determinadas condiciones y se puede plantear:
𝑝1 𝑣2 𝑝2 𝑣2
+ 𝑍1 + 𝛼 2𝑔1 + ℎ𝑚 + 𝑄𝑐 = + 𝑍2 + 𝛼 2𝑔2 + (𝐸𝐼2 − 𝐸𝐼1 ) (3.10)
𝛾 𝛾
Donde:
En un sistema donde el flujo es incompresible 1=2=
Para flujos turbulentos 1, pero muy cercano a 1 y se puede considerar =1
31
Si no hay transferencia de calor Qc=0
La ecuación 3.10 queda entonces expresada como:
𝑝1 𝑣2 𝑝2 𝑣2
+ 2𝑔1 + 𝑍1 + ℎ𝑚 = + 2𝑔2 + 𝑍2 + (𝐸𝐼2 − 𝐸𝐼1 ) (3.11)
𝛾 𝛾
𝑝1 𝑣2 𝑝2 𝑣2
+ 2𝑔1 + 𝑍1 + ℎ𝑚 = + 2𝑔2 + 𝑍2 + ℎ𝑝 (3.11a)
𝛾 𝛾
La ecuación 3.12 es la Ecuación de Bernoulli en honora Daniel Bernoulli (1700 – 1782). El teorema
data de 1738 y es una forma de manifestar la Ecuación de Conservación de la Energía.
3.2. Problemas resueltos
Problema 1
Para medir el gasto de gasolina que fluye por un tubo de Diámetro d=14 mm se halla instalada una
tobera (Figura 3.1) de diámetro d=9 mm y van acoplados los piezómetros (ver figura).
Determinar el gasto de gasolina Q en L/s, si la diferencia de niveles de gasolina en los piezómetros
es h=1,5 m. ¿Cómo cambiara la diferencia de alturas H, si la gasolina se sustituye por agua con el
mismo gasto Q?
Kg
Nota: La densidad para la gasolina es ρ = 680 m3
Datos
Diámetro del tubo: D= 14 mm
Tobera d= 9 mm
h= 1,5 m
Aplicando la Ecuación de Bernoulli 3.12
P1 V12 P2 V22
+ + Z1 = + + Z2
ρg 2g ρg 2g
Debido a que ambos niveles del líquido son iguales dentro del
tubo Z1=Z2=0
P1 V2 P V2
ρg
+ 2g1 = ρg2 + 2g2 Figura 3.1 Tobera.
32
De la ecuación de continuidad Q = A1 V1 = A2 V2
πD2
A1 4
V2 = V = V1
A2 1 πd2
4
3,14∗0,014m2
4
V2 = 3,14∗0,009m2
V1
4
V2 = 2,42V1
Entonces
P1 − P2 V22 − V12
=
ρg 2g
P1 − P2 (2,42V1 )2 − V12
=
ρg 2g
P1 − P2 4,856V12
=
ρg 2g
kg m
Determinación de P1 − P2 = ρgh = 680 m3 ∗ 9,81 ∗ 1,5 m = 10006,2 Pa
s2
2 ∗ 10006,2 Pa
4,856V12 = kg
680
m3
V1 = 2,46 m/s
Determinación del gasto Q en L/s
3,14 ∗ 0,014m2 m −4
m3 l
Q = A1 V1 = ∗ 2,46 = 3,786 x10 = 0,378
4 s s 𝑠
Problema 2
Al medidor de Venturi cuyas dimensiones son: D=50 mm y d=30 mm, esta acoplado un manómetro
diferencial en U de mercurio (Figura 3.2). Determinar el gasto de agua Q, si la indicación del
manómetro es h= 300 mm.
33
Figura 3.2. Medidor de Venturi.
Datos
D=50 mm
D=30 mm
Aplicando la ecuación de Bernoulli
P1 V12 P2 V22
+ + Z1 = + + Z2
ρg 2g ρg 2g
Debido a que ambos niveles del líquido son iguales dentro del tubo Z1=Z2=0
P1 V12 P2 V22
+ = +
ρg 2g ρg 2g
De la ecuación de continuidad Q = A1 V1 = A2 V2
πD2
A1 4
V2 = V = V1
A2 1 πd2
4
3,14∗0,05m2
4
V2 = 3,14∗0,03m2
V1
4
V2 = 2,78V1
Entonces
kg m
P1 − P2 = ρhg gh = 13600 ∗ 9,81 ∗ 0,3 m = 40024,8 Pa
m3 s2
P1 − P2 (2,78V1 )2 − V12
=
ρg 2g
34
40024,8 7,73V12
=
1000 ∗ 9,81 2 ∗ 9,81
Despejando V1 tenemos:
m
V1 = 3,449
s
Determinación del gasto Q
3,14 ∗ 0,05m2 m m3 l
Q = A1 V1 = ∗ 3,45 = 6,77 x10−3 = 6,77
4 s s 𝑠
Problema 3
El ventilador centrifugo (Figura 3.3) aspira aire de la atmósfera a través de una tobera (ver figura). A
la parte cilíndrica de la tobera cuyo diámetro es D= 200 mm va acoplado un tubo de cristal cuyo
extremo inferior está sumergido en un recipiente con agua. El agua en el tubo se elevó hasta la
kgf
altura h=250 mm. Determinar la cantidad de aire que se aspira por segundo 𝛾 = 1,29 m3
35
PA V12 Pagua V22
+ + Z1 = + + Z2
ρg 2g ρg 2g
Debido a que ambos niveles del líquido son iguales dentro del tubo Z1=Z2=0
A su vez debido al aire en calma a la entrada del tubo V1 = 0
kg m
Pagua = 1000 3
∗ 9,81 2 ∗ 0,25 m = 2452,5 Pa
m s
V22 −2452,5 Pa
=
2g 1,29 kg ∗ 9,81 m
m3 s2
m
V2 = 61,7
s
Determinación de la cantidad de aire que se aspira.
πD2 π ∗ 0,22 m3
Q = A2 V2 = ∗ V2 = ∗ 61,7 = 1,94
4 4 s
Problema 4
Hallar la velocidad de la corriente de agua en el tubo si la lectura del manómetro de mercurio unido
al tubo de Pitot y a los orificios de presión estática es h= 600 mm (Figura 3.4)
PE = 24 930 Pa
2. En una tubería de diámetro D=50 mm se ha colocado delante de la válvula un manómetro (Figura
3.6). Estando cerrada la válvula, el manómetro indica una presión igual a 6 at. Cuando la válvula
está abierta la lectura disminuye hasta 2 at. Determinar el gasto de agua en la tubería.
37
Figura 3.6. Esquema de análisis.
Solución
Escribiendo la ecuación de Bernoulli para la superficie libre en el depósito donde la presión es
constante y para la sección en que está puesto el manómetro, obtenemos que Q=55 l/s.
3. Determinar la rarificación de aire (en mm c.m) en la sección estrecha de la boquilla adicional de
aire del carburador del motor del avión (Figura 3.7), si el aire se aspira directamente de la
atmósfera a la altura cero (Q=0,15 m3/s). El diámetro de la sección estrecha de la boquilla
adicional de aire es D= 40 mm, el diámetro del pulverizador de gasolina es d = 10 mm y 𝛾 =
kgf
1,29 m3
38
Capítulo IV. Teoría de Semejanzas
39
Datos
Escala: 1:15
Agua salada
vp: 12 mph
Igualando los números adimensionales y despejando la velocidad del modelo tenemos:
Rep = Rem
Fi Fi
( ) =( )
Fv p Fv m
v∗l v∗l
( ) =( )
ν p ν m
v ∗ l = v ∗ lm
lm
120 mph l = vm ∗
15
vm = 180 mph
Si las viscosidades fueran diferentes
v∗l v∗l
( ) =( )
νagua 15 νagua
p m
vm l νagua = 15 v l νagua m
15 νagua v l 15 v νagua
m
vm = =
l νagua νagua
𝑝 p
Problema 2
A través de una tubería de 60 cm de diámetro está circulando aire a 20 0C a una velocidad media de
2,0 m/s. ¿Cuál debe ser el diámetro de la tubería que a transportar agua a 15 0C y a una velocidad
de 1,22 m/s de lugar a un flujo dinámicamente semejante?
Datos
Øt : 60 cm
Aire 20 0C υ= 1,448 .10-5 m2/s
vaire = 2,0 m/s
Agua 15 0C υ= 1,142 .10-6 m2/s
vagua= 1,22 m/s
Rep = Rem
40
v∗d v∗d
( ) =( )
ν p ν m
m m
2 ∗ 0,6m ∗d 1, 22
s s
( m 2 ) = ( m2
)
1,488 . 10−5 𝑠 1,142 . 10−6 𝑠
p m
d = 0,075 m = 7,5 cm
Problema 3
Un barco cuyo casco tiene una longitud de 140 m ha de moverse a 7,50 m/s.
a) Calcular el número de Froude.
b) Para la semejanza dinámica ¿A qué velocidad debe remolcarse en agua un modelo
construido a una escala 1:30?
Cálculo del número de Froude
v 7, 5 m/s
ℱ= = = 0,2
√l ∗ g √140 m ∗ 9,81
m
s2
v v
( ) =( )
√l ∗ g p √l ∗ g m
7,5 v
( ) =
√140 ∗ 9,81 p 140
√ ∗ 9,81
( 30 )m
v = 1, 37 m/s
4.3. Problemas propuestos
1. Un aceite de viscosidad cinemática 4,70 *10-5 m2/s va a utilizarse en un prototipo en el que
son fuerzas predominantes las debidas a la viscosidad y a la gravedad. También se desea
experimentar sobre un modelo a escala 1:5. ¿Qué valor debe tener la viscosidad del líquido
del modelo para que tanto el número de Froude como el de Reynolds sean iguales en
modelo y prototipo?
2. El modelo de un recipiente se vacía en 4 minutos al abrir una compuerta de tajadera. El
modelo está construido a una escala 1:225. ¿Cuánto tiempo tardara en vaciarse el
prototipo?
41
3. La resistencia medida en agua dulce, presentada a un modelo de barco de 2, 50 m,
moviéndose a una velocidad de 2,0 m/s, fue de 4,40 kg.
¿Cuál será la velocidad del prototipo de 40 m?
¿Cuál será la fuerza necesaria para mover a esta velocidad el barco en agua salada?
42
Capítulo V. Teoría de la capa límite
43
Como ponen de manifiesto las fórmulas dadas, el espesor de la capa límite es directamente
proporcional a la raíz cuadrada de la longitud x y la raíz cuadrada de la viscosidad cinemática e
inversamente proporcional a la raíz cuadrada de la velocidad. Análogamente, la tensión cortante en
la superficie del contorno τo es directamente proporcional a la raíz cuadrada del producto de ρ y µ a
la potencia de tres medios de V e inversamente proporcional a la raíz cuadrada de x.
2. Capa límite turbulenta (contorno liso)
a. Coeficiente de resistencia medio
1,328 (V.4)
CD = 𝑅𝑒 0,20 para 2.105< Re< 107
45
Capítulo VI. Flujo de Fluidos en tubos
46
V22
Hb = Z2 + Z1 + + ∑ hf (Ec. 1)
2g
Donde
2 2
∑ hf = ∑ hf T.asp + ∑ hf T.desc
1 1
2
L V2
∑ hf tramo_recto = λ ∗
d 2g
1
47
λ = 0,021
Dejando la expresión de la perdida en función del caudal tenemos:
2
7m Q2
∑ hf tramo_recto_asp = 0,021 ∗ 2∗9,81∗π2 ∗0,154
0,15 m
1 16
2
∑ hf tramo_recto_asp = 163,3 Q2
1
2
V2
∑ hf accesorios_asp = ξ ∗
2g
1
∑ 𝜉accesorios = 2,69
∑ hf accesorios_asp = 439 Q2
1
48
De la ecuación de continuidad
Q= A∗V
Q 0,0277 m3 /𝑠
V= = 𝜋∗0,1252
= 2,25 m/s
A
4
kg m
ρvd 995,7 m3 ∗ 2,25 s ∗ 0,125 m
Re = =
μ 0,8 ∗ 10−3 Pa. s
R e = 350 050,8 Estamos en presencia de régimen turbulento
Según el material de la tubería la rugosidad superficial es de ε=0,2 mm
Determinación de la rugosidad relativa
ε 0,2 mm
ξ= = = 0,0016
d 125 mm
Por tanto entrando al diagrama de Moody con R e = 350 050,8 y ξ = 0,0016 para un régimen
turbulento
λ = 0,023
Dejando la expresión de la pérdida en función del caudal tenemos:
2
215 m 16 Q2
∑ hf tramo_recto_imp = 0,023 ∗ 2∗9,81∗π2 ∗0,1254
0,125 m
1
2
∑ hf tramo_recto_imp = 13 388,6 Q2
1
∑ 𝜉accesorios = 3,4
49
2
16Q2 16Q2
∑ hf accesorios_imp = =
2 ∗ 9,81 ∗ 𝜋 2 ∗ d4 2 ∗ 9,81 ∗ 𝜋 2 ∗ 0,1254
1
2
∑ hf accesorios_imp = 338Q2 ∑ 𝜉
1
2
∑ hf accesorios_imp = 1150, 6 Q2
1
Entonces de la Ec.1
V22
Hb = Z2 + Z1 + + ∑ hf (Ec. 1)
2g
Hb = 32 m + 5 m + 338 Q2 + 13 388,6 Q2 + 1150,6 Q2 + 163,3 Q2 + 439 Q2
Hb = 37 + 15 472,2 Q2
Construyendo la característica del sistema con ayuda del Excel dado los valores en la Tabla 1.
Tabla 1. Valores para la construcción de la curva característica del sistema de tuberías.
Caudal (m3/h) Carga H(m)
20 37,478
40 38,910
60 41,298
80 44,641
100 48,938
120 54,191
La figura 6.2 muestra la representación de la carga estática y dinámica del sistema de tuberías del
problema anterior.
50
Figura 6.2. Curva característica del sistema de tuberías.
Fuente: Elaboración propia.
6.2. Problemas propuestos
1. En la Figura 6.3 muestra un esquema de un sector de aspiración del conducto de aceite en un
avión. Determinar la presión absoluta del aceite a la entrada de la bomba, siendo: la altura de
vuelo H=8 km (presión atmosférica: 266,8 mm c.m), el gasto del aceite Q=1,25 l/s, el diámetro del
tubo d=40 mm, la longitud l=1,5 m, la temperatura del aceite 60 ℃ , el coeficiente de resistencia
kgf
de la llave ξ=2,0 γac = 860 y la altura h=1,0 m. La presión del aire en el depósito es igual a
m3
la atmosférica.
51
2. La bomba aspira el agua del depósito A por un tubo de diámetro d=20 mm en una cantidad Q=1
l/s. La longitud total del tubo de aspiración es l=10 m. La temperatura del agua es 60 ℃ lo que
corresponde a la presión de formación del vapor p t=0,2 at. Despreciando las resistencias
locales y aceptando el coeficiente λ=0,02 determinar la altura máxima admisible de aspiración
Haspa condición de que la presión absoluta del agua a la entrada en la bomba supere en 1 m c.a
a la presión de formación de vapor.
Solución
Hasp =1 m
52
Capítulo VII: Flujo de fluidos compresibles. Sustentación
VII.1Fuerza de empuje ascensional
La fuerza de empuje ascensional o de sustentación es una fuerza que actúa sobre un cuerpo
sumergido en dirección perpendicular al movimiento relativo entre el fluido y el cuerpo. (Smith, 1980)
(Mataix, 1982)(Munson, 1999)(Potter, 2002)
A b
Figura VII.1. Perfiles de velocidad (a) y presión (b) en un ala moviéndose en un fluido en reposo.
Aplicando el teorema de Bernoulli a las líneas de corriente resulta que hay una presión menor en la
parte superior del ala y una presión mayor en la parte inferior Figura I.1 b) lo que da lugar a una
sustentación neta hacia arriba. (Elaboración propia)
La mayor velocidad en la parte superior del ala Figura VII.1 y la menor velocidad en la parte inferior
se pueden explicar, sin dejar de tener en cuenta que se induce una circulación Γ (Gamma)
VII.2 Circulación
Para poder explicar el concepto de Circulación es necesario conocer la ecuación que rige en el
campo fluido de un torbellino o vórtice libre.
En el vórtice libre no hay consumo de energía procedente de una fuente exterior. El fluido gira
gracias a una rotación previamente dada al mismo o debido a alguna acción interna (Fay, 1996)
(Shames, 1967). Ejemplo de ello son: El torbellino en un río, el flujo giratorio que se produce en
53
recipientes de poca profundidad cuando el líquido se desaloja a través de un agujero en el fondo
(lavamanos, bañaderas etc.) y el flujo en la carcasa de una bomba centrífuga fuera del rodete o
impelente (en la voluta).
d (m. V)s
∑F = 𝐄𝐜. 𝐕𝐈𝐈. 𝟏
dt
La sumatoria de todas las fuerzas exteriores que actúan sobre un sistema fluido es igual a la
velocidad de cambio de la cantidad de movimiento ( (m. V)s )
El momento de las fuerzas de impulso con respecto al centro de giro de la partícula de masa m
moviéndose en una trayectoria circular con una velocidad V t. viene dada por la segunda ley de
Newton la cual queda expresada a partir de la Ec. VII.3:
d (m. Vs). r
𝜏 = F. r = 𝐄𝐜. 𝐕𝐈𝐈. 𝟑
dt
54
d (m. Vu. r)
𝜏= 𝐄𝐜. 𝐕𝐈𝐈. 𝟑. 𝟏
dt
Para el caso del Vórtice Libre (Fluido sin fricción) no hay un Par o Torque exterior aplicado por lo que
la Ecuación VII.3.1 queda:
d (m. Vu. r)
𝜏= =0
dt
Por lo que el término m. Vu. r = Constante y por ello para el vórtice libre se cumple que:
m V1 r1 = m V2 r2 𝐄𝐜. 𝐕𝐈𝐈. 𝟑. 𝟐
La ecuación VII.3.2 a) es la ecuación que rige en el campo fluido de un vórtice o torbellino libre.
La circulación se puede definir haciendo uso de la Figura VII.3 donde las líneas de corriente
representan un campo de flujo bidimensional y L representa cualquier trayectoria cerrada dentro del
campo.
La circulación (Γ) se define matemáticamente como una integral de línea de la velocidad a lo largo
de una trayectoria cerrada, por tanto.
55
Γ = ∫ V. dl = ∫ V cos β. dl 𝐄𝐜. 𝐕𝐈𝐈. 𝟒
Aclarar que esta integral de línea es cerrada donde:
Retomando a la Figura I.1 la mayor velocidad en la parte superior del ala y la menor velocidad en la
parte inferior se pueden explicar al tener en cuenta que se induce una circulación Γ. Cuando el ala
se mueve con relación al campo de flujo Figura VII.4
Existe una estrecha relación entre la sustentación y la circulación, la cual será estudiada en el
empuje ascensional de un cilindro donde se superpone la circulación alrededor del cilindro al campo
de flujo. Se puede calcular fácilmente la circulación si se elige como trayectoria cerrada la línea de
corriente circular L1 concéntrica en el centro del torbellino libre.
La circulación Γ es entonces:
56
Γ= ∫ V cos β. dl 𝐄𝐜. 𝐕𝐈𝐈. 𝟒como V es tangente a L el ángulo es 00 y por tanto cos β = 1 por lo
que:
Γ = ∫ V. dl 𝐄𝐜. 𝐕𝐈𝐈. 𝟒
Dado que V= Constante, queda la integral cerrada de la trayectoria circular, que es la longitud de la
circunferencia quedando:
Γ = V ∫ dl = V1 (2πr1 ) = V2 (2πr2 )
V1 r1 = V2 r2 = C
Se puede demostrar con mayor rigor que la circulación alrededor de cualquier trayectoria que
encierra el centro del torbellino viene dada por la ecuación VII. 𝟓
Como se ve de la ecuación I.5 la circulación depende únicamente del término (C = Vr) que es la
Constante del Torbellino y se denomina Fuerza del Torbellino
57
VII.3 Empuje ascensional de un cilindro circular
Considerando un flujo uniforme de un fluido ideal alrededor de un cilindro que tiene longitud infinita.
(Figura VII.6)
La distribución de presión sobre el cilindro se puede calcular aplicando el teorema de Bernoulli entre
un punto en el infinito de la corriente principal y un punto en la superficie del cilindro. Dado que la
distribución de presión alrededor del cilindro es completamente simétrica no se producen fuerzas de
sustentación, ni resistencia para este caso ideal.
58
Estudiando un flujo estacionario alrededor del cilindro (Figura VII.7) y asignando un signo positivo (+)
a la circulación en el sentido de las agujas del reloj se deduce a partir de la Ecuación VII.5 que la
velocidad periférica Vt2 (en el sentido de las agujas del reloj) sobre la superficie del cilindro debido a
la circulación se puede expresar como:
Γ
Vt2= 𝐄𝐜. 𝐈. 𝟕
2π. R
La velocidad en el sentido de las agujas del reloj en la periferia del cilindro es la suma de las dos
contribuciones ya presentadas en las Figuras VII.6 y VII.7 quedando:
Γ
Vt= 2Usen θ + 𝐄𝐜. 𝐕𝐈𝐈. 𝟖
2πR
La ecuación general para la presión (p) en cualquier punto sobre la superficie del cilindro se obtiene
a partir de la ecuación de la energía (Teorema de Bernoulli).
Pu U 2 P V𝑡2
+ = + 𝐄𝐜. 𝐕𝐈𝐈. 𝟗
γ 2g γ 2g
Donde:
59
Pu Presión en un punto alejado donde la velocidad U es uniforme.
P Pu U 2 V𝑡2
− = − 𝐄𝐜. 𝐕𝐈𝐈. 𝟗. 𝟏
γ γ 2g 2g
𝛾
P − Pu = (U 2 − V𝑡2 )
2g
Donde
ρg 2 ρ
P − Pu = (U − V𝑡2 ) = (U 2 − V𝑡2 )
2g 2
ρ 2
P − Pu = (U − V𝑡2 ) 𝐄𝐜. 𝐕𝐈𝐈. 𝟗. 𝟐
2
ρ 2 Γ 2
P − Pu = [U − (2U sen θ + ] 𝐄𝐜. 𝐕𝐈𝐈. 𝟏𝟎
2 2πR
La Ecuación VII.10 ofrece el valor de la presión que actúa sobre la superficie del cilindro respecto a
la presión en un punto alejado en el campo de flujo.
Considerando la diferencia de presión dada por VII.10 y el área de la superficie del cilindro, se puede
determinar la fuerza de sustentación o empuje ascensional (Fl ), referirse a la Figura VII.8. El área
elemental descrita por el ángulo dθ es
60
dF𝑙 = Pr . dA
Donde
Pr . Presión resultante
Pr = P − Pu
Integrando la Ecuación VII.12 a lo largo del perímetro exterior del cilindro queda:
2𝜋
F𝑙 = B ∫ (P − Pu )R. sen θ. dθ 𝐄𝐜. 𝐕𝐈𝐈. 𝟏𝟑
0
2𝜋
ρ 2 Γ 2
F𝑙 = B ∫ [U − (2Usenθ + ) ] R. sen θ. dθ 𝐄𝐜. 𝐕𝐈𝐈. 𝟏𝟒
0 2 2πR
La fuerza de sustentación total viene dada cuando se resuelve la integral de la Ecuación VII.14 la
que se reduce a:
F𝑙 = ρ. B. U. Γ𝐄𝐜. 𝐈. 𝟏𝟓
Esta fuerza transversal que actúa sobre el cilindro que gira está originada por el Efecto Magnus
(Granados, 2009; Streeter, 1963) denominado así en honor de Heinrich G. Magnus (1802-1870)
científico alemán que lo observó por primera vez en 1852.
61
La Ecuación I.15 se denomina además el Teorema de Kutta-Joukowwski en honor de W.M. Kutta y
N.E. Joukowski pioneros por separado en la investigación cuantitativa de la fuerza de sustentación
en 1902 y 1906 respectivamente.
M L L2 M. L
F𝑙 = 3 . L. . = 2
L T T T
M. L kg. m
F𝑙 = = 2 = Newton
T2 s
El Teorema de Kutta -Joukowski es de gran importancia porque se puede aplicar no solo al cilindro
circular, sino también a un cilindro de cualquier forma, incluyendo un perfil de ala de avión como será
estudiado con posterioridad.
En la Figura VII.8 se ve claramente que los puntos de remanso se mueven hacia abajo respecto a la
horizontal, pero siguen situándose de forma simétrica alrededor del eje vertical.
En el punto de remanso sobre el cilindro la velocidad V t dada por la ecuación I.8 tendrá un valor
cero. Por el punto de remanso se obtiene la siguiente expresión.
Γ
−2U sen θ = 𝐄𝐜. 𝐈. 𝟏𝟔
2πR
62
Si se puede medir el ángulo del punto de remanso y se conoce la velocidad de la corriente principal,
se puede obtener la circulación Γ de la Ecuación VII.16
Donde
-Cuando el perfil de ala es de longitud finita, se producen condiciones que influyen tanto en la
sustentación como en la resistencia.
- Como la presión sobre la superficie interior es mayor que la existente en la superficie superior, el
fluido escapará alrededor de los extremos produciéndose:
Un flujo general hacia afuera desde el centro hasta los extremos, sobre la superficie interior
del ala.
Un flujo hacia dentro desde los extremos hasta el centro sobre la superficie superior.
63
Estos movimientos combinados originaron pequeños torbellinos de borde marginal que se
desprenderán de los puntos de las alas, que constituyen un gran anillo de torbellinos, dentro del cual
se produce una velocidad hacia abajo inducida por los propios torbellinos.
La velocidad inducida cambia la dirección del flujo cerca del ala del sentido U a U 0, por lo que se
reduce el ángulo de ataque efectivo de α a α0 tal como se muestra en la Figura VII.10
Figura VII.10.Perfil de ala de avión con las fuerzas actuantes. (James, 1981)
Donde:
V2
FDi = CDi ρ A 𝐄𝐜. 𝐕𝐈𝐈. 𝟏𝟖
2
Donde
64
ρ Densidad del fluido
La fuerza de resistencia total FD que actúa sobre el perfil de ala se determina según la expresión
(VII.19)
Donde:
FD0 Resistencia del perfil de ala (Incluye la resistencia de presión +la resistencia de fricción del
propio perfil).
La Fuerza de sustentación FL se determina según la expresión dada por la ecuación I.15 pero
resulta además conveniente expresarla como:
V2
FDi = CL ρ A 𝐄𝐜. 𝐕𝐈𝐈. 𝟐𝟎
2
Donde:
A = B. C
Donde:
65
Los Coeficientes de resistencia CD , resistencia inducida CDi y de sustentación CL son obtenidos por
ensayos de los perfiles de alas en los túneles de vientos donde se realizan los ensayos a diferentes
perfiles.
siguientes expresiones:
CL
αi = B2
𝐄𝐜. 𝐕𝐈𝐈. 𝟐𝟏
π A
Y que
CL2
CDi = B2
𝐄𝐜. 𝐕𝐈𝐈. 𝟐𝟐
π A
CD
FD = .F 𝐄𝐜. 𝐕𝐈𝐈. 𝟐𝟑
CL L
Pv = RT 𝐄𝐜. 𝐕𝐈𝐈. 𝟐𝟒
1 P
Entonces si v = ρ obtenemos que: ρ = RT
Donde:
P- Presión absoluta.
N.m
*Para el aire el valor de R = 287 kg.K
P
Por lo que si = RT y ρ = γ/g (siendo 𝛾 peso específico)
ρ
Pg
γ= 𝐄𝐜. 𝐕𝐈𝐈. 𝟐𝟓
RT
De lo que se deduce que se puede calcular el peso específico (𝛾) de cada gas si se conoce R y g.
Dado que las ecuaciones VII − 24 e VII − 25 relacionan varias propiedades de los gases en un
estado dado estas ecuaciones se denominan Ecuaciones de Estado.
Dado que la ley de Avogadro establece que “Volúmenes iguales de todos los gases perfectos a una
presión y temperatura particulares, contienen el mismo número de moléculas “se deduce que el peso
específico de un gas especifico es proporcional a su peso o masa molecular la cual se representa
por M. De la ecuación 𝑉𝐼𝐼 − 24 se puede plantear:
γ M
Dado que: γ2 = M2
1 1
67
Si sustituimos 𝛾 por el término de la ecuación 𝑉𝐼𝐼 − 25 obtenemos:
P2 g
R 2 T2 M R M
P1 g = M2 Por tanto: R1 = M2
1 2 1
R 1 T1
De lo que se despeja:
̅𝐄𝐜. 𝐕𝐈𝐈. 𝟐𝟔
R1 ∗ M1 = R 2 ∗ M2 = R
N. m
̅ = 8312
R
kg. K
Ecuación de Proceso
Donde:
P- Presión absoluta
v- Volumen especifico
-Para un proceso sin transferencia de calor hacia o desde el gas se denomina adiabático (Q=0); pero
Cp
si el además el proceso es reversible, se denomina Isotrópico donde n = k = donde k es la
Cv
68
Otras relaciones útiles que surgen de las combinaciones anteriores son:
n−1
T2 v1 n−1 ρ2 n−1 P2 n
=( ) =( ) =( ) 𝐄𝐜. 𝐕𝐈𝐈. 𝟐𝟖
T1 v2 ρ1 P1
Para los Gases Reales a bajas presiones y temperaturas moderadas se puede considerar que
cumplen las leyes de los Gases Perfectos y son tratados como tales.
P1 v12 P2 v22
( + Z1 + + I1 ) + Wm + Qc = ( + Z2 + + I2 ) 𝐄𝐜. 𝐕𝐈𝐈. 𝟐𝟗
γ 2g γ 2g
I-Energía Interna.
Muchos procesos reales donde es baja la irreversibilidad se aproximan a los procesos reversibles,
ejemplos son las Toberas convergentes –(la fricción es pequeña y resulta despreciable la
transferencia de calor), así como el flujo en tuberías es un proceso irreversible debido a la fuerza de
fricción.
69
En un sistema cerrado la entropía siempre aumenta, por lo que la energía disponible disminuye.
p
h = i + pv o h = i + 𝐄𝐜. 𝐕𝐈𝐈. 𝟑𝟎
ρ
Donde:
𝑖 - Energía interna por unidad de masa. Esta es debida a la energía cinética del movimiento
molecular y a las fuerzas de interacción molecular. Es una función que depende de la temperatura.
Por lo tanto, la Entalpia representa la energía que posee un gas en virtud de su temperatura.
Se define como el aumento de entalpia por unidad de masa cuando la temperatura del gas aumenta
un grado mientras la presión se mantiene constante.
∂h
Por tanto: Cp = (∂T) a P = cte 𝐄𝐜. 𝐕𝐈𝐈. 𝟑𝟏
Se define como el aumento de la energía interna por unidad de masa cuando la temperatura del gas
aumenta un grado mientras el volumen se mantiene constante.
∂i
Por tanto: Cv = (∂T) a v = cte 𝐄𝐜. 𝐕𝐈𝐈. 𝟑𝟐
Ecuación de continuidad
Donde:
ṁ-Flujo másico
70
Q -Flujo volumétrico (Q = A ∗ V)
Ecuación de energía
P1 v12 P2 v22
( + Z1 + α + I1 ) + Wm + Qc = ( + Z2 + α + I2 ) 𝐄𝐜. 𝐕𝐈𝐈. 𝟑𝟒
γ 2g γ 2g
P
Pero el término( γ1 + I) = h (Entalpia)
Dado que las columnas de gases 𝑍1 y 𝑍2 pueden despreciarse por su bajo valor se puede escribir la
𝐸𝑐. 𝑉𝐼𝐼 − 33 como:
v12 v22
(h1 + α ) + Wm + Qc = (h2 + α ) 𝐄𝐜. 𝐕𝐈𝐈. 𝟑𝟓
2g 2g
v2
El coeficiente 𝛼 en el términoα 2g1 es el denominado Coeficiente de Coriolis o Factor de Conexión de
la Energía Cinética. Para flujo turbulento α=1,01 ÷1,15 En la generalidad de los problemas
de ingeniería α =1 para fluidos turbulentos. En fluidos laminares α =2.
v2
En ocasiones el termino h1 + α 2g1 = h0 (Entalpia de remanso)
1
v = f( )
A
Se logran incrementos de la velocidad con la reducción del área de la sección transversal del
conducto.
71
Para el fluido ideal
ṁ = ρ ∗ A ∗ V = ρ ∗ Q = Constante
dρ dA dV
+ + =0
ρ A V
Considerando que:
dp E
-c 2 = dρ = √ρ (velocidad sónica o acústica, velocidad de una perturbación de la presión que viaja
a través de un fluido)
dρ
+ VdV = 0
ρ
dρ
Multiplicando el primer término de la ecuación de Euler por dρ esta no se alterará, la misma se
dρ dρ
∗ + VdV = 0
dρ ρ
dp
Sustituyendo el términodρ =c 2 la ecuación queda:
dρ
c2 ∗ + VdV = 0
ρ
72
dρ
Despejando el término ρ se obtiene
dρ VdV
=− 2
ρ c
VdV dA dV
− + + =0
c2 A V
𝐯
Conociendo que M = donde M es el número de Mach
𝐜
v2
Por tanto c 2 = M2 y sustituyendo en la ecuación queda:
vdv dA dv
− v2
+ + =0
A v
M2
Reordenando la ecuación:
𝑀2 𝑣𝑑𝑣 𝑑𝐴 𝑑𝑣
− + + =0
𝑣2 𝐴 𝑣
M 2 dv dA dv
− + + =0
v A v
M 2 dv dv dA
− + =−
v v A
dv dA
(1 − M 2 ) = −
v A
dv
Despejando v
dA dA dA
dv −A −A A
= = =
v (1 − M 2 ) −(M 2 − 1) (M 2 − 1)
73
dv 1 dA
= 2
∗ 𝐄𝐜. 𝐕𝐈𝐈. 𝟑𝟔
v Ma − 1 A
De la que se concluye:
𝑑𝐴 𝑑𝑣
Si < 0; > 0 (la disminución del área hace que la velocidad aumente)
𝐴 𝑣
𝑑𝐴 𝑑𝑣
Si > 0; < 0 (El aumento del área hace que la velocidad disminuya)
𝐴 𝑣
𝑑𝐴 𝑑𝑣
Si < 0; < 0 (la disminución del área hace que la velocidad disminuya)
𝐴 𝑣
𝑑𝐴 𝑑𝑣
Si > 0; > 0 (El aumento del área hace que la velocidad aumente)
𝐴 𝑣
1
= ∞ → Excesivamente grande (Indeterminado)
𝑀2 −1
𝑑𝑣
𝑑𝐴
Por tanto 𝑣
= 0 y por ello =0
∞ 𝐴
-Para acelerar un fluido a velocidad subsónica (𝑀𝑎 < 1) se requiere un conducto convergente.
(Tobera del tipo convergente)
-Para acelerar un fluido a velocidad supersónica ( 𝑀𝑎 > 1) se requiere un conducto del tipo
divergente (Tobera1 del tipo divergente)
1
nozzle
74
Figura VII.11.Efecto de la variación del área en el flujo compresible. (White, 2003)
La velocidad de la salida del fluido en una tobera convergente se calcula a partir de la integración de
la ecuación de Euler para flujos isotrópicos de gases perfectos.
k−1
v22 P2 k P0 k
= ∗ [( ) − 1] 𝐄𝐜. 𝐈. 𝟑𝟕
2 ρ2 k − 1 P2
kP
Pero c2 = √ ρ 2 donde 𝑐2 es la velocidad sónica para un gas perfecto.
2
k−1
v2 2 2 P0 k
( ) = Ma22 = [( ) − 1] 𝐄𝐜. 𝐈. 𝟑𝟖
c2 k − 1 P2
75
k
P2∗ 2 k−1
=( ) 𝐄𝐜. 𝐕𝐈𝐈. 𝟑𝟗
P0 k+1
2 k+1
2k P2 k P2 k
ṁ = A2 √ ∗ P0 ρ0 [( ) − ( ) ] 𝐄𝐜. 𝐕𝐈𝐈. 𝟒𝟎
k−1 P0 P0
P2 P∗2
Siempre que >
P0 P0
k
P2 P∗2 2 (k−1)
Pero cuando =P = [(k+1)]
P0 0
(k+1)
A2 P0k 2 ((k−1))
ṁmax = √[ ( ) ] 𝐄𝐜. 𝐕𝐈𝐈. 𝟒𝟏
√T0 R k+1
Solución
Datos
-Diámetro: 10 mm
76
-Temperatura: 15℃
v2
FD = CD ∗ ρ ∗ ∗ A (I)
2
Para determinar el valor del coeficiente de resistencia (𝐶𝐷 ) es necesario conocer el Nro. de Reynolds
dado por la ecuación (𝐼𝐼)
v∗d
Re = (II)
ν
μ
ν = ρ (III)
Dondeμ es la viscosidad absoluta, que para presión atmosférica a nivel del mar, Temperatura igual a
15℃ y densidad 1.225 kg/m3 es:
N kg∗m
μ 1.789 ∗ 10−5 Pa ∗ s ∗s s2 m2
m2
ν= = → kg = kg
= kg
= m2 /s
ρ 1.225
m3 m3 m3
ν = 1.46 ∗ 10−5 m2 /s
v∗d 27 ∗ 0.010
Re = = = 18 493.15
ν 1.46 ∗ 10−5
Análisis dimensional
m m2
v∗d ∗m
s s
Re = = m2
= m2
ν
s s
77
Determinación del coeficiente de resistencia (CD )
CD = 1.2
v2 (27)2
FD = CD ∗ ρ ∗ ∗ A = 1.2 ∗ 1.225 ∗ ∗ (1 ∗ 0.01)
2 2
FD = 5.36 N⁄m
Solución
Datos
C = 2m
B = 10 m
α = 5.4°
v = 90 m/s
hvuelo = 2 600 m
Del diagrama Polar para este tipo de ala (Y de Clark) con el ángulo α se obtiene:
78
-Para una altura de 2600 m, la densidad del aire es:
ρ = 0.9843 kg/m3
Esquema de cálculo
R a-)
W = FL Dado que:
v2
FL = CL ∗ ρ ∗ ∗ A
2
(90)2
FL = 0.8 ∗ 0.9843 ∗ ∗ (10 ∗ 2)
2
FL = 63 782.64 N
Análisis dimensional
kg m2 kg ∗ m
FL = 3 ∗ 2 ∗ m2 = =N
m s s2
63782.64 N
FL = N = 6501.7 kgf
9.81 kgf
79
R b-)
Potencia = FD ∗ v
kg
ρ = 0.9843
m3
C
FD = ( CD ) ∗ FL por lo que:
L
0.047
FD = ( ) ∗ 63782.64
0.8
FD = 3 747.23 N
La potencia es:
Pot = FD ∗ v
Pot = 3 747,23 ∗ 90
Problema 3) Un flujo de aire a 1 MPa y 600 0C entra a una tobera convergente como se muestra en
la Figura I.1 con una velocidad de 150 m/s. Determine el flujo másico por la tobera, si la garganta de
esta posee un área de 50 cm2 cuando la presión en la salida es a) 0,7 MPa b) 0,4 MPa.
Esquema de cálculo
80
Solución
Suposiciones: 1. Aire es un gas ideal con calor específico constante. 2. El flujo es estacionario,
unidimensional e isoentrópico.
Propiedades.
El calor específico a presión constante CP, así como la constante politrópicak tiene valor de Cp =
kJ
1,005 kg.K k = 1,4, la densidad es ρ = 3,093 kg/m3
Para el análisis del problema los subíndices i y t representan las propiedades en la tobera a la
entrada y en la garganta de esta respectivamente. La presión y temperatura de estancamiento a la
entrada de la tobera puede ser determinada por:
kJ
Vi2 (150 m/s2 )2 1 kg
Toi = Ti + = 873 K + kJ
( m2
) = 884 K
2Cp 2(1,005 ) 1 000
kg.K s2
To = Toi = 884 K
k
P∗ 2 k−1
=( ) = 0,5283
P k+1
81
Pb 0,7 MPa
= = 0,670
Po 1,045 MPa
Pb P∗
Como la relación de presiones reales es mayor que la relación de presiones criticas > . La
Po P
2 k−1
2k P2 k P2 k
ṁ = A2 √ ∗ P0 ρ0 [( ) − ( ) ]
k−1 P0 P0
2 1,4+1
2 ∗ 1,4 kg 700 kPa 1,4 700 kPa 1,4
ṁ = 50x 10−4 √ ∗ 1,045 x 106 Pax 3,093 3 [( ) −( ) ]
1,4 − 1 m 1 045 kPa 1 045 kPa
ṁ = 5,88 kg/s
Pb 0,4 MPa
= = 0,383
Po 1,045 MPa
Pb P∗
Como la relación de presiones reales es menor que la relación de presiones criticas < . existe
Po P
condición sónica en la garganta de la tobera Ma=1, el flujo másico en este caso puede ser
determinado por la ecuación:
k+1
k 2 [2(k−1)]
ṁ = A2 P0 √ ( )
RTo k + 1
1,4+1
1,4 2 [2(1,4−1)]
ṁ = 50x 10−4 . 1,045. 106 Pa√ J ( )
287 kgK . 884 K 1,4 + 1
ṁ = 7,1 kg/s
82
VII.7 Problemas Propuestos
Problema 1) Un perfil aerodinámico rectangular. Tiene 9 m de envergadura y 1.8 m de cuerda se
mueve horizontalmente bajo cierto ángulo de ataque a través de aire en calma a 240 km/h .Calcule
la fuerza de sustentación, la de resistencia al avance y la potencia necesaria para impulsar el perfil a
esta velocidad a través de:
1.1-) Calcule la velocidad y la potencia requerida para obtener la condición (b) y la fuerza de
sustentación de la condición (a).
Soluciones
FL = 20 286 N
FD = 1 543,5 N
Pot = 102,9 kW
FL = 17 934,48 N
FD = 1 364,58 N
Pot = 90,97 kW
Se requiere en las condiciones (b) volar a 255,24 km/h y consumir una potencia de 96.7 kW para
lograr la FL de (a).
83
Problema 2) Un perfil aerodinámico rectangular. Tiene 12 m de envergadura y 2 m de cuerda y
coeficiente de sustentación y de resistencia al avance de 0.5 y 0.04 respectivamente bajo un
ángulo de ataque de 6°.
Calcule la fuerza de resistencia al avance y la potencia necesaria para impulsar este perfil a 80,
160,240 km/h en dirección horizontal a través de aire en calma (2℃ y 79.5 kPa) ¿Qué fuerza de
sustentación se obtiene a esas velocidades? Exprese sus condiciones.
Soluciones
b) El momento de flexión en la parte inferior de su apoyo cuya altura desde el suelo hasta la parte
inferior es de 1,5 m. No se tenga en consideración la resistencia del poste.
84
Figura VII.13 Señalización circula sometida a viento normal. (Cengel, 2006)
a) Determine la fuerza de arrastre que actúa sobre las superficies superior y lateral y la potencia
requerida para superar esta resistencia.
Problema 5) Un tronco cilíndrico de pino de 2 m de longitud y 0,2 m de diámetro con una densidad =
513 kg / m3 es suspendido por una grúa en posición horizontal. El tronco está expuesto a una
corriente de aire normal con una magnitud de 40 km/h a 5 0C y 88 kPa. Sin tener en cuenta el peso
del cable y su resistencia, determine:
85
a) El ángulo θ que el cable hará con la horizontal.
b) La tensión en el cable.
Solución a) θ= 840b) T=318 N
Problema 6) Desde un depósito de grandes dimensiones se está descargando aire a 28℃, a través
de una tobera con diámetro de salida de 10mm. La descarga se efectúa a la atmósfera donde la
presión es de 96,5 kPa. La presión del aire en el interior del depósito es de 40,0 kPa (manométrica).
Determine:
El caudal másico máximo de aire que se puede aspirar a través de este tubo.
kg
Solución ṁ = 0.0136 s
86
Capítulo VIII: Equipos de bombeo. Teoría de funcionamiento
Existen dos categorías fundamentales para las turbomáquinas, Bombas y Turbinas, la palabra
bomba es un término general para cualquier máquina de flujo que añade energía al fluido. Algunos
autores llaman a los equipos de bombeo como dispositivos absorbedores de energía. El incremento
de la energía del fluido se traduce en un aumento de la presión. Las turbinas por otro lado son
dispositivos que producen energía, estos extraen energía del fluido y transforman la mayor parte de
la energía en forma mecánica, el fluido a la salida de la turbina sufre pérdidas energéticas,
típicamente en forma de presión. (Mott, 2006; Mott, 1996; Zubicaray, 2005) (Ramírez, 2012)
Por tanto, las máquinas de flujo que mueven líquido son llamadas bombas, aunque existe por
ejemplo otro nombre para las máquinas que manejan gases (Ventiladores) caracterizándose por su
baja presión y grandes volúmenes de flujo manejado. (Robles, 1997) (Pérez, 2004) (García, 2003)
Donde:
kg
ρ. Densidad de la sustancia (m3)
Carga de impulsión (H). Cantidad de energía entregada al fluido por unidad de masa (kJ/kg; m).
El rendimiento de una bomba es caracterizado también adicionalmente por la carga neta (H) definida
como el cambio en la carga de Bernoulli entre la entrada y la salida de la bomba.
P V2 P V2
H=( + + 𝑧) −( + + 𝑧) 𝐄𝐜. 𝐕𝐈𝐈𝐈. 𝟐
ρg 2g salida
ρg 2g entrada
Potencia (N o P). Cantidad de energía por unidad de tiempo entregada al fluido o comunicada de la
red (Putil ; Pinterna Pfreno Pelec )
Donde
g. gravedad (m/s2)
H. Carga de impulsión
88
Curvas de rendimiento de las bombas y punto de operación.
Donde
T. Torque
ρgVH
𝜂= 𝐄𝐜. 𝐕𝐈𝐈𝐈. 𝟓
w. T
E1 + Hb = E2 + ∑ hf 𝐄𝐜. 𝐈𝐈. 𝟔
1
2
Reservoir: Depósito, Valve: Válvula, Pump: Bomba
89
2
P1 V1 2 P2 V2 2
+ Z1 + + Hb = + Z2 + + ∑ hf 𝐄𝐜. 𝐈𝐈. 𝟕
ρg 2g ρg 2g
1
2
P2 − P1 V2 2 − V1 2
Hb = + Z2 − Z1 + + ∑ hf 𝐄𝐜. 𝐕𝐈𝐈𝐈. 𝟖
ρg 2g
1
Las bombas de desplazamiento positivo son ideales para aplicaciones de alta presión, para el
bombeo de líquidos viscosos, así como aplicaciones donde se precisan una cierta cantidad de flujo
90
A continuación, se describe por fases el funcionamiento de una bomba de desplazamiento positivo
de dos lóbulos.
Figura. VIII.4.Cuatro fases de una bomba de desplazamiento positivo de dos lóbulos, la parte clara
representa la entrada del fluido por el rotor superior, así como la parte oscura la entrada por el
inferior. (Cengel, 2006)
91
Los procesos son:
Donde
n. rpm de la bomba
Donde
S S
1-D En general = 0.8 − 2(Relación carrera-diámetro)
D
2Sn Sn
Cmp = = C = 0,5 − 0,9 m/𝑠
60 30 mp
S
Mientras más rápida es la bomba, menor es la relación D
Haciendo uso del diagrama de indicador se aprecia que existe un trabajo en la aspiración y otro en la
descarga
Wcas = Ac . p1 . S
Wcd = Ac . p2 . S
Pi ∗ Ac S ∗ n
Ni = (kW)𝐄𝐜. 𝐕𝐈𝐈𝐈. 𝟏𝟑
60 ∗ 1000
93
Para las bombas de doble efecto y de émbolos múltiples la potencia interior se calcula como la suma
de las potencias interiores de los cilindros por separado.
La potencia real (Potencia al freno o al eje) aplicada desde el motor a la bomba es mayor que la
potencia interna por lo que parte de ella se gasta en vencer el rozamiento mecánico.
Ni Ni
ηmec = Neje =
Neje ηmec
Rendimiento interno 𝜂𝑖
Nu Nu
ηi = Ni =
Ni ηi
Nu ρgQH
Neje = = 𝐄𝐜. 𝐕𝐈𝐈𝐈. 𝟏𝟒
ηi ∗ ηmec ηi ∗ ηmec
Pero
ηi = ηv ∗ ηh Quedando la Ecuación
ρgQH
Neje = 𝐄𝐜. 𝐕𝐈𝐈𝐈. 𝟏𝟓
ηmec ∗ ηv ∗ ηh
Bombas de engranes
94
Figura. VIII.7. Bomba de engranes (Cengel, 2006)
Q (m3/min)
𝑍2
ξ=
𝑍1
95
u1 = u2 = u
w1 ∗ r1 = w2 ∗ r2
2πn1 D1 2πn2 D2
=
60 2 60 2
n1 D1 = n2 D2
Queda entonces:
n1 D2
=
n2 D1
f ∗ l ∗ 𝑍2 ∗ n2
Q = f ∗ l ∗ 𝑍1 ∗ n1 (1 + ) ∗ ηv
f ∗ l ∗ 𝑍1 ∗ n1
D2
Q = f ∗ l ∗ 𝑍1 ∗ n1 (1 + ξ ) η 𝐄𝐜. 𝐕𝐈𝐈𝐈. 𝟏𝟕
D1 v
D1
ξ=1 =1
D2
Q = 2 ∗ f ∗ l ∗ 𝑍1 ∗ n1 ∗ ηv
96
Generalmente el ηv = 0,7 − 0,9
m3
Q = 0,233 l ∗ R ∗ h ∗ n ( ) 𝐄𝐜. 𝐕𝐈𝐈𝐈. 𝟏𝟖
min
Donde:
R. Radio primitivo.
Bombas de Paletas
Donde
n. rpm
97
Figura. VIII.11. Bomba de paletas (Elaboración propia)
La bomba centrifuga es una máquina hidráulica compuesta en esencia por un impulsor con alabes,
que accionado desde el exterior transmite al líquido la energía necesaria para obtener una presión
determinada (Abella, 2003) (Gómez, 1991) (Martínez, 2002) (Santos, 2014). El cuerpo de la bomba o
voluta recibe el líquido salido del impulsor y por su construcción especial transforma su energía
cinética en presión, dirigiéndolo al mismo tiempo hacia el exterior por un canal de descarga. Existe
una gran variedad de bombas en cuanto se refiere a sus formas de construcción, de acuerdo con las
aplicaciones particulares o imperativos mecánicos.
d(m · v)
F=
dt
98
-La acción de las fuerzas externas es igual a la variación en el tiempo de la cantidad de movimiento.
Isaac Newton3 planteó exactamente ``El cambio del movimiento en el tiempo es proporcional a la
fuerza matriz aplicada y se efectúa en la dirección de la línea recta en que se aplica esa fuerza´´).
El término usado por Newton como ``movimiento´´ hoy corresponde al término cantidad de
movimiento.
F dt = d (m . v)
Fdt = mdv
2 2
∫ F · dt = m ∫ dv = m(V2 − V1 )
1 1
F∆t = m(V2 − V1 )
𝑚
N · s = Kg · 𝑠
Kg · m m
2
· s = kg ·
s s
3
Isaac Newton (1643-1727). Físico, filósofo, teólogo, inventor, alquimista y matemático inglés.
99
Suposiciones:
Donde:
w1 : Velocidad angular.
100
Figura VIII.12. Paralelogramos de velocidad a la entrada y la salida de los canales interálabes.
(White, 2003)
F · r · ∆t = m(v2 r2 − v1 r1 )
Mt ∞ · ∆t = mv2 r2 − mv1 r1
Mt ∞ · ∆t = (ρQC2 l2 − ρQC1 l1 )
Mt ∞ = ρQ(C2 l2 − C1 l1 )
l2 = R 2 cos ∝2
l1 = R1 cos ∝1
Por tanto.
101
Mt ∞ = ρQ(C2 R 2 cos ∝2 − C1 R1 cos ∝1 )
Mt ∞ = ρQ(C2u R 2 − C1u R1 )
Nt ∞ = Mt ∞ · w
Por tanto.
Nt ∞ = ρQ(C2u R 2 − C1u R1 ) · w
u
Como: u = w · R queda entonces w=R
Nt ∞ = ρgQ · Ht ∞
Por lo que:
𝐮𝟐 𝐂𝟐𝐮 − 𝐮𝟏 𝐂𝟏𝐮
𝐇𝐭 ∞ = 𝐄𝐜. 𝐕𝐈𝐈𝐈. 𝟐𝟎
𝐠
Ht = μ · Ht ∞ 𝐄𝐜. 𝐕𝐈𝐈𝐈. 𝟐𝟏
102
Donde 𝜇 : Coeficiente de corrección que tiene en cuenta el número finito de álabes.
𝝁 < 1 Generalmente 𝛍 = 𝟎, 𝟖
Hr = ηh · Ht = ηh · μ · Ht ∞ 𝐄𝐜. 𝐕𝐈𝐈𝐈. 𝟐𝟐
u2 C2u − u1 C1u
Ht ∞ =
g
u2 −C2u
Del triángulo de velocidades cot β2 = C2r
Por tanto.
u2 2 − u2 C2r cot β2
Ht ∞ =
g
u2 2 u2 2
Cuando β2 = 90°Ht ∞ = -Para β2 = 90°Ht ∞ =
g g
103
La carga estática por:
u2 − (C2r cot β2 )2
(Hest)t ∞ = 𝐄𝐜. 𝐕𝐈𝐈𝐈. 𝟐𝟒
2g
Figura VIII.13. Dependencia de 𝐻𝑡∞ en función del ángulo 𝛽2. (Elaboración propia)
Figura VIII.14Representación gráfica de los distintos tipos de carga en función del ángulo 𝛽2 .
(Elaboración propia)
Curva característica de H vs Q.
u2 2 − u2 C2r cot β2
Ht ∞ =
g
104
Pero Q = C2r · Área
Q = C2r · π · D2 · b2
Q
u2 2 − u2 · π·D · cot β2
2 ·b2
Ht ∞ =
g
Q
C2r =
π · D2 · b2
ρgQH
N= 𝐄𝐜. 𝐕𝐈𝐈𝐈. 𝟐𝟓
K
Representación gráfica.
Figura VIII.15 Características teóricas de las máquinas centrífugas para distintos ángulos
β2 .(Elaboración propia)
105
Potencia y Rendimiento.
ρgQH
Nu = (kW) 𝐄𝐜. 𝐕𝐈𝐈𝐈. 𝟐𝟔
1000
Ni
ρg(Q + ∆Q) · (H + ∆H)
= 𝐄𝐜. 𝐕𝐈𝐈𝐈. 𝟐𝟕
1000
Rendimiento interno. ( 𝜂𝑖 )
ρgQH
Nu 1000
ηi = =
Ni ρg(Q+∆Q)·(H+∆H)
1000
ηi = ηv · ηh
106
Ni
ηm =
Neje
Ni
Neje =
ηm
Nu Nu
Pero Ni = =
ηi ηv · ηh
Nu Nu
Ne = =
ηv · ηh · ηm ηmaq
ρgQH
Ne = 𝐄𝐜. 𝐕𝐈𝐈𝐈. 𝟐𝟖
1000 · ηmaq
VIII.5Teoría de Semejanza
El movimiento de los líquidos y los gases en los pasajes interálabes de las maquinas centrífugas es
complicado. Es por esta razón que el cálculo de los elementos de trabajo de las máquinas
centrifugas presentan grandes dificultades y se acude por ello a los datos experimentales obtenidos
del diseño y prueba de máquinas anteriores. No obstante, la posibilidad de usar los datos
experimentales está dada únicamente porque se hayan observado las Leyes de Semejanza.
En las máquinas de flujo hay que observar tres tipos o formas de semejanza:
Semejanza geométrica.
Semejanza cinemática.
Semejanza dinámica.
107
Semejanza geométrica:
β1a = β1b = β1
β2a = β2b = β2
βia = βib = βi
De donde:
Semejanza cinemática:
α1a = α1b
α2a = α2b
αia = αib
108
Semejanza dinámica:
Se mantiene constante la relación entre las fuerzas de igual naturaleza que actúan en
puntos homólogos de las máquinas geométricas y cinemáticamente semejantes.
En la semejanza dinámica
Los criterios generales de semejanza de los flujos estudiados en la Mecánica de los Fluidos se
aplican también para los flujos en las máquinas centrifugas. Por lo cual al estudiar la semejanza de
estas máquinas debe observarse la igualdad de estos criterios adimensionales.
V∙l
Número de Reynolds -------------------R e = υ
V2
Número de Froude ----------------------Fr = g∙l
P
Número de Euler-------------------------Eu = ρ∙V2
l∙n
Número de Strouhal---------------------Sh = V
La semejanza de las corrientes de flujo se caracteriza por la igualdad de los criterios adimensionales
en puntos homólogos de las maquinas:
109
De la semejanza cinemática se puede escribir la relación de velocidad C2ra ⁄C2rb en función de las
dimensiones físicas del impelente y la velocidad de rotación.
Qa D22a b2a na
= ∙ ∙
Qb D22b b2b nb
ηva
∙ 𝐄𝐜. 𝐕𝐈𝐈𝐈. 𝟐𝟗
ηvb
b2a D2a
=
b2b D2b
Qa D32a na ηva
= ∙ ∙ 𝐄𝐜. 𝐕𝐈𝐈𝐈. 𝟑𝟎
Qb D32b nb ηvb
Para el análisis de una misma máquina D2a = D2b y a diferentes rpm la expresión (𝐕𝐈𝐈𝐈. 𝟑𝟎)
queda:
Qa na ηva
= ∙ 𝐄𝐜. 𝐕𝐈𝐈𝐈. 𝟑𝟎. 𝟏
Qb nb ηvb
La 𝑉𝐈𝐈𝐈. 𝟑𝟎. 𝟏 nos indica que para una máquina dada el caudal varía con la primera potencia de las
frecuencias de rotación del árbol y los rendimientos volumétricos.
Como se sabe la ecuación de Euler la Carga de Impulsión creada por la máquina centrifuga es:
U2∙ C2u
H= ∙ ηh 𝐄𝐜. 𝐕𝐈𝐈𝐈. 𝟑𝟏
g
110
Ha U2a C2ua ηha
= ∙ ∙ 𝐄𝐜. 𝐕𝐈𝐈𝐈. 𝟑𝟏. 𝟏
Hb U2b C2ub ηhb
U2a C U n ∙D
= C2ua yU2a = n a ∙D2a
U2b 2ub 2b b 2b
Para la misma máquina funcionando en diferente régimen; dado que D2a = D2b queda:
Ha n2a ηha
= ∙ 𝐄𝐜. 𝐕𝐈𝐈𝐈. 𝟑𝟑
Hb n2b ηhb
Para el caso de los ventiladores se acostumbra a expresar esta relación en función de presión y no
de las cargas de impulsión. Para ello multiplíquese la ecuación 𝐕𝐈𝐈𝐈. 𝟑𝟐 por los términos de ρg.
Para la misma máquina, siendo diferente la frecuencia de rotación del árbol, se tiene:
Pa n2a ρa ηha
= ∙ ∙ 𝐄𝐜. 𝐕𝐈𝐈𝐈. 𝟑𝟒. 𝟏
Pb n2b ρb ηhb
Como ya se había estudiado la potencia de la máquina de flujo centrífuga está dada por la ecuación:
ρ∙g∙Q∙H
N= 𝐄𝐜. 𝐕𝐈𝐈𝐈. 𝟑𝟓
ηmaq
Estableciendo la relación entre los consumos de potencia de las máquinas semejantes trabajando en
regímenes semejantes se llega a la ecuación:
Na ρa ∙ Qa ∙ Ha ∙ ηb
= 𝐄𝐜. 𝐈𝐈. 𝟑𝟔
Nb ρb ∙ Qb ∙ Hb ∙ ηa
111
Como se tiene en las ecuaciones 𝐕𝐈𝐈𝐈. 𝟐𝟗 y 𝐕𝐈𝐈𝐈. 𝟑𝟐 las relaciones de Q y H en función de las
dimensiones físicas y la frecuencia de rotación la ecuación 𝐕𝐈𝐈𝐈. 𝟑𝟔 se expresa por:
Na ρa D52a n3a ηa
= ∙ ∙ ∙ 𝐄𝐜. 𝐕𝐈𝐈𝐈. 𝟑𝟕
Nb ρb D52b n3b ηb
Na ρa D52a Na3 ρa ηa
= ∙ ∙ ∙ ∙ 𝐄𝐜. 𝐕𝐈𝐈𝐈. 𝟑𝟖
Nb ρb D52b D3b ρb ηb
Para fluidos incompresibles, la misma máquina trabajando a diferentes frecuencias de rotación tiene
variaciones en el consumo de potencia según la ecuación:
Na n3a ηb
= ∙ 𝐄𝐜. 𝐕𝐈𝐈𝐈. 𝟑𝟗
Nb n3b ηa
Una magnitud muy importante, que determina la semejanza de las corrientes de flujo en las bombas
centrífugas es la Velocidad Específica.
Qa D32a Na
= ∙ 𝐄𝐜. 𝐕𝐈𝐈𝐈. 𝟒𝟎
Qb D32b Nb
Ha n2a D22a
= ∙ 𝐄𝐜. 𝐕𝐈𝐈𝐈. 𝟒𝟏
Hb n2b D22b
Si a la máquina (a) en escala reducida o modelo, se define que entrega la potencia útil (Nu) de 1CV
(caballo de vapor) y desarrolla una carga de impulsión de 1m, el flujo entregado por ella es:
112
Nua
Qa = 𝐄𝐜. 𝐕𝐈𝐈𝐈. 𝟒𝟑
ρ ∙ g ∙ Ha
N∙m
1CV ∙ 735.5 m3
s∙cv
Qa = Kg m =
103 m3 ∙ 9.81 s2 ∙ 1 s
3
Qa = 0.07497 ≈ 0.075 m ⁄s
Si se denomina (ns ) “Velocidad Específica” a la frecuencia de rotación a que debe guiar la máquina
3
(a) para entregar el flujo o caudal (Q) = 0.075 m ⁄s, a la carga Ha = 1m. Sustituyendo en la
D
ecuación 𝐕𝐈𝐈𝐈. 𝟑𝟐y despejando la relación de diámetro Da se obtiene:
b
Ha
D22a Hb η
= n2
Considerando ηha = 1. ∴ ηℎ𝑎 ≈ ηℎ𝑏
D22b a ∙ ηha hb
n2 η
b hb
1⁄
Ha H 2
D2a (Ha )
H
= √ n2b = b
1⁄ 𝐄𝐜. 𝐕𝐈𝐈𝐈. 𝟒𝟒
D2b a
n2a 2
n2b (n2 )
b
1
1⁄
D2a Hb 2
= ns 𝐄𝐜. 𝐕𝐈𝐈𝐈. 𝟒𝟓
D2b
nb
113
Qa D32a na ηva
= ∙ ∙ 𝐄𝐜. 𝐕𝐈𝐈𝐈. 𝟒𝟔
Qb D32b nb ηvb
1 3
1⁄
0.075 Hb 2 ns
=( ns ) ∙ 𝐄𝐜. 𝐕𝐈𝐈𝐈. 𝟒𝟕
Qb nb
nb
Desarrollando:
1 1
3⁄ 3⁄
0.075 Hb 2 ns Hb 2
= ∙ =
Qb n3s nb n2s
n3b n2b
0.075 n2b
= 3⁄
Qb
Hb 2 ∙ n2s
Despejando (ns ):
n2b ∙ Qb
n2s = 3⁄
0.075 ∙ Hb 2
1⁄
2
n2b
∙ Qb
Ns = ( 3⁄ )
0.075 ∙ Hb 2
1⁄
nb ∙ Q b 2
Ns = 1⁄
1⁄ 3⁄ 2
0.075 2 ∙ (H 2)
1⁄
1 nb ∙ Q b 2
Ns = 1⁄ ∙ 3⁄
0.075 2 H 4
nb ∙√Qb
Ns = 3.65 ∙ 3⁄ (Ecuación de la Velocidad Específica)
Hb 4
114
La velocidad específica caracteriza la forma constructiva de la máquina.
Por el contrario, para ↓ Q y ↑ H → Ns ↓ – Típico en máquinas radiales con canales estrechos y alta
D2 ⁄D1 .
En la Tabla VIII.1 se ofrecen a modo de ejemplo, los intervalos de velocidades específicas para
distintos tipos de bombas.
Tabla VIII.1. Intervalos de velocidades específicas para distintos tipos de bombas4. (Ramos, 1986)
De remolino 10-40
Centrifugas 40-300
Mixtas 300-600
Axiales 600-1200
Nota: Estos valores dependen de las unidades empleadas, los dados corresponden al sistema
internacional: Q (m3/s), H (m), N (kW)
a) La fuerza que se requiere para mover el émbolo en las carreras de ida y vuelta.
4
Los valores dependen de las unidades empleadas, los dados corresponden al sistema internacional
115
b) El caudal entregado por la bomba
c) La potencia absorbida por la bomba.
En el problema desprecie las pérdidas y el rozamiento.
Solución
Datos
Démbolo= 250 mm
dvástago=50 mm
S=375 mm (carrera)
n=60 rpm.
Psuc
= −4,5 m c. a
ρg
Pdesc
= 18 m c. a
ρg
πD2 π. 0.2502
Ap = = = 0,049 m2
4 4
πd2 π. 0.0502
Av = = = 0,0019 m2
4 4
P
Si hPa = −4,5 m c. a y hPa = ρga se tiene que:
kg
Pa = ρghPa = 1000 . 9,81. (−4,5) = −44 145 Pa
m3
116
Pd
Si hPd = 18 m c. a y hPd = se tiene que:
ρg
kg
Pd = ρghPd = 1000 . 9,81. 18 m c. a = 176 580 Pa
m3
Aunque la presión es negativa hay que hacer una fuerza para mover el pistón.
m3
Qi = 0,01767
s
117
𝜋 2 𝜋 60 m3
Qd = Area. S = D . S. n = 0,2502 . 0,375 . s = 0,01841
4 4 60 s
min
Caudal total
m3 m3 l
Qt = Qi + Qd = 0,01767 + 0,01841 = 0,03608 = 36,08
s s s
ρ. g. Q. H
Nb =
1000
Cálculo de H (Carga)
Pdesc − Psuc
H= = 18 − (−4,5) = 22,5 m
ρg
kg m3
1000 . 9,81 x 0,03608 x 22,5 m
m3 s
Nb = = 7,96 kW
1000
Problema 2) Una bomba de émbolo de simple efecto tiene un diámetro de pistón de 180 mm y una
carrera de 200 mm. El líquido a bombear tiene una densidad relativa (ρrel = 0,96). La bomba toma
líquido de un colector a presión atmosférica y los descarga en un recipiente a una presión de 4 at. La
altura geométrica de impulsión o descarga es de 15 m. La pérdida total en la tubería de succión es
de 2 m y en la descarga es de 7m.
¿A cuántas rpm debe operar la bomba y de que potencia debe instalarse el motor para suministrar
un caudal de 600 l/min?
Datos
Démbolo= 180 mm
118
S=200 mm (Carrera)
ρrel= 0,96
Pd= 4 at.
Zd= 15 m.
∑ hf succ = 2 m
∑ hf des = 7 m
Q=600 l/min.
ηv = 0,85
ηins = 0,76
Solución
𝜋 2
Q= D . S. n ηv
4
4Q 4. 600
n= = = 138,69 rpm~ 140 rpm
π D . S. ηv 1000 . π . 0,182 . 0,2 . 0.85
2
Cálculo de la potencia
ρ. g. Q. H
Nb =
1000 ηins
Determinación de la carga H
2
∆P ∆V 2
H= + ∆Z + + ∑ hf
ρ. g 2g
1
(4 − 0). 104
H= + 15 + 2 + 7 = 65,66 m
960
119
La potencia es:
Problema 3) un rodete que gira a 1 150 rpm tiene las siguientes dimensiones: b 1=31,75 mm;
b2=19,05 mm; D1=177,8 mm; D2=381,0 mm; β1=180 ; β2=200 (b1 y b2 son respectivamente, la
anchura de los canales de entrada y a la salida; la superficie de sección transversal es A = πDb si
se desprecia el espesor de los álabes). Suponiendo que la corriente de entrada es radial y
despreciando el espesor de los álabes, dibujar los triángulos de velocidades y calcular el caudal en
litros por minutos, así como la carga teórica para z=ɷ, despreciando la corriente circulatoria.
Datos
n= 1150 rpm
D2=381 mm
D1=177,8 mm
β2=200
β1=180
b1=31,75 mm
b2=19,05 mm
Z=ɷ
Solución
120
Cálculo de U1
2 π n D1
U1 = w . r1 =
60 2
π n D1 π. 1150 . 0,1778 m
U1 = = = 10,7
60 60 s
Cr1 m m
tan β1 = Cr1 = U1 . tan β1 = 10,70 tan 180 = 3,47
U1 s s
Cálculo de Q1
m3
Q1 = A1 . Cr1 = π D1 b1 . Cr1 = π . 0,1788 . 0,03175 x3,47 = 0,0615
s
Q1 = 3690 l/min
Cálculo de U2
2 π n D2
U2 = w . r2 =
60 2
π n D2 π. 1150 . 0,381 m
U2 = = = 22,94
60 60 s
m3
Q Q 0,0615 s
Q = A2 . Cr2 Cr2 = = =
A2 π D2 b2 π . 0,381 . 0,01905
121
m
Cr2 = 2,69
s
Cálculo de C2u
U2 − C2u
cot β2 =
Cr2
m
C2u = 22,94 − 2,69 cot 200
s
m
C2u = 15,6
s
Ht∞ = 36,57 m
Problema 4) Se desea bombear 3 785 litros por hora de gasolina (densidad 0,85) a la temperatura
de 10 0C (viscosidad= 0,8 centipoise) desde la parte inferior de un depósito al nivel del suelo a la
parte superior de una columna de fraccionamiento de 27,43 m de altura. La longitud total de la
tubería (tubería+longitud equivalente debida a los accesorios) es de 38,1 m. Se trata de una tubería
de acero corriente de 2 pulgadas (50,8 mm) de diámetro. Calcular la potencia del motor requerido, si
el rendimiento del conjunto motor-bomba se supone del 55 %. Considere iguales las presiones en el
tanque y la parte superior de la torre de fraccionamiento.
Datos
Sustancia: Gasolina
t= 10 0C
122
μ=0,8 cp.
ΔZ=27,43 m
Ltotal=LT.RECTA+LEQ=38,10 m
ηmotor+ηbomba=55%
Solución
2
P1 v12 P2 v22
+ Z1 + + Hb = + Z2 + + Hb + ∑ hf
ρg 2g ρg 2g
1
Hb = Z2 − Z1 + ∑ hf
1
2
L v2
∑ hf = λ
d 2g
1
Hay que determinar el factor de fricción λ por lo que es necesario conocer el número de Reynolds
v. d ρ. v. d
Re = =
υ μ
Q = A. v
Q 4Q 4. 1,051 x 10−3 m
v= = 2
= = 0,026
A π. d π. 0,0508 s
El coeficiente de fricción λ es
64 64
λ= = = 0,0456
R e 1 403,35
2
L v2 38,10 0,0262
∑ hf = λ = 0,0456
d 2g 0,0508 2 . 9,81
1
∑ hf = 8,11 10−7
1
Hb = Δ𝑍 + ∑ hf = 27,43 m
1
Problema 5) Se necesita una bomba para un caudal de 4 540 lpm para una carga de impulsión de
91,5 m, acoplada a un motor que gira 1150 rpm. Probablemente será más indicada una bomba de 2
rodetes, que una de rodete único. ¿Por qué? Argumente su respuesta.
Datos:
Q = 4540 lpm.
124
H = 91,5 m.
N = 1150 rpm.
Solución:
n√Q
ᶯs =
H 3/4
1150√4540/60
ᶯs =
(91,5)3/4
ᶯs = 338 rpm
Para este valor de velocidad específica el rendimiento estaría en el orden del 63%, extrapolando la
línea de caudal en el gráfico de Worthington Pump & Machinery Corp. (Vennard, 1979)
HT 91,5
H= = = 45,75 m
2 2
n√Q
ᶯs =
H 3/4
1150√4540/60
ᶯs =
(45,75)3/4
ᶯs = 570 rpm
ɳb = 70%
Si es más indicada una bomba de dos etapas, el impelente para un radial de menor relación D2/D1
que resulta más eficiente y la máquina alcanza un mayor rendimiento durante la explotación.
125
VIII.7 Problemas propuestos
Problema 1) Se expone el diseño de una bomba de engranaje (Figura VIII.17) que se emplea en los
sistemas hidráulicos de desplegar y replegar el chasis del avión. Determinar el suministro de la
bomba para la diferencia de presiones igual a cero y n= 2 200 rpm, si se sabe que el área máxima
de sección de cada diente limitada por la circunferencia exterior del engranaje vecino es igual a
S=0,2 cm2. La anchura del engranaje (Longitud del diente) es b=12 mm.
Problema 2) Una bomba de émbolo que trabaja a 150 rpm debe impulsar agua calentada hasta 60
0C. Los cálculos previos demostraron que el consumo de energía para crear la velocidad, para
separar las pérdidas de inercia y las resistencias hidráulicas de la línea de succión constituyen en
total 6,5 m de agua, la presión atmosférica media en el lugar de instalación de la bomba es de 736
mm de Hg. ¿A qué altura sobre el nivel del mar se debe de instalar la bomba?
SoluciónZg=1,66 m
Problema 3) Una bomba centrífuga tiene un rodete de trabajo de diámetro D 2=250 mm con seis
álabes radiales (β2= 900), el diámetro de la circunferencia formado por los bordes de entrada de los
álabes es D1=125 mm. ¿Qué número de rpm hay que comunicar para obtener una diferencia de
presiones ΔP= 3 kgf/cm2? El rendimiento hidráulico de la bomba es ηh=0,75.
126
Problema 4) Si el rodete del problema resuelto 3 trabaja en vacío (C2r=0) ¿Cuál es la carga teórica
infinita? ¿Qué le ocurre al ángulo α2? Tiene ello que ver con algunos de los rendimientos de la
máquina. Explique.
SoluciónHt= 53,92 m
Problema 6) Se desea una bomba para un caudal de agua de 10 600 l/min y una altura de elevación
53,34 m. a) ¿Cuál sería la velocidad aproximada de funcionamiento que nos daría el mejor
rendimiento? ¿Con que tipo de impelente?
Problema 7) Se necesita una bomba para un caudal de 4 540 l/min para una altura manométrica de
91,5 m acoplada a un motor que gira a 1 150 rpm. Probablemente será más indicada una bomba
de dos rodetes (impelentes) que una de rodete único. ¿Por qué?
Problema 8) De la curva de ensayo de una bomba del tipo de flujo mixto de un solo rodete con
aspiración unilateral funcionando a 575 rpm y calculada para dar un caudal de 378,5 litros de
agua por segundo, se han obtenido los siguientes resultados.
Magnitud 1 2 3 4 5 6
127
Η (%) 0 23,0 45,0 75,5 88,0 81,0
a. Trace las gráficas de las características hidráulica y dibuje la curva de potencia absorbida.
b. En la misma hoja dibuje las tres curvas con líneas de trazos discontinuos para el caso en que el líquido
bombeado sea de características análogas a las del agua, pero de densidad relativa 0,9. Razone de
manera resumida cada una de las curvas.
c. Trace en la misma hoja las tres curvas de funcionamiento de la bomba trabajando a una
velocidad de 500 rpm.
Problema 9) ¿Qué potencia requerirá una bomba para entregar un caudal de 1 200 l/min venciendo
una carga de 30 metros? Cuando gira a 1750 rpm.
Problema 10) Hay que suministrar 1 225 l/min contra una carga de 126 m a 3 600 rpm. Suponiendo
un rendimiento aceptable de la bomba a valores de velocidad específica establecida por usted
¿Cuántas etapas de bombeo se necesitarán?
Solución4 Etapas
Problema 11) Una bomba centrífuga tiene un impelente de dimensiones: r1 = 75 mm; r2 = 200 mm;
b1 = 50º; b2 = 40º. La anchura del rodete a la entrada es, b 1 = 40 mm y a la salida, b2 = 20 mm. Se
puede suponer que funciona en condiciones de rendimiento máximo (c1m = c1). Considere el
rendimiento en un 78 %. Determinar, para un caudal Q = 0,1 m3/s lo siguiente:
Solución: b) Htmax=8,62 m
128
equivalen a 5 veces la energía cinética relativa, a la salida del agua de la bomba, y las pérdidas en la
tubería equivalen a 15 q2. El diámetro a la salida de la bomba es d2 = 0,2 m, y la sección útil de
salida del rodete es A2 = 0,2 d22. El rendimiento puede ser considerado como un 75 %. Determinar:
d) La potencia útil.
Solución a) β2=34,40b) u2=31,9 m/s, n=3 044 rpm, c) Nu= 32 CV, d) ns=119, 07 rpm
Problema 13) Una bomba centrífuga tiene las siguientes características: q = 50 litros/s;Hm = 100 m;
n = 1500 rpm; nm= 0,67; N = 100 CV. Se quiere bajar para utilizarla en una mina en donde va a
funcionar a un mayor número de revoluciones.
Determinar:
Solución Puede auxiliarse para la solución de los criterios de semejanza a) Hmaq1=200 mb)
N1=282,8 CV, c) Q1=70,7 litros/s
Problema 14) Dado un modelo de bomba centrífuga de 1 000 CV y n = 1200 r.p.m. cuya curva
característica es de la forma, Hm = 180 - 375 q2, se acopla a una tubería de impulsión que presenta
una curva característica Δe =15 q2.
Determinar:
129
b) El número de revoluciones por minuto necesarias, si las pérdidas de carga en la tubería aumentan
a 8 veces la inicial.
d) El radio r2
e) Los triángulos de velocidades a 1 200 rpm, sabiendo que: α1 = 90º; β2 = 40º; b2 = 0,05 r2; r1 = 0,3
r2
Problema 15) Una bomba centrífuga trabaja a n = 1500 rpm, e impulsa agua de un pozo a un
depósito por medio de una tubería, cuya curva característica viene dada por:
Hm (metros) 8 9 10 11 12
Q (litros/s) 0 6 9,3 11,4 12,8
y la curva característica de la bomba, para las revoluciones de trabajo por:
Hm (metros) 0 2 4 6 8 10 9,4
Q (litros/s) 12 10 7,2 2,5 0
a) Determine el número de rpm que hay que comunicar a esta bomba para aumentar el gasto, en la
tubería mencionada, al doble.
Problema 16) En el ensayo de una bomba centrífuga con agua, que tiene iguales las cotas y
diámetros de aspiración e impulsión, se tomaron los siguientes resultados:
Presión de impulsión, 3,5 Kg/cm2; Presión de aspiración, 294 mm de columna de mercurio; Caudal,
q= 6,5litros/seg; Torque, C= 4,65 m. Kg; Número de revoluciones por minuto, n= 800
Determinar:
130
a) La potencia efectiva en CV
Problema 17) Se instala una bomba centrífuga para elevar agua de manera que las bridas de
aspiración e impulsión son horizontales y la de impulsión está medio metro más elevada que la de
aspiración. Un vacuómetro de mercurio conectado a la aspiración marca una depresión de 300 mm
mientras que el manómetro de impulsión marca 19 m. c.a. La tubería de aspiración es de 225 mm, y
la de impulsión de 200 mm. El rodete o impelente tiene un diámetro exterior de 300 mm y un ancho a
la salida de 25 mm; b2 = 22º; n = 1.320 rpm. En los cálculos se supondrán álabes afilados a la salida;
ηh= 0,80; ηm = 0,85 y entrada de la corriente en los álabes sin circulación.
Calcular:
c) El caudal
d) La potencia útil
e) La potencia de accionamiento
Conclusiones.
1. La obra cuenta con ocho capítulos donde se presenta la teoría de cada uno de los temas
tratados como base para la solución de los problemas resueltos y propuestos.
2. Cada uno de los problemas desarrollados parte de la metodología seguida en clases con el
propósito de crear en los estudiantes la habilidad en la solución de problemas aplicados a la
mecánica de los fluidos.
131
3. Los problemas propuestos presentados complementan el autoestudio en función de un
mejor aprovechamiento por parte de los estudiantes.
132
Apéndices.
Simbología
Γ- Circulación gpm.- Galones/minutos
𝜏- Torque lpm (l/min) .- litros/minutos
ρ.- Densidad de la sustancia (kg/m ) 3 Tv.- Tensión del vapor
M0.- Masa de la sustancia (kg) ns.- Velocidad específica (rpm)
V.- Volumen (m ) 3 η.- Rendimiento
W.- Peso de la sustancia (N) ηh .- Rendimiento hidráulico
γ.- Peso específico (N/m ) 3 ηm .- Rendimiento mecánico
Hb.- Carga de la bomba (m) CV .- Caballos de Vapor
υ.- Viscosidad cinemática (m /s) 2 mph.- Velocidad en Millas por hora
P.- Presión (Pa)
F.- Fuerza (N)
A.- Área (m )
2
REFERENCIAS BIBLIOGRÁFICAS
CRANE Co. (1989). Flujo de fluidos en válvulas, accesorios y tuberías. México: Mc- Graw-Hill.
Giles, R. V., Evett, J. B., & Liu, C. (1994). Mecánica de los Fluidos e Hidráulica. Madrid: Mc- Graw-
Hill.
Howell R, J. & Buckius, O. R. (1990). Principios de termodinámica para Ingeniería. México: Mc-
Graw-Hill.
Mataix, C. (1986). Mecánica de los Fluidos y Máquinas Hidráulicas (segunda). Madrid: Ediciones del
Castillo.
Mott, R. L. (2006). Mecánica de los Fluidos (sexta edición). México: Pearson Education.
Streeter, V. L., Wylie, E. B., & Bedford, K.W. (2008). Mecánica de Fluidos (Novena Edición, Vol. 1).
México: Mc- Graw-Hill.
Vennard, J. K., &Street, R. L. (1983). Elementos de Mecánica de Fluidos. México: CECSA.
133
Abella, M. A., Chenlo, F., Blanco, J., & España, A. U. L. L. T. (2003). Optimización del uso de
convertidores de frecuencia con bombas centrífugas y motores trifásicos en sistemas de bombeo
fotovoltaico. Madrid, España.
Çengel, Y. A., Cimbala, J. M., & Skarina, S. F. (2006). Mecánica de fluidos: fundamentos y aplicaciones
(Vol. 1). McGraw-Hill.
de las Heras Jiménez, S. (2011). Fluidos, bombas e instalaciones hidráulicas. Universidad Politécnica
de Catalunya. Iniciativa Digital Politécnica.
García, C., Robaina, A. D., & Calgaro, M. (2003). Selección de bombas para riego. INIA Serie Técnica.
Giles, R. V., Evett, J. B., & Liu, C. (1999).Mecánica de los fluidos e hidráulica. McGraw Hill.
Gómez, J. S., & Zamora, R. E. (1991). Estimación de costos de intercambiadores de calor, tanques y
bombas centrífugas. Tecnología y Ciencia Ed. (IMIQ), 7(1), 10-19.
James, D., & Donald, H. Dinámica de los Fluidos con Aplicaciones a la Ingeniería, Ed. Trillas, México,
l981
Martínez, L, E. R., & Pérez Robalino, B. A. (2002). Mejoras al funcionamiento de bombas centrifugas
cuando se combinan con conversores de frecuencia (Tesis de Grado).
Mataix, C. M., Giles, R. V. E., Liu, J. B., Cheng, M. J (1982). Mecánica de fluidos y máquinas
hidráulicas. Oxford University Press:.
Munson, B. R., Young, D. F., & Okiishi, T. H. (1999).Fundamentos de mecánica de fluidos. Limusa.
Pérez-Barreto, R. (2004). Procedimiento para evaluar y seleccionar bombas centrífugas. Minería &
Geología, 20(1-2), 84-88.
Potter, M. C., Wiggert, D. C., Hondzo, M., & Shih, T. I. (2002).Mecánica de fluidos. Thomson.
Ramírez, E. A., & Mayorga, A. (2012). Estudio de bombas centrífugas y su factibilidad de aplicación en
la Facultad de Ingeniería Mecánica con la finalidad de obtener parámetros técnicos con variación
134
de caudal (Tesis de Grado, Universidad Técnica de Ambato. Facultad de Ingeniería Civil y
Mecánica. Carrera de Ingeniería Mecánica).
Royo, C., HASTA, E., Royo, C. C., Hasta, C. M., & Daugherty, R. L. (1996). Bombas centrífugas.
Organización Mundial de la Salud.
Santos, R., & Guerrero, J. A. M. (2014). Procedimiento para la elaboración de la curva de NPSHr en
bombas centrífugas a partir de parámetros adimensionales. Ciencia e Ingeniería, 35(2), 77-84.
Shames, I. H., Moneva, J. M., & Crusells, S. P. (1967).La mecánica de los fluidos (Vol. 1). McGraw-Hill.
Smith, J. M., Van Ness, H. C., & Abbott, M. M. (1980).Introducción a la termodinámica en ingeniería
química (No. 660.296 9 S724i 2003.). McGraw-Hill.
Soler, M. A. (2009). Fundamentos de clasificación de las bombas centrífugas. Tecnología del Agua,
(310), 22-34.
Streeter, V. L., Wylie, E. B., Bedford, K. W., & Saldarriaga, J. G. (1963). Mecánica de los fluidos.
Ediciones del Castillo.
Vennard, J. K., Fox, R. W. M., Alan, T., Brun, E. A., & A Mathieu, J. (1979).Elementos de mecánica de
fluidos: versión si. Compañía Editorial Continental,.
135
Anexos
Simbología
ρ.- Densidad de la sustancia (kg/m3) h1-2.- Pérdida de carga (m)
M0.- Masa de la sustancia (kg) Z.- Altura geométrica de presión (m)
V.- Volumen (m3) L1.- Longitud de tubería de impulsión
(m)
W.- Peso de la sustancia (N) La.- Longitud de tubería de aspiración
(m)
𝛾.- Peso específico (N/m3) Hb .- Carga de la bomba (m)
δ.- Densidad relativa ε.- Rugosidad relativa (mm)
β.-Coeficiente de compresión volumétrica (m2/N) ξ .- Coeficiente de Resistencia en los
accesorios
∆V λ.-Coeficiente de fricción en tramo recto
136
Sistemas de Unidades
Al iniciar el estudio de la Mecánica de los Fluidos es necesario realizar un recordatorio de los
Sistemas de Unidades. En los estudios de Física y otras asignaturas básicas específicas el
estudiante ha trabajado con dimensiones y unidades, pero la experiencia demuestra que,
insuficiencias en su conocimiento y aplicación, lleva a errores conceptuales inaceptables en el
campo de la ingeniería.
El desarrollo histórico de la ciencia en diferentes países, ha hecho que en la actualidad coexistan
varios sistemas de unidades, ellos son:
Sistema Inglés (Usado en Inglaterra, Estado Unidos y países en la esfera de influencia de la
ciencia y la ingeniería).
Sistema Giorgi5 o Sistema MKS (metro, kilogramo, segundo) (Usado en Países Europeos
tales como Alemania, Rusia, España, Francia, y países en la esfera de influencia de la
ciencia y la ingeniería). Una variedad de este sistema usado desde la antigüedad por los
físicos es el Sistema Cegesimal (centímetro, gramo, segundo), donde sus unidades son
submúltiplos del sistema MKS.
Sistema Técnico, similar al sistema MKS y usado por los mismos países, donde la
diferencia fundamental está dada en que la fuerza es una unidad básica y la masa es una
unidad derivada, contrario a lo que ocurre en el sistema Giorgi o MKS.
Sistema Internacional (SI). Es una extensión y perfeccionamiento del sistema Giorgi o MKS
y en la actualidad es adoptado por la mayoría de los países a nivel mundial. En Cubaes el
sistema vigente según “decreto-ley 62 adoptado en 1960 y se estable con carácter
obligatorio el uso del Sistema Internacional de Unidades”.
Dado que en nuestra región los sistemas más utilizados son el Sistema Técnico y el Sistema
Internacional, el análisis de las unidades será referido a ambos sistemas.Todos los sistemas de
unidades a los que se ha hecho referencia se basan en la segunda ley de Newton.
Realizando el análisis para el sistema técnico la fuerza del mismo queda dada por la ecuación (1.1)
F=m∗a
IGiovanni Giorgi. 27/11/1871 19/08/1950. Ingeniero electricista italiano. En 1901 propuso a la Asociación Electrotécnica Italiana el sistema de
unidades que lleva su nombre. Fue adoptado en 1935 por la Comisión Electrotécnica Internacional (IEC) como sistema MKSΩ (metro, kilogramo,
segundo, ohmio) y es el origen del actual Sistema Internacional de Unidades. Fue profesor de Física y Matemática en la Universidad de Roma y desde
1935 profesor de Ingeniería Eléctrica.
137
La definición del kilogramo fuerza (kilopondio), es la fuerza ejercida sobre una masa de 1 kg masa
por la gravedad estándar en la superficie terrestre. La aceleración gravitatoria se considera de
9,80665 m/s2.
Sustituyendo valores en la ecuación (1.1)
m kg masa∗m
1kgf (kp ) = 1kgmasa ∗ 9,80665 = 9,80665
s2 s2
En honor a Sir Isaac Newton6la unidad de fuerza en el Sistema Internacional (SI) es nombrado
Newtony es representa por(N). Ecuación (1.7)
6 Isaac Newton. 4/01/1643 31/03/1727. Físico, filósofo, teólogo, inventor, alquimista y matemático inglés. Autor de los Philosophiaenaturalis
principia mathematica, más conocidos como los Principia, donde describió la ley de la gravitación universal y estableció las bases de la mecánica
clásica mediante las leyes que llevan su nombre. Entre sus otros descubrimientos científicos destacan los trabajos sobre la naturaleza de la luz y la
óptica y el desarrollo del cálculo matemático.
138
m kg masa∗m
1 Newton = 1kgmasa ∗ 1 = 1
s2 s2
Dado que ambos sistemas de unidades trabajan con una unidad denominada kilogramo, es
necesario tener muy claro cuando dicha unidad corresponde a unidades de fuerza (Sistema Técnico)
y cuando a unidades de masa (Sistema Internacional). Es frecuente cometer el error de trabajar la
fuerza en kilogramos fuerza o kilopondio y considerar la masa en kilogramos, cuando realmente le
corresponde como unidad de masa la Unidad Técnica de Masa (UTM), la cual es una masa 9,81
veces mayor que la del kilogramo.
Prefijos usados en el Sistema Internacional
Prefijo Símbolo Valor (10n)
yotta Y 1024
zetta Z 1021
exa E 1018
peta P 1015
tera T 1012
giga G 109
mega M 106
kilo k 103
hecto h 102
deca da 101
unidad ninguno 1
deci d 10-1
centi c 10-2
mili m 10-3
micro µ 10-6
nano n 10-9
pico p 10-12
femto f 10-15
atto a 10-18
zepto z 10-21
yocto y 10-24
139
Conversiones de unidades
La existencia de diferentes sistemas de unidades obliga a convertir de un sistema a otro. Es
aconsejable para su vida estudiantil y profesional, disponer de alguno de los convertidores digitales
existente en la red de la facultad de ingeniería.
Para facilitar la solución de los problemas propuestos se añaden algunas conversiones de unidades
frecuentes en el campo de la mecánica de los fluidos.
1 kgfuerza (kilopondio kp) = 9,80665N
1 UTM (unidad técnica de masa) = 9,80665 Kgmasa
1 kgm (kilográmetro) (kgfuerza*metro) = 9,80665 Joule (J) = (N*m)
1 kgfuerza/m2= 9,80665 Pa (Pascal= N/m2)
1 kgfuerza/cm2= 98,0665 kPa (kilo Pascal)
1 kgfuerza/cm2= 14,2233 libras/pulgada cuadrada.
140
Anexo 1: Factores de conversión
141
Anexo 2: Propiedades del agua
142
Anexo 4: Propiedades del aire
143
Anexo 6: Coeficiente de resistencia CD de varios cuerpos en dos dimensiones para Re>104 basado
en el área frontal A= b*D donde b es la longitud en dirección normal a la hoja. (Cengel, 2006)
144
Anexo 7: Coeficiente de resistencia CD de varios cuerpos en tres dimensiones para Re>104 basado
en el área frontal.(Cengel, 2006)
145
Anexo 8: Criterio de selección a partir de la velocidad específica ns
146