Entrega 5 Mecánica de Fluidos Q1. Curso 2016/2017 M.A.B. , A.I.S. i A.P.E. 20/12/2016
1
ENTREGA 5. MECÁNICA DE FLUIDOS
Problema 1. El esquema de la figura 1 presenta dos depósitos para un grupo de bombeado.
La diferencia de cotas entre el nivel del líquido de los dos depósitos es de 65 metros, la
longitud del conducto que une los dos depósitos es de 2500 metros, el conducto es de PVC y
tienen un diámetros de 0,12 metros.
Sabiendo que las características de la bomba se detallan con el diagrama siguiente (la gráfica
está dada en H (m) y Q (l/min), curva de 40-250-220-. Determinar:
a.- El caudal circulante. (Punto de funcionamiento de la bomba).
b.- El diámetro del conducto para que la velocidad del fluido en el interior del mismo
sea de 1 m/s. ¿Cuál es el punto de funcionamiento en este caso?
c.- Qué diferencia de cotas entre niveles del depósito sería necesaria para que
manteniendo el diámetro de conducto inicial se tenga una velocidad del fluido en el
interior del conducto de 1 m/s.
a) El cabal circulante. (Punto de funcionamiento de la bomba).
Primero de todo aplicaremos la ecuación de Bernoulli entre los puntos 1 y 2, el punto que está
en la superficie del líquido del primer depósito y el punto de la superficie del líquido del
segundo.
𝑃1𝜌𝑔
+ 𝑧1 +𝑉12
2𝑔+ 𝐻 =
𝑃2𝜌𝑔
+ 𝑧2 +𝑉22
2𝑔+ ∆𝑌12
Ambos términos de la presión quedan descartados ya que las consideramos presiones
relativas. Y también anularemos las velocidades porque la teoría de los grandes depósitos nos
dice que la velocidad de un punto de la superficie del líquido en un depósito grande que se
vacía por un agujero pequeño, es cero.
𝐻 = (𝑧2 − 𝑧1) + ∆𝑌12
𝐻 = 65 + ∆𝑌12
Representación de los depósitos
Entrega 5 Mecánica de Fluidos Q1. Curso 2016/2017 M.A.B. , A.I.S. i A.P.E. 20/12/2016
2
𝐻 = 65 + 𝑓8𝐿𝑄2
𝐷5𝑔𝜋2
Como tenemos de incógnitas los términos Q, H y f, resolveremos la ecuación mediante el
método de iteración. Cogemos un valor de f cualquiera y lo verificaremos con la gráfica de
Moody, sino volveremos a coger otro valor de este que nosotros consideremos que se
encuentra más cerca de la opción correcta. Tenemos la siguiente información:
𝐻 = −83297𝑄2 + 127,95𝑄 + 99,985
Comenzamos con un factor de fricción de valor 0,02:
𝐻 = 65 + 0,02 ·8 · 2500 · 𝑄2
0,125 · 9,81 · 𝜋2
𝐻 = 65 + 166029,4038 𝑄2
𝐻 = −83297𝑄2 + 127,95𝑄 + 99,985
Tenemos un sistema de ecuaciones con dos incógnitas, Q y H:
−83297𝑄2 + 127,95𝑄 + 99,985 = 65 + 166029,4038 𝑄2
249326,4038𝑄2 − 127,95𝑄 − 34,985 = 0
𝑄 = 0,01210496297𝑚3
𝑠 𝑖 − 0,01159178026
𝑚3
𝑠
Obviamente, descartamos el valor negativo del cabal, entonces obtendremos un valor de H de
89,32831 m
Mediante la hipótesis que tratamos con el agua como fluido, con una viscosidad de 1,02·106
m2/s, buscamos el valor de Re para comprobar el valor conseguido.
𝑅𝑒 =4𝑄
𝐷𝜋𝜐 → 𝑅𝑒 =
4 · 0,01210496297
0,12 · 𝜋 · 1,02 · 10−6= 125919,261
Y ahora calculamos la rugosidad relativa, sabiendo que la rugosidad del PVC es de 0,0015 mm:
𝜀𝑟 =0,0015
120= 1,25 · 10−5
Con el valor de Re y el de la rugosidad relativa podemos encontrar el valor del factor de
fricción f’=0,0175. Como este valor no es exactamente el que habíamos supuesto al inicio
(0,02), volvemos a aplicar el mismo método pero ahora con f=0,017:
Entrega 5 Mecánica de Fluidos Q1. Curso 2016/2017 M.A.B. , A.I.S. i A.P.E. 20/12/2016
3
𝐻 = 65 + 0,017 ·8 · 2500 · 𝑄2
0,125 · 9,81 · 𝜋2
𝐻 = 65 + 141124,9932 𝑄2
𝐻 = −83297𝑄2 + 127,95𝑄 + 99,985
Encontramos H y Q:
−83297𝑄2 + 127,95𝑄 + 99,985 = 65 + 141124,9932 𝑄2
224421,9932𝑄2 − 127,95𝑄 − 34,985 = 0
𝑄 = 0,01277388538𝑚3
𝑠 𝑖 − 0,01220375409
𝑚3
𝑠
Volviendo a descartar el valor negativo obtenemos una H=88,027668 m.
𝑅𝑒 =4𝑄
𝐷𝜋𝜐 → 𝑅𝑒 =
4 · 0,01277388538
0,12 · 𝜋 · 1,02 · 10−6= 132877,5817
Buscamos el valor de Reynolds y para la rugosidad relativa nos da un valor de fricción de
f’=0,0172. Para ser más exactos iteramos con un valor entre 0,0172 y 0,0175.
𝐻 = 65 + 0,0172 ·8 · 2500 · 𝑄2
0,125 · 9,81 · 𝜋2
𝐻 = 65 + 142785,2872 𝑄2
𝐻 = −83297𝑄2 + 127,95𝑄 + 99,985
Volvemos a encontrar H y Q:
−83297𝑄2 + 127,95𝑄 + 99,985 = 65 + 142785,2872 𝑄2
226082,2872𝑄2 − 127,95𝑄 − 34,985 = 0
𝑄 = 0,01272582618𝑚3
𝑠 𝑖 − 0,0121598818
𝑚3
𝑠
Descartando cabal negativo, tenemos H=88,12360 m. Con un número Re:
𝑅𝑒 =4𝑄
𝐷𝜋𝜐 → 𝑅𝑒 =
4 · 0,01272582618
0,12 · 𝜋 · 1,02 · 10−6= 132377,6563
Y finalmente obtenemos un factor de fricción f’=0,0174. Volvemos a iterar por última vez con
el valor de f=0,0174, que más adelante veremos que se trata del valor exacto porque nos
daremos cuenta que f=f’.
𝐻 = 65 + 0,0174 ·8 · 2500 · 𝑄2
0,125 · 9,81 · 𝜋2
Entrega 5 Mecánica de Fluidos Q1. Curso 2016/2017 M.A.B. , A.I.S. i A.P.E. 20/12/2016
4
𝐻 = 65 + 144445,5813 𝑄2
𝐻 = −83297𝑄2 + 127,95𝑄 + 99,985
Calculamos H y Q:
−83297𝑄2 + 127,95𝑄 + 99,985 = 65 + 144445,5813 𝑄2
227742,5813𝑄2 − 127,95𝑄 − 34,985 = 0
𝑄 = 0,01267830132𝑚3
𝑠 𝑖 − 0,0121164828
𝑚3
𝑠
Descartamos el valor negativo y obtenemos una H=88,21809 m. Con un número de Reynolds
que veremos que nos da un valor de f’=0,0174, tal y como hemos mencionado antes, sucede
que f=f’. Ya tenemos el valor exacto del coeficiente de fricción.
Entonces realizamos los cálculos de Q y H con este justo factor de fricción, que son:
𝑄 = 0,01267830132𝑚3
𝑠
𝐻 = 88,21808514 m
b) El diámetro del conducto para que la velocidad del fluido en el interior del mismo
sea de 1 m/s. ¿Cuál es el punto de funcionamiento en este caso?
Volvemos a aplicar Bernoulli igual que anteriormente pero en este caso conocemos la
velocidad:
𝑄 = 𝑣 · 𝑆 = 1 ·𝜋 · 𝐷2
4=𝜋 · 𝐷2
4
Sustituyendo a la ecuación de la bomba:
𝐻 = −83297(𝜋 · 𝐷2
4)
2
+ 127,95(𝜋 · 𝐷2
4) + 99,985
𝐻 = −51381,777𝐷4 + 100,491695 𝐷2 + 99,985
Otra vez Bernoulli
𝐻 = 65 + 𝑓8𝐿 (
𝜋 · 𝐷2
4 )2
𝐷5𝑔𝜋2= 65 + 𝑓
𝐿
2𝐷𝑔
Ahora haremos lo mismo que en el apartado anterior, suponemos una semilla de f, como por
ejemplos f=0,01.
𝐻 = 65 + 0,01 ·2500
19,62 · 𝐷
𝐻 = −51381,77736 𝐷4 + 100,491695 𝐷2 + 99,985
Entrega 5 Mecánica de Fluidos Q1. Curso 2016/2017 M.A.B. , A.I.S. i A.P.E. 20/12/2016
5
Para resolver este sistema de dos ecuaciones precisaremos del programa matemático Wolfram
Alpha, que no da las posibles soluciones de:
𝐻1 = 73,2585 𝑚 𝑖 𝐷1 = 0,154291 𝑚
𝐻2 = 100,028 𝑚 𝑖 𝐷2 = 0,0363769 𝑚
Se puede ver que con un diámetro pequeño no nos acercaremos a la solución, ya que nos da
un factor de fricción f’ el cual no se encuentra entre el intervalo [0.1-0.2], un intervalo de
factores de fricción posibles para conductos comerciales. Nos quedamos con la restante, que
nos dará un número Re de:
𝑅𝑒1 =𝑉 · 𝐷1𝑣
=1 · 0,154291
1,02 · 10−6= 151265,6863
𝜀𝑟1 =0,0015
154,291 = 9,7218 · 10−6
Con estos valores, en el diagrama encontraremos un valor de factor de fricción de f’=0,0167.
Aún seguiremos intentando acercarnos más a la solución y optaremos por usar una iteración
de f=0,017.
𝐻 = 65 + 0,017 ·2500
19,62 · 𝐷
𝐻 = −51381,77736 𝐷4 + 100,491695 𝐷2 + 99,985
Resolvemos el sistema con el programa:
𝐻1 = 80,0574 𝑚 𝑖 𝐷1 = 0,14386 𝑚
𝐻2 = 99,5888 𝑚 𝑖 𝐷2 = 0,062626 𝑚
Cogemos el primer caso, ya que el segundo tiene un diámetro demasiado pequeño.
Obtenemos un número de Re de:
𝑅𝑒1 =𝑉 · 𝐷1𝑣
=1 · 0,14386
1,02 · 10−6= 141039,2157
𝜀𝑟1 =0,0015
143,86 = 1,0427 · 10−5
Y finalmente conseguimos un factor de fricción f’=0,0168. Por última vez probaremos con este
mismo factor de fricción f=0.0168 y veremos que será el exacto, con el cual f=f’.
𝐻 = 65 + 0,0168 ·2500
19,62 · 𝐷
𝐻 = −51381,77736 𝐷4 + 100,491695 𝐷2 + 99,985
Resolvemos el sistema con el programa:
𝐻1 = 79,8424 𝑚 𝑖 𝐷1 = 0,144227 𝑚
Entrega 5 Mecánica de Fluidos Q1. Curso 2016/2017 M.A.B. , A.I.S. i A.P.E. 20/12/2016
6
𝐻2 = 99,618 𝑚 𝑖 𝐷2 = 0,061837 𝑚
Cogemos la primera solución ya que la segunda tiene un diámetro demasiado pequeño.
Obteniendo entonces un número de Re:
𝑅𝑒1 =𝑉 · 𝐷1𝑣
=1 · 0,144227
1,02 · 10−6= 141399,0196
𝜀𝑟1 =0,0015
144,227= 1,0400 · 10−5
Y ahora sí, con estos valores encontramos que f=f’=0.0168, el valor exacto del factor de
fricción.
Ya podemos realizar el cálculo del diámetro con nuestro factor de fricción exacto.
𝐷 = 0,144227 𝑚
El cabal consecuente a este diámetro tiene un valor de:
𝑄 = 𝑣 · 𝑆 = 1 ·𝜋 · 𝐷2
4=𝜋 · 𝐷2
4=𝜋 · 0,1442272
4= 0,01634
𝑚3
𝑠
Pudiendo calcular el valor final de H, que es:
𝐻 = 79,8424 𝑚
c) ¿Qué diferencia de cotas entre niveles del depósito sería necesaria para que
manteniendo el diámetro del conducto inicial se tenga una velocidad del fluido en
el interior del conducto de 1 m/s?
En este último apartado debemos calcular la diferencia de z ya que nos dan la velocidad del
fluido y su diámetro. Siempre aplicando Beroulli:
𝑄 = 𝑣 · 𝑆 = 1 ·𝜋 · 𝐷2
4=𝜋 · 𝐷2
4=𝜋 · 0,122
4= 0,01131
𝑚3
𝑠
𝐻 = ∆𝑧 + 𝑓8 · 𝐿 · 𝑄2
𝐷5 · 𝑔 · 𝜋2 → 𝐻 = 𝑧 + 𝑓
8 · 2500 · 0,011312
0,125 · 9,81 · 𝜋2
𝐻 = 𝑧 + 𝑓1061,891691
Por otro lado podemos calcular la H con la ecuación de la bomba, sustuyendo el valor del cabal
que hemos calculado al inicio del apartado
𝐻 = −83297𝑄2 + 127,95𝑄 + 99,985
𝐻 = −83297(0,01131)2 + 127,95(0,01131) + 99,985
𝐻 = 99,7771 𝑚
Entonces ya podemos encontrar el valor del número de Reynolds y su rugosidad relativa
Entrega 5 Mecánica de Fluidos Q1. Curso 2016/2017 M.A.B. , A.I.S. i A.P.E. 20/12/2016
7
𝑅𝑒 =𝑉 · 𝐷
𝑣=
1 · 0,12
1,02 · 10−6= 117647,0588
𝜀𝑟 =0,0015
120= 1,25 · 10−5
En este caso nos resultará un valor de f=0.018, valor del diagrama de Moody para conductos
comerciales. Y podremos encontrar la diferencia de las cotas z, de los dos depósitos.
Sustituyendo el valor de f en la ecuación de la H
∆𝑧 = 𝐻 − 𝑓1061,891691
∆𝑧 = 99,7771 − 0,018 · 1061,891691
∆𝑧 = 80,66 𝑚
Entrega 5 Mecánica de Fluidos Q1. Curso 2016/2017 M.A.B. , A.I.S. i A.P.E. 20/12/2016
8
Problema 2. Sea el conjunto depósitos y conductos que se expone en los dos esquemas
siguientes. Si se conoce: Que la curva característica de la bomba está dada por la ecuación
𝒀𝒃𝒐𝒎𝒃𝒂 = 𝒀𝟎 − 𝒂𝑸𝟐, donde 𝒀𝟎 y 𝒂, son constantes conocidas, conociendo además las
longitudes, diámetros y rugosidades absolutas de todos los tramos, se pide determinar el
caudal que circula por las dos instalaciones y por cada uno de los tramos. Realizar el cálculo
mediante la determinación de las constantes equivalentes de los conductos. Supóngase que
los tramos situados a la entrada y salida de la bomba son muy cortos y se puede despreciar
su efecto.
La ecuación que deberá utilizarse para determinar las pérdidas de carga en cada tramo es la
de Darcy Weisbach.
∆𝒀𝒊 = 𝒇𝒊𝑳𝒊
𝑫𝒊𝟓
𝟖𝑸𝒊𝟐
𝝅𝟐
Tomar como primera aproximación, el factor de fricción “f” función de la rugosidad relativa
(ɛ
𝑫).
Resolución del primer caso:
Para empezar, debemos aplicar la ecuación de la energía entre las superficies libres de los dos
depósitos:
𝑃1𝜌+𝑣12
2+ 𝑔𝑧1 + 𝑌𝐵 =
𝑃2𝜌+𝑣22
2+ 𝑔𝑧2 + ∆𝑌12
Como trabajamos con presiones relativas i consideramos que la superficie de los dos depósitos
es muy grande, obtenemos las expresiones siguientes:
Entrega 5 Mecánica de Fluidos Q1. Curso 2016/2017 M.A.B. , A.I.S. i A.P.E. 20/12/2016
9
𝑊
𝜌𝑄= 𝑌𝐵 = 𝑔(𝑧2 − 𝑧1) + ∆𝑌12
∆𝑌12 = ∆𝑌𝑝𝑎𝑟𝑎𝑙·𝑙𝑒𝑙 + ∆𝑌4 → ∆𝑌𝑝𝑎𝑟𝑎𝑙·𝑙𝑒𝑙 = ∆𝑌1 + ∆𝑌2 + ∆𝑌3
∆𝑌1 = 𝐾1𝑄12
∆𝑌2 = 𝐾2𝑄22
∆𝑌3 = 𝐾3𝑄32
Si estas tres últimas las comparamos con la ecuación de Darcy Weisbach, vemos que la 𝐾
representa todos los términos a excepción del cabal.
Ésta constante depende del factor de fricción, que es función de la rugosidad relativa y del
número de Reynolds. Estas constantes se suponen conocidas, ya que supondremos como
primera aproximación que el régimen en todo momento sea turbulento desarrollado.
De la ecuación de continuidad, tenemos lo siguiente:
𝑄4 = 𝑄1 + 𝑄2 + 𝑄3 = √∆𝑌𝑝𝑎𝑟𝑎𝑙·𝑙𝑒𝑙
𝐾1+√
∆𝑌𝑝𝑎𝑟𝑎𝑙·𝑙𝑒𝑙𝐾2
+√∆𝑌𝑝𝑎𝑟𝑎𝑙·𝑙𝑒𝑙
𝐾3
Si definimos la constante equivalente de las pérdidas del conjunto de las tres brancas en
paralelo como 𝐾𝑒𝑓, tenemos la expresión siguiente:
∆𝑌𝑝𝑎𝑟𝑎𝑙·𝑙𝑒𝑙 = 𝐾𝑒𝑓𝑄42 = 𝐾𝑒𝑓 (√
∆𝑌𝑝𝑎𝑟𝑎𝑙·𝑙𝑒𝑙
𝐾1+√
∆𝑌𝑝𝑎𝑟𝑎𝑙·𝑙𝑒𝑙
𝐾2+√
∆𝑌𝑝𝑎𝑟𝑎𝑙·𝑙𝑒𝑙
𝐾3)
2
∆𝑌𝑝𝑎𝑟𝑎𝑙·𝑙𝑒𝑙 = 𝐾𝑒𝑓∆𝑌𝑝𝑎𝑟𝑎𝑙·𝑙𝑒𝑙 (1
√𝐾1+
1
√𝐾2+
1
√𝐾3)
2
→ 𝐾𝑒𝑓 =1
(1
√𝐾1+
1
√𝐾2+
1
√𝐾3)
2
Nos falta definir la constante para la tobera 4:
∆𝑌4 = 𝐾4𝑄42
Así pues, la ecuación de la energía planteada inicialmente quedará expresada de la siguiente
manera:
[𝑚2
𝑠2]
𝑊
𝜌𝑄4= 𝑔(𝑧2 − 𝑧1) + 𝐾𝑒𝑓𝑄4
2 + 𝐾4𝑄42
Donde, multiplicando por 𝑄4, podemos encontrar la expresión siguiente:
𝑄4 =
𝑊𝜌
(𝑔(𝑧2 − 𝑧1) + 𝐾𝑒𝑓𝑄42 + 𝐾4𝑄4
2)
Entrega 5 Mecánica de Fluidos Q1. Curso 2016/2017 M.A.B. , A.I.S. i A.P.E. 20/12/2016
10
Finalmente, las pérdidas de carga en el tramo 4 y en las tres brancas en paralelo serán:
∆𝑌4 = 𝐾4𝑄42
∆𝑌𝑝𝑎𝑟𝑎𝑙·𝑙𝑒𝑙 = 𝐾𝑒𝑓𝑄42
Por lo tanto, el cabal en cada uno de los tramos se obtendrá a partir de la siguiente expresión:
∆𝑌𝑝𝑎𝑟𝑎𝑙·𝑙𝑒𝑙 = 𝐾𝑒𝑓𝑄42 = 𝐾1𝑄1
2 = 𝐾2𝑄22 = 𝐾3𝑄3
2
Donde tenemos:
𝑄1 = √𝐾𝑒𝑓𝑄4
2
𝐾1 ; 𝑄2 = √
𝐾𝑒𝑓𝑄42
𝐾2 ; 𝑄3 = √
𝐾𝑒𝑓𝑄42
𝐾3
Resolución del segundo caso:
Al igual que con el primer apartado, aplicando la ecuación de la energía entre las superficies
libres de los dos depósitos, obtenemos que:
𝑊
𝜌𝑄𝑡𝑜𝑡𝑎𝑙= 𝑌𝐵 = 𝑔(𝑧𝐵 − 𝑧𝐴) + ∆𝑌𝐴𝐵
La diferencia es que, en este caso, las pérdidas de energía se definirán como:
∆𝑌𝐴𝐵 = ∆𝑌𝑝𝑎𝑟𝑎𝑙·𝑙𝑒𝑙 1−3 + ∆𝑌𝐷
Donde tenemos que:
∆𝑌𝑝𝑎𝑟𝑎𝑙·𝑙𝑒𝑙 1−3 = ∆𝑌𝑝𝑎𝑟𝑎𝑙·𝑙𝑒𝑙 1−2 + ∆𝑌2−3 = 𝐾𝑒𝑓(1−3)𝑄𝑡𝑜𝑡𝑎𝑙2
𝑄𝑡𝑜𝑡𝑎𝑙 = 𝑄(1−3)𝐴 + 𝑄23 = 𝑄(1−3)𝐴 +𝑄(1−2)𝐵 +𝑄(1−2)𝐶
Por otro lado tenemos las siguientes expresiones:
∆𝑌12 = ∆𝑌(1−2)𝐵 = ∆𝑌(1−2)𝐶
∆𝑌(1−2)𝐵 = 𝐾(1−2)𝐵𝑄(1−2)𝐵2
∆𝑌(1−2)𝐶 = 𝐾(1−2)𝐶𝑄(1−2)𝐶2
Si lo sustituimos en la ecuación entre los puntos 1 y 2, obtenemos que:
𝑄23 = √∆𝑌(1−2)𝐵𝐾(1−2)𝐵
+√∆𝑌(1−2)𝐶𝐾(1−2)𝐶
∆𝑌12 = 𝐾(1−2)𝑒𝑞𝑢𝑖𝑣𝑎𝑙𝑒𝑛𝑡𝑄2−32 = 𝐾(1−2)𝐵𝑄(1−2)𝐵
2 = 𝐾(1−2)𝐶𝑄(1−2)𝐶2
∆𝑌12 = 𝐾(1−2)𝑒𝑞𝑢𝑖𝑣𝑎𝑙𝑒𝑛𝑡𝑄2−32 = 𝐾(1−2)𝑒𝑞𝑢𝑖𝑣𝑎𝑙𝑒𝑛𝑡∆𝑌(1−2) (
1
√𝐾(1−2)𝐵+
1
√𝐾(1−2)𝐶)
2
Entrega 5 Mecánica de Fluidos Q1. Curso 2016/2017 M.A.B. , A.I.S. i A.P.E. 20/12/2016
11
𝐾(1−2)𝑒𝑞𝑢𝑖𝑣𝑎𝑙𝑒𝑛𝑡 =1
(1
√𝐾(1−2)𝐵+
1
√𝐾(1−2)𝐶)
2
Por otro lado:
∆𝑌(1−3)𝑝𝑎𝑟𝑎𝑙·𝑙𝑒𝑙 = 𝐾(1−3)𝑒𝑞𝑢𝑖𝑣𝑎𝑙𝑒𝑛𝑡𝑄𝑡𝑜𝑡𝑎𝑙2 = ∆𝑌(1−3)𝐴 = ∆𝑌(1−2)𝑝𝑎𝑟𝑎𝑙·𝑙𝑒𝑙 + ∆𝑌(2−3)
𝐾(1−3)𝑒𝑞𝑢𝑖𝑣𝑎𝑙𝑒𝑛𝑡𝑄𝑡𝑜𝑡𝑎𝑙2 = 𝐾(1−3)𝐴𝑄(1−3)𝐴
2 = 𝐾(1−2)𝑒𝑞𝑢𝑖𝑣𝑎𝑙𝑒𝑛𝑡𝑄2−32 + 𝐾2−3𝑄2−3
2
Así pues:
𝐾(1−3)𝑒𝑞𝑢𝑖𝑣𝑎𝑙𝑒𝑛𝑡𝑄𝑡𝑜𝑡𝑎𝑙2 = 𝐾(1−3)𝐴𝑄(1−3)𝐴
2 = 𝑄2−32 (𝐾(1−2)𝑒𝑞𝑢𝑖𝑣𝑎𝑙𝑒𝑛𝑡 + 𝐾2−3)
= ∆𝑌(1−3)𝑡𝑟𝑎𝑚 𝑠𝑢𝑝𝑒𝑟𝑖𝑜𝑟 = ∆𝑌(1−3)𝑡𝑟𝑎𝑚 𝑖𝑛𝑓𝑒𝑟𝑖𝑜𝑟
𝐾(1−2)𝑒𝑞𝑢𝑖𝑣𝑎𝑙𝑒𝑛𝑡 +𝐾2−3 =1
(1
√𝐾(1−2)𝐵+
1
√𝐾(1−2)𝐶)
2 + 𝐾2−3 = 𝐾(1−3)𝑡𝑟𝑎𝑚 𝑖𝑛𝑓𝑒𝑟𝑖𝑜𝑟
Si recordamos que:
𝐾(1−3)𝑒𝑞𝑢𝑖𝑣𝑎𝑙𝑒𝑛𝑡𝑄𝑡𝑜𝑡𝑎𝑙2 = 𝐾(1−3)𝑡𝑟𝑎𝑚 𝑖𝑛𝑓𝑒𝑟𝑖𝑜𝑟𝑄2−3
2 = 𝐾(1−3)𝐴𝑄(1−3)𝐴2
𝑄𝑡𝑜𝑡𝑎𝑙 = 𝑄2−3 + 𝑄(1−3)𝐴 = √∆𝑌(1−3)𝑡𝑟𝑎𝑚 𝑖𝑛𝑓𝑒𝑟𝑖𝑜𝑟
𝐾(1−2)𝑒𝑞𝑢𝑖𝑣𝑎𝑙𝑒𝑛𝑡 + 𝐾2−3+√
∆𝑌(1−3)𝐴𝐾(1−3)𝐴
𝐾(1−3)𝑒𝑞𝑢𝑖𝑣𝑎𝑙𝑒𝑛𝑡𝑄𝑡𝑜𝑡𝑎𝑙2 = ∆𝑌(1−3)𝑒𝑞𝑢𝑖𝑣𝑎𝑙𝑒𝑛𝑡
= 𝐾(1−3)𝑒𝑞𝑢𝑖𝑣𝑎𝑙𝑒𝑛𝑡 [√∆𝑌(1−3)𝑡𝑟𝑎𝑚 𝑖𝑛𝑓𝑒𝑟𝑖𝑜𝑟
𝐾(1−2)𝑒𝑞𝑢𝑖𝑣𝑎𝑙𝑒𝑛𝑡 + 𝐾2−3+√
∆𝑌(1−3)𝐴𝐾(1−3)𝐴
]
2
Donde tenemos que:
∆𝑌(1−3)𝑒𝑞𝑢𝑖𝑣𝑎𝑙𝑒𝑛𝑡
= 𝐾(1−3)𝑒𝑞𝑢𝑖𝑣𝑎𝑙𝑒𝑛𝑡∆𝑌(1−3)𝑒𝑞𝑢𝑖𝑣𝑎𝑙𝑒𝑛𝑡 [1
√𝐾(1−2)𝑒𝑞𝑢𝑖𝑣𝑎𝑙𝑒𝑛𝑡 + 𝐾2−3+
1
√𝐾(1−3)𝐴]
2
𝐾(1−3)𝑒𝑞𝑢𝑖𝑣𝑎𝑙𝑒𝑛𝑡 =1
[1
√𝐾(1−2)𝑒𝑞𝑢𝑖𝑣𝑎𝑙𝑒𝑛𝑡 + 𝐾2−3+
1
√𝐾(1−3)𝐴]
2
∆𝑌𝐴𝐵 = 𝐾(1−3)𝑒𝑞𝑢𝑖𝑣𝑎𝑙𝑒𝑛𝑡𝑄𝑡𝑜𝑡𝑎𝑙2 + 𝐾𝐷𝑄𝑡𝑜𝑡𝑎𝑙
2
Por lo tanto, la ecuación de la energía que habíamos planteado inicialmente nos queda
expresada de la siguiente manera:
𝑊
𝜌𝑄𝑡𝑜𝑡𝑎𝑙= 𝑔(𝑧𝐵 − 𝑧𝐴) + 𝑄𝑡𝑜𝑡𝑎𝑙
2 [𝐾(1−3)𝑒𝑞𝑢𝑖𝑣𝑎𝑙𝑒𝑛𝑡 + 𝐾𝐷]
Entrega 5 Mecánica de Fluidos Q1. Curso 2016/2017 M.A.B. , A.I.S. i A.P.E. 20/12/2016
12
Si multiplicamos la ecuación por 𝑄𝑡𝑜𝑡𝑎𝑙, podremos aislar su valor:
𝑄𝑡𝑜𝑡𝑎𝑙 =
𝑊𝜌
𝑔(𝑧𝐵 − 𝑧𝐴) + 𝑄𝑡𝑜𝑡𝑎𝑙2 [𝐾(1−3)𝑒𝑞𝑢𝑖𝑣𝑎𝑙𝑒𝑛𝑡 +𝐾𝐷]
Seguidamente, tenemos que la pérdida de la carga en cada uno de los tramos 1-3 será de:
𝑄𝑡𝑜𝑡𝑎𝑙2 = ∆𝑌1−3 [
1
√𝐾(1−2)𝑒𝑞𝑢𝑖𝑣𝑎𝑙𝑒𝑛𝑡 + 𝐾2−3+
1
√𝐾(1−3)𝐴]
2
Y por lo tanto:
∆𝑌1−3 =𝑄𝑡𝑜𝑡𝑎𝑙2
[1
√𝐾(1−2)𝑒𝑞𝑢𝑖𝑣𝑎𝑙𝑒𝑛𝑡 + 𝐾2−3+
1
√𝐾(1−3)𝐴]
2
Así pues, los cabales en cada uno de los tramos 1-3 serán:
∆𝑌1−3 = 𝐾(1−3)𝐴𝑄(1−3)𝐴2 → 𝑄(1−3)𝐴 = √
∆𝑌1−3𝐾(1−3)𝐴
∆𝑌1−3 = (𝐾(1−2)𝑒𝑞𝑢𝑖𝑣𝑎𝑙𝑒𝑛𝑡 + 𝐾2−3)𝑄2−32 → 𝑄2−3 = √
∆𝑌1−3
(𝐾(1−2)𝑒𝑞𝑢𝑖𝑣𝑎𝑙𝑒𝑛𝑡 + 𝐾2−3)
Finalmente, nos damos cuenta de que tenemos que tener en cuenta las expresiones siguientes
para encontrar el cabal entre los puntos 1-2 de la tobera B y el cabal que circula entre estos
dos puntos en la tobera C:
𝑄2−3 = 𝑄(1−2)𝐵 + 𝑄(1−2)𝐶
∆𝑌1−3 = ∆𝑌2−3 + ∆𝑌(1−2)𝑒𝑞𝑢𝑖𝑣𝑎𝑙𝑒𝑛𝑡
Por lo tanto, llegamos a la conclusión de que las expresiones finales que nos permiten
encontrar todas las incógnitas que teníamos inicialmente son:
∆𝑌2−3 = 𝐾2−3𝑄2−3
∆𝑌(1−2)𝑒𝑞𝑢𝑖𝑣𝑎𝑙𝑒𝑛𝑡 = ∆𝑌1−3 − ∆𝑌2−3
∆𝑌(1−2)𝑒𝑞𝑢𝑖𝑣𝑎𝑙𝑒𝑛𝑡 = ∆𝑌(1−2)𝐵 = ∆𝑌(1−2)𝐶
∆𝑌(1−2)𝐵 = 𝐾(1−2)𝐵𝑄(1−2)𝐵2 → 𝑄(1−2)𝐵 = √
∆𝑌(1−2)𝐵𝐾(1−2)𝐵
∆𝑌(1−2)𝐶 = 𝐾(1−2)𝐶𝑄(1−2)𝐶2 → 𝑄(1−2)𝐶 = √
∆𝑌(1−2)𝐶
𝐾(1−2)𝐶
Entrega 5 Mecánica de Fluidos Q1. Curso 2016/2017 M.A.B. , A.I.S. i A.P.E. 20/12/2016
13
Problema 3. El esquema definido en la figura muestra una instalación de bombeo de agua.
Las cotas de los puntos 1, 2, 3, 4, 5 y 6 son respectivamente: 𝒁𝟏 = 𝟎𝒎, 𝒁𝟐 = 𝟐𝟎𝒎, 𝒁𝟑 =
𝟓𝟎𝒎, 𝒁𝟒 = 𝟐𝟎𝒎, 𝒁𝟓 = 𝟒𝟎𝒎 i 𝒁𝟔 = 𝟖𝟎𝒎.
Se conoce que las pérdidas de carga en cada uno de los tramos están dadas por:
∆𝒉𝟏→𝟐 = 𝟓𝟎𝟎𝟎 · 𝑸𝟏−𝟐𝟐
∆𝒉𝟐→𝟑 = 𝟖𝟎𝟎𝟎 · 𝑸𝟐−𝟑𝟐
∆𝒉𝟒→𝟓 = 𝟏𝟎𝟎𝟎𝟎 · 𝑸𝟒−𝟓𝟐
∆𝒉𝟐→𝟒 = 𝟑𝟎𝟎𝟎 · 𝑸𝟐−𝟒𝟐
∆𝒉𝟒→𝟔 = 𝟔𝟎𝟎𝟎 · 𝑸𝟒−𝟔𝟐
Recuérdese que la pérdida de carga genérica está definida por:
∆𝒉𝒊 = 𝐟𝐢𝑳𝒊
𝑫𝒊𝟓
𝟖𝑸𝒊𝟐
𝝅𝟐 · 𝒈
Se conoce además que la curva característica de cada una de las bombas esta definida por:
𝑯𝑩𝟏 = 𝟏𝟎𝟎 − 𝟑𝟎𝟎𝟎 · 𝑸𝟏−𝟐𝟐
𝑯𝑩𝟐 = 𝟒𝟎 − 𝟔𝟎𝟎𝟎 · 𝑸𝟒−𝟔𝟐
Determinar:
1. El caudal circulante en cada tramo, así como las presiones en cada uno de los nudos. En el
esquema siguiente se indica el sentido de circulación del fluido.
2. Si se conoce que la longitud de cada uno de los 5 tramos de conducto es de 500 m. y se
estima que f=0,02. Hallar el diámetro de cada uno de los conductos.
3. Comentar como se podría mejorar la instalación partiendo de los datos obtenidos en el
caso. Determinar los nuevos caudales con la mejora establecida. (Realizar las hipótesis que
se crean oportunas).
Esquemas de los conjuntos de depósitos y conductos a estudiar
Entrega 5 Mecánica de Fluidos Q1. Curso 2016/2017 M.A.B. , A.I.S. i A.P.E. 20/12/2016
14
Apartado 1:
Para encontrar las ecuaciones que caracterizan el sistema, aplicaremos la ecuación de
Bernoulli por un fluido incompresible, entre los extremos de los diferentes tramos que
conforman el sistema. Así pues, aplicando Bernoulli entre los puntos 1 y 2, tenemos que:
𝑃1𝜌 · 𝑔
+𝑉12
2 · 𝑔+ 𝑧1 +𝐻𝑏1→2 =
𝑃2𝜌 · 𝑔
+𝑉1→22
2 · 𝑔+ 𝑧2 + ∆ℎ1→2
Trabajamos con presiones relativas, y por lo tanto sabemos que:
𝑃1 = 𝑃3 = 𝑃5 = 𝑃6 = 0
Como no se conocen las dimensiones del depósito 1, se establece que la velocidad en este
depósito es menospreciable (𝑉1 ≅ 0). De manera que:
𝐻𝑏1→2 =𝑃2𝜌 · 𝑔
+𝑉1→22
2 · 𝑔+ 𝑧2 + ∆ℎ1→2
Sustituyendo las dadas del enunciado:
𝐻𝑏1→2 =𝑃2𝜌 · 𝑔
+𝑉1→22
2 · 𝑔+ 𝑧2 + ∆ℎ1→2
100 − 3000 · 𝑄1−22 =
𝑃2𝜌 · 𝑔
+𝑉1→22
2 · 𝑔+ 20 + 5000 · 𝑄1−2
2 [1]
Ahora, aplicando Bernoulli entre los puntos 2 y 3:
𝑃2𝜌 · 𝑔
+𝑉2→32
2 · 𝑔+ 𝑧2 =
𝑃3𝜌 · 𝑔
+𝑉32
2 · 𝑔+ 𝑧3 + ∆ℎ2→3
Como consideramos que la sección del depósito 3 es muy grande comparada con la sección de
los conductos, entonces se puede considerar que la velocidad en este depósito es de nuevo
menospreciable (𝑉3 ≅ 0). De forma que:
𝑃2𝜌 · 𝑔
+𝑉2→32
2 · 𝑔+ 20 = 50 + 8000 · 𝑄2−3
2 [2]
Con la ecuación de Bernoulli nuevamente entre los puntos 2 y 4:
𝑃2𝜌 · 𝑔
+𝑉2→42
2 · 𝑔+ 𝑧2 =
𝑃4𝜌 · 𝑔
+𝑉42
2 · 𝑔+ 𝑧4 + ∆ℎ2→4
Tal y como se observa en la figura del enunciado, los puntos 2 y 4 se encuentran a la misma
altitud (Z), y por lo tanto podemos afirmar que 𝑉2 = 𝑉4. Y entonces:
𝑃2𝜌 · 𝑔
=𝑃4𝜌 · 𝑔
+ 3000 · 𝑄2−42 [3]
Entrega 5 Mecánica de Fluidos Q1. Curso 2016/2017 M.A.B. , A.I.S. i A.P.E. 20/12/2016
15
Aplicando otra vez Bernoulli pero entre los puntos 4 y 5:
𝑃4𝜌 · 𝑔
+𝑉4→52
2 · 𝑔+ 𝑧4 =
𝑃5𝜌 · 𝑔
+𝑉52
2 · 𝑔+ 𝑧5 + ∆ℎ4→5
Siguiendo el mismo razonamiento que en el depósito 3, afirmamos que el depósito 5 tiene la
sección muy grande comparada con la sección de los conductos, y por eso consideramos que la
velocidad en el depósito es menospreciable (𝑉5 ≅ 0). De forma que:
𝑃4𝜌 · 𝑔
+𝑉4→52
2 · 𝑔+ 20 = 40 + 10000 · 𝑄4−5
2 [4]
Aplicando Bernoulli entre los puntos 4 y 6
𝑃4𝜌 · 𝑔
+𝑉4→62
2 · 𝑔+ 𝑧4 +𝐻𝑏4→6 =
𝑃6𝜌 · 𝑔
+𝑉62
2 · 𝑔+ 𝑧6 + ∆ℎ4→6
Siguiendo el mismo razonamiento que en el depósito 3 y 5, podemos igualmente afirmar que
el depósito 6 tiene la sección muy grande comparada con la sección de los conductos, y por lo
tanto se puede considerar que la velocidad del depósito es menospreciable (𝑉6 ≅ 0). Y de
manera que como en los casos anteriores tenemos:
𝑃4𝜌 · 𝑔
+𝑉4→62
2 · 𝑔+ 20 + 40 − 6000 · 𝑄4−6
2 = 80 + 6000 · 𝑄4−62 [5]
Como de entrada no se conocen los diámetros de los conductos, los términos
𝑉1→22
2·𝑔,𝑉2→32
2·𝑔,𝑉2→42
2·𝑔,𝑉4→52
2·𝑔,𝑉4→62
2·𝑔 son muy pequeños en comparación con los otros términos, y se
pueden despreciar. Teniendo en cuenta que las incógnitas de las ecuaciones anteriores son los
cabales de cada tramo i las presiones en las intersecciones 2 y 4, debemos aplicar la ecuación
de continuidad en estos puntos y así obtener un sistema del mismo número tanto en
ecuaciones como en incógnitas.
Aplicamos esta ecuación de continuidad en la intersección 2, y obtenemos la siguiente
expresión:
𝑄1→2 = 𝑄2→3 + 𝑄2→4 [6]
En la intersección 4, tenemos:
𝑄2→4 = 𝑄4→5 + 𝑄4→6 [7]
Y resolviendo el sistema de ecuaciones correspondiente a todas las ecuaciones encontradas a
lo largo del desarrollo, [1] a [7], determinadas en el siguiente paréntesis:
Entrega 5 Mecánica de Fluidos Q1. Curso 2016/2017 M.A.B. , A.I.S. i A.P.E. 20/12/2016
16
{
100 − 3000 · 𝑄1−2
2 =𝑃2𝜌 · 𝑔
+𝑉1→22
2 · 𝑔+ 20 + 5000 · 𝑄1−2
2
𝑃2𝜌 · 𝑔
+𝑉2→32
2 · 𝑔+ 20 = 50 + 8000 · 𝑄2−3
2
𝑃2𝜌 · 𝑔
=𝑃4𝜌 · 𝑔
+ 3000 · 𝑄2−42
𝑃4𝜌 · 𝑔
+𝑉4→52
2 · 𝑔+ 20 = 40 + 10000 · 𝑄4−5
2
𝑃4𝜌 · 𝑔
+𝑉4→62
2 · 𝑔+ 20 + 40 − 6000 · 𝑄4−6
2 = 80 + 6000 · 𝑄4−62
𝑄1→2 = 𝑄2→3 + 𝑄2→4𝑄2→4 = 𝑄4→5 + 𝑄4→6 }
Obtenemos los siguientes valores de Q i P:
{
𝑄12 = 0.07519 𝑚
3/𝑠
𝑄23 = 0.02443 𝑚3/𝑠
𝑄24 = 0.05076 𝑚3/𝑠
𝑄45 = 0.02654 𝑚3/𝑠
𝑄46 = 0.02452 𝑚3/𝑠
𝑃2 = 3.4112 · 105 𝑃𝑎
𝑃4 = 2.6529 · 105 𝑃𝑎}
Apartado 2:
Como en el primer apartado del problema hemos determinado los cabales que circulan en
cada tramo, ya podemos calcular directamente el diámetro de cada conducto, acordándonos
de que las pérdidas de energía debidas a la fricción entre el conducto y el fluido vienen
determinadas por la ecuación de Darcy-Weisbach:
∆ℎ𝑖→𝑗 = 𝐾𝑖→𝑗 · 𝑄𝑖→𝑗2 = fi ·
𝐿𝑖→𝑗
𝐷𝑖→𝑗5 ·
8𝑄𝑖→𝑗2
𝜋2 · 𝑔
𝐷𝑖→𝑗 = √fi ·𝐿𝑖→𝑗
𝐾𝑖→𝑗·
8
𝜋2 · 𝑔
5
En el enunciado del problema se especifican los valores de 𝐾𝑖→𝑗; la longitud de todos los
tramos que es 𝐿 = 500 𝑚; y f = 0.02 en todos los tramos.
- Tramo 1-2:
𝐷𝑖→𝑗 = √0.02 ·500
5000·
8
𝜋2 · 9.81
5
= 0.17524 𝑚
- Tramo 2-3:
Entrega 5 Mecánica de Fluidos Q1. Curso 2016/2017 M.A.B. , A.I.S. i A.P.E. 20/12/2016
17
𝐷𝑖→𝑗 = √0.02 ·500
8000·
8
𝜋2 · 9.81
5
= 0.15752 𝑚
- Tramo 2-4:
𝐷𝑖→𝑗 = √0.02 ·500
3000·
8
𝜋2 · 9.81
5
= 0.19409 𝑚
- Tramo 4-5:
𝐷𝑖→𝑗 = √0.02 ·500
10000·
8
𝜋2 · 9.81
5
= 0.15255 𝑚
- Tramo 4-6:
𝐷𝑖→𝑗 = √0.02 ·500
6000·
8
𝜋2 · 9.81
5
= 0.16896 𝑚
También queremos remarcar que los diámetros obtenidos no son posibles en realidad, ya que
no existen conductos con estas dimensiones.
Apartado 3:
Para poder mejorar el sistema de distribución del agua, calculamos inicialmente cuál es la
velocidad del fluido en cada tramo. De los apartados anteriores, conocemos el cabal y el
diámetro correspondiente en cada tramo, y por esta razón podemos calcular la velocidad
siguiendo lo siguiente:
𝑄𝑖→𝑗 = 𝑉𝑖→𝑗 · 𝑆 = 𝑉𝑖→𝑗 · 𝜋 · (𝐷𝑖→𝑗2
4) ; 𝑉𝑖→𝑗 =
4 · 𝑄𝑖→𝑗
𝜋 · 𝐷𝑖→𝑗2
- Tramo 1-2:
𝑉𝑖→𝑗 =4 · 0.075189
𝜋 · 0.175242= 3.11743 𝑚/𝑠
- Tramo 2-3:
𝑉𝑖→𝑗 =4 · 0.024426
𝜋 · 0.157522 = 1.25340 𝑚/𝑠
- Tramo 2-4:
𝑉𝑖→𝑗 =4 · 0.050763
𝜋 · 0.194092= 1.71574 𝑚/𝑠
- Tramo 4-5:
Entrega 5 Mecánica de Fluidos Q1. Curso 2016/2017 M.A.B. , A.I.S. i A.P.E. 20/12/2016
18
𝑉𝑖→𝑗 =4 · 0.026538
𝜋 · 0.15255 2= 1.45193 𝑚/𝑠
- Tramo 4-6:
𝑉𝑖→𝑗 =4 · 0.0245225
𝜋 · 0.168962= 1.08047 𝑚/𝑠
Observamos entonces que hay dos posibles mejoras:
1. Velocidad del fluido cercana a 1m/s en todos los tramos. 2. Diámetros de los conductos que existen en realidad.
Para conseguirlas, debemos determinar cuál es el diámetro necesario en los conductos para
que la velocidad del fluido sea 1m/s.
- Aplicando Bernoulli entre los puntos 1 y 2:
100 − 3000 · 𝑄1−22 =
𝑃2𝜌 · 𝑔
+𝑉1→22
2 · 𝑔+ 20 + 5000 · 𝑄1−2
2
100 − 3000 · (𝑉12 · 𝜋 ·𝐷122
4)
2
=𝑃2𝜌 · 𝑔
+1
2 · 𝑔+ 20 + 5000 · (𝑉12 · 𝜋 ·
𝐷122
4)
2
100 − 3000 · (𝜋 ·𝐷122
4)
2
=𝑃2𝜌 · 𝑔
+1
2 · 𝑔+ 20 + 5000 · (𝜋 ·
𝐷122
4)
2
- Aplicando Bernoulli entre los puntos 2 y 3:
𝑃2𝜌 · 𝑔
+𝑉2→32
2 · 𝑔+ 20 = 50 + 8000 · 𝑄2−3
2
𝑃2𝜌 · 𝑔
+1
2 · 𝑔+ 20 = 50 + 8000 · (𝑉23 · 𝜋 ·
𝐷232
4)
2
𝑃2𝜌 · 𝑔
+1
2 · 𝑔+ 20 = 50 + 8000 · (𝜋 ·
𝐷232
4)
2
- Aplicando Bernoulli entre los puntos 2 y 4:
𝑃2𝜌 · 𝑔
=𝑃4𝜌 · 𝑔
+ 3000 · 𝑄2−42
𝑃2𝜌 · 𝑔
=𝑃4𝜌 · 𝑔
+ 3000 · (𝑉24 · 𝜋 ·𝐷242
4)
2
𝑃2𝜌 · 𝑔
=𝑃4𝜌 · 𝑔
+ 3000 · (𝜋 ·𝐷242
4)
2
Entrega 5 Mecánica de Fluidos Q1. Curso 2016/2017 M.A.B. , A.I.S. i A.P.E. 20/12/2016
19
- Aplicando Bernoulli entre los puntos 4 y 5:
𝑃4𝜌 · 𝑔
+𝑉4→52
2 · 𝑔+ 20 = 40 + 10000 · 𝑄4−5
2
𝑃4𝜌 · 𝑔
+1
2 · 𝑔+ 20 = 40 + 10000 · (𝑉45 · 𝜋 ·
𝐷452
4)
2
𝑃4𝜌 · 𝑔
+1
2 · 𝑔+ 20 = 40 + 10000 · (𝜋 ·
𝐷452
4)
2
- Aplicando Bernoulli entre los puntos 4 y 6:
𝑃4𝜌 · 𝑔
+𝑉4→62
2 · 𝑔+ 20 + 40 − 6000 · 𝑄4−6
2 = 80 + 6000 · 𝑄4−62
𝑃4𝜌 · 𝑔
+1
2 · 𝑔+ 20 + 40 − 6000 · (𝑉46 · 𝜋 ·
𝐷462
4)
2
= 80 + 6000 · (𝑉46 · 𝜋 ·𝐷462
4)
2
𝑃4𝜌 · 𝑔
+1
2 · 𝑔+ 20 + 40 − 6000 · (𝜋 ·
𝐷462
4)
2
= 80 + 6000 · (𝜋 ·𝐷462
4)
2
Y ahora aplicando las ecuaciones de continuidad a las intersecciones 2 y 4:
- Intersección 2:
𝑄1→2 = 𝑄2→3 + 𝑄2→4
𝑉12 · 𝜋 ·𝐷122
4= 𝑉23 · 𝜋 ·
𝐷232
4+ 𝑉24 · 𝜋 ·
𝐷242
4
𝐷122 = 𝐷23
2 +𝐷242
- Intersección 4:
𝑄2→4 = 𝑄4→5 + 𝑄4→6
𝑉24 · 𝜋 ·𝐷242
4= 𝑉45 · 𝜋 ·
𝐷452
4+ 𝑉46 · 𝜋 ·
𝐷462
4
𝐷242 = 𝐷45
2 + 𝐷462
Y resolviendo el sistema definido en el siguiente paréntesis:
Entrega 5 Mecánica de Fluidos Q1. Curso 2016/2017 M.A.B. , A.I.S. i A.P.E. 20/12/2016
20
{
100 − 3000 · (𝜋 ·
𝐷122
4)
2
=𝑃2𝜌 · 𝑔
+1
2 · 𝑔+ 20 + 5000 · (𝜋 ·
𝐷122
4)
2
𝑃2𝜌 · 𝑔
+1
2 · 𝑔+ 20 = 50 + 8000 · (𝜋 ·
𝐷232
4)
2
𝑃2𝜌 · 𝑔
=𝑃4𝜌 · 𝑔
+ 3000 · (𝜋 ·𝐷242
4)
2
𝑃4𝜌 · 𝑔
+1
2 · 𝑔+ 20 = 40 + 10000 · (𝜋 ·
𝐷452
4)
2
𝑃4𝜌 · 𝑔
+1
2 · 𝑔+ 20 + 40 − 6000 · (𝜋 ·
𝐷462
4)
2
= 80 + 6000 · (𝜋 ·𝐷462
4)
2
𝐷122 = 𝐷23
2 +𝐷242
𝐷242 = 𝐷45
2 + 𝐷462 }
Obtenemos que:
{
𝐷12 = 0.40176 𝑚𝐷23 = 0.25769 𝑚𝐷24 = 0.30840 𝑚𝐷45 = 0.22272 𝑚𝐷46 = 0.21280 𝑚
𝑃2 = 5.0804 · 105 𝑃𝑎
𝑃4 = 3.4460 · 105 𝑃𝑎}
Obviamente estos conductos no existen en la realidad, y aproximando los diámetros hemos
obtenido los siguientes valores:
𝐷12(𝑟𝑒𝑎𝑙) = 0.40 𝑚
𝐷23(𝑟𝑒𝑎𝑙) = 0.25 𝑚
𝐷24(𝑟𝑒𝑎𝑙) = 0.30 𝑚
𝐷45(𝑟𝑒𝑎𝑙) = 0.25 𝑚
𝐷46(𝑟𝑒𝑎𝑙) = 0.20 𝑚
Una vez definidos los valores de los diámetros reales, calculamos el valor de los cabales de
cada tramo, utilizando las expresiones siguientes:
Entrega 5 Mecánica de Fluidos Q1. Curso 2016/2017 M.A.B. , A.I.S. i A.P.E. 20/12/2016
21
{
100 − 3000 · 𝑄1−2
2 =𝑃2𝜌 · 𝑔
+8 · 𝑄1−2
2
𝜋2 · 𝑔 · 0.44+ 20 + 0.02 ·
500
0.45·8 · 𝑄1−2
2
𝜋2 · 𝑔
𝑃2𝜌 · 𝑔
+𝑄2−32
𝜋2 · 𝑔 · 0.254+ 20 = 50 + 0.02 ·
500
0.255·8 · 𝑄2−3
2
𝜋2 · 𝑔
𝑃2𝜌 · 𝑔
=𝑃4𝜌 · 𝑔
+ 0.02 ·500
0.305·8 · 𝑄2−4
2
𝜋2 · 𝑔
𝑃4𝜌 · 𝑔
+𝑄4−52
𝜋2 · 𝑔 · 0.254+ 20 = 40 + 0.02 ·
500
0.255·8 · 𝑄4−5
2
𝜋2 · 𝑔
𝑃4𝜌 · 𝑔
+𝑄4−62
𝜋2 · 𝑔 · 0.204+ 20 + 40 − 6000 · 𝑄4−6
2 = 80 + 0.02 ·500
0.205·8 · 𝑄4−6
2
𝜋2 · 𝑔𝑄1→2 = 𝑄2→3 + 𝑄2→4𝑄2→4 = 𝑄4→5 + 𝑄4→6 }
Obtenemos:
{
𝑄12 = 0.12708 𝑚
3/𝑠
𝑄23 = 0.01555 𝑚3/𝑠
𝑄24 = 0.11152 𝑚3/𝑠
𝑄45 = 0.08507 𝑚3/𝑠
𝑄46 = 0.02645 𝑚3/𝑠
𝑃2 = 2.9626 · 105 𝑃𝑎
𝑃4 = 2.5477 · 105 𝑃𝑎}
; 𝑉𝑖→𝑗 =4 · 𝑄𝑖→𝑗
𝜋 · 𝐷𝑖→𝑗2
Y así, en cada uno de los tramos tenemos que:
- Tramo 1-2:
𝑉𝑖→𝑗 =4 · 0.127076
𝜋 · 0.42= 1.0112 𝑚/𝑠
- Tramo 2-3:
𝑉𝑖→𝑗 =4 · 0.0155529
𝜋 · 0.252 = 0.3168 𝑚/𝑠
- Tramo 2-4:
𝑉𝑖→𝑗 =4 · 0.111523
𝜋 · 0.32= 1.577 𝑚/𝑠
- Tramo 4-5:
𝑉𝑖→𝑗 =4 · 0.085068
𝜋 · 0.25 2= 1.733 𝑚/𝑠
- Tramo 4-6:
𝑉𝑖→𝑗 =4 · 0.0264548
𝜋 · 0.22= 0.84208 𝑚/𝑠
Entrega 5 Mecánica de Fluidos Q1. Curso 2016/2017 M.A.B. , A.I.S. i A.P.E. 20/12/2016
22
Problema 4. En una central térmica de producción de energía eléctrica se tiene una
instalación cuyo esquema se muestra a continuación, siendo el fluido de trabajo agua.
Se sabe que la presión y la temperatura en la caldera de vapor son de P = 2 bar (absoluta), T
= 400ºC, y que en la turbina se produce una expansión adiabático-isentrópica con un salto
entálpico de 777 KJ/Kg. (Considérese este salto entre los puntos 3 y 6.)
Si las pérdidas de carga en la tubería de aspiración son Δh = 104 Q2 y en la tubería de
impulsión Δh = 312 Q2, siendo Q [m3/s], Δh [m columna de agua], y sabiendo que la bomba
que se utiliza es el modelo 150/315, con un diámetro de rodete de 270 mm, (y se considera
que la cota del nivel del líquido del condensador está 1m por encima de la cota del nivel del
líquido de la caldera,) se pide hallar:
1. El punto de funcionamiento de la bomba.
2. La cota Z (respecto al nivel del líquido del
condensador) a la que hay que colocar la bomba
para que no se produzca cavitación.
3. Debido a que se ha hecho un reajuste en el
proceso, se precisa aumentar el caudal en un 20%. Si
al motor de accionamiento de la bomba se le acopla
un variador de frecuencia, determine a qué
revoluciones debería girar para que la bomba
suministre el nuevo caudal. ¿Con qué rendimiento
trabaja ahora la bomba?
Antes de empezar debemos encontrar las condiciones termodinámicas que nos faltan a partir
de los diagramas adjuntos en el enunciado.
En la caldera, que es el punto 3, el enunciado
nos dice que:
𝑃3 = 2𝑏𝑎𝑟
𝑇3 = 400º𝐶
dónde la presión es absoluta.
En el segundo diagrama (el que muestra
entalpía-entropía para el vapor de agua)
buscamos el punto que se corresponde a 2bar
y 400ºC:
ℎ3 = 3250𝐾𝑗/𝐾𝑔
En la figura de la derecha mostramos el punto
de intersección encontrado con un punto azul.
Entrega 5 Mecánica de Fluidos Q1. Curso 2016/2017 M.A.B. , A.I.S. i A.P.E. 20/12/2016
23
Hemos estimado el valor de ℎ3, pues es muy difícil encontrar el valor exacto.
Por otro lado, el enunciado indica que el salto entálpico entre 3 y 6 es de 777Kj/Kg con una
expansión adiabático-isentrópica.
ℎ6 = (3250 − 777)𝐾𝑗/𝐾𝑔 ≅ 2500𝐾𝐽/𝐾𝑔
Introduciendo esta entalpía en el diagrama anterior encontramos (marcado en punto rojo):
𝑃6 = 0,05𝑏𝑎𝑟
𝑇6 = 35º𝐶
Habiendo encontrado esto ya podemos empezar con el ejercicio.
1. El punto de funcionamiento de la bomba.
Primero de todo aplicamos Bernoulli entre las superficies libres del condensador (subíndice co)
y la caldera (subíndice ca):
𝑃𝑐𝑜𝜌𝑔
+ 𝑍𝑐𝑜 +𝑉𝑐𝑜
2
2𝑔+ 𝐻 =
𝑃𝑐𝑎𝜌𝑔
+ 𝑍𝑐𝑎 +𝑉𝑐𝑎
2
2𝑔+ ∆ℎ7−2
dónde ∆ℎ7−2 es la variación de entalpía entre el punto 7 y el 2.
De aquí aislamos H:
𝐻 =𝑃𝑐𝑎 − 𝑃𝑐𝑜𝜌𝑔
+ 𝑍𝑐𝑎 − 𝑍𝑐𝑜 +𝑉𝑐𝑎
2 − 𝑉𝑐𝑜2
2𝑔+ ∆ℎ7−2
Considerando despreciables las energías cinéticas en las superficies libres de los depósitos:
𝐻 =𝑃𝑐𝑎 − 𝑃𝑐𝑜𝜌𝑔
+ 𝑍𝑐𝑎 − 𝑍𝑐𝑜 + ∆ℎ7−2
Sustituyendo los valores obtenidos anteriormente tenemos:
𝐻 =(2 − 0,05)105
1000 · 9,81− 1 + 416𝑄2
𝐻 = 18,877 + 416𝑄2 (1)
La intersección entre esta curva y la curva característica de la bomba da lugar al punto de
funcionamiento. Para obtenerla primero cogemos una serie de puntos de la gráfica de la
bomba, pasando el cabal a m3/s:
Q (m^3/h) Q (m^3/s) H
0 0 24
50 0,01388889 24
100 0,02777778 24
150 0,04166667 23,5
200 0,05555556 23
Entrega 5 Mecánica de Fluidos Q1. Curso 2016/2017 M.A.B. , A.I.S. i A.P.E. 20/12/2016
24
250 0,06944444 22
300 0,08333333 18
350 0,09722222 14
Ahora trazamos en una misma gráfica estos puntos unidos por una línea continua y la ecuación
(1) obtenida:
𝑄 = 0,0726𝑚3/𝑠 = 261,36𝑚3/ℎ
𝐻 = 21,0912𝑚
2. La cota Z (respecto al nivel del líquido del condensador) a la que hay que colocar la bomba
para que no se produzca cavitación.
Aplicamos Bernoulli entre las superficies libres del condensador (subíndice co) y la bomba
(subíndice 8):
𝑃𝑐𝑜𝜌𝑔
+ 𝑍𝑐𝑜 +𝑉𝑐𝑜
2
2𝑔=𝑃8𝜌𝑔
+ 𝑍8 +𝑉82
2𝑔+ ∆ℎ7−8
𝑍𝑐𝑜 − 𝑍8 =𝑃8 − 𝑃𝑐𝑜𝜌𝑔
+𝑉82 − 𝑉𝑐𝑜
2
2𝑔+ ∆ℎ7−8
Como no se conoce el diámetro del conducto y el término de energía cinética de la brida de
aspiración de la bomba es mucho mayor que en el condensador se desprecian ambos
términos. Por tanto NPSDHd queda:
𝑁𝑃𝑆𝐻𝑑 =𝑃8𝜌𝑔
−𝑃𝑣𝑎𝑝𝑜𝑟𝜌𝑔
Sustituyendo en Bernoulli:
𝑍𝑐𝑜 − 𝑍8 = 𝑁𝑃𝑆𝐻𝑑 +𝑃𝑣𝑎𝑝𝑜𝑟𝜌𝑔
−𝑃𝑐𝑜𝜌𝑔
+ ∆ℎ7−8
Entrega 5 Mecánica de Fluidos Q1. Curso 2016/2017 M.A.B. , A.I.S. i A.P.E. 20/12/2016
25
En el punto de cavitación incipiente se cumple que:
𝑁𝑃𝑆𝐻𝑑 = 𝑁𝑃𝑆𝐻𝑅
Buscamos en la gráfica de la curva característica de la bomba el valor de 𝑁𝑃𝑆𝐻𝑅 para un
caudal de 𝑄 = 261,36𝑚3/ℎ:
Se cumple aproximadamente 𝑁𝑃𝑆𝐻𝑅 = 3𝑚.
Ahora suponiendo que la temperatura del líquido en el condensador es la misma que en la
bomba (35ºC), la presión del vapor también será la misma que la del condensador (0,05bar):
𝑍𝑐𝑜 − 𝑍8 = 3 +0,05 · 105
𝜌𝑔−0,05 · 105
𝜌𝑔+ 104(
261,36
3600)2
= 3,55𝑚
3. Debido a que se ha hecho un reajuste en el proceso, se precisa aumentar el caudal en un
20%. Si al motor de accionamiento de la bomba se le acopla un variador de frecuencia,
determine a qué revoluciones debería girar para que la bomba suministre el nuevo caudal.
¿Con qué rendimiento trabaja ahora la bomba?
Nuevo caudal de trabajo: 𝑄′ = 261,36 ∗ 1,2𝑚3/ℎ = 313,63𝑚3/ℎ
Aplicando la ecuación (1) obtenida en el apartado 1 encontramos la nueva altura de
funcionamiento:
𝐻 = 18,877 + 416(313,63
3600)2
= 22,03𝑚
La constante de una curva que pase por el nuevo punto de trabajo se denomina curva de
afinidad y se halla como:
∆ℎ = 𝐾𝑄2
Entrega 5 Mecánica de Fluidos Q1. Curso 2016/2017 M.A.B. , A.I.S. i A.P.E. 20/12/2016
26
𝐾 =22,03
(313,633600
)2 = 2903,123
Trazamos en la gráfica que ya habíamos usado antes una 2ª ecuación, ∆ℎ = 2903,123𝑄2, y
hallamos el punto de intersección con la curva inicial:
𝑄 = 0,0805𝑚3/𝑠 = 289,8𝑚3/ℎ
𝐻 = 18,813𝑚
Aplicamos entre los puntos {𝐻 = 22,03𝑚, 𝑄 = 313,63𝑚3/ℎ} y {𝐻 = 18,813𝑚, 𝑄 =
289,8𝑚3/ℎ} los grupos adimensionales, cifra característica, altura de elevación y cifra
característica de caudal:
𝐻1
𝑤12𝐷1
2 =𝐻2
𝑤22𝐷2
2
Consideramos los diámetros iguales:
22,03
𝑤12=18,813
14502 → 𝑤1 = 1569,08𝑟𝑝𝑚
El rendimiento al que trabaja la bomba será prácticamente el mismo que antes, según
observamos en la gráfica característica.