estudio de pérdidas de carga en tuberías

LABORATORIO DE TERMOFLUIDOSPRÁCTICA 4: Estudio de pérdidas de carga en tuberías

Pág. 1

INTRODUCCIÓN TEÓRICA.INTRODUCCIÓN TEÓRICA.INTRODUCCIÓN TEÓRICA.INTRODUCCIÓN TEÓRICA.----

Las pérdidas de carga en las tuberías son de dos clases: primarias y secundariasdos clases: primarias y secundariasdos clases: primarias y secundariasdos clases: primarias y secundarias. Las pérdidas primarias se definen como las pérdidas de superficie en el contacto del fluido con la tubería, rozamiento de unas capas del fluido con otras (régimen laminar) o de las partículas del fluido entre sí (régimen turbulento). Tienen lugar en flujo uniforme, por lo que principalmente suceden en los tramos de tubería de sección constante.

Las pérdidas secundarias o locales se definen como las pérdidas de forma, que tienen lugar en las transiciones (estrechamientos o expansiones de la corriente), codos, válvulas y en toda clase de accesorios de tubería.

A continuación estudiamos ambos tipos de pérdidas:

I.- Pérdidas Primarias:

Supongamos una tubería horizontal de diámetro constante por la que circula un fluido cualquiera. Aplicando la ecuación de Bernouilli entre dos puntos 1 y 2:

P1/ρ⋅g+z1+v1

2/2⋅g= P2/ρ⋅g+z2+v2

2/2⋅g+∆h,

donde ∆h representa las pérdidas primarias entre 1 y 2.

Existen muchas ecuaciones para calcular estas pérdidas. Una de ellas es la ecuación de Darcy-Weisbach, que se desarrolló para tuberías rellenas de agua con un diámetro constante:

∆h=f⋅L⋅v2/(2⋅g⋅D),

donde f es el coeficiente de fricción, L la longitud de la tubería, D ó ∅ el diámetro de la tubería y v la velocidad media del fluido.

El coeficiente f es adimensional, y depende de la velocidad (v), del diámetro (D), de la densidad (ρ), de la viscosidad (µ) y de la rugosidad (ε).

Es decir:

f=h(v, D, ρ, µ, ε)

Mediante análisis dimensional obtenemos:

f=h(v⋅D⋅ρ/µ, ε/D)

Al primer término de la relación anterior se le conoce como número de Reynolds :

Re=v⋅D⋅ρ/µ

El segundo término se denomina rugosidad relativa. Ambos juegan un papel fundamental en el cálculo de las pérdidas de carga primarias, puesto que la f se calcula mediante estos coeficientes en el “diagrama de Moody”.

Este diagrama es un ábaco que permite calcular el coeficiente de fricción conociendo la rugosidad relativa y el nº de Reynolds.

El coeficiente de fricción (f) puede calcularse mediante un amplio grupo de ecuaciones, aparte de la aplicación del “diagrama de Moody”. Muchas de estas funciones sirvieron incluso para dibujar el diagrama.

En esta práctica se emplean dos de estas ecuaciones:

1.- Ecuación de Poiseuille. Aplicable en fluidos bajo régimen laminar en tuberías rugosas o lisas, puesto que en dicho régimen el coeficiente de fricción no es función de la rugosidad relativa.

f=64/Re


Pág. 2

2.- Ecuación de Blasius. Aplicable en fluidos bajo régimen turbulento y con Re<100000. La tubería ha de ser lisa. (rugosidad ε=0).

f=0.316⋅Re

0.25

NOTA: Generalmente el coeficiente de fricción (f) se calcula mediante “diagrama de Moody”.

II.- Pérdidas Secundarias.

En este caso se aplica la ecuación de Bernouilli entre dos puntos entre los cuales existen distintos accesorios de tubería.

El factor ∆h se dividirá entonces en dos: hf (pérdidas primarias) y he (pérdidas secundarias), ocasionadas por los accesiorios de las tuberías.

� Cálculo de he. Aplicamos la ecuación:

he=K⋅v1

2/2⋅g,

donde v1 es la velocidad antes del accesorio y K es un coeficiente determinado experimentalmente. Este coeficiente es necesario excepto en el caso debido a una expansión brusca de la tubería.

En este caso:

he=v1

2/2⋅g ,

siempre que el diámetro de la tubería sea despreciable frente al ensanchamiento de la misma.

Las pérdidas menores también pueden expresarse en términos de longitud equivalente, que es la longitud de tubo que haría falta para ocasionar una pérdida de carga similar a la que ocasiona el accesorio de la tubería.

� Cálculo de la longitud equivalente.

f(Le/D)⋅(v2/2⋅g)=K⋅v2/2⋅g,

donde K puede referirse a una sola pérdida o a la suma de varias pérdidas. Al despejar llegamos a la expresión definitiva de la longitud equivalente:

Le=K⋅D/f


Pág. 3


Pág. 4

PROCEDIMIENTO EXPERIMENTAL.PROCEDIMIENTO EXPERIMENTAL.PROCEDIMIENTO EXPERIMENTAL.PROCEDIMIENTO EXPERIMENTAL.----

Para poner en funcionamiento el equipo se abren la válvula de flujo y la válvula de control para permitir que circule el agua por el circuito. Una vez que el aire existente en el interior del mismo ha sido expulsado, se conecta la válvula antirretorno y se presuriza el sistema.

A continuación tomamos las lecturas del manómetro de agua, y se mide el caudal mediante una probeta. Este proceso se realiza para distintas posiciones de la válvula antirretorno.

Los datos obtenidos se muestran en la tabla 1.

TABLA 1:

Lecturas Manómetro de H2O (mm.c.a.)

MEDIDA Volumen (l) Tiempo (s) Caudal (m3/s)

h1 h2

Pérdida de carga ∆h (m.c.a.)

1 0.08 30 2.67*10-6 284 248 362 0.10 30 3.33*10-6 283 243 403 0.23 60 3.83*10-6 141 92 494 0.13 30 4.33*10-6 287 233 545 0.15 30 5.00*10-6 292 224 686 0.31 60 5.17*10-6 151 70 817 0.16 30 5.33*10-6 297 220 778 0.17 30 5.67*10-6 295 221 749 0.19 30 6.33*10-6 301 201 10010 0.39 60 6.50*10-6 158 49 10911 0.21 30 7.00*10-6 310 189 12112 0.24 30 8.00*10-6 313 185 12813 0.48 60 8.00*10-6 174 20 15414 0.25 30 8.33*10-6 314 162 15215 0.51 60 8.50*10-6 327 145 18216 0.26 30 8.67*10-6 317 178 13917 0.27 30 9.00*10-6 321 167 15418 0.56 60 9.33*10-6 340 131 20919 0.30 30 1.00*10-5 336 142 19420 0.60 60 1.00*10-5 354 107 24721 0.32 30 1.07*10-5 364 88 27622 0.33 30 1.10*10-5 369 78 291

((((NOTA:NOTA:NOTA:NOTA: LAS MEDIDAS SE HAN REALIZADO DE TAL MANERA QUE ESTAN ORDENA LAS MEDIDAS SE HAN REALIZADO DE TAL MANERA QUE ESTAN ORDENA LAS MEDIDAS SE HAN REALIZADO DE TAL MANERA QUE ESTAN ORDENA LAS MEDIDAS SE HAN REALIZADO DE TAL MANERA QUE ESTAN ORDENADAS DE DAS DE DAS DE DAS DE MENOR A MAYOR VALOR DE CAUDAL)MENOR A MAYOR VALOR DE CAUDAL)MENOR A MAYOR VALOR DE CAUDAL)MENOR A MAYOR VALOR DE CAUDAL)


Pág. 5

Con las anteriores medidas se ha realizado la siguiente tabla:

TABLA 2:

Medida Q (m3/s) Velocidad (m/s) V2 (m2/s2) ∆h (m.c.a.) f Re Log(V) Log(∆h) Log(f) Log(Re)1 2.67*10-6 0.378 0.143 36 0.0285 1128 -0.423 1.556 -1.545 3.0522 3.33*10-6 0.471 0.222 40 0.0204 1406 -0.327 1.602 -1.690 3.1483 3.83*10-6 0.542 0.294 49 0.0188 1618 -0.267 1.690 -1.726 3.2094 4.33*10-6 0.613 0.376 54 0.0162 1830 -0.213 1.732 -1.790 3.2625 5.00*10-6 0.707 0.500 68 0.0154 2110 -0.150 1.833 -1.812 3.3246 5.17*10-6 0.731 0.535 81 0.0171 2182 -0.136 1.908 -1.767 3.3397 5.33*10-6 0.754 0.569 77 0.0153 2251 -0.123 1.886 -1.815 3.3528 5.67*10-6 0.802 0.643 74 0.0130 2394 -0.096 1.869 -1.886 3.3799 6.33*10-6 0.896 0.802 100 0.0141 2675 -0.048 2.000 -1.851 3.42710 6.50*10-6 0.920 0.846 109 0.0145 2746 -0.036 2.037 -1.839 3.43911 7.00*10-6 0.990 0.981 121 0.0139 2955 -0.004 2.083 -1.857 3.47112 8.00*10-6 1.132 1.281 128 0.0113 3379 0.054 2.107 -1.947 3.52913 8.00*10-6 1.132 1.281 154 0.0135 3379 0.054 2.188 -1.866 3.52914 8.33*10-6 1.178 1.389 152 0.0124 3516 0.071 2.182 -1.907 3.54615 8.50*10-6 1.203 1.446 182 0.0142 3591 0.080 2.260 -1.848 3.55516 8.67*10-6 1.227 1.504 139 0.0105 3663 0.089 2.143 -1.979 3.56317 9.00*10-6 1.273 1.621 154 0.0107 3800 0.105 2.188 -1.971 3.58018 9.33*10-6 1.320 1.742 209 0.0136 3940 0.121 2.320 -1.866 3.59519 1.00*10-5 1.415 2.001 194 0.0109 4224 0.151 2.288 -1.963 3.62620 1.00*10-5 1.415 2.001 247 0.0140 4224 0.151 2.393 -1.855 3.62621 1.07*10-5 1.514 2.291 276 0.0136 4519 0.180 2.441 -1.866 3.65522 1.10*10-5 1.556 2.422 291 0.0136 4645 0.192 2.464 -1.867 3.667

Se ha decidido no despreciar ninguna de las medidas por diversas razones:

1. Consideramos que, a pesar de que individualmente existan medidas que no se ajustan a las curvas teóricas (en algunas medidas la pérdida de carga disminuye al aumentar la velocidad), en su conjunto sí lo hacen.


Pág. 7

2.- Dado el gran error de precisión del aparato (seguramente producido porque éste no se encontraba en equilibrio, como ejemplo ver las parejas de medidas 12-13 y 19-20, donde a igualdad de velocidad se registran pérdidas de carga muy distintas), no puede asegurarse que ninguna medida sea errónea, por lo cual ninguna es despreciable.

CUESTIONES.CUESTIONES.CUESTIONES.CUESTIONES.

1.1.1.1. Realizar, a escala conveniente, los siguientes gráficos:Realizar, a escala conveniente, los siguientes gráficos:Realizar, a escala conveniente, los siguientes gráficos:Realizar, a escala conveniente, los siguientes gráficos:

Gráfico 1: log(Gráfico 1: log(Gráfico 1: log(Gráfico 1: log(∆∆∆∆h) en función de log(v)h) en función de log(v)h) en función de log(v)h) en función de log(v)Gráfico 2: log(f) en función de log(Re)Gráfico 2: log(f) en función de log(Re)Gráfico 2: log(f) en función de log(Re)Gráfico 2: log(f) en función de log(Re)

Con ayuda de estos gráficos tratar de determinar el valor de la velocidad vc por debajo de la cual Con ayuda de estos gráficos tratar de determinar el valor de la velocidad vc por debajo de la cual Con ayuda de estos gráficos tratar de determinar el valor de la velocidad vc por debajo de la cual Con ayuda de estos gráficos tratar de determinar el valor de la velocidad vc por debajo de la cual el el el el régimen es laminar.régimen es laminar.régimen es laminar.régimen es laminar.

De los gráficos (se presentan al final de la cuestión), extraemos las siguientes conclusiones:

GRÁFICO 1: Se sabe que, mientras el régimen sea laminar, la pérdida de carga es proporcional a la primera potencia de la velocidad. En régimen turbulento, en cambio, la pérdida de carga es proporcional a la segunda potencia de la velocidad. Es decir:

� Régimen laminar:

∆h=k*v

� Régimen turbulento:

∆h=k*v2

Tomando logarítmos en estas expresiones:

� Régimen laminar:

log(∆h)=log(k)+log(v)

� Régimen turbulento:

log(∆h)=log(k)+2log(v)

Es decir, en el gráfico 1 los datos se ajustarán a una recta de orden 1 en régimen laminar, y de orden 2 en régimen turbulento.

GRÁFICO 2: Gracias al diagrama de Moody sabemos que, mientras el régimen sea laminar, el coeficiente de fricción en función del número de Reynolds se ajusta a una recta de la forma f=64/Re. Por ello, el logaritmo del coeficiente de fricción en función del logaritmo del número de Reynolds se ajustará a una recta de pendiente –1.

En cambio, en régimen turbulento, con número de Reynolds menor de 10000 y para un tubo liso como el nuestro, el coeficiente de fricción en función del log(Re ) se ajusta a una curva de la forma f≈0.316/Re

0.25. Por ello, el log del coeficiente de fricción en función de log(Re ) se ajustará a una recta de pendiente –0.25.

Por todo ello, la vc será aquella velocidad a partir de la cual los gráficos 1 y 2 se ajustan a rectas con pendiente mayor que 1 (en el caso del gráfico 1), y mayor que –1 (gráfico 2).

En el gráfico 1, los datos de la medida 5 se ajustan a la siguiente recta:

y=1.016x+1.963,


Pág. 8

con r=0.97810.

En el gráfico 2, en la misma medida, los datos se ajustan a la recta:

y=-0.98453x+1.473,

con r=-0.98453.

Si consideramos un mayor número de medidas, los datos se ajustan a rectas que ya no tienen la pendiente buscada, acercándose a las pendientes descritas para régimen turbulento. Por ello, consideramos que la vc es la velocidad correspondiente a la medida 5, es decir, vc=0.707 [m/s].

¿Coincide c¿Coincide c¿Coincide c¿Coincide con el valor teórico para Ron el valor teórico para Ron el valor teórico para Ron el valor teórico para Reeee<2300?<2300?<2300?<2300?

Calculamos el valor teórico para Re<2300 despejando la velocidad de Re.

v=Re⋅µ/(ρ⋅D) ⇒ vvvvcccc=0.7705=0.7705=0.7705=0.7705 [m/s]

donde: Re=2300; µ=1.005⋅10-3 [Ns/m2]; ρ=1000 [Kg/m3] y D=0.003 [m].

El valor de la medida no coincide con el valor teórico, aunque la desviación es mínima.

Gráfico 3: Con los valores vGráfico 3: Con los valores vGráfico 3: Con los valores vGráfico 3: Con los valores v≥≥≥≥vvvvcccc, dibujar el gráfico de , dibujar el gráfico de , dibujar el gráfico de , dibujar el gráfico de ∆∆∆∆h en función de v.h en función de v.h en función de v.h en función de v.Gráfico 4: Con los valores v Gráfico 4: Con los valores v Gráfico 4: Con los valores v Gráfico 4: Con los valores v ≤≤≤≤vvvvcccc, dibujar el gráfico de , dibujar el gráfico de , dibujar el gráfico de , dibujar el gráfico de ∆∆∆∆h en función de v. Unir a continuación los gráficos h en función de v. Unir a continuación los gráficos h en función de v. Unir a continuación los gráficos h en función de v. Unir a continuación los gráficos

3 y 4 en uno solo.3 y 4 en uno solo.3 y 4 en uno solo.3 y 4 en uno solo.

LOS GRÁFICOS SE ENCUENTRAN AL FINAL DE LA CUESTIÓN.

Del gráfico 1, determinar las relaciones empíricas de la forma Del gráfico 1, determinar las relaciones empíricas de la forma Del gráfico 1, determinar las relaciones empíricas de la forma Del gráfico 1, determinar las relaciones empíricas de la forma ∆∆∆∆h=Kh=Kh=Kh=K⋅⋅⋅⋅vvvvnnnn que pueden expresar las dos zonas que pueden expresar las dos zonas que pueden expresar las dos zonas que pueden expresar las dos zonas del gráfico.del gráfico.del gráfico.del gráfico.

Tomando log en la ecuación anterior, obtenemos:

log(∆h)=log(k)+n⋅log(v)

que es la expresión a la cual se ajusta el gráfico 1.

Como se ha dicho anteriormente, si el régimen es laminar, el gráfico 1 se ajusta a una recta de pendiente 1, y si el régimen es turbulento a una recta de pendiente 2. Tomando como régimen laminar las medidas de 1 a 5 y como régimen turbulento las medidas 6 a 22, obtenemos que el gráfico se ajusta a las siguientes rectas:

� Régimen laminar:log(∆h)=1.016⋅log((v)+1.936, con r=0.97810

� Régimen turbulento:log(∆h)=1.72⋅log((v)+2.080, con r=0.97626

No se ajusta a una recta de pendiente 2 debido a que en la zona de transición entre régimen laminar y régimen turbulento, los datos se ajustan a una recta con pendiente de valor entre 1 y 2.

Tomando exponenciales queda:

� Régimen laminar:∆h=86.30⋅v1.016

� Régimen turbulento:∆h=120.23⋅v1.721


Pág. 9

Del gráfico 2, buscar las relaciones empíricas de la forma f=KDel gráfico 2, buscar las relaciones empíricas de la forma f=KDel gráfico 2, buscar las relaciones empíricas de la forma f=KDel gráfico 2, buscar las relaciones empíricas de la forma f=K⋅⋅⋅⋅RRRReeee

nnnn para las dos zonas del gráfico. para las dos zonas del gráfico. para las dos zonas del gráfico. para las dos zonas del gráfico.

Tomando logaritmos en la ecuación anterior, obtenemos:

log(f)=n⋅log(Re)+log(K),

que es la expresión a la que se ajusta la gráfica 2.

Como se dijo anteriormente, si el régimen es laminar el gráfico 2 se ajusta a una recta de pendiente –1, y si el régimen es turbulento se ajusta a una recta de pendiente –0.25. Tomando como régimen laminar las medidas de 1 a 5 y como régimen turbulento las medidas 6 a 22, obtenemos que el gráfico se ajusta a las siguientes rectas:

� Régimen laminar:log(f)=-0.98453⋅Re+1.437, con r=-0.98453

� Régimen turbulento:log(f)=-0.2805⋅Re-0.8913, con r=0.51

Tomando exponenciales queda:

� Régimen laminar:f=27.35⋅Re

-0.98453

� Régimen turbulento:f=0.1284⋅Re

-0.2805

Del gráfico 3, determinar un valor medio de f para el régimen turbulento.Del gráfico 3, determinar un valor medio de f para el régimen turbulento.Del gráfico 3, determinar un valor medio de f para el régimen turbulento.Del gráfico 3, determinar un valor medio de f para el régimen turbulento.

Se sabe que la pérdida de carga en función de la velocidad en régimen turbulento se ajusta a una ecuación de la forma:

∆h=k⋅v2,

la cual es la ecuación de una parábola. Esta fórmula es la ecuación de Darcy-Weisbach, la cual se difine como:

∆h=f⋅L⋅v2/(2⋅g⋅D) ⇒ k= f⋅L/(2⋅g⋅D).

Tomando logarítmos en esta expresión queda:

log(∆h)=2⋅log(v)+log[f⋅L/(2⋅g⋅D)],

expression que puede ajustarse a una recta mediante regresión lineal. La ordenada en el origen de la recta será log[f⋅L/(2⋅g⋅D)].

La recta que obtenemos es:

log(∆h)=1.72⋅log(v)+2.080 ⇒ log[f⋅L/(2⋅g⋅D)]=2.080

Llegados a este punto y, sabiendo que L=0.52[m] y D=0.003[m]:

f=0.1136

Del gráfico 4, determinar el valor de Del gráfico 4, determinar el valor de Del gráfico 4, determinar el valor de Del gráfico 4, determinar el valor de µµµµpara el agua a la temperatura del ensayo.para el agua a la temperatura del ensayo.para el agua a la temperatura del ensayo.para el agua a la temperatura del ensayo.


Pág. 10

Se sabe que para el régimen laminar:

∆h=8⋅µ⋅L⋅v/(ρ⋅D2⋅g) ⇒ ∆h=k⋅v, donde k=8⋅µ⋅L/(ρ⋅D2⋅g)

Esta ecuación puede ajustarse a una regresión lineal con los datos de la gráfica 4. La recta obtenida me diante la regresión es:

∆h=0.0975⋅v-3.4715⋅10-3 ⇒ k=8⋅µ⋅L/(ρ⋅D2⋅g)=0.0975,

siendo L=0.52[m], ρ=1000 [Kg/m3] y D=0.003 [m].

Despejando µ , calculamos su valor:

µ=2.067⋅10-3 [Ns/m2]


Pág. 11

- GRÁFICOS CUESTIÓN 1 -


Pág. 12

GRÁFICO 1GRÁFICO 1GRÁFICO 1GRÁFICO 1y = 1.0162x + 1.9631R2 = 0.9567

y = 1.7179x + 2.0808R2 = 0.9296

0

0.5

1

1.5

2

2.5

3

-0.5 -0.4 -0.3 -0.2 -0.1 0 0.1 0.2 0.3

log(v)

log PÉ

RDIDA DE CARGA

Log(V) Log(∆h)-0,423 1,556-0,327 1,602-0,267 1,69-0,213 1,732-0,15 1,833

-0,136 1,908-0,123 1,886-0,096 1,869-0,048 2-0,036 2,037-0,004 2,0830,054 2,1070,054 2,1880,071 2,1820,08 2,26

0,089 2,1430,105 2,1880,121 2,320,151 2,2880,151 2,3930,18 2,441

0,192 2,464


Pág. 13

GRÁFICO 2GRÁFICO 2GRÁFICO 2GRÁFICO 2y = -0.9845x + 1.4369R2 = 0.9531

y = -0.2805x - 0.8913R2 = 0.2598

-2.5

-2

-1.5

-1

-0.5

0

0 0.5 1 1.5 2 2.5 3 3.5 4

log Re

log f

Log(Re) Log(f)3,052 -1,5453,148 -1,693,209 -1,7263,262 -1,793,324 -1,8123,339 -1,7673,352 -1,8153,379 -1,8863,427 -1,8513,439 -1,8393,471 -1,8573,529 -1,9473,529 -1,8663,546 -1,9073,555 -1,8483,563 -1,9793,58 -1,971

3,595 -1,8663,626 -1,9633,626 -1,8553,655 -1,8663,667 -1,867


Pág. 14

AMPLIACIÓN GRÁFICO 2 (A)AMPLIACIÓN GRÁFICO 2 (A)AMPLIACIÓN GRÁFICO 2 (A)AMPLIACIÓN GRÁFICO 2 (A) y = -0.9845x + 1.4369R2 = 0.9531

-1.9

-1.85

-1.8

-1.75

-1.7

-1.65

-1.6

-1.55

-1.5

3 3.05 3.1 3.15 3.2 3.25 3.3 3.35

log(Re)

log(f)


Pág. 15

AMPLIACIÓN GRÁFICO 2 (B)AMPLIACIÓN GRÁFICO 2 (B)AMPLIACIÓN GRÁFICO 2 (B)AMPLIACIÓN GRÁFICO 2 (B) y = -0.3343x - 0.7056R2 = 0.3227

-2

-1.95

-1.9

-1.85

-1.8

-1.75

3.3 3.35 3.4 3.45 3.5 3.55 3.6 3.65 3.7

log(Re)

log(f)


Pág. 16

GRÁFICO 3GRÁFICO 3GRÁFICO 3GRÁFICO 3

0

50

100

150

200

250

300

350

0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8

VELOCIDAD (m/s)

PÉRDIDA DE CARGA (m.c.a.)

Velocidad (m/s) ∆h (m.c.a.)0.707 680.731 810.754 770.802 740.896 1000.92 1090.99 121

1.132 1281.132 1541.178 1521.203 1821.227 1391.273 1541.32 209

1.415 1941.415 2471.514 2761.556 291


Pág. 17

GRÁFICO 4GRÁFICO 4GRÁFICO 4GRÁFICO 4

0

10

20

30

40

50

60

70

80

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8

VELOCIDAD (m/s)

PÉRDIDA DE CARGA (m.c.a.)

Velocidad (m/s) ∆h (m.c.a.)0,378 360,471 400,542 490,613 540,707 68


Pág. 18

0

50

100

150

200

250

300

350

0 0.5 1 1.5 2

Velocidad (m/s) ∆h (m.c.a.)0,378 360,471 400,542 490,613 540,707 680,731 810,754 770,802 740,896 1000,92 1090,99 121

1,132 1281,132 1541,178 1521,203 1821,227 1391,273 1541,32 209

1,415 1941,415 2471,514 2761,556 291


Pág. 19

2.2.2.2.---- ¿Se pone de manifiesto la existencia de dos tipos o regímenes en el movimie ¿Se pone de manifiesto la existencia de dos tipos o regímenes en el movimie ¿Se pone de manifiesto la existencia de dos tipos o regímenes en el movimie ¿Se pone de manifiesto la existencia de dos tipos o regímenes en el movimiento?nto?nto?nto?

Si, se pone de manifiesto la existencia de régimen laminar y turbulento, puesto que se cumple que la pérdida de carga en función de la velocidad se ajusta primero a una recta (régimen laminar) y luego a una curva de orden cercano a dos (régimen turbulento).

Por otra parte, también se confirma que el factor de fricción en función del número de Reynolds se ajusta primero a una recta (régimen laminar) y luego a una curva de orden –0.25 (régimen turbulento).

3.3.3.3.---- ¿Se confirman las relaciones f=64/R ¿Se confirman las relaciones f=64/R ¿Se confirman las relaciones f=64/R ¿Se confirman las relaciones f=64/Reeee cuando el régimen es laminar, y f=0.316/R cuando el régimen es laminar, y f=0.316/R cuando el régimen es laminar, y f=0.316/R cuando el régimen es laminar, y f=0.316/Reeee

0.250.250.250.25 cuando el régimen cuando el régimen cuando el régimen cuando el régimen es turbulento?es turbulento?es turbulento?es turbulento?

Los datos experimentales si que se ajustan a curvas de la forma de las teóricas, excepto por un desajuste en los coeficientes numéricos. El régimen laminar se ajusta a una curva de la forma f=27.35⋅Re

-0.98 y el régimen turbulento a una curva de la forma f=0.1284⋅Re

-0.28.

Como se ve, el desajuste estriba en que los coeficientes numéricos de nuestras curvas son aproximadamente 0.4 veces los coeficientes numéricos de las relaciones teóricas

4.4.4.4.---- ¿Coinciden los valores medios experimentales aceptados de f y ¿Coinciden los valores medios experimentales aceptados de f y ¿Coinciden los valores medios experimentales aceptados de f y ¿Coinciden los valores medios experimentales aceptados de f y µµµµ con los valores determinados con ayuda con los valores determinados con ayuda con los valores determinados con ayuda con los valores determinados con ayuda de los datos recogidos?. Si no es asi, señalar las razones de la discrepancia.de los datos recogidos?. Si no es asi, señalar las razones de la discrepancia.de los datos recogidos?. Si no es asi, señalar las razones de la discrepancia.de los datos recogidos?. Si no es asi, señalar las razones de la discrepancia.

Los valores experimentales aceptados de f serán los obtenidos aplicando las fórmulas f=64/Re para el régimen laminar y f=0.316/Re

0.25 para el régimen turbulento, ya que estos valores se construye el “diagrama de Moody”, considerando régimen turbulento para Re>3000.

Tomando como viscosidad 1.005 cp, los valores experimentales aceptados para nuestro experimento serían:

TABLA 3:

Medida f1 0.05672 0.04553 0.03964 0.03505 0.03036 0.02937 0.02848 0.02679 0.0239

10 0.023311 0.021712 0.041513 0.041514 0.041015 0.040816 0.040617 0.040218 0.039919 0.039220 0.039221 0.038522 0.0383


Pág. 20

Realizando el cálculo del valor medio obtenemos un valor muy aproximado de f:

f=0.0364

El valor medio determinado con la ayuda de los datos recogidos es f=0.1136, el cual es aproximadamente tres veces mayor.

Como se aprecia, nuestro valor medio del coeficiente de fricción es mucho mayor que el valor medio experimental. Esto puede ser debido a varios factores:

1.- Errores producidos en el cálculo del número de Reynolds: el nº de Reynolds presenta la ecuación Re=ρ⋅v⋅D/µ. Para el cálculo de los valores teóricos, se ha supuesto un valor de viscosidad de 1.004 cp, correspondiente a una temperatura de 20 ºC. En realidad, el agua sale a una temperatura menor, por lo cual la viscosidad aumentaría. Ambos factores determinan que el nº de Reynolds sería menor.Otro valor en el cual podemos arrastrar error puede ser el del diámetro del tubo. Se ha supuesto que ∅=0.003 m, pero pueden existir precipitados en el interior del tubo que hagan ese diámetro menor. Esto también llevaría a un valor del nº de Reynolds menor.Si el nº de Reynolds es menor, los coeficientes de fricción calculados mediante las fórmulas f=64/Re

y f=0.316/Re

0.25 serían mayores, por lo que se aproximarían al valor calculado dado por los datos recogidos.

2.- Errores al considerar la rugosidad: si existen los precipitados nombrados anteriormente, la tubería ya no será lisa, sino que presentará cierta rugosidad. Por ello, el coeficiente de fricción será mayor, de forma queel valor se aproxima al calculado mediante los datos recogidos.

3.- Errores en el cálculo de las pérdidas de carga: si las pérdidas de carga reales son menores que las tomadas por nosotros, el coeficiente de fricción calculado mediante los datos recogidos sería mayor que el valor real. Asi, el valor calculado mediante los datos recogidos se aproximaría al valor experimental aceptado. Este error al calcular las pérdidas de carga podría deberse a que no se ha dejado suficiente tiempo entre unas medidas y otras ó una inadecuada puesta en marcha del aparato.

POSIBLEMENTE EL ERROR SE DEBA A UNA MEZCLA DE ESTOS TRES FACTORES.

En lo que difiere a la viscosidad, el valor teórico tomado es el correspondiente a 20 ºC, es decir, 1.004 cp. El valor obtenido utilizando nuestros datos es 2.067 cp.

El error se debe a la mezcla de dos factores:

1.- El valor de viscosidad considerado es el correspondiente a una temperatura de 20 ºC, siendo realmente menor la temperatura del agua. Por ello, la viscosidad del agua utilizada es mayor de 1.004 cp, aunque es imposible que tome el valor de 2.067 cp, ya que la viscosidad del agua a una temperatura de 0 ºC es menor.

2.- Las pérdidas de carga reales pueden ser menores que las tomadas por nosotros, de forma que la viscosidad calculada sería menor, puesto que viscosidad y pérdida de carga son directamente proporcionales.


Pág. 21

5.5.5.5.---- Para un tubo del mismo diámetro y material que el del dispositivo de la práctica y de 5 m de longitud, Para un tubo del mismo diámetro y material que el del dispositivo de la práctica y de 5 m de longitud, Para un tubo del mismo diámetro y material que el del dispositivo de la práctica y de 5 m de longitud, Para un tubo del mismo diámetro y material que el del dispositivo de la práctica y de 5 m de longitud, cuya presión a la entrada es de 2 atm, calcular la cuya presión a la entrada es de 2 atm, calcular la cuya presión a la entrada es de 2 atm, calcular la cuya presión a la entrada es de 2 atm, calcular la caída de presión en dicho tubo.caída de presión en dicho tubo.caída de presión en dicho tubo.caída de presión en dicho tubo.

La caída de presión en un tubo se calcula aplicando la ecuación de Bernouilli entre los puntos de entrada y salida de la tubería:

P1/ρ+v1

2/2+z1= P2/ρ+v2

2/2+z2+∆h

Dado que la densidad en los puntos 1 y 2 es la misma (se trata del mismo fluido), y la velocidad también coincide (el caudal y la sección son constantes), la ecuación anterior se reduce a:

P1+z1= P2+z2+∆z ⇒ P1-P2=z2-z1+∆z

Esta es la ecuación utilizada para el cálculo de la caída de presión en el tubo. Como z2-z1=-5 [m], ya que el tubo está en posición vertical:

P1-P2=-5+∆z

∆z se calcula mediante la ecuación de Darcy-Weisbach: ∆h=f⋅L⋅v2/(2⋅g⋅D), donde desconocemos f y v. Podemos calcular f mediante el nº de Reynolds: Re=ρ⋅v⋅D/µ., donde desconocemos v.

Sustituyendo los datos conocidos en las ecuaciones anteriores queda:

∆h=∆z=0.8503⋅f⋅v2

Re=2988⋅v

Sin embargo, teniendo en cuenta que rugosidad, diámetro, viscosidad y densidad son iguales en nuestro experimento y en esta cuestión, el factor rugosidad/diámetro será el mismo en ambos casos y el nº de Reynolds será también igual siempre que se evalúe para las mismas velocidades. Por ello, el factor f será el mismo en ambos casos.

Consideramos como valor de f el obtenido aplicando las ecuaciones f=64/Re y f=0.316/Re

0.25, puesto que se asemeja mucho mas que el calculado con la ayuda de los datos recogidos al valor adecuado para las características de la tubería de la cuestión.

Este valor es f=0.0364.

Por lo cual, las anteriores ecuaciones quedarán:

∆h=0.031⋅v2

Re=2988⋅v

Para calcular la pérdida de carga, como se conoce el factor ε/∅ =0, dado que la rugosidad es 0 y conocemos el valor de f, calculamos el valor de Re y con éste la velocidad.

El valor de Re que obtenemos en el “diagrama de Moody” es Re=5000 ⇒ v=1.67 [m/s] ⇒ ∆h=0.086 [m].

Luego:

P1-P2=-5+0.086=-4.91 [atm]

Esta es la pérdida de presión en la tubería. Como se ve no existe pérdida sino ganancia, lo cual es lógico dado que el fluido desciende y la tubería es lisa.

La solución será válida siempre que el coeficiente de fricción se ajuste al utilizado, para lo cual la velocidad ha de ser similar a la correspondiente en nuestras medidas experimentales a un coeficiente de fricción como el utilizado.


Pág. 22

BIBLIOGRAFÍA

C. Mataix, Mecánica de fluidos y máquinas hidráulicas.

White, Mecánica de fluidos.

Streeter-Wylie, Mecánica de los fluidos.

Costa, Curso de Ingeniería Química.

McCabe, Operaciones Unitarias en Ingeniería Química.

Realizado por:

SERGIO DE LAS PEÑAS LÓPEZ (NP 7460)

GRUPO IND-3102

estudio de pérdidas de carga en tuberías

Documents