Laboratorio de Mecánica de Fluidos: Estudio y patronamiento de medidores de caudal en conductos a presión.

SEXTO INFORME DEL LABORATORIO DE MECÁNICA DE FLUIDOS

PRÁCTICA X: ESTUDIO Y PATRONAMIENTO DE MEDIDORES DE CAUDAL EN

CONDUCTOS A PRESIÓN

PRESENTADO POR:

Kelly Johana Castillo Calambas

[email protected]
Carlos Mario Manzano Garzón
[email protected]
Cristian Camilo Solarte Urbano
[email protected]

PRESENTADO A:
Ing. Andrea Marcela Bolaños Rodríguez

UNIVERSIDAD DEL CAUCA

FACULTAD DE INGENIERÍA CIVIL
POPAYÁN, CAUCA
2022

1
 Laboratorio de Mecánica de Fluidos: Estudio y patronamiento de medidores de caudal en conductos a presión.

PRÁCTICA X: ESTUDIO Y PATRONAMIENTO DE MEDIDORES DE CAUDAL EN

CONDUCTOS A PRESIÓN

Introducción

El presente informe que trata sobre el “estudio y patronamiento de medidores de caudal en

conductos a presión”, está determinado por la utilización de un tubo de Venturi y de un
diafragma. La parte experimental de esta práctica consiste en obtener datos como el tiempo para
determinar el valor del caudal en el rotámetro y de lecturas de alturas (h n) tanto del tubo de
Venturi como del diafragma, para determinar finalmente el caudal que pasa a través de este
conducto a presión y además de los coeficientes de descarga, de velocidad, de contracción y
amplitud, correspondientes.

Objetivo

 Conocer varios sistemas de medición de caudal en conductos a presión.

 Determinar las ecuaciones de patronamiento de distintos dispositivos para medición de
caudal en conductos a presión.

Montaje experimental

Ilustración 1. Aparato para el estudio de medidores de caudal.

2
 Laboratorio de Mecánica de Fluidos: Estudio y patronamiento de medidores de caudal en conductos a presión.

Desarrollo procedimental
Datos:

 T = 20 ºC
2
cm
 v=1.004∗10−2
s

Tabla No.1. Datos experimentales de laboratorio.

Tubo de Venturi Diafragma

No. Obs Tiempo (s) V (Litros) Q (Litros/s) (Litros/s)
(Litros/min) (Litros/s) h1 h2 h3 h6 h7 h8
14.03 0.35638
1 13.56 5 0.36873 0.36461 22 0.36667 38.5 14.9 31.5 23.8 7.8 12.7
13.56 0.36873
15.35 0.32573
2 16.06 5 0.31133 0.32164 19 0.31667 34.5 16.9 28.6 21.3 9.4 13.1
15.25 0.32787
17.6 0.28409
3 17.88 5 0.27964 0.27914 17 0.28333 31.4 17.7 26.4 19.6 10.5 13.3
18.27 0.27367

El caudal hallado para el aforo estuvo determinado por el volumen en litros y el tiempo
correspondiente, siendo la ecuación del caudal la siguiente:

∀ Volumen
Q= =
t Tiempo
Ecuación 1. Caudal (Q)

Luego de esto, con ayuda del promedio aritmético para cada observación de los 3 tiempos
determinados, se halló un caudal medio (Q), el cual se registró en la anterior tabla.
n

∑ Qi
i=1
Q=
n
Ecuación 2. Caudal Promedio (Q ¿

Se puede ver claramente como el caudal del aforo obtenido, comparado con el que nos brinda el
mismo equipo con el rotámetro son muy parecidos, permitiéndonos así determinar que el error
entre ellos es muy mínimo.
En este caso, se va a trabajar con el valor dado por el equipo, es decir: Q rotámetro.
Seguidamente, procedemos a hallar el valor de las áreas tanto de la tubería como de la garganta
del diafragma y del Venturi, ya sabemos que está determinada por la siguiente ecuación:
2
π∗D
A=
4

3
Laboratorio de Mecánica de Fluidos: Estudio y patronamiento de medidores de caudal en conductos a presión.

Ecuación 3. Área (A1, A2)

Tabla No.2. Parámetros del diafragma y tubo de Venturi.

Venturi Diafragma
Diametro de tuberia (cm) 3.175 3.175
Diametro de garganta (cm) 1.5 2
2
A1: Área de tuberia (cm ) 7.917 7.917
2
A2: Área de garganta (cm ) 1.767 3.142

A. CALCULOS PARA EL MEDIDOR DE VENTURI

1. Velocidad real, velocidad teórica y el coeficiente de velocidad para la garganta del

tubo de Venturi.

 Velocidad real (V ¿¿ r )¿. Se determina para cada caudal la velocidad real en la garganta (
V r 2 ) utilizando la siguiente ecuación:

Q
V r=
A2

Ecuación 4. Velocidad real (V ¿¿ r )¿.

Donde:
- Q=Caudal
- A2=area de la garganta

Remplazando los valores de Q rotámetro dados en la Tabla No.1 y los valores de área
( A¿¿ 2)¿de la Tabla No.2, en la anterior ecuación, se calculan los respectivos valores de
cm3
las velocidades reales. Es necesario trabajar el caudal en .
s
Observación 1:
Litro 3
∗1000 cm
s
0.36667
1 Litro cm
V r= =207.508
1.767 cm 2
s

Observación 2:
Litro 3
∗1000 cm
s
0.31667
1 Litro cm
V r= =179.2114
1.767 cm 2
s

Observación 3:

4
Laboratorio de Mecánica de Fluidos: Estudio y patronamiento de medidores de caudal en conductos a presión.

Litro
∗1000 cm3
s
0.28333
1 Litro cm
V r= =160.3471
1.767 cm
2
s

 Velocidad teórica (V t ): Para la determinación de la velocidad teórica en el caso de la

garganta (V t ) del Venturi, es necesario primero, calcular la diferencia de presiones entre
la entrada y la garganta. Se debe utilizar la siguiente expresión:

∆ h=h 1−h2
Ecuación 5. Diferencia de presiones (∆ h) .

Tabla No.3. Datos correspondientes a las diferencias de presiones entre la entrada y la garganta para
el medidor de Venturi

Tubo de Venturi (cm)

No. Obs
h1 h2 h3 ∆h
1 38.5 14.9 31.5 23.6
2 34.5 16.9 28.6 17.6
3 31.4 17.7 26.4 13.7

Luego, con los datos registrados en las Tablas No.2 y 3, se procede a calcular las
velocidades teóricas en el medidor de Venturi para cada caudal (u observación), utilizando
la siguiente expresión.

√()
2∗g∗∆ h
V t=
A2 2
1−
A1

Ecuación 6. Velocidad teórica (V ¿¿ t) ¿.

Observación 1:

√
cm
2∗981 2
∗23.6 cm
s cm
V t= =220.75025

( )
1.767 cm
2 2 s
1−
7.917 cm2

Observación 2:

5
Laboratorio de Mecánica de Fluidos: Estudio y patronamiento de medidores de caudal en conductos a presión.

√
cm
2∗981 2
∗17.6 cm
s cm
V t= =190.63451

( )
1.767 cm
2 2 s
1−
7.917 cm2

Observación 3:

√
cm
2∗981 2
∗13.7 cm
s cm
V t= =220.75

( )
1.767 cm
2 2 s
1−
7.917 cm 2

 Coeficiente de velocidad (Cv). Ahora se procede a calcular el coeficiente de velocidad,

utilizando la expresión:
V r2
C v=
Vt2
Ecuación 7. Coeficiente de velocidad (Cv).

Donde, los valores de las velocidades reales y teóricas para cada observación o caudal, ya
fueron obtenidas en los anteriores puntos, permitiéndonos así obtener los siguientes
resultados:
Observación 1:
cm
207.508
s
C v= =0.94001
cm
220.75
s
Observación 2:
cm
179.211
s
C v= =0.94008
cm
190.634
s
Observación 3:
cm
160.347
s
C v= =0.95336
cm
168.192
s

Tabla No. 4. Velocidades reales, teóricas y los coeficientes de velocidad en el medidor de Venturi para
cada caudal en el caso de la garganta.

6
 Laboratorio de Mecánica de Fluidos: Estudio y patronamiento de medidores de caudal en conductos a presión.

Tu bo de Venturi
No. Obs
Q (Litros/s) ∆h (cm) V real V teórica
1 22 23.6 207.50802 220.75025 0.94001
2 19 17.6 179.21147 190.63451 0.94008
3 17 13.7 160.34710 168.19206 0.95336

2. Número de Reynolds (ℜ) en la garganta, para la velocidad real. Este valor de número
de Reynolds está determinado por la siguiente ecuación:

V r∗D g
ℜ=
v
Ecuación 8. Número de Reynolds en el tubo Venturi para la garganta (Re2).

Donde:
- V r =Velocidad real
- D g=diametro de la garganta
- v=viscosidad cinematica

Finalmente, reemplazando los datos de la Tabla No. 1 y 2, junto con los valores de las
velocidades reales obtenidas para cada caudal, tenemos que:
Observación 1:
cm
207.508
∗1.5 cm
s
ℜ= =31002.19
1.004∗10−2 cm2 /s
Observación 2:
cm
179.211∗1.5 cm
s
ℜ= −2 2
=26774.62
1.004∗10 cm / s
Observación 3:
cm
160.347
∗1.5 cm
s
ℜ= 2
=23956.24
−2 cm
1.004∗10
s

Realizando la gráfica (Cv vs Re) de los valores obtenidos experimentalmente en el tubo de

Venturi, tenemos:

7
Laboratorio de Mecánica de Fluidos: Estudio y patronamiento de medidores de caudal en conductos a presión.

Tubo de Venturi
No. Obs
Q (Litros/min) Re2
1 22 0.94001 31002.19383
2 19 0.94008 26774.62194
3 17 0.95336 23956.24069

Gráfica No.1. Coeficientes de Velocidades (Cv) vs Reynolds (Re) para el medidor Venturi.
Coeficiente de velocidad (Cv)

Medidor Venturi (Cv vs Re)

0.955
0.953
0.951
0.949
0.947
0.945
0.943
0.941
0.939
0.937
0.935
23000 Reynolds (Re): en escala logaritmica

3. La gráfica obtenida de Cv vs Re, sobre la Figura X.4. Este sería el resultado final:

4. Determinación de las perdidas por el Venturi metro (hp), la contracción (hc) y el

coeficiente Kc.

 Cálculo de las perdidas por el Venturi metro (h p): Donde se tiene que la pérdida total
de carga en el dispositivo (h p ) es igual a la diferencia de energía entre la entrada (h¿ ¿1)¿
y la salida (h3 ) del Venturi metro.
h p =h1−h3

8
Laboratorio de Mecánica de Fluidos: Estudio y patronamiento de medidores de caudal en conductos a presión.

Ecuación 9. Perdidas del Venturi metro (hp)

Con ayuda de los datos anexados en la tabla No.1 de las alturas presentes en el Venturi
metro, podemos obtener los siguientes resultados:
Observación 1:
h p =38.5 cm−31.5 cm=7 cm

Observación 2:
h p =34.5 cm−28.6 cm=5.9 cm

Observación 3:
h p =31.4 cm−26.4 cm=5 cm

 Cálculo del coeficiente: KC. El cual se obtiene para cada caudal a través de la siguiente
expresión:

) ( ( ))
2

(1 A2
K c = 2 −1 ∗ 1−
cv A1

Ecuación 10. Coeficiente Kc

Con ayuda de los datos como los coeficientes de velocidad (Tabla No.4) y los de área
(Tabla No.2), presentes en el Venturi metro, podemos reemplazar los valores
correspondientes para obtener así, los siguientes resultados:
Observación 1:

( )( ( ))
2 2
1 1.767 cm
K c= 2
−1 ∗ 1− 2
=0.1251
(0.94001) 7.917 cm

Observación 2:

( )( ( ))
2 2
1 1.767 cm
K c= 2
−1 ∗ 1− 2
=0.1249
(0.94008) 7.917 cm

Observación 3:

( )( ( ))
2 2
1 1.767 cm
K c= 2
−1 ∗ 1− 2
=0.0952
(0.95336) 7.917 cm

 Cálculo de la contracción (hc). Calculados los coeficientes KC, se puede proceder a

determinar la contracción hc a partir de la siguiente ecuación:
2
K c ∗V t 2
hc=
2∗g

9
 Laboratorio de Mecánica de Fluidos: Estudio y patronamiento de medidores de caudal en conductos a presión.

Ecuación 10. Coeficiente Kc

Por lo tanto, al reemplazar los datos correspondientes se tiene:

Observación 1:

( )
2
cm
0.1251∗ 220.7502
s
hc= =3.1082 cm
cm
2∗981 2
s
Observación 2:

( ) =2.3152 cm
2
cm
0.1249∗ 190.6345
s
hc=
cm
2∗981
s2
Observación 3:

( ) =1.3733 cm
2
cm
0.0952∗ 168.192
s
hc=
cm
2∗981 2
s
Tabla No.5. Valores de perdidas por el Venturi metro (hp), la contracción (hc) y el coeficiente Kc para cada
caudal.

Tubo de Venturi
No. Obs
h p (cm) h c (cm) Kc
1 7 3.1082013 0.125142763
2 5.9 2.3151751 0.124991301
3 5 1.3733328 0.095249747

5. Cálculo del caudal teórico (Qt ). Determinación del caudal teórico para cada observación,
mediante la siguiente expresión.
Q t = A2∗V t

Ecuación 11. Caudal teórico (Q t ¿

Donde:
- A2=área de garganta
- V t =velocidad teorica

Reemplazando los respectivos valores en la anterior ecuación, se tiene que:

Observación 1:

10
 Laboratorio de Mecánica de Fluidos: Estudio y patronamiento de medidores de caudal en conductos a presión.

2 cm 3
Qt =1.767 cm ∗220.75 =390.0657 cm / s
s
Observación 2:
2 cm 3
Qt =1.767 cm ∗190.63 =336.8512 cm /s
s
Observación 3:
2 cm 3
Qt =1.767 cm ∗168.19 =297.1954 cm /s
s

6. Cálculo del coeficiente de descarga (Cd) del medidor Venturi. A continuación, se

procede realizar el cálculo del coeficiente de descarga para cada caudal, utilizando la
siguiente ecuación:
Qr
C d=
A2 √ 2∗g∗∆ h

Ecuación 12. Coeficiente de descarga (C d ¿

Donde:
- Q r =caudal del rotametro
- A2=area de la garganta
- ∆ h=diferencia de presiones en la entrada y garganta
Gracias a los datos brindados en la Tabla No. 1, 2 y 3, podemos hallar el correspondiente valor de
coeficiente para cada observación, siendo así estos:
Observación 1:
Litro 3
∗1000 cm
s
0.36667
1 Litro
C d= =0.9643
2

√ cm
1.767 cm ∗ 2∗981 2 ∗23.6 cm
s
Observación 2:
Litro 3
∗1000 cm
s
0.31667
1 Litro
C d= =0.9644
2

√ cm
1.767 cm ∗ 2∗981 2 ∗17.6 cm
s
Observación 3:

11
 Laboratorio de Mecánica de Fluidos: Estudio y patronamiento de medidores de caudal en conductos a presión.

Litro 3
∗1000 cm
s
0.28333
1 Litro
C d= =0.9780
2

√ cm
1.767 cm ∗ 2∗981 2 ∗13.7 cm
s
Tabla No. 6. Coeficiente de descarga en el medidor de Venturi para cada caudal, junto con el promedio y el
número de Reynolds hallado anteriormente.

Tubo de Venturi
No. Obs
Q (Litros/min) Cd Re2
1 22 0.96434 31002.19383
2 19 0.96441 26774.62194
3 17 0.97803 23956.24069
Promedio 19.33333 0.96892 27244.35215

A continuación, se grafica el coeficiente de descarga vs Reynolds (Cd vs Re).

Gráfica No. 2. Coeficientes de descarga vs Reynolds para el medidor Venturi.
Coeficiente de descarga (Cd)

0.980

0.975

0.970

0.965

0.960

0.955
23500
Reynolds (Re) : en escala logaritmica

7. La gráfica obtenida de Cd vs Re, sobre la Figura X.5. Este sería el resultado final:

12
 Laboratorio de Mecánica de Fluidos: Estudio y patronamiento de medidores de caudal en conductos a presión.

8. Cálculo del coeficiente Ci.

Se procede a calcular el coeficiente Ci para cada caudal, que depende del tipo de Venturi y de la
relación d/D dado por el fabricante, dicho coeficiente se calcula mediante la siguiente expresión:
hp
C i=
h1−h 2
Ecuación 13. Coeficiente Ci del medidor (C i ¿

Con los valores anteriormente calculados y registrados en las tablas No. 4 y 5 se realiza el
correspondiente reemplazo, teniendo así, los siguientes resultados:
7 cm
Observación 1: C i= =0.2966
23.6 cm
5.9 cm
Observación 2: C i= =0.3352
17.6 cm
5 cm
Observación 3: C i= =0.3649
13.7 cm
Tabla No.7. Coeficiente Ci para cada caudal registrado en la práctica para el medidor de Venturi.

13
 Laboratorio de Mecánica de Fluidos: Estudio y patronamiento de medidores de caudal en conductos a presión.

Tubo de Venturi
No. Obs
Ci
1 0.29661
2 0.33523
3 0.36496
Promedio 0.33227

9. Cálculo de la ecuación de patrona miento del medidor de Venturi

Para establecer la ecuación de patrona miento del medidor Venturi se debe calcular antes el
coeficiente K, el cual está dado por la siguiente expresión:
K=Cd∗A2∗√ 2∗g
Ecuación 14. Valor de la constante K

Por lo tanto, para determinar dicho coeficiente K, se trabajará con el valor más probable del
coeficiente de descarga “Cd”, es decir con el promedio del valor de Cd registrado en la Tabla
No.6.
Así, el valor del coeficiente K para el medidor de Venturi esta dado por:

K=0.96892∗1.767 cm ∗ 2∗981
2

√ cm
s
2
=75.83575

Así, la ecuación de patrona miento para el medidor Venturi estará definida como:
m
Q=K∗∆ h
Ecuación 15. Ec. De Patronamiento K

con :m ≈ 0.5
0.5
Q=75.83575∗∆ h
Gráfica No.3. Primera curva de patrona miento para el medidor de Venturi, obtenido encontrando el valor de
K con la ecuación 14 y 15.

14
 Laboratorio de Mecánica de Fluidos: Estudio y patronamiento de medidores de caudal en conductos a presión.

Podemos determinar otra ecuación de patrona miento a partir de la utilización del método de
mínimos cuadrados. Para ello es necesario realizar linealización a la ecuación 15:
Se tiene que, aplicando logaritmo y sus propiedades a ambos lados:
m
Q=K∗∆ h
ln ( Qr )=ln ( K∗∆ h )
m

ln ( Qr )=ln ¿ ¿
ln ( Qr )=ln ¿ ¿

Así finalmente se tiene una ecuación lineal tipo:

y=ax+ b
Donde, para nuestro caso:
- y =ln ( Qr )
- a=m
- x=ln ⁡(∆ h)
- b=ln ( K )

Ahora bien, se procede a calcular los valores necesarios para la realización de mínimos
cuadrados.
Tabla No.8. Datos para aplicar mínimos cuadrados y determinar parámetros de la ecuación de patrona
miento para el medidor de Venturi.

15
 Laboratorio de Mecánica de Fluidos: Estudio y patronamiento de medidores de caudal en conductos a presión.

No. Obs ∆h (cm) Qr (cm^3/s) Xi (Ln (∆h)) Yi (Ln (Qr)) Xi * Yi Xi * Xi

1 23.6 366.6667 3.1612 5.9045 18.6654 9.9935
2 17.6 316.6667 2.8679 5.7578 16.5129 8.2248
3 13.7 283.3333 2.6174 5.6466 14.7795 6.8508
Sumatoria (∑) 8.6465 17.3089 49.9578 25.0691

Para hallar el valor de la pendiente (m) o el valor de a, que en este caso es el exponente, viene
dado por la siguiente formula:
m=n∗¿ ¿

3∗( 49.9578 )−8.6465∗17.3089

m=
3∗25.0691−¿¿
m=0.4748
Ahora bien, calculamos el intercepto “b”, con la siguiente ecuación:
b=¿ ¿
17.3089∗25.0691−8.6465∗49.9578
b=
3∗25.0691−¿ ¿
b=4.40116
De la linealización se tiene que: b=ln ( K )
Por lo tanto, al despejar K se tiene que:
b
K=e
4.40116
K=e
K=81.5458

Finalmente, la ecuación de patrona miento para el medidor de Venturi que se obtiene por el
método de mínimos cuadrados es:
m
Q=K∗∆ h
0.4748
Q=81.5458∗∆ h
Gráfica No.4. Segunda curva de patrona miento para el medidor de Venturi, obtenida encontrando el valor de
K y m, a través de mínimos cuadrados.

16
 Laboratorio de Mecánica de Fluidos: Estudio y patronamiento de medidores de caudal en conductos a presión.

10. Comparación de los coeficientes de descarga obtenidos por diferentes métodos.

Ahora, a partir del valor de K obtenido por mínimos cuadrados, se puede calcular el valor del
coeficiente de descarga característico para el medidor de Venturi, de la siguiente manera:
K
C d=
A2∗√ 2∗g
Ecuación 16. Cd obtenida con el valor de K

Reemplazando los valores correspondientes, se obtiene:

81.5458
C d= =1.04187
2

√ cm
1.767 m ∗ 2∗981 2
s
Haciendo la correspondiente comparación de los dos Cd, se tiene que:
- Cdp (obtenida con el promedio) = 0.9689
- Cdk (obtenida con el valor de K, hallado) = 1.04187
La desviación entre los valores es de:
ε %=¿ C dp−C dk ∨ ¿ ∗100 % ¿
Cdk
Ecuación 17. Error relativo porcentual del coeficiente de descarga

ε %=¿ 0.9689−1.04187∨ ¿ ∗100 % ¿

1.04187

ε %=7.004

11. Finalmente tenemos los siguientes resultados para el medidor de Venturi:

Tabla No.9. Datos y resultados para el medidor de Venturi.

17
 Laboratorio de Mecánica de Fluidos: Estudio y patronamiento de medidores de caudal en conductos a presión.

Datos
N° Obs
(cm^ 3/s) (cm^ 3/s) h1 (cm) h2 (cm) h3 (cm)
1 364.614105 366.6667 38.5 14.9 31.5
2 321.644752 316.6667 34.5 16.9 28.6
3 279.135218 283.3333 31.4 17.7 26.4

Cálculos
N° Obs
∆h
1 23.6 7 207.5080 220.7502 390.0657 0.9400 0.9643 0.2966 3.1082 0.1251 31002.1938
2 17.6 5.9 179.2115 190.6345 336.8512 0.9401 0.9644 75.8358 0.3352 2.3152 0.1250 26774.6219
3 13.7 5 160.3471 168.1921 297.1954 0.9534 0.9780 0.3650 1.3733 0.0952 23956.2407

B. CALCULOS PARA EL DIAFRAGMA.

1. Diferencia de presiones y las perdidas por carga en el diafragma. En este caso es

necesario conocer la diferencia de presiones entre la entrada y la garganta del diafragma,
para esto, se debe utilizar la misma estructura que la ecuación 5, solo que esta vez, h 1 es h6
y h2 es h7, y para el caso de las perdidas por carga, esta se realiza utilizando igualmente la
estructura de la ecuación 9, solo que h 1 es h6 y h3 es h8. Finalmente, el resultado es el
siguiente:

Tabla No.10. Datos correspondientes a las diferencias de presiones y los valores calculados de perdida
por carga para cada caudal del diafragma.

Diafragma (cm)
No. Obs
h6 h7 h8 ∆h hp
1 23.8 7.8 12.7 16 11.1
2 21.3 9.4 13.1 11.9 8.2
3 19.6 10.5 13.3 9.1 6.3

2. Cálculo del coeficiente de descarga (Cd) del diafragma. La ecuación necesaria para el
cálculo del coeficiente es la siguiente (donde el valor de A0 =A 2, que en este caso es el de
la garganta del diafragma, recordemos que ese valor está en la tabla No.2):
Qr
C d=
A0 √ 2∗g∗∆ h
Ecuación 18. Coeficiente de descarga para el diafragma

Observación 1:

18
Laboratorio de Mecánica de Fluidos: Estudio y patronamiento de medidores de caudal en conductos a presión.

3
cm
366 .66 67
s
C d= =0.6587
2

√ cm
3.142 cm ∗ 2∗981 2 ∗16 cm
s
Observación 2:
3
cm
316 .66 67
s
C d= =0.6596
2

√ cm
3.142 cm ∗ 2∗981 2 ∗11.9 cm
s
Observación 3:
3
cm
283.33 33
s
C d= =0.6749
2

√ cm
3.142 cm ∗ 2∗981 2 ∗9.1 cm
s

Tabla No.11. Coeficiente de descarga en el Diafragma para cada caudal y su valor medio.

Diafragma (cm)
No. Obs
Qr (cm^3/s) ∆h A0 Cd
1 366.6667 16 3.142 0.6587
2 316.6667 11.9 3.142 0.6596
3 283.3333 9.1 3.142 0.6749
Promedio de Cd 0.6644
.
3. Cálculo de las velocidades reales en el diafragma ( V r ). En este caso para el diafragma
se utilizará la siguiente expresión, recordando que A0 será igual a A2.
Qr
V r=
A0

Observación 1:
3
cm
366.6667
s cm
V r= =116.6985
3.142 cm 2
s

Observación 2:

cm3
316.6667
s cm
V r= =100.7851
3.142 cm
2
s

19
 Laboratorio de Mecánica de Fluidos: Estudio y patronamiento de medidores de caudal en conductos a presión.

Observación 3:

cm3
283.3333
s cm
V r= 2
=90.1761
3.142 cm s

4. Cálculo del número de Reynolds en el diafragma ( ℜ). Con los valores anteriormente
calculados tanto de velocidad como del diámetro en el diafragma, se procede a determinar
el número de Reynolds:
V r∗D o
ℜ=
v
Donde:
- Do =diametro de la garganta del diafragma
- v=viscosidad cinematica

Observación 1:
cm
116.6985 ∗2 cm
s
ℜ= −2 2
=23246.712
1.004∗10 cm / s
Observación 2:
cm
100.7851 ∗2cm
s
ℜ= −2 2
=20076.706
1.004∗10 cm / s
Observación 3:
cm
90.1761 ∗2 cm
s
ℜ= 2
=17963.368
−2 cm
1.004∗10
s

5. Se nos permite obtener la gráfica de Cd vs Re, del diafragma. En este caso, se

realizará encima de la Figura X.10.
Para esto, utilizamos la siguiente tabla:
Tabla No. 12. Numero de Reynolds para las velocidades reales de cada caudal en el diafragma y los valores del
coeficiente de descarga.

20
 Laboratorio de Mecánica de Fluidos: Estudio y patronamiento de medidores de caudal en conductos a presión.

No. Obs Cd Re
1 0.6587 23246.7119
2 0.6596 20076.7057
3 0.6749 17963.3683

Gráfica No. 3. Coeficientes de descarga vs Reynolds para el diafragma.

Ahora, ubicando estos puntos obtenidos experimentalmente de Cd y Re, respectivamente en la

figura IV.10 de la guía, tenemos que:

6. Obtención del coeficiente Ci para el diafragma

21
 Laboratorio de Mecánica de Fluidos: Estudio y patronamiento de medidores de caudal en conductos a presión.

Se procede a determinar el coeficiente de perdida de carga para el diafragma “Ci”, relativo a la

caída de carga por efecto de la contracción, utilizando la siguiente expresión:
hp
C i=
h6−h 7

Para el cálculo de este coeficiente, se deben trabajar con los valores ya determinados de ∆ h y hp
para el diafragma y registrados en la tabla No.10.
11.1 cm
Observación 1: C i= =0. 69375
1 6 cm
8.2 cm
Observación 2: C i= =0. 68907
1 1.9 cm
6.3 cm
Observación 3: C i= =0. 69231
9.1 cm
Tabla No. 13. Coeficiente Ci para cada caudal registrado en la práctica (diafragma).

No. Obs Ci
1 0.69375
2 0.68908
3 0.69231
Promedio 0.69171

7. Cálculo de la ecuación de patrona miento del diafragma

Ahora bien, para establecer la ecuación de patrona miento del diafragma se debe calcular antes el
coeficiente K, el cual está dado por la ecuación 14:
K=Cd∗A0∗√ 2∗g

Por lo tanto, para determinar dicho coeficiente K, se trabajará con el valor más probable del
coeficiente de descarga “Cd”, es decir con el promedio de los Cd: 0.6644.
Así, el valor del coeficiente K para el diafragma, está dado por:

K=0.66 44∗3.142 cm2∗ 2∗981

√ cm
s
2
=92. 462611

Así, la ecuación de patrona miento para el diafragma estará definida como:

m
Q=K∗∆ h
con :m ≈ 0.5
0.5
Q=92.462611∗∆ h

22
 Laboratorio de Mecánica de Fluidos: Estudio y patronamiento de medidores de caudal en conductos a presión.

Gráfica No.5. Primera curva de patrona miento para el diafragma, obtenido encontrando el valor de K con la
ecuación 14 y 15.

Podemos determinar otra ecuación de patrona miento a partir de la utilización del método de
mínimos cuadrados. Para ello es necesario realizar linealización a la ecuación 15:
Se tiene que, aplicando logaritmo y sus propiedades a ambos lados:

Q=K∗∆ hm
ln ( Qr )=ln ( K∗∆ h m )
ln ( Qr )=ln ¿ ¿
ln ( Qr )=ln ¿ ¿

Así finalmente se tiene una ecuación lineal tipo:

y=a x+b
Donde, para nuestro caso:
- y =ln ( Qr )
- a=m
- x=ln ⁡(∆ h)
- b=ln ( K )

Ahora bien, se procede a calcular los valores necesarios para la realización de mínimos
cuadrados.
Tabla No. 14. Datos para aplicar mínimos cuadrados y determinar parámetros de la ecuación de patro
namiento para el diafragma.

23
 Laboratorio de Mecánica de Fluidos: Estudio y patronamiento de medidores de caudal en conductos a presión.

No. Obs ∆h (cm) Qr (cm^3/s) Xi (Ln (∆h)) Yi (Ln (Qr)) Xi * Yi Xi * Xi

1 16 366.6667 2.7726 5.9045 16.3706 7.6872
2 11.9 316.6667 2.4765 5.7578 14.2595 6.1332
3 9.1 283.3333 2.2083 5.6466 12.4693 4.8765
Sumatoria (∑) 7.4574 17.3089 43.0995 18.6970

Para hallar el valor de la pendiente (m) o el valor de a, que en este caso es el exponente, viene
dado por la siguiente formula:
m=n∗¿ ¿

3∗( 4 3.0995 ) −7.4574∗17.3089

m=
3∗18.697−¿ ¿
Ahora bien, calculamos el intercepto “b”, con la siguiente ecuación:
b=¿ ¿
17.3089∗18.697−7.4574∗4 3.0995
b=
3∗18.697−¿¿
De la linealización se tiene que: b=ln ( K ).
Por lo tanto, al despejar K se tiene que:

K=e b
K=e 4.6323
K=102.7468607

Finalmente, la ecuación de patrona miento para el diafragma que se obtiene por el método de
mínimos cuadrados es:
0.4575
Q=102.7468607∗∆ h
Gráfica No.5. Segunda curva de patrona miento para el diafragma, obtenida encontrando el valor de K y m, a
través de mínimos cuadrados.

24
 Laboratorio de Mecánica de Fluidos: Estudio y patronamiento de medidores de caudal en conductos a presión.

8. Comparación de los coeficientes de descarga obtenidos por diferentes métodos.

A partir del valor de K obtenido por mínimos cuadrados, se puede calcular el valor del
coeficiente de descarga característico para el diafragma, de la siguiente manera:
K
C d=
A o∗√ 2∗g

Reemplazando los valores correspondientes, se obtiene:

102.7468607
C d= =2.56924
3.142 cm ∗√ 2∗981 cm/ s
2 2

Haciendo la correspondiente comparación de los dos Cd, se tiene que:

- Cdp (obtenida con el promedio) = 0.6644
- Cdk (obtenida con el valor de K, hallado) = 2.56924
La desviación entre los valores es de:
ε %=¿ C dp−C dk ∨ ¿ ∗100 % ¿
Cdk
Ecuación 17. Error relativo porcentual del coeficiente de descarga

ε %=¿ 0.6644−2.56924∨ ¿ ∗100 % ¿

2.56924

ε %=74,14

9. Finalmente tenemos los siguientes resultados para el diafragma:

Tabla No. 15. Datos y resultados para el diafragma.

25
 Laboratorio de Mecánica de Fluidos: Estudio y patronamiento de medidores de caudal en conductos a presión.

Datos Cálculos
No. Obs
(cm^ 3/s) (cm^ 3/s) h6 (cm) h7 (cm) h8 (cm) ∆h (cm) hp (cm) Cd K Ci Re
1 364.6141 366.6667 23.8 7.8 12.7 16 11.1 0.0065 0.6938 23246.7119
2 321.6448 316.6667 21.3 9.4 13.1 11.9 8.2 0.0105 102.74686 0.6891 20076.7057
3 279.1352 283.3333 19.6 10.5 13.3 9.1 6.3 0.0158 0.6923 17963.3683

C. COMPARACIÓN DE LOS COEFICIENTES DE DESCARGA Y DE LAS

CURVAS DE PATRONA MIENTO ENTRE EL VENTURI Y EL DIAFRAGMA.

 Comparación de los coeficientes de descarga del Venturi y el diafragma.

Tabla No. 16. Datos para comparar los coeficientes de descarga del Venturi y el diafragma.

Tipo Venturi Diafragma

por promedio 0.9689 0.6644
por mínimos cuadrados 1.04187 2.56924
Desviación % entre 7.004% 74.14%

 Comparación de las curvas de patrona miento del Venturi y diafragma, a través de

la obtención de K por ecuación.

Tabla No.17. Datos para comparar curva de patrona miento del Venturi y diafragma.

Venturi Ecuación de Patronamiento

∆h (cm) Qr (cm^3/s)
23.6 368.4088
17.6 318.1488
13.7 280.6947
Diafragma Ecuación de Patronamiento
∆h (cm) Qr (cm^3/s)
16 369.8504
11.9 318.9625
9.1 278.9246

Gráfica No. 6. Curvas de patrona miento para el diafragma y el Venturi, obtenida encontrando el valor de K
con ecuación.

26
 Laboratorio de Mecánica de Fluidos: Estudio y patronamiento de medidores de caudal en conductos a presión.

Curvas de patronamiento
390

370
f(x) = 92.4626109101025
f(x) = x^0.500000000000002
75.8357499999995 x^0.500000000000002
350
Q (cm^3/s)

330

310

290

270

250
5 7 9 11 13 15 17 19 21 23 25

∆h (cm)

 Comparación de las curvas de patrona miento del Venturi y diafragma, a través de

la obtención de K por método de mínimos cuadrados.

Tabla No.18. Datos para comparar curva de patrona miento del Venturi y diafragma.

Venturi Ecuación de Patronamiento

∆h (cm) Qr (cm^3/s)
23.6 366.6667
17.6 316.6667
13.7 283.3333
Diafragma Ecuación de Patronamiento
∆h (cm) Qr (cm^3/s)
16 366.6667
11.9 316.6667
9.1 283.3333

Gráfica No.7. Segunda curva de patrona miento para el diafragma, obtenida encontrando el valor de K y m, a
través de mínimos cuadrados.

27
 Laboratorio de Mecánica de Fluidos: Estudio y patronamiento de medidores de caudal en conductos a presión.

ANALISIS DE RESULTADOS
En lo que respecta a la práctica podemos encontrar que primeramente gracias a la obtención de
los diferentes valores de lo que respecta al desarrollo de esta como los coeficientes y las curvas
de patrona miento, justo en estos datos va a verse reflejado entre lo que se obtiene
experimentalmente y lo que en verdad se espera de la teoría.
Primero que todo, dando un desarrollo circunstancial de lo anterior, si analizamos y comparamos
los valores de velocidades para cada caudal registrados en la tabla no. 4, tanto teóricos como
reales en el caso del medidor de Venturi, se evidencia que las velocidades teóricas obtenidas son
mayores que las velocidades reales registradas en la práctica. Con lo anterior, es posible
corroborar que lo obtenido por teoría, en general sus resultados van a ser mayores a los resultados
reales, esto debido a que, en este caso, las velocidades teóricas consideran un fluido y un
conducto ideal (sin considerar las perdidas), mientras que en la realidad las perdidas si existen y
se deben tener en cuenta, ya que afectan directamente la velocidad de fluido.
En otro caso, si observamos los valores del Cd de descarga obtenidos por el método de promedios
y el de mínimos cuadrados, podemos ver como este en ambos casos, varía en un porcentaje que
supera al admitido, es así como en la tabla no.16 miramos como en el Venturi este se presenta de
un 7% y en el caso del diafragma, es de un 74%. Una posible explicación es por lo mencionado
anteriormente, el valor hallado por promedios es básicamente un valor teórico, obtenido por
valores como el caudal, el área y la diferencia de presiones, mientras que el obtenido hallando K,
por mínimos cuadrados, es un valor real o experimental que fue deducido gracias a la curva que
acabamos de hallar, este al ser un valor inferido, va a presentar errores o diferencias, pero que en

28
 Laboratorio de Mecánica de Fluidos: Estudio y patronamiento de medidores de caudal en conductos a presión.

cierta manera, va a ser más representativo porque nos va a permitir inferir mas datos de la
experimentación.
Lo mencionado anteriormente, también aplicaría para la obtención de K y de las ecuaciones de
patrona miento. Es por esto, que los valores varían en gran proporción y que, en sí, las más
ajustadas a los datos siempre va a ser la experimental y la que en sí, nos sirva para el posterior
estudio que deseemos realizar a futuro.
Es preciso mencionar que los procedimientos aplicados en la práctica, presentaron problemas
puesto, que ambos métodos si logran mostrar una gran diferencia en lo que respecta a obtención
de curvas de patrona miento, pues sus valores en K varían en 20 unidades, siendo así, no posible
hacer una comparación a margen de cuál sería la más beneficiosa para la obtención de caudales.
En cierto caso, es posible mencionar que dicha variación se pudo haber ocasionado debido a las
lecturas que se realizaron en el laboratorio, ya que estas dependían mucho de la agudeza visual
que se tenía para realizar las respectivas lecturas de las alturas piezométricas.
En lo que respecta a la obtención de caudales por aforo y por rotámetro, en sí, al comparar estos
valores ellos son demasiado mínimos, pero, aun así, presenta diferencias, las cuales se pueden
justificar o asociar, directamente en la toma de datos específicamente en los tiempos. O en otro
caso, en las observaciones al leer el caudal en el rotámetro, ya que la vista puede fallar en
ocasiones realizando en todo caso lecturas no tan exactas.
Casi finalizando, tenemos que al observar la gráfica no.1 (coeficiente de velocidad vs Reynolds
del medidor de Venturi), se tiene que en cuanto mayor es el número de Reynolds, el coeficiente
de velocidad tiende a 1, es decir que en este aumento, va a llevar a una mejor relación entre la
velocidad real y la velocidad teórica que se están presentando.
Al analizar y comparar los resultados obtenidos de los coeficientes de descarga para el medidor
de Venturi y el diafragma podemos evidenciar en base a la tabla no. 16, que el diafragma posee
un menor coeficiente de perdida por descarga, en lo que respecta al obtenido por promedio, pero
que el obtenido por mínimos cuadrados es mayor a todos los demás; es por esto, que tal vez, sería
más aprovechable trabajar con el Venturi, en el momento de medición y análisis de caudales,
puesto que precisamente dicho coeficiente nos va indicar las posibles pérdidas presentes en el
recorrido del fluido en el conducto, y al no haber tanta certeza en el del diafragma, sería más
aprovechable el Venturi, o en otro caso, realizar nuevamente la toma de datos para corroborar lo
obtenido.

CONCLUSIONES

 Se concluye que cada medidor de caudal, en su ecuación de patrona miento tendrá un

coeficiente K característico y que, por esto, circunstancialmente se tendrá una ecuación de
patrona miento de caudal característico para cada uno de los medidores.

29
Laboratorio de Mecánica de Fluidos: Estudio y patronamiento de medidores de caudal en conductos a presión.

 Concluimos que las velocidades teóricas que se calculan son mayores a las velocidades
reales que se tienen, pues en el cálculo de velocidades teóricas se trabaja con un fluido
ideal y en condiciones ideales, mientras que en la realidad existen factores que generan
perdida de velocidades.

 Se concluye que al comparar los coeficientes de descarga del Venturi y del diafragma, se
nota que el diafragma posee un menor coeficiente de descarga en lo que respecta a lo
teórico, mientras que experimentalmente, en este caso fue mucho mayor a comparación
con el del Venturi.

 Se concluye que, al medir caudales por el método volumétrico con tiempo, presenta
variaciones con el obtenido con el rotámetro, debido a que existen variaciones debido a la
eficiencia y precisión del aforo, mas que todo a lo que respecta a la obtención de los
tiempos necesarios.

30