Práctica 3. Medición de Flujo en Ductos
Práctica 3. Medición de Flujo en Ductos
Práctica 3. Medición de Flujo en Ductos
INTRODUCCIÓN
En nuestra cotidianidad, para el transporte de fluidos o flujos de un lugar a otro es común emplear
ductos que faciliten dicho transporte, por ejemplo en una casa, en una industria, en una ciudad el
agua es transportada y distribuida por medio de redes de tuberías.
Ahora bien, el flujo de fluidos se puede clasificar como externo e interno, dependiendo de la
manera en la que el fluido sea obligado o “puesto”, por decirlo así, a fluir. Esto es, sobre una
superficie o dentro de un ducto. Cabe mencionar que es importante realizar dicha clasificación
puesto que los flujos internos y externos tienen características diferentes.
En este reporte se considera tanto el flujo interno o en ductos, así como el flujo en canales abierto,
por lo que es importante entonces aclarar que diferencia caracteriza a estos dos tipos de flujos.
El flujo interno, se considera como el flujo en donde el conducto está totalmente lleno con el
fluido, el cual se impulsa principalmente mediante diferencia de presión; mientras que en el caso
de flujos en canales abiertos, el conducto está parcialmente lleno con el fluido y por lo tanto está
parcialmente acotado por superficies sólidas, como una zanja de irrigación o cauce y el flujo se
impulsa sólo con la gravedad, en otras palabras, cuando hablamos de flujo en canales abiertos
implica que el flujo en el canal está abierto a la atmósfera, también es importante mencionar que
un flujo en canal abierto implica sólo líquido (usualmente agua o agua residual) expuesto a un gas
(por lo general aire, el cual se encuentra a la presión atmosférica).
El tema de medición de flujos por ductos tanto por canales abiertos gana relevancia en lo que
respecta al quehacer del ingeniero químico, ya que el flujo de un líquido o de un gas a través de
tuberías o ductos se usa comúnmente en sistemas de calefacción y enfriamiento y en redes de
distribución de fluido poniendo particular atención a la fricción, que se relaciona directamente con
la caída de presión y las pérdidas de carga durante el flujo a través de tuberías y ductos, contexto
con el cual el ingeniero se ve totalmente familiarizado. En lo que respecta a los sistemas de flujo
en canal abierto hechos por el hombre incluyen sistemas de irrigación, alcantarillas, desagües y
cunetas, y el diseño de estos sistemas es un área importante de aplicación de la ingeniería.
Ahora bien, comprendiendo estos temas, se procede a reportar los modelos y métodos
implementados para el desarrollo de un modelo lineal para estimar el coeficiente de velocidad de
una placa de orificio, y el tubo Pitot, incluyendo el coeficiente de correlación y los intervalos de
confianza para estos coeficientes, calculados a partir del intervalo de confianza de la pendiente, así
como la calibración del vertedero en diferentes formas, utilizando regresión lineal, reportando
coeficiente de correlación de la regresión para cada ajuste, el intervalo de confianza del
coeficiente Cd, a partir del intervalo de confianza para la pendiente encontrado. Finalmente se
discute y analiza sobre los resultados obtenidos.
MATERIALES Y MÉTODOS:
En las prácticas de laboratorio de Medición de fujo en ductos y Medición de flujo en canales
abiertos, se utilizaron los siguientes insumos y materiales:
• Agua
• Tubo Pitot
• Medidor Venturi
• Anemómetro
• Vertedero triangular
Donde se presenta para cada uno de ellos una descripción clara y detallada, de las principales
propiedades fisicoquímicas de la sustancia empleada, como viscosidad, densidad, presión, así
como también temperatura crítica, peso molecular, etc.
Agua:
El agua es un compuesto que se forma a partir de la unión, mediante enlaces covalentes, de dos
átomos de hidrógeno y uno de oxígeno; quedando así su fórmula molecular H2O y se trata de una
molécula muy estable. En la estructura de la molécula los dos átomos de hidrógeno y el de oxígeno
están dispuestos en un ángulo de 105°, lo cual le confiere características relevantes.
Adicionalmente el agua, es una molécula dipolar en la que el átomo de oxígeno central comparte
un par de electrones con cada uno de los dos átomos de hidrógeno con un exceso de carga
negativa junto al oxígeno, compensada por otra positiva repartida entre los dos átomos de
hidrógeno (A.C., 2010).
El agua presenta otras propiedades, como su facilidad para disolver una gran cantidad de
sustancias, su fuerte tendencia a empapar las paredes de los tubos capilares y ascender por ellos.
Además, tiene una capacidad calorífica alta. Su calor específico, por definición, es de una caloría
por gramo (Sánchez, y otros, 2016).
Medidor Venturi:
El medidor Venturi pertenece al grupo de los medidores deprimógenos los cuales consisten
básicamente en una reducción gradual o brusca de la sección transversal de la tubería,
ocasionando con esto un aumento de velocidad y una disminución de la presión en el fluido.
Además, es uno de los dispositivos más precisos para medir el gasto en tuberías y tiene la
desventaja de tener un costo elevado. Adicionalmente, causa una muy baja pérdida de carga y,
con las precauciones debidas, se puede usar para líquidos con determinadas concentraciones de
sólidos. El tubo Venturi se compone de tres secciones.
1. Entrada
2. Garganta
3. Salida
Donde, la sección de entrada tiene un diámetro inicial igual al diámetro de la tubería y una sección
cónica convergente que termina con un diámetro igual al de la garganta; y la salida consiste en una
sección cónica divergente que concluye con el diámetro de la tubería.
En adición, es esencial que el flujo entrando al Venturi sea uniforme, por lo que, es necesario un
largo tramo continuo con tubería recta aguas arriba y aguas abajo de la localización del Venturi
para mejorar la precisión en la medición del flujo. El tubo Venturi clásico puede construirse de
cualquier material, incluso de plástico, el material más empleado es acero al carbono, también se
utiliza el latón, bronce, acero inoxidable, cemento, y revestimientos de elastómeros para paliar los
efectos de la corrosión (Rivas Acosta, 2001).
El tubo Venturi ofrece varias ventajas con respecto a otros medidores, como son:
• Menor pérdida de carga permanente, que la producida por del diafragma y la tobera de
flujo, gracias a los conos de entrada y salida.
• Medición de caudales superiores a un 60% a los obtenidos por el diafragma para la misma
presión diferencial e igual diámetro de tubería.
• El Venturi requiere un tramo recto de entrada más corto que otros elementos primarios.
• Facilidad para la medición de flujo de líquidos con sólidos en suspensión
El tubo Venturi se recomienda en casos donde el flujo es grande y que se requiera una baja caída
de presión, o bien, el fluido sea altamente viscoso, se utiliza donde se requiera el máximo de
exactitud, en la medición de fluidos altamente viscosos, y cuando se necesite una mínima caída de
presión permanente (Rivas Acosta, 2001).
Pitot compuesto:
El tubo Pitot es un instrumento destinado a la medición del caudal a través de la cuantificación de
la velocidad del flujo. Es un medidor indirecto de caudal, y puede utilizarse para la medida de
caudales de gas en grandes conducciones, como chimeneas de industrias pesadas, generando una
pequeña diferencia de presión.
Este tubo contiene dos orificios, los cuales están ubicados diametralmente opuestos. En uno de
ellos, el que se encuentra perpendicular y en dirección opuesta al flujo de agua, se mide la presión
total y en el otro, el que se encuentra perpendicular y en la misma dirección del flujo de agua, se
mide la presión estática. El elemento secundario, el cual puede ser una columna de agua o bien
transductores de presión, debe conectarse a dichos orificios, para obtener la diferencia entre la
presión total y la estática, es decir, la presión dinámica (Figueroa Ibarra, Rodríguez Martínez, &
Santabárbara Botello, 2007).
Anemómetro:
Un anemómetro es un instrumento usado para medir la velocidad o rapidez de los gases ya sea en
un flujo contenido, como el flujo de aire en un conducto, o en flujos no confinados, como un
viento atmosférico. Un anemómetro detecta el cambio en alguna propiedad física del fluido o el
efecto del fluido en un dispositivo mecánico insertado en el flujo para así determinar la velocidad
de dicho fluido (company, 2020). Además, dependiendo del modelo es también posible medir no
solo la velocidad del aire sino también el caudal del aire, dirección del viento, temperatura, altitud
y presión atmosférica.
El tipo de anemómetro usado fue el anemómetro de rotación, el cual está dotado de cazoletas
(Robinson) o hélices unidas a un eje central cuyo giro, proporcional a la velocidad del viento, es
registrado convenientemente; en los anemómetros magnéticos, dicho giro activa un diminuto
generador eléctrico que facilita una medida precisa (Drogallega, 2011).
Para una mayor claridad del marco matemático, se procede a explicar clara y detalladamente cada
una de las ecuaciones utilizadas tanto para la práctica #3 como para la práctica #4.
En este sentido para la práctica #3. Medición de flujos en ductos se implementó el siguiente
modelo matemático.
Parte a. Tubo Pitot. Determinación del coeficiente de velocidad (Cp)
1. Modelo lineal para estimar el coeficiente de velocidad
- Corrección de presiones
√ Po−Pc (2)
Luego, con los volúmenes (datos) y con raíz de P o – Pc, se gráfica V vs √ (P o-Pc) con la finalidad de
obtener la pendiente de la línea recta.
2. Cálculo de la densidad del aire
Otro cálculo importante para determinar el Cp, es el cálculo de la densidad del aire, la cual tiene
una relación directa con éste.
P M waire
ρ= (3)
RT
Donde ρ: densidad (kg/m3)
P: Presión (kPa)
Mw aire: Peso molecular (kg/kmol)
R: Constante universal de los gases ideales (kPa*m 3/kmol*K)
T: Temperatura (K)
y N 2=( 1− y H 2 O ) 0,78(5)
y O 2=( 1− y H 2 O ) 0,21(6)
y Ar= ( 1− y H 2 O ) 0,01(7)
Con las condiciones de trabajo y con la ecuación (8) se determina la composición del agua
V
Q= (10)
t
Siendo V: Volumen (m3)
t: tiempo (s)
V = A∗h (11)
Para la determinación del área
π
Atanque = ∗d 2 (12)
4
Con los caudales determinados así como con las raíces de las diferencias de presiones, se procede
a graficar Q vs √ (P2-P1) con el fin de determinar la pendiente de la línea recta
2. Determinación de Co
Utilizando y manipulando la ecuación (13), se puede despejar Co (ecuación 14).
C0 A2 2
m=
√ 1−β 4
√ ρ
(13)
m √ 1−β 4
=C 0 (14)
2
A2
ρ √
Ahora bien, es necesario el cálculo de
D2
β= (15)
D1
Siendo D2: do, diámetro del orificio.
D1: D, Diámetro de la garganta.
Es importante mencionar que para el cálculo del área se emplea igual que para el caso interior, la
ecuación (12) teniendo cuidado en el diámetro a emplear.
A continuación, se presentan las ecuaciones que se utilizaron para hacer el cálculo del coeficiente
de velocidad de un tubo Venturi, y el tubo Pitot, así como, el C p y el Cv de cada medidor
respectivamente, incluyendo el coeficiente de correlación y los intervalos de confianza para estos
coeficientes, calculados a partir del intervalo de confianza de la pendiente.
Tubo Pitot: De acuerdo a los datos medidos en la práctica de velocidad de la corriente de aire
y de la presión dinámica, además, teniendo en cuenta el factor del manómetro inclinado que
es 0.05, se calculó la raíz del diferencial de presión para así poder realizar una regresión lineal
con los datos obtenidos
Vertedero triangular:
Una vez usado el método balde-reloj se obtuvieron los datos de altura, volumen y tiempo,
permitiendo así hallar el caudal para cada dato y realizar el diagrama Q vs h 5/2. Posterior a esto fue
necesario calcular cada uno de los volúmenes con la masa del balde y la densidad del agua así:
V =m∗ρ
Donde,
Así, con cada volumen calculado se pudo encontrar el caudal para cada uno de los datos.
V
Q=
t
Donde,
Una vez realizada la gráfica Q vs h5/2 se obtiene la pendiente con el modelo de regresión lineal ya
que la finalidad es hallar el Cd (estimado) del vertedero triangular, este se encuentra usando la
ecuación ().
8 θ
m=Cd
15
√ 2 g tan
2()
Donde,
σ2 σ2
^β 1−t α
2
∗
, ( n−2 ) √ S xx
< β 1 < ^β 1+ t α
2 √
, ( n−2)
∗
S xx
1 x́ 2 1 x́ 2
^β 0−t α
2
∗ σ∗
, ( n−2 )
√ 2
( +
n S xx )
< β 0 < ^β 0+ t α
2 √
, ( n−2) (2
∗ σ∗ +
n S xx )
Finalmente se calculó el porcentaje de error absoluto para el Cd con la ecuación ().
ANÁLISIS Y DISCUSIÓN
Para ductos cerrados se debe calibrar un tubo pitot y una placa orificio, teniendo en cuenta cada
una de las mediciones de las variables que se deben de tomar, como se expuso en la metodología;
Con la finalidad de hallar el Cp y C0 de cada medidor de flujo.
Para vertederos abiertos se debe calibrar un vertedero triangular, teniendo en cuenta las variables
que se deben medir según la metodología para hallar el Cd de este.
1. Ductos Cerrados
Tubo pitot:
Esta calibración se realizó con un túnel de viento, en donde se midió la velocidad de la corriente de
aire con ayuda de un anemómetro y se midió la presión dinámica con el instrumento a calibrar,
además se tuvo en cuenta el factor del manómetro inclinado que es 0.05, dicho factor se
multiplicó por la presión dinámica leída en el dispositivo. Luego se halló la raíz del diferencial de
presión para así poder realizar una regresión lineal con los datos obtenidos con la finalidad de
hallar el Cp (estimado) del tubo pitot, los datos obtenidos anteriormente se tabulan en la Tabla 1.
Por consiguiente, se procede a graficar V vs (Po-Pc)1/2 puesto que con base al modelo lineal
propuesto y al modelo matemático se tiene la ecuación (), la cual nos indica que la pendiente es:
Dato v (m/s)
Calibración
(Po-Pc) kPa
ρ √
m=Cp de tubo pitot
(Po-Pc) Pa (Po-Pc)1/2
Así, se obtiene 1 3.5 0.007 7 2.646 la gráfica 1.
2 2.9 0.007 6.5 2.550
3 2.5 0.006 5.5 2.345
4 1.9 0.004 4 2.000
5 0.8 0.003 2.5 1.581
6 0.2 0.002 2 1.414
2.0
R² = 0.99
1.5
1.0
0.5
0.0
1.2 1.4 1.6 1.8 2.0 2.2 2.4 2.6 2.8
(Po-Pc)1/2
Condiciones de trabajo
Humedad del aire 0.588
T aire (K) 296.85
P (kPa) 84.694
Es importante destacar que con el valor de la humedad relativa se calcula la fracción del agua, para
así poder calcular fracciones de las demás sustancias, tales como: Nitrógeno, Oxigeno y Argón. Las
ecuaciones implementadas fueron las siguientes:
M waireSeco =∑ x i M w i
PM w Aire
ρ=
RT
Tabla 3. Datos para el cálculo de la densidad del aire.
Luego, ya es posible hallar el valor del estimado de Cp, implementando la ecuación (). Este tiene
un valor de 1.7227. Para tener una mejor aproximación a el valor de la pendiente se procedió a
identificar los intervalos de confianza con ayuda del programa Statgraphics, para así identificar
cuál sería el límite superior e inferior de la pendiente, intercepto y del Cp. Donde se obtuvo los
siguientes resultados tabulados en las tablas 4 y 5.
Tabla 4. Intervalos de confianza del 95% para los parámetros del modelo ajustado.
Intervalos de confianza del 95 % para los estimados de los coeficientes
Error
Parámetro Estimado Estándar Límite inferior Límite Superior
Constante -3.15647 0.320549 -4.04646 -2.26648
Pendiente 2.45204 0.149723 2.0363 2.86774
Seguidamente se obtuvo del programa Statgraphics la gráfica del modelo ajustado, incorporando
así el límite superior e inferior.
Finalmente, por medio de la gráfica 1 se evidencia que el coeficiente de correlación indica que el
modelo explica el 98.53% de la variabilidad en la velocidad, teniendo en cuenta que este valor es
cercano a 100% se puede afirmar que los puntos están próximos a la línea de tendencia, dado que
teóricamente esto muestran que las variables están fuertemente asociadas, pero como se observa
hay puntos que se alejan del modelo lineal, evidenciando así desviaciones entre los datos tomados
y afectando de manera significativa este coeficiente. Así mismo, se puede afirmar que el valor de
Cp estimado es un valor alejado al teórico que es 1[ CITATION Cen20 \l 9226 ], con un porcentaje
de error absoluto de 78.459%, y esto se pudo deber a diversos errores, tales como: errores
humanos, como por ejemplo la posición del anemómetro y la del tubo Pitot, también pudo haber
errores de lectura en la presión en el manómetro de tubo inclinado y falta de calibración en los
instrumentos implementados. Aunque, es importante destacar que por medio de la gráfica 2,
evidenciamos que los datos obtenidos se encuentran dentro del intervalo de confianza.
Por consiguiente, se procede a calcular el volumen del tanque teniendo en cuenta los datos de
este, tabulados en la tabla 7 y las ecuaciones () y ().
π 2
Atanque = D
4
V tanque = Atanque h
Se procede a calcular el caudal de cada dato con la ecuación () y los datos obtenidos se encuentran
en la tabla 8.
V tanque
Q=
t
Tabla 8. Caudales para la placa de orificio.
Caudales
Dato Q (m3/s)
Se grafica Q vs (P1-P2)1/2 puesto 1 0.00115 que con base al modelo lineal
propuesto y al modelo 2 0.00107 matemático se tiene la
ecuación (), la cual nos indica 3 0.00098 que la pendiente es
4 0.00088
C0 A2 2
m=
√ 1−β 4
√ ρ
5
6
7
0.00074
0.00061
0.00041
Así, se obtiene la gráfica 3.
R² = 0.95
0.0008
0.0006
0.0004
0.0002
0.0000
30 40 50 60 70 80 90 100 110
(P2-P1)1/2
Gráfico 3. Caudal Vs (P1-P2)1/2
Donde se obtuvo que le valor de la pendiente es: 1E-05, para hallar el C0 de la placa de orificio es
necesario calcular el valor de y el área de la vena contracta (A 2), que se calcula con la ecuación ()
y (), respectivamente, además, se tiene en cuenta los datos tabulados en la tabla 9.
D2
β=
D1
π 2
A2= D
4
Tabla 9. Datos para calcular C0, y A2
Cálculo de Co
D1 (m) 0.03145
D2 (m) 0.0206
3
H2O (kg/m ) 1000
Donde, se obtuvieron estas variables con un valor de = 0.1841 y A2= 3E-04. Luego, con la
ecuación () se pudo calcular el C 0 estimado de la placa de orificio, teniendo un valor de 0.606. Para
tener una mejor aproximación a el valor de la pendiente se procedió a identificar los intervalos de
confianza con ayuda del programa Statgraphics, para así identificar cuál sería el límite superior e
inferior de la pendiente, intercepto y del C 0. Donde se obtuvo los siguientes resultados tabulados
en las tablas 10 y 11.
Tabla 10. Intervalos de confianza del 95% para los parámetros del modelo ajustado.
Seguidamente se obtuvo del programa Statgraphics la gráfica del modelo ajustado, incorporando
así el límite superior e inferior.
A partir de los intervalos de confianza para los dos medidores de flujo, es posible afirmar que el
valor de C0 tiene un grado alto de confiabilidad, puesto que los intervalos son estrechos lo que
hace que el margen de error en la toma de datos sea menor. Lo anterior puede corroborarse a
partir del coeficiente de correlación que se muestra en la gráfica de caudal vs el diferencial de
presión.
2. Canales abiertos:
Vertedero triangular: Esta calibración se efectuó empleando la técnica de balde- reloj para así
hallar el caudal, por ende, se varió el flujo y se midió la altura en el vertedero. Así, se
implementó un balde donde se vertió el agua recaudada y luego se pesó con una balanza; con
la densidad del agua se pudo calcular el volumen de agua recogida. En este orden de ideas,
con el volumen y el tiempo se puede hallar el caudal para cada dato. Se determinó la altura 5/2
para así poder realizar una regresión lineal con los datos obtenidos con la finalidad de hallar el
Cd (estimado) del vertedero triangular. Los datos teóricos se tabulan en la tabla 13 y los datos
obtenidos en la tabla 14.
Dat mH2O+balde(g) mH2O(g) VH2O (mL) VH2O (m3) t (s) h (m) h5/2 (m5/2)
o
30.1 0.01
1 1091 1005 1007.397606 0.001007398 3.23817E-05
7 6
0.02
2 1105 1019 1021.431006 0.001021431 12.1 8.92335E-05
4
0.03
3 1145 1059 1061.526433 0.001061526 6.61 0.000213156
4
0.05
4 3961 3480 3488.302159 0.003488302 8.44 0.000587389
1
0.05
5 5146 4665 4676.129187 0.004676129 8.11 0.000845532
9
0.06
6 5033 4552 4562.859606 0.00456286 5.98 0.001161949
7
0.07
7 3090 2609 2615.224234 0.002615224 2.42 0.001592333
6
Por consiguiente, se procede a calcular el caudal para cada dato con la ecuación () y se obtienen
los datos de la tabla 15.
V agua
Q=
t
Tabla 15. Caudales de agua para cada dato.
Dato Q (m3/s)
1 3.3E-05
2 8.4E-05
3 0.00016
4 0.00041
5 0.00058
6 0.00076
7 0.00108
Se grafica Q vs h5/2 puesto que con base al modelo lineal propuesto y al modelo matemático se
tiene la ecuación (), la cual nos indica que la pendiente es
8 θ
m=Cd
15 ()
√ 2 g tan
2
8 θ
Q=Cd
15 ()
√ 2 g tan H 5 /2
2
Se pudo calcular el Cd estimado de vertederos triangulares, teniendo un valor de 0.8334. Para
tener una mejor aproximación a el valor de la pendiente se procedió a identificar los intervalos de
confianza con ayuda del programa Statgraphics, para así identificar cual sería el límite superior e
inferior de la pendiente, intercepto y del Cd. Donde se obtuvo los siguientes resultados tabulados
en las tablas 16 y 17.
Tabla 16. Intervalos de confianza del 95% para los parámetros del modelo ajustado.
Seguidamente se obtuvo del programa Statgraphics la gráfica del modelo ajustado, incorporando
así el límite superior e inferior.
El valor de Cd tiene un grado alto de confiabilidad, puesto que los intervalos son estrechos lo que
hace que el margen de error en la toma de datos sea menor. Lo anterior puede corroborarse a
partir del coeficiente de correlación que se muestra en la gráfica de caudal vs la altura a la 5/2.
Conclusiones:
Referencias
Cenam. (10 de 03 de 2020). Calibración de tubos pitot tipo S en el labinthap. Obtenido de
http://www.cenam.mx/memorias/descarga/Memorias%20Simposio/documentos/ta-
or061.pdf