Hidráulica de tuberías

Redes ramificadas

Métodos numéricos
aplicados al cálculo del
coeficiente de fricción
Holger Benavides Muñoz
hmbenavides@utpl.edu.ec
1. Diseño de redes ramificadas o
abiertas
1.2 Modelos matemáticos para hidráulica a
presión.
 Conferencia:
“Métodos numéricos aplicados al cálculo del
coeficiente de fricción”
 Contenidos:
 Punto fijo.

 Newton - Raphson.

 Ejercicios de aplicación.

 Software de aplicación.

11/04/2017 hmbenavides@utpl.edu.ec 2
Exordio del tema de estudio

Resultados de aprendizaje (RA)

 El estudiante podrá identificar los métodos

numéricos para el cálculo del coeficiente de
fricción (f) de Darcy-Weisbach

 El estudiante podrá analizar tuberías simples

con nudo de altura conocido.

11/04/2017 hmbenavides@utpl.edu.ec 3
Objetos de aprendizaje.- Bibliografía
 SALDARRIAGA, J. 2008. Hidráulica de tuberías. Bogotá, Col.,
Alfaomega. 690 p.

 Gould, E. Sherman (Edward Sherman), b. 1837. Practical hydraulic

formulae for the distribution of water through long pipes.

 SOTELO, G. 1987. Hidráulica General, vol. I, fundamentos. México,

Limusa. 561 p.

 MATAIX, C. 1982. Mecánica de fluidos y máquinas hidráulicas. 2ª

de. México, HARLA.660 P.

 PEREZ, R. 1993. Dimensionado óptimo de redes de distribución de

agua ramificadas. Universidad Politécnica de Valencia. Tesis
doctoral.
 Reconocimiento y agradecimiento: Dr. Ing. Rafael Pérez García. GMMF.

11/04/2017 hmbenavides@utpl.edu.ec 4
 ANTECEDENTES

En el interior de las tuberías existen

protuberancias o irregularidades de diferentes
formas y tamaños que se conocen como
rugosidad.

Observaciones microscópicas en
Pineda & Benavides – UTPL - 2008 5
 Flujo hidráulicamente liso

La rugosidad () queda cubierta por la subcapa

laminar (). Entonces, la rugosidad no influye en
el valor del factor de fricción (ƒ), comportándose
la tubería como un material liso.

  0.305 

6
 Flujo hidráulicamente semirrugoso o
zona de transición
El espesor de la subcapa laminar () se aproxima
al valor medio de rugosidad absoluta (), el
coeficiente de fricción (ƒ) depende tanto del
número de Reynolds (Re) como de la rugosidad
relativa (/D).

7
 Flujo hidráulicamente rugoso

El espesor de la capa límite () es menor que la

rugosidad absoluta (), por lo tanto, las
irregularidades internas de la conducción
sobrepasan la subcapa laminar, produciendo
turbulencia completa.

  6.1 

8
 RUGOSIDAD
Irregularidades o asperezas de diferentes
alturas y profundidades que se forman en el
interior de un conducto.

Tuberías de hierro dúctil

9
 Origen de la rugosidad en una tubería

El mal terminado en la fabricación, longitud corta

entre uniones, características físico - químicas y
bacteriológicas del agua, el fenómeno de corrosión
no es controlado y el tiempo de servicio.

Tuberías de hierro dúctil 10

 Influencia del tiempo en el envejecimiento
de conductos

El transporte de fluidos por tubos puede producir

corrosión o adherencia de materia orgánica o
sedimentos a las paredes del mismo; (formándose
irregularidades en su superficie interna).

Tuberías de hierro
dúctil; investigación
UTPL -2008

11
 Rugosidad absoluta y rugosidad relativa

Las irregularidades del interior de un conducto se

expresan como la altura media () de las
asperezas (rugosidad absoluta).

Rugosidad relativa =  / D

12
 PÉRDIDAS CONTÍNUAS

Las pérdidas por rozamiento son función de la

rugosidad del conducto, de la viscosidad
cinemática del fluido, del régimen de
funcionamiento (flujo laminar o turbulento) y del
caudal circulante, es decir, de la velocidad; (a
más velocidad, más pérdidas).

13
DIRECCIÓN DEL FLUJO EN EL INTERIOR DE
UNA TUBERÍA

14
Ecuaciones para el cálculo de pérdidas
de carga por longitud:
 DARCY - WEISBACH
L  V   8  λ L 
2

hf  λ    
 
 2 
5 
Q 2

D  2  g   π  g D 

0.2
 λi 
 0.6073   *   Qi 
0.4
Di  teorico 
J 

15
 HAZEN - WILLIAMS
1.85185
10.646  L  Q 
hf  
4.87037 

D  C HW 
0.38
1.6252  Qi 
Di  teorico    0.21  
J*  C HW 
16
 CHEZY - MANNING
 L 
h f  10.29  n   16   Q
2 2

D 3

Di  teorico  
1.5482
* 0.1875
n  Qi 
0.375

J
17
 Ecuación de Darcy – Weisbach
para el cálculo de pérdidas por longitud
L  V 2   8  f  L 
hf  f     Q 2
D  2  g    2  g  D 5 
Donde:
hf _ pérdida de carga.
f _ coeficiente de pérdida de carga.
L _ longitud de la tubería.
D _ diámetro de la tubería.
V _ velocidad media del fluido.
g _ aceleración de la gravedad.
Q _ caudal

18
 El factor f ó
Depende de la velocidad (V), del diámetro de la
tubería (D), de la densidad (ρ), de la viscosidad
cinemática (ν) y de la rugosidad (ε)

V D  
f  F , 
  D
Número de Reynolds Rugosidad relativa

19
 64
 Para flujo laminar: f 
Re

 En flujo laminar (Re < 2000), la

rugosidad relativa del ducto no influye
en las pérdidas por fricción.

20
 Según Blassius:
 Para flujo turbulento, hidráulicamente liso ( 
0.305 ´):
0.25
f  0.316  Re
  
 ´ 32.81 
 f  V 

 Según Prandtl-von Kármán:

 Para flujo turbulento, hidráulicamente rugoso (
 6.1 ´):
1 D
 2  log   1.14
f  
21
 Modelos matemáticos
Cálculo del factor de fricción.

 Para determinar el factor de fricción f en tubos lisos

o rugosos, para flujos en la zona de transición (2200
< Re < 4000) o turbulenta (Re > 4000), se utiliza la
ecuación de Colebrook-White

 Forma implícita.

1   2.51 

 2 log10   

f  3.7 D Re f 

22
23
 Factor de fricción (  ó f )

 Aproximación de Barr,
cuando Re > 4000

1
 2
  ε 5.1286 
 2  log10   0.89 
  3.71 D Re 

24
 Free surface flow
Another form of the Colebrook-White equation
exists for free surfaces. Such a condition may
exist in a pipe that is flowing partially full of fluid.
For free surface flow:

25
 Haaland equation
 The Haaland equation is used to solve directly
for the Darcy–Weisbach friction factor f for a
full-flowing circular pipe.
 It is an approximation of the implicit Colebrook–White
equation. It was developed by S. E. Haaland in 1983.

26
Brkić solution
 Brkić shows one approximation of the
Colebrook equation based on the Lambert
W-function

27
 Ecuación de Barr:

 Otra forma:

28
 Factor de fricción (  ó f )

 Aproximación de Swamee y Jane,

cuando Re > 4000

1
λ 2
  ε 5.74 
4  Log  0.9 
  3.71 D Re 

29
 Otras ecuaciones para calcular f

 Sousa-Cunha-Marques, 1999 (error = 0.123%):

 Haaland, 1983 (error = 0.220%):

30
 Barr, 1972 (error = 0.375%):

 Swamee-Jain, 1976 (error = 0.386%):

31
Malafaya - Baptista

11/04/2017 hmbenavides@utpl.edu.ec 32
Chen

Manadilli

33
Streeter

Miller

Pavlov

34
Camaraza y Landa

35
 Avci, Kargoz

Zigrang and Sylvester

36
 Evangleids, Papaevangelou, Tzimopoulos

Moody
  
1

   10 
6 3

  0.0055 1   20000    
   D  Re  

37
Interpolación
cúbica al
diagrama de
Moody para
calcular f,
cuando
2000< Re<4000

38
39
Rugosidad absolutas ()

11/04/2017 hmbenavides@utpl.edu.ec 40
 V - Velocidad (m/s) Q 4Q
 V 
 πD 
2
π  D 2
 
 4 
  - Coef. de viscosidad cinemática (m2/s)
0.000497
υ
Tº42.51.5

VD
 Re - Número de Reynolds Re 


hf
J S
 J - Gradiente hidráulico (m/m) L

41
Re crítico
 El valor numérico del número adimensional de
Reynolds, a partir del cual el factor de fricción comienza
a ser constante es determinado mediante la siguiente
expresión, recomendada por [Kreit et al, 1999], [Crane,
2000], [Fernández Diez, 2000], [Mark, 2001], [Nayyar,
2003], [Shames, 2001], [Wallas, 1999], [White, 2001],
[Streeter, 2000], [Sámano, 2003] y [Perry et al, 1999]:

Coeficiente de Viscosidad Cinemática (agua fría) = 1.010 e -6 m2/s

Coeficiente de Viscosidad Cinemática (agua caliente) = 0.478 e -6 m2/s

11/04/2017 hmbenavides@utpl.edu.ec 42
 CÁLCULO DE f POR TANTEOS

t 12 ºC

v 1.235E-06
m2/s
1   2.51 
 2 log10   

f  3.7 D Re f 
Q D 
(L/s) (m) (m)
185 0.3518 0.000015
Q 4Q
V 
 π  D2  π  D2
 
x
1    4  2.51  f i 1 
1
g ( x)  2 log10    g ( x )2
f
 3.7 D Re f 

43
 CÁLCULO DE f POR TANTEOS

x
1   2.51  f i 1 
1
g ( x)  2 log10   

f
 3.7 D Re f  g ( x )2

fi x g(x) f i+1
0 0.0100000 10.0000000 8.4756781 0.013920380
1 0.0139204 8.4756781 8.5887461 0.013556278
2 0.0135563 8.5887461 8.5798347 0.013584453
3 0.0135845 8.5798347 8.5805337 0.013582240
4 0.0135822 8.5805337 8.5804789 0.013582413
5 0.0135824 8.5804789 8.5804832 0.013582400
6 0.0135824 8.5804832 8.5804828 0.013582401
7 0.0135824 8.5804828 8.5804829 0.013582401
8 0.0135824 8.5804829 8.5804829 0.013582401
9 0.0135824 8.5804829 8.5804829 0.013582401
10 0.0135824 8.5804829 8.5804829 0.013582401
44
 Modelos matemáticos
CÁLCULO DE f POR NEWTON
RAPHSON
g xi   xi
xi 1  xi 
g ' xi   1

1   2.51 
x 
g ( x)  2 log10   
f 
 3.7 D Re f 

 2.51 
 
 2   Re  1
g ' ( x)    
 ln 10     2.51  fi1 
 3.7 D  (xi1)2
 Re f 

45
 Cálculo de f por Newton - Raphson
t 12 ºC

v 1.235E-06
m2/s

Q D 
(L/s) (m) (m)
185 0.3518 0.000015
Q 4Q
V 
 π  D2  π  D2
fi   Xi+1 f i+1
x g(x) g'(x)
 4 
0 0.0100000 10.0000000 8.4756781 -0.0695509 8.5748019 0.0136004040616
1 0.0136004 8.5748019 8.5809288 -0.0785106 8.5804828 0.0135824011871
2 0.0135824 8.5804828 8.5804829 -0.0784703 8.5804829 0.0135824008510
3 0.0135824 8.5804829 8.5804829 -0.0784703 8.5804829 0.0135824008510
4 0.0135824 8.5804829 8.5804829 -0.0784703 8.5804829 0.0135824008510

f i - f i+1
fi x g(x) g'(x) Xi+1 f i+1 D
0.0135824082 8.5804805 8.5804830 -0.0784703 8.5804829 0.013582400851 0.0000

46
 Cálculo de f por Newton - Raphson

t 14 ºC Q  D visc. Cine Re
u 1.1703E-06 m2/s 2
(L/s) (m) (m) (m /s)
V 1.80035 m/s 175 0.0000015 0.3518 1.1703E-06 541213.49

47
 RELACIÓN ENTRE ECUACIONES DE PÉRDIDAS
CONTÍNUAS

Darcy – Weisbach y Hazen - Williams

hf Darcy - Weis bach  hf Hazen Williams

 8 f  L 
1.85185
10.646  L  Q 
h f   2   Q2
 hf  
C


 gD
5 4.87037
 D  HW 

128.84  D 0.12963
f  1.85185
CHW  Q 0.14815

1.000057  D 0.070
CHW 
f 0.540  Q 0.080
48
RELACIÓN ENTRE ECUACIONES DE PÉRDIDAS
CONTÍNUAS

Darcy – Weisbach y Chezy - Manning

hf Darcy - Weisbach  hf Chezy- Manning

 8 f  L   L  2
h f   2   Q2
 h f  10.29  n   16 2
  Q
   g  D 5
 D 3 

n2
f  124.536  1
D 3

1 1
n  0.090  D  f 6 2
49
 RELACIÓN ENTRE LA RUGOSIDAD ABSOLUTA Y
EL TIPO DE AGUA QUE TRANSPORTA LA TUBERÍA

t  o
Genijew  t   o  (a  t ) a
t
Donde:
o _ rugosidad del tubo (nuevo), en mm.

a _ coeficiente que depende del grupo en el que

se clasifique el agua.
t _ número de años de servicio de la tubería.
t _ rugosidad absoluta del conducto, después
de t años de servicios en mm.
50
 Identificación del tipo de agua que circula por la
tubería

GRUPO COEFICIENTE a CARACTERÌSTICAS DEL AGUA

I varia de 0.005 a 0.055; Agua con poco contenido mineral que no

valor medio; 0.025 origina corrosión.
II varia de 0.055 a 0.18 ; Agua con poco contenido mineral que
valor medio; 0.07 origina corrosión.
III varia de 0.18 a 0.40; valor Agua que origina fuerte corrosión y con
medio; 0.20 escaso contenido de cloruros y sulfatos.

IV varia de 0.40 a 0.60; Agua que origina corrosión, con un gran

valor medio; 0.51 contenido de cloruros, con una gran
cantidad de materia orgánica.
V varia de 0.6 a mas de 1 Agua con cantidades importantes de
carbonatos, pero de dureza pequeña
permanente.
51
 Concentración total de Calcio y Magnesio
AGUAS (mg/L)
Blandas 1. - 75
Dureza Moderadamente duras 75- 150
Duras 150 - 300
Muy duras pasa de 300

Escala de pH entre 0 (muy ácido) a

Potencial
14 (muy básico). Agua a 25 C es
Hidrógeno (pH) neutral si tiene pH = 7

•Bicarbonatos.
Alcalinidad •Carbonatos.
•Hidróxidos
52
 Neutraliza bases, y reaccionar con
iones hidróxido, cede protones como
Acidez
medida de su contenido total de
sustancias ácidas

•Valores de pH bajo.
•Oxígeno disuelto alto.
Corrosividad •Cloro residual libre alto.
•Alcalinidad baja.
•Concentración de sólidos disueltos alta.
•Dureza por Calcio baja.
•Temperatura alta.

Calcio y Magnesio Ocasionan la dureza en el agua

53
MEDICIÓN DE LA RUGOSIDAD ABSOLUTA EN (mm)

Escala (mm)

Pineda & Benavides – UTPL - 2008 54

 PERFÍL LONGITUDINAL DE LA RUGOSIDAD

TUBERÍA DE HIERRO GALVANIZADO

5.00
4.50
4.00
ALTURA (mm)

3.50
3.00
2.50
2.00
1.50
1.00
0.50
0.00
0.10

10.10

20.10

30.10

40.10

50.10

60.10

70.10
LONGITUD (mm)

Pineda & Benavides – UTPL - 2008 55

 CÁLCULO DEL FACTOR DE FRICCIÓN
Varios métodos

DIAMETRO RUGOSIDAD RUGOSIDAD TEMPERTURA VISCOSIDAD VISCOSIDAD

mm m ABSOLUTA RELATIVA ºc cm2/s m2/s

46 0.046 0.00011797 0.0025645 14 0.011756 1.1756E-06

Newton Punto Swamme Moody Colebrook Swamme Colebrook Swamme

Jain White Jain White Jain
Raphson fijo

VELOCIDAD REYNOLDS f f f f C HW C HW n n
0.60 23477.37 0.0065 0.0065 0.0133 0.0305 101.43 68.86 0.0043 0.0062
1.00 39128.96 0.0059 0.0059 0.0132 0.0289 97.94 63.69 0.0041 0.0062
2.00 78257.91 0.0054 0.0054 0.0131 0.0275 91.94 57.14 0.0040 0.0061
3.00 117386.8 0.0052 0.0052 0.0131 0.0270 87.90 53.58 0.0039 0.0061
5.00 195644.7 0.0050 0.0050 0.0130 0.0266 82.49 49.38 0.0038 0.0061
6.00 234773.7 0.0050 0.0050 0.0130 0.0265 80.51 47.96 0.0038 0.0061
8.00 313031.6 0.0049 0.0049 0.0130 0.0264 77.37 45.79 0.0038 0.0061
10.00 391289.5 0.0049 0.0049 0.0130 0.0263 74.94 44.18 0.0038 0.0061
56
Pineda & Benavides – UTPL - 2008
TUBERÍA DE HIERRO FUNDICIÓN DÚCTIL

TUBERÍA DE ASBESTO
CEMENTO

57
Pineda & Benavides – UTPL - 2008
 DETERMINACIONES PRÁCTICAS

•PVC, diámetro interno de 46 mm, valor promedio de

rugosidad ( =0.001193 m) ó ( =1.1930 mm), con un
periodo de servicio de 18 años.

•Polietileno, diámetro interno de 50 mm, rugosidad

promedio de ( =0.00065557 m) ó ( =0.65557 mm),
con un periodo de servicio de 20 años.

58
Pineda & Benavides – UTPL - 2008
Ejercicio de aplicación
Calcular la energía en el nudo 1

 Datos:

38.38

1 2

59
 Formulario
4 Q
V
  D2

0.000497 VD
υ Re  J S
hf
tº42.51.5  L

1
  L  V 2 
hf   λ  
2
  ε 5.1286 
 2  log10   
 D  2g 
 3.71 D Re
0.89
 

60
Resultados

61.40

38.38

1 2
Flujo turbulento

848.52 m3/h

61
62
 Velocidad en función de pérdidas
2 g  D hf
f 
 Si: L V 2

V D
 Y: Re 
v
 Entonces:

  2 2  g  D  h f    2.51 v  L 
V    log   
 L  10 
3.7  D D  2  g  D  h f 
   

63
Taller de clase: Diseño de tuberías simples
Texto básico ~ Ejem: 2.1
 Calcular el Q de agua que puede ser movido a través de
una tubería de PVC de 327.6 mm de diámetro. Longitud:
800 m, y con una carga gravitacional de: 52.5 m. Se
asume que la rugosidad absoluta () en m es: 0.00000152.
Por coefc. global de pérdidas menores se espera un Kfm:
15. La densidad del agua: 998.2 kg/m³. Viscosidad
dinámica (m): 0.001005 Kg/(m.s).

 Viscosidad cinemática [m2/s] (): 0.001005 / 998.2 =

Viscocidad cinemática: 1.00681E-06 m2/s

64
   2 2  g  D  hf   
V    log    2.51 v  L 
 L  10
 3.7 D D 2  g  D  h f 
   

H /d hf i V hf i+1
(m) (m) (m/s) (m)
V 2 
52.5
52.5
4.6398E-06
4.6398E-06
52.500000
21.758545
6.341122
3.934865
21.758545
40.662721
h f i 1  H  K fm  
52.5 4.6398E-06 40.662721 5.522929 29.179859  2g 
52.5 4.6398E-06 29.179859 4.614277 36.222056
52.5 4.6398E-06 36.222056 5.187785 31.924226
52.5 4.6398E-06 31.924226 4.844689 34.555801
52.5 4.6398E-06 34.555801 5.057139 32.947510 m3/s
Q = 0.41955
52.5 4.6398E-06 32.947510 4.928232 33.931595
52.5 4.6398E-06 33.931595 5.007446 33.329881
52.5 4.6398E-06 33.329881 4.959140 33.697959
52.5 4.6398E-06 33.697959 4.988737 33.472860
52.5 4.6398E-06 33.472860 4.970655 33.610542
52.5 4.6398E-06 33.610542 4.981722 33.526337
NOTA: Este procedimiento
52.5 4.6398E-06 33.526337 4.974956 33.577839 es válido siempre que las
52.5 4.6398E-06 33.577839 4.979095 33.546340 pérdidas menores no sean
52.5 4.6398E-06 33.546340 4.976564 33.565606
52.5 4.6398E-06 33.565606 4.978112 33.553823 mayores al 30% que las de
52.5 4.6398E-06 33.553823 4.977165 33.561029 fricción.
52.5 4.6398E-06 33.561029 4.977744 33.556622
52.5 4.6398E-06 33.556622 4.977390 33.559317
65
Tarea

 Calcular el factor f (Colebrook – White),

(aplicar el método de Newton – Raphson).

 Comparar resultados para cada método listado

 Comparar resultados por cada material
 Conclusiones y recomendaciones

11/04/2017 hmbenavides@utpl.edu.ec 66