REDES DE TUBERÍAS

Prof. Andrés Granados (Nov/2020-Corregido Abr/2023)

ECUACIÓN DE TUBERÍA

QQ
F (Q ) = PB − PA + ρ K f + K g L senα = 0
2 A2

Kg = ρ g Flujo Incompresible

K f = ∑ K i + f (L + ∑ L j )
$
4A

K ′f = cK f Unidades de Presión

K g′ = cK g c = 10-6 P (MPa)
-5
c = 10 P (Bares)
K ′f′ = 10 −6 K ′f c=1 P (Pa)
Q (Litros/seg) c = (ρg)-1 P (m)

Tubería curvada L senα = ZB - ZA Dh = 4A/ $ = Diámetro Hidráulico $ = Perímetro Mojado

Tubería sección circular ( $ = π D A = π D / 4)
2

ρKf 8ρ K f 1 $ ρ Q D 4ρ Q
Cf = = C ′f = c C f = Re D = =
2A 2
π D
2 4
D 4A µA µ$

[Granados,2021]
Granados, A.L. Mecánica y Termodinámica de Sistemas Materiales Continuos. Universidad Simón
Bolívar, Agosto, 2021.
Para el caso de flujo incompresible revisar la Sección IX.1.4.3, ec.(36)-(38), pág.202.
Para el caso de flujo compresible revisar la Sección IX.3.2.2, ec.(45)-(47), pág. 230.

Ecuación de Colebrook [1938-39]

 2,51X K e / Dh  f =
1
X K + 1 = −2 log + 
 Re D 3,7  X2
1
 3,7 
2
ft = = lim f
X 0 = log  X 02 Re→∞
 ε / Dh 
ECUACIÓN DE NODO

F (P ) = ∑Q − ∑Q i j =0
Entrando Saliendo

Tubería Entrando Qi > 0 se suma el caudal

Tubería Saliendo − Q j < 0 se resta el caudal

ECUACIÓN FUENTE / SUMIDERO

Caudal Conocido Presión Conocida

F (P ) = ∑ Qi −∑ Q j + Q = 0 F (Q ) = Pm − P = 0
Entr Sal n

∑P i

Fuente Sumidero
Pm = i =1
ó Pm = P1
n
Q>0 Q<0

ESTIMADOS INICIALES

 Todas las presiones iguales a la presión media de las fuentes de presión (debe haber al
menos una fuente de presión).

 Los caudales siguen la ley de proporcionalidad D5/2 supuesta bajo la hipótesis de que todos
los gradientes de presión (sin influir en f) son todos iguales.

∆Q = ∑ Qi −∑ Q j Qi′ = Qi − r ∆Q Tubería entrando Dhi5 / 2

r=
Q ′j = Q j + r ∆Q Tubería saliendo ∑ Dhi5 / 2
Entr Sal
OTRAS CORRELACIONES EXPLICITAS

1  7  0,9 e / Dh 
Swamee & Jain [1976] and X = = −2 log   + 
f  Re  3,7 
Churchill [1973]
Zigrang & Sylvester [1982] 5x103 < Re < 107
4x10 < e / Dh < 0,05
-5

1  6,9  e / D 1,11 
Haaland [1983] X = = −1,8 log  + h
  (2% error)
f  Re  3,7  

BOMBAS

F (PA ) = QB − Q A = 0
c W η
F (PB ) = PB − PA − =0 Bombas Reciprocantes η<1
QA
F (PB ) = PB − PA − K g H B (Q A ) = 0 Bombas Centrífugas
H B (Q) = H o − C B Q 2
Curva de la bomba
c' = c /( ρ g ) c' ' = 10 3 c' C B′ = 10 −6 C B metros Litros/seg

TURBINAS

F (PA ) = QB − Q A = 0
c W
F (PB ) = PB − PA + =0 Potencia η<1
QB η
F (PB ) = PB − PA + K g H T (Q A ) = 0 Curva
c' = c /( ρ g ) c' ' = 10 c'
3
metros Litros/seg
TANQUES Y RECIPIENTES

Tanque Recipiente

F (PB ) = PB − Patm − K g H T = 0 F (PB ) = PB − PV − K g H V = 0

TANQUES Y RECIPIENTES DE COMPENSACIÓN

Tanque Recipiente

F (PA ) = PA − Patm = 0 F (PA ) = PA − PV = 0

F (PB ) = PB − Patm − K g H T = 0 F (PB ) = PB − PV − K g H V = 0

FLUJO COMPRESIBLE

Para redes de tuberías el tratamiento más apropiado es el de flujo isotérmico. El flujo

isotérmico evoluciona con la temperatura constante. El flujo tiende hacia un flujo crítico
determinado por el número de Mach 𝑀𝑀𝑀𝑀+ = 1/√𝑘𝑘 condición indicada en las variables
con el superíndice “+”. Para el cálculo del flujo másico del flujo isotérmico se puede usar la
siguiente ecuaci\'on integral
Esta ecuación se puede colocar en función del caudal volumétrico 𝑄𝑄� a condiciones
estándares, donde 𝑆𝑆𝑔𝑔 = 𝜌𝜌�/𝜌𝜌
�𝑎𝑎 es la densidad relativa a la del aire 𝜌𝜌𝑎𝑎 a condiciones
estándares de presión 𝑃𝑃𝑎𝑎 (101325 pa) y temperatura 𝑇𝑇𝑎𝑎 (20°C).

La expresión para el caudal estándar en el flujo isotérmico se puede colocar en una forma
similar a la expresión para el flujo incompresible [Granados,2021;IX.3.2.2]

donde los coeficientes de p\'erdida 𝐶𝐶𝑓𝑓 y del energía potencial 𝐾𝐾𝑔𝑔 involucrados son
distintos a los anteriores.

El primer término de 𝐾𝐾𝑓𝑓 , donde involucra las 𝐾𝐾𝑖𝑖 , se denomina pérdidas localizadas e
incluye las pérdidas de presión ocasionadas por accesorios (codos, tees, yees, reducciones,
expansiones, válvulas, etc.) a lo largo de la tubería, afectado por un factor 𝑌𝑌𝑖𝑖 de expansión
neta que tiene en cuenta el efecto de la compresibilidad en los accesorios [Crane,1986]. El
segundo término de 𝐾𝐾𝑓𝑓 son las pérdidas distribuidas ocasionadas por la fricción viscosa
con la pared de la tubería a lo largo de la longitud L real de la misma tubería.

Las longitudes 𝐿𝐿𝑖𝑖 son pérdidas menores (localizadas), pero expresadas como longitud
equivalente de tubería que produciría la misma pérdida de la presión, que si se reemplaza
por dicha longitud equivalente adicionada a la longitud original de la tubería. El tercer
término incluye el efecto de la diferencia de presión en la reducción del flujo másico y se
ha encerrado entre valores absolutos por si el flujo se revierte, ( 𝑃𝑃𝐴𝐴 debe ser siempre mayor
que 𝑃𝑃𝐵𝐵 cuando el sentido del flujo es el correcto ).

El factor de fricci\'on de Darcy-Weisbach f interviene como una función del tipo

dependiente del número de Reynolds Re y la rugosidad relativa 𝜀𝜀/𝐷𝐷 mediante, por

ejemplo, la correlación de Colebrook. La longitud D es el denominado Diámetro
Hidraúlico y en el caso de una tubería de sección circular coincide con su propio diámetro.
El esfuerzo cortante en la pared de la tubería 𝜏𝜏𝑤𝑤 expresa básicamente la ecuación de
Darcy-Weisbach, que para el caso incompresible coincide con ℎ𝑓𝑓 = 𝐶𝐶𝑓𝑓 𝑄𝑄� 2 /(𝜌𝜌𝜌𝜌) . El
coeficiente clásico 𝑐𝑐𝑓𝑓 = 𝑓𝑓/4 , denominado coeficiente de fricción de Fanning, es un
número adimensional de Euler. Cuando el tramo de la tubería es recto y forma un ángulo 𝛼𝛼
con la horizontal, entonces 𝑍𝑍𝐵𝐵 − 𝑍𝑍𝐴𝐴 = 𝑙𝑙 sin 𝛼𝛼.

La fracción 𝑅𝑅ℎ se denomina Relación Hidraúlica y expresa el balance entre las fuerzas
inerciales de flujo en la sección A y las fuerzas viscosas en la pared del conducto en el
perímetro S por la condición de no delizamiento. Por eso dicho perímetro se denomina
“mojado”, porque si no lo está (mojado), no produce esfuerzo viscoso en dicha pared.

Particularmente, cuando la tubería es de sección circular las variables se reducen a las

últimas ecuaciones de arriba.

El valor absoluto que aparece en Q en las ecuaciones de arriba es para contemplar que, si el
flujo no va del punto A al Punto B, sino al contrario, las ecuaciones conservan de igual
manera su signo correcto. Si Q es positivo entonces la presión 𝑃𝑃𝐵𝐵 es menor que 𝑃𝑃𝐴𝐴 en una
tubería horizontal. Si Q es negativo, entonces ocurre lo contrario. El coeficiente 𝐾𝐾𝑔𝑔 en la
ecuación interviene para considerar el efecto de la gravedad como se hizo en el caso
incompresible. Dicha ecuación se ha colocado igualada a cero para formar, en caso de una
red de tuberías, un sistema de ecuaciones algebraicas no lineales del tipo homogénea
f(x)=0.

PROCEDIMIENTO

Identificados cada uno de los elementos de la red de tuberías, cada uno con una ecuación
homogénea, se crea un sistema de ecuaciones no lineales homogénea f(x)=0. Este sistema
se resuelve con algún método, e.g. Método de Newton-Raphson, aunque el ideal sería un
método de segundo orden (ver: Granados A. L., Método Numérico para Redes), ya que las
tuberías y turbomáquinas tienen una ecuación aproximadamente cuadrática (las otras son
lineales). Cada elemento de la red agrega una ecuación y una incógnita al sistema, que es
diagonalmente dominante en la matriz jacobiana, la cual es dispersa, si la numeración se
hace consecutiva en elementos próximos. El sistema de unidades escogido debe hacerse
para que las variables (presiones y caudales) sean del orden de la unidad, de manera que la
matriz jacobiana no esté mal condicionada. Los iterados iniciales se escogen como arriba se
indica (ESTIMADOS INICIALES), tal que se esté próximo a una cuenca de aproximación
convergente.
BIBLIOGRAFIA:

Juan G. Saldarriaga V. Hidráulica de Tuberías: Abastecimiento de Agua, Redes y Riegos.

Alfaomega (Bogotá) & Universidad de Los Andes, 2007.

Manuel Vicente Méndez. Tuberías a Presión en Los Sistemas de Abastecimiento de

Agua. Fundación Polar & Universidad Católica Andrés Bello, 1995.

Hardy Cross. “Analysis of Flow in Networks of Conduits or Conductors”. Univ. Illinois

Bull., Engineering Experiment Station, No.286, (1936).

Crane Co. Flow of Fluids Through Valves, Fittings, and Pipe. Metric Edition-S.I. Units.
Technical Paper No. 410M, Sixth Printing, 1986.

Granados, A. L. “Métodos Numéricos Para Redes”. Universidad Simón Bolívar, Ene.,

(2016).

Granados, A.L. Mecánica y Termodinámica de Sistemas Materiales Continuos.

Universidad Simón Bolívar, Agosto, 2021.

[Link]