Hidráulica Básica [2015961]

Tema # 1: Flujo real y disipación de energı́a

Luis Alejandro Morales (Ph.D)
Profesor Asistente

Universidad Nacional de Colombia-Bogotá

Facultad de Ingenierı́a
Departamento de Ingenieria Civil y Agrı́cola

Contents
1 Fluido ideal y fluido real 2
1.1 Flujo ideal . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
1.2 Flujo real . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3

2 Capa limite en flujo a presión 3

3 Esfuerzo de corte y perdidas de cabeza de energı́a 4

4 Experimentos de Reynolds 7

5 Flujo laminar 7
5.1 Flujo laminar en conductos y coronas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
5.2 Flujo laminar en tuberı́as circulares: ecuación de Hagen-Poiseuille . . . . . . . . . . . . . . 11

6 Flujo turbulento 13
6.1 Esfuerzo de corte turbulento . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
6.2 Longitud de mezcla de Prandtl . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
6.3 Distribuciones de velocidad . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
6.4 Pérdidas de energı́a . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
6.5 Cálculo del factor de fricción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
6.6 Otras ecuaciones para calcular las pérdidas de energı́a . . . . . . . . . . . . . . . . . . . . . 24

7 Pérdidas menores 25
7.1 Longitud equivalente . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30

8 Flujo normal a superficie libre 31

8.1 Estimación de la profundidad normal . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
8.2 Sección óptima de un canal . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
8.3 Cálculo de la profundidad normal en cauces naturales . . . . . . . . . . . . . . . . . . . . . 36

1
1 Fluido ideal y fluido real
1.1 Flujo ideal
Un fluido ideal es un fluido hipotético en donde se asume que el fluido no tiene viscosidad por lo tanto la
ley de viscosidad de newton

du
τ =µ (1)
dy
en donde τ es el esfuerzo de corte, µ es la viscosidad dinámica y u = f (y) es la velocidad del flujo, no
es aplicable. Esto quiere decir que la fricción en el flujo es despreciable por lo tanto no existen esfuerzos de
corte entre capas ni con los contornos, lo que implica que no hay disipación de energı́a debido a la fricción
ni formación de remolinos. En un fluido ideal las partı́culas se mueven unas sobre otras sin ningún tipo
resistencia, sometidas a fuerzas hidroestáticas aplicadas sobre su superficie. El movimiento y la aceleración
de dichas partı́culas se presenta gracias al desbalance de fuerzas actuantes de acuerdo con la segunda ley
de Newton. La suposisión de fluido ideal es de gran ayuda para el análisis de problemas prácticos en
ingenierı́a en donde las fuerzas viscosas son despreciables dando resultados precisos. Por ejemplo si se
quiere determinar la fuerza de levantamiento del ala de un avión es posible asumir un fluido ideal, sin
embargo, dicha suposisión no serı́a correcta si se quisiera determinar la fuerza de arrastre sobre el ala de un
avión. Asumiendo el flujo de partı́culas de fluido ideal e incompresible (en donde la densidad no cambia)
y de acuerdo con la segunda ley de Newton, se deduce la ecuación de Bernoulli :

p V2
+ + z = H = Constante (2)
γ 2g
donde p es la presión (absoluta o manométrica), V es la velocidad media del flujo, z es la altura del
sistema con respecto a un nivel de referencia y H es la cabeza de energı́a total en una sección del flujo
la cual es constante (H1 = H2 ) y equivale a la suma de la cabeza de energı́a de presión (p/γ), cabeza de
2
energı́a cinética (V 2 /2g) y cabeza de energı́a potencial (z). Note que al termino γp + V2g se le conoce como
cabeza de presión dinámica la cual se puede medir usando un tubo Pitot. La ecuación de Bernoulli, puede
ser expresada gráficamente a través de la lı́nea de energı́a (LE=H) y lı́nea de gradiente hidráulico (LGH
= p/γ + z).
La inclusión (a travéz de una bomba) o la extracción (a través de una turbina) de energı́a a un flujo de
un fluido ideal da lugar a una forma mas completa de la ecuación de Bernoulli conocida también como la
ecuación de trabajo-energia:

p1 V12 p2 V22
+ + z1 + hB = + + z 2 + hT (3)
γ 2g γ 2g
donde hB es la energı́a suministrada por una bomba y hT es la enegı́a sustraı́da por una turbina. La
potencia hidráulica (PH ) suministrada (bomba) al flujo o extraı́da (turbina) del flujo, se calcula como:

PH = γhQ (4)
donde h es la cabeza de energı́a mecánica (hB o hT ) y Q es el caudal. La potencia nominal (o mecánica)
(Pn ) es:
PH
Pn = (5)
η
donde η es la eficiencia de la bomba.

2
1.2 Flujo real
En un flujo real su movimiento es controlado por las fuerzas de fricción y las fuerzas turbulentas. Esto
quiere decir que para mover un flujo real, es necesario realizar trabajo sobre el flujo para vencer estos
esfuerzos y dicha energı́a se convierte en calor. Es por esto que en un flujo laminar las capas de fluido
adyancentes se mueven a velocidades diferentes en función de la transmisión de esfuerzos de corte en la
interface. Lo mismo ocurre en las fronteras solidas en donde la fricción de las paredes son transmitidas
a las capas de flujo haciendo que su velocidad aumente a medida que se alejan de las paredes. El grado
de ”pegajosidad” depende de la viscosidad del fluido. Los fluidos reales también se conocen como flujos
Newtonianos por que siguen la ley de visosidad de Newton (ver Ecuación 1).
En el caso de flujos turbulentos, flujos a velocidades altas generalmente, los esfuerzos viscosos generan
vortices en el flujo. Si a las ecuaciones de Euler se le adicionan los términos debido a los esfuerzos de corte,
se obtienen las ecuaciones de Navier-Stokes las cuales son un sistema de ecuaciones diferenciales parciales
no lineales y de segundo grado que describen el movimiento de flujos reales, compresibles o incompresibles
y permanentes o no permanentes. Las ecuaciones de Navier-Stokes para fluidos incompresibles son:

1 ∂ du
− (p + γh) + ν∇2 u = (6)
ρ ∂x dt
1 ∂ dv
− (p + γh) + ν∇2 v = (7)
ρ ∂y dt
1 ∂ dw
− (p + γh) + ν∇2 w = (8)
ρ ∂z dt
d ∂ ∂ ∂ ∂
donde ν es la viscosidad cinemática (constante para este caso), dt = u ∂x + v ∂y + w ∂z + ∂t es la derivada
∂2 ∂2 ∂2
total y ∇2 = ∂x2
+ ∂y 2
+ ∂z 2
(laplaciano). Para un fluido no viscoso, las ecuaciones de Navier-Stokes se
convierten en las ecuaciones de Euler.
Los efectos de los esfuerzos viscosos, son mas notorios en cercanı́as a las fronteras solidas (e.g fondo del
canal o paredes de una tuberı́a); dicha región es conocidad como capa lı́mite.

2 Capa limite en flujo a presión

En un flujo real la capa lı́mite es una porción de la sección de flujo en donde los esfuerzos debido a la
fricción (viscosos) están confinados o cobran gran importancia y en donde el flujo es rotacional ∇ ⃗ ×U ⃗ ̸= 0.
Esto quiere decir por fuera de la capa lı́mite la viscosidad del fluido es inoperativa y el flujo es irrotacional
⃗ ×U
∇ ⃗ = 0.
Un flujo a presión o flujo interno es aquel que viaja por un conducto y ocupa toda su sección transversal
(ver figura 1). El movimiento del flujo en el conducto se da por el gradiente de presión entre dos puntos
en el conducto separadas una distancia L. Dicho gradiente se presenta gracias a la perdida de energı́a a lo
largo de L debido a los esfuerzos de fricción y a la separación del flujo.

Figure 1: Tipos comunes de secciones transversales de tuberı́as (tomado de [2]).

Si se tiene un tanque grande del cual se conecta en la parte baja una tuberı́a (ver figura 2), los esfuerzos
viscosos empiezan a crecer una vez el fluido ingresa a la tuberı́a y por tanto la capa lı́mite (δ) también

3
empieza a crecer a lo largo de la tuberı́a. La zona inicial es una zona de flujo ideal en donde los esfuerzos
viscosos son despreciales y por lo tanto la velocidades son uniformes. Una vez, el flujo sale de zona inicial,
la capa lı́mite crece en una zona de flujo no establecido en donde se desarrollan los esfuerzos de corte.
Es posible que cuando la entrada a la tuberı́a no se hace a través de una transición suave, se presente
separación del flujo de las paredes de la tuberı́a generandose remolinos que viajan y desaparecen a lo largo
de la zona de flujo no establecido y presiones negativas a velocidades muy altas. Una vez las capas lı́mites
alrededor y crecientes en dirección del flujo se encuentran es cuando se tiene un flujo establecido o flujo
real gobernado por los esfuerzos de corte con una distribución no uniforme de velocidades.

Figure 2: Desarrollo de la capa lı́mite a lo largo de una tuberı́a (tomado de [4]).

El mecanismo de crecimiento de la capa lı́mite se puede describir como sigue. Cuando el fluido entra a
la tuberı́a se desarrollan altos valores de dv/dy, donde y es el eje vertical. Estos altos gradientes ocurren
dentro de la capa lı́mite y son debido a los esfuerzos debido a la fricción los cuales tratan de frenar el flujo.
Dicha capa crece en la dirección del flujo hasta el punto en el que se encuentran. A partir de este punto
de encuentro la acción de las fuerzas de friccion influencian el flujo y toda la sección es rotacional.
El flujo dentro de la capa lı́mite puede ser laminar o turbulento. Si el numero de Reynolds Re = Vνd ,
donde V es la velocidad media del flujo, d es el diámetro de la tuberı́a y ν es la viscosidad cinemática,
es Re < 2100, se puede inferir que el flujo laminar establecido resulta del crecimiento de la capa lı́mite
laminar. En este caso la longitud que toma el establecimiento de este flujo es xd ≈ Re 20 . Si el Re aumenta
levemente el flujo será laminar a lo largo de xd ≈ Re 20 y luego sera transisional antes de que el flujo este
establecido. Si Re >> 2100 la capa lı́mite será turbulenta. Para Re altos, en casos prácticos se puede
decir que la longitud x de la zona de flujo puede ser hasta x ≈ 100d. Sin embargo, el flujo es establecido
para valores mayores a xd ≈ 20. De acuerdo con esto, es posible notar que la energı́a en un flujo establecido
disminuye a lo largo de la tuberı́a. Note que la energı́a debido a la cabeza de presión disminuye a lo largo
de la tuberı́a debido a los esfuerzos de corte generados por la fricción dentro del flujo establecido (ver
figura 3).

Figure 3: Lı́nea de gradiente hidraáulico (LGH) y de energı́a (LE) en un a) fluido ideal y en un b) fluido
real (tomado de [2]).

3 Esfuerzo de corte y perdidas de cabeza de energı́a

Los esfuerzos de corte son producidos debido a la turbulencia del flujo o la viscosidad del fluido lo que
con lleva a una resistencia al flujo que se traduce en perdidas de energı́a. Una pregunta clave, que se
derivarı́a es ¿Cuales son los efectos de las fuerzas de fricción sobre la superficie de un volumen de control,

4
por ejemplo, en una tuberı́a?. Para esto, analizaremos los esfuerzos de corte (τ ) en un flujo 1D compresible
y permanente a través de la tuberı́a inclinada de la figura 4.

Figure 4: Fuerzas actuantes sobre el volumen de control en la tuberı́a inclinada (tomado de [4]).

Aplicando la ecuación de conservación de candidad de movimiento lineal para las fuerzas actuantes
en la dirección del flujo sobre el volumen de control entre las secciones 1 y 2 en la figura 4, tenemos que
las fuerzas fundamentales que actúan son las fuerzas de presión, las fuerzas gravitacionales y las fuerzas
viscosas. Por esto se tiene:
 
dγ dz
pA − (p + dp)A − τo P dl − γ + Adl = (V + dV )2 A(ρ + dρ) − V 2 Aρ (9)
2 dl
donde V es la velocidad del flujo, p es la presión en la sección, P es el perı́metro de la sección, A es
el área de la sección transversal, τo es el esfuerzo de corte en la superficie de control, dl es la longitud del
volumen de control, dV es un cambio en la velocidad a través de dl, dz es la diferencia de alturas entre las
secciones 1 y 2, ρ es la densidad del fluido, γ es el peso especı́fico del fluido,
 dρ esun cambio de ρ a través
de dl y dγ es el cambio del peso especı́fico a lo largo de dl. Note que γ + dγ 2 Adl dz dz
dl = dW dl , donde
dW dz dz
dl es el peso del fluido en el volumen de control en la dirección contraria del flujo, donde dl = sin θ
y θ es el angulo de inclinación de la tuberı́a. Teniendo en cuenta que entre 1 y 2 los efectos de turbinas
y bombas son despreciables, dividiendo por Aγ, donde γ = ρg y despreciando los terminos que contengan
productos de diferenciales, la ecuación 9 queda:
 2
dp V τo dl
+d + dz = − (10)
γ 2g γRh
donde Rh = PA es el radio hidráulico de la sección. Para flujo incompresible en la tuberı́a de la figura 4
y suponiendoque
 la tuberı́a es de sección constante, significa que τo no es función de l y que γ es constante
1
por lo que d γ = 0. Por lo tanto la ecuación 10 queda:

p V2
 
τo dl
d + +z =− (11)
γ 2g γRh
Integrando la ecuación 11 entre las secciones 1 y 2 (note que al integrar de esa manera el signo de la
función cambia), queda:

p1 V12 p2 V22
   
τo (l2 − l1 )
+ + z1 − + + z2 = (12)
γ 2g γ 2g γRh
Note que la diferencia de energı́a entre las secciones 1 y 2 (termino izquierdo de la ecuación 12) es la
caı́da de energı́a entre las dos secciones, por lo que la ecuación 12 se expresa como:

5
 p1 V12 p2 V22
   
+ + z1 − + + z2 = ∆(EL) = hL1−2 (13)
γ 2g γ 2g
donde ∆(LE) es la caı́da de la lı́nea de energı́a o la pérdida de energı́a entre 1 y 2 (hL1−2 ). Las pérdidas
de energı́a se pueden expresar como:

τo (l2 − l1 )
hL1−2 = (14)
γRh

Note que en la ecuación 14 quiere decir que las pérdidas de energı́a en el volumen de control son
directamente proporcionales a la longitud del volumen de control y a los esfuerzos cortantes ejercidos por
las paredes de la tuberı́a sobre las paredes del volumen de control. Las pérdidas de energı́a son además
inversamente proporcionales al radio hidráulico del volumen de control. Si la tuberı́a es de sección circular
πr2
de radio R, Rh = 2πr = 2r , la ecuación 14 quedarı́a:

2τo (l2 − l1 )
hL1−2 = (15)
γr

En términos generales, τ se puede expresar a partir de la ecuación 15:

 
γhL
τ= r (16)
2l

donde r es una distancia radial y l es la longitud de la tuberı́a. Note que τ varia linealmente con r
(ver figura 5) donde τmax = τo se logra cuando r = R (paredes de la tuberı́a). Note que las ecuaciones
anteriores fueron deducidas independiente si el flujo es laminar o turbulento.

Figure 5: Perfil de velocidades y esfuerzos de corte (tomado de [2]).

Example 1. Agua fluye en un conducto rectangular de sección 0.9 m de ancho por 0.6 m de alto. La
pérdida de cabeza de energı́a en este conducto de 60 m de longitud fue determinada experimentalmente
e igual 10 m. a) Calcular el esfuerzo de corte en las paredes del conducto. Si el conducto es de sección
circular de diámetro D = 0.6m, b) ¿cual es el esfuerzo cortante en las paredes? y c) ¿dentro del flujo en
un punto a 200 mm de las paredes?

6
4 Experimentos de Reynolds
Osborne Reynolds en 1883 mediante un experimento el cual consistió en establecer un flujo de agua a través
de una tuberı́a de vidrio en el que la velocidad era controlada por una valvula a la salida de la tuberı́a
(ver figura 6). A la entrada de la tubeı́a se inyecta una tinta que tiene un peso especı́fico igual al del agua.
Reynolds encontró que cuando la válvula está ligeramente abierta, las partı́culas de tinta se mueven de
forma ordenada formando un filamento y a manera de capas que se deslizan una sobre otra sin mezclarse.
Sin embargo, a medida que la válvula se va abriendo, se alcanza una condición en la cual la tinta presenta
un movimiento fluctuante a medida que avanza en la tuberı́a, en donde las partı́culas de la tinta se mueven
caóticamente mezclandose. Al primer tipo de flujo se le llamo laminar y al segundo turbulento.

Figure 6: Experimento de Reynolds (tomado de [3]).

Figure 7: Velocidad para flujo laminar y flujo turbulento (tomado de [3]).

Reynolds encontró que el comportamiento de flujo se podı́a correlacionar con un parametro adimensional
que relacionaba las fuerzas de inercia y las fuerzas viscosas, el cual es conocido como el numero de Reynolds
(Re):

VD 1.273Q
Re = = (17)
ν νD

donde D es el diámetro de la tuberı́a.

Reynolds encontró que flujo laminar se obtenı́a para valores de Re < 12000, mientras que el flujo turbu-
lento se lograba con Re > 50000. Sin embargo estos valores encontrados por Reynolds se obtuvieron para
condiciones alejadas de lo que es un sistema de conducción real y que ocurren comunmente en ingenierı́a.
Para propositos prácticos en tuberı́as comerciales se ha encontrado que:

Re < 2100 → flujo laminar

2100 < Re < 4000 → flujo de transición

Re > 4000 → flujo de turbulento

5 Flujo laminar
Se analiza primero el caso general de flujo permanente e incompresible entre placas paralelas e inclinadas
a un angulo θ, en donde la placa superior se mueve con una velocidad constante U (ver figura 8) y la

7
Figure 8: Flujo en movimiento entre placas paralelas con placa superior en movimiento (tomado de [5]).

velocidad del flujo es u = f (y). Si analizamos las fuerzas actuantes sobre el elemento diferencial de fluido
de ancho unitario, espesor δy y longitud δl de la figura 8, tenemos:
   
∂p ∂τ
pδy − pδy + δlδy − τ δl + τ δl + δyδl + γδlδy sin θ = 0 (18)
∂l ∂y
donde p es la presión, γ es el peso especı́fico del fluido y τ es el esfuerzo de corte. Simplificando la
ecuación 18, dividiendo por el diferencial de volumen δlδy y reemplazando sin θ = − ∂h∂l , se tiene:

∂τ ∂
= (p + γh) (19)
∂y ∂l
donde h es una distancia vertical positiva. Teniendo en cuenta que u y por lo tanto τ varian con
respecto a y únicamente, ∂τ dτ
∂y = dy . De manera similar, si (p + γh) no cambia en y (no hay aceleracion en
∂(p+γh) d(p+γh)
y), esta cambia en dirección de l, por lo tanto ∂l = dl . La ecuación queda:

dτ d
= (p + γh) (20)
dy dl
dτ 2
De la ley de viscosidad de Newton, derivando τ con respecto a y, tenemos dy = µ ddyu2 , e igualando con
la ecuacion 20, se tiene:

dτ d2 u d
= µ 2 = (p + γh) (21)
dy dy dl
Integrando la ecuación 21 con respecto a y, se tiene:

du d
µ = y (p + γh) + A (22)
dy dl
Inegrando la ecuación 22, queda:

1 d A
u= (p + γh)y 2 + y + B (23)
2µ dl µ
donde A y B son constantes de integración. Si u = 0 (velocidad en la placa fija) para y = 0 y u = U
(velocidad de la placa movil) para y = a, las constantes pueden ser encontradas resultando que B = 0 y
que:

µ ad
A=U − (p + γh) (24)
a 2 dl
Reemplazando las constantes en la ecuación 23:

8
 Uy 1 d
u= − (p + γh)(ay − y 2 ) (25)
a 2µ dl

Para el caso de placas horizontales, el gradiente debido a la presión o altura es constante p + γh = C,

por lo que, de la ecuación 25, la distribucion de velocidades es u = Uay (lineal). Para el caso en el que la
1 d 2
placa superior es fija (U = 0), u = − 2µ dl (p + γh)(ay − y ). Notese que la velocidad maxima es cuando
a µU
y= 2 − d
a dl (p+γh)
.
El caudal (Q) que pasa a través de una sección transversal (por unidad de ancho), se obtiene como:

Z a  2  a
y 1 d 2 1 d 3 Ua 1 d
Q= udy = U − (p + γh)y a + (p + γh)y = − (p + γh)a3 (26)
0 2a 4µ dl 6µ dl 0 2 12µ dl

5.1 Flujo laminar en conductos y coronas

Supongamos que consideramos una tuberı́a circular inclinada con flujo permanente (ver figura 9). Si se
toma un anillo de flujo de espesor infinitesimal δr, las fuerzas actuantes en la dirección de flujo son:

Figure 9: Diagrama de cuerpo de libre de un elemento cilińdrico para flujo laminar en un tubo circular
inclinado (tomado de [5]).

   
dp d
2πrδrp − 2πrδrp + 2πrδr δl + 2πrδlτ − 2πrδlτ + (2πrδlτ )δr + γ2πrδrδl sin θ = 0 (27)
dl dr

Teniendo en cuenta que sin θ = − dh

dl y simplificando la ecuación 27, se tiene:

dp d dh
−2πrδr
δl − (2πrδlτ )δr − γ2πrδrδl =0 (28)
dl dr dl
dividiendo por el volumen del anillo 2πrδrδl, se tiene que:

d 1 d
(p + γh) + (τ r) = 0 (29)
dl r dr
Teniendo en cuenta que d(p + γh)/dl no es una función de r y separando términos en la ecuación 29
para integrar con respecto a r:

9
 d
(p + γh)rdr + d(τ r) = 0 (30)
dl
Integrando la ecuación 30, tenemos:

r2 d
(p + γh) + τ r = A (31)
2 dl
donde A es una constante de integración. Reemplazando por τ = −µ du
dr (note que el signo menos indica
que du/dy es negativo ya que la velocidad decrece cuando r aumenta), reemplazando en la ecuación 31
tenemos:

r2 d du
(p + γh) − µ r = A (32)
2 dl dr
Separando términos para integración en la ecuación 32, se tiene:

1 d A dr
du = (p + γh)rdr − (33)
2µ dl µ r
integrando la ecuación 33, tenemos:

r2 d A
u= (p + γh) − ln r + B (34)
4µ dl µ
donde B es una constante de integración. Si se tiene una corona circular como la que se muestra en la
figura 10, en donde la velocidad para el radio interior r = b es u = 0 y donde la velocidad para el radio
exterior r = a es u = 0, las constantes de integración A y B en la ecuación 34 quedarı́an:

Figure 10: Flujo a través de un cilindro concentrico (tomado de [5]).

A a2 d
B= ln a − (p + γh) (35)
µ 4µ dl
1 d
A= (p + γh)(a2 − b2 ) (36)
4 ln(a/b) dl
Reemplazando las constantes en la ecuación 34, se tiene una expresión para el perfil de velocidades en
una corona circular:

10
 a2 − b2 a
 
1 d 2 2
u=− (p + γh) a − r + ln (37)
4µ dl ln(b/a) r

El caudal Q se calcula como:

a a
a2 − b2 a
Z Z  
πr d 2 2
Q= 2πrudr = − (p + γh) a − r + ln dr
b b 2µ dl ln(b/a) r
Z a
(a2 − b2 )r a

π d 2 3
=− (p + γh) ra − r + ln dr (38)
2µ dl b ln(b/a) r
Z a
(a2 − b2 )

π d 2 3
=− (p + γh) ra − r + (r ln a − r ln r) dr
2µ dl b ln(b/a)
R
Note que el termino r ln r en la integral de la ecuación 38 se debe integrar por partes como uv − vdu
(u = ln r y dv = r). Integrando la ecuación 38 y simplificando, se tiene:

(a2 − b2 )2
 
π d 4 4
Q=− (p + γh) a − b − (39)
8µ dl ln(a/b)

5.2 Flujo laminar en tuberı́as circulares: ecuación de Hagen-Poiseuille

En un tubo circular, en la ecuación 34, A = 0 cuando r = 0. Por lo tanto, la ecuación 34 queda:

r2 d
u= (p + γh) + B (40)
4µ dl
a d 2
Si u = 0 cuando r = a, se tiene que B = − 4µ dl (p + γh), reemplazando en la ecuación 40, se tiene:

r 2 − a2 d
u= (p + γh) (41)
4µ dl

La ecuación 41 representa la distribución de velocidades en una tuberı́a circular, cuya forma es un

paraboloide de revolución (ver figura 11).

Figure 11: Distribución de velocidades y esfuerzos de corte en una tuberı́a circular (tomado de [5]).

La velocidad máxima umax ocurre cuando r = 0, por lo tanto de la ecuación 41 se tiene:

a2 d
umax = − (p + γh) (42)
4µ dl

El caudal Q se calcula como:

11
 Z a Z a
π d
Q= 2πrudr = (p + γh) (r3 − ra2 )dr (43)
0 2µ dl 0
Integrando la ecuación 43

πa4 d
Q=− (p + γh) (44)
8µ dl

La velocidad media se define como V = Q Q

A = πa2 , que ademas es la mitad del paraboloide de revolución
de velocidades (mitad de la velocidad máxima), de acuerdo con esto:

a2 d
V =− (p + γh) (45)
8µ dl

Para un tubo horizontal de longitud L, h es constante, por lo tanto:

∆p dp
=− (46)
L dl
donde ∆p es la caı́da de presión a lo largo de L. Para el caso de un tubo horizontal, Q se expresa en
términos del diámetro D como:

∆pπD4
Q= (47)
128µL

La ecuación anterior es conocida como la ecuación de Hagen-Poiseuille. La velocidad media a partir

de la ecuación 47:

∆pD2
V = (48)
32µL

Example 2. Determinese la dirección de flujo del tubo mostrado en la figura en donde γ = 8000 N/
m3 y µ = 0.04 kg/m.s. Además, encontrar el caudal en litros por segundo y el número de Reynolds.

De la ecuación 47, despejando la caı́da de presión ∆p para una tuberı́a horizontal, la cual representa
las pérdidas por unidad de volumen he , se tiene:

128µLQ
∆p = he = (49)
πD4

12
 El factor de corrección de la energı́a cinética α para flujo laminar en una tuberı́a se calcula como:
Z  
1 u 3
α= dA (50)
A V
Reemplazando las ecuaciones 41 y 45 en la ecuación 50, tenemos:

1
Z a   r 2 3
α= 2 1− 2πrdr = 2 (51)
πa2 0 a
Notese que la energı́a cinética es el doble de la energı́a si se considera una distribución real (u = f (r))
de velocidades. El factor de correccı́on de la cantidad de movimiento β, se expresa como:
Z  
1 u 2
β= dA (52)
A V
Reemplazando las ecuaciones 41 y 45 en la ecuación 52, tenemos que β = 4/3.

6 Flujo turbulento
6.1 Esfuerzo de corte turbulento
En un flujo turbulento, las partı́culas de agua viajan de manera desordenada siguiente trayectorias caóticas.
Esto implica que las velocidades y las presiones fluctuen rápidamente en el tiempo. Esto implica que el
análisis de las ecuaciones de Navier-Stokes sea dificil analı́ticamente, e incluso, numéricamente. Para esto
es necesario expresar dichas cantidades en sus valores medios en el tiempo y en sus cantidades fluctuantes.
Para la componente de la velocidad en x, u, esta se puede representar:

Figure 12: Velocidad en x para flujo turbulento (tomado de [5]).

u = u + u′ (53)
donde u es la velocidad media o promedio y u′ es la fluctuación. Para un tiempo T , u se expresa:
Z T
1
u= udt (54)
T 0
y el promedio de las fluctuaciones de u en T , se expresa como:
Z T
1
u′ = (u − u)dt = 0 (55)
T 0
Nótese que el promedio de u′ es cero para T . Sin embargo el cuadrado de las fluctuaciones:

13
 Z T
1
u′ 2 = (u − u)2 dt ̸= 0 (56)
T 0
es diferente de cero. Reynolds descompuso cada propiedad del flujo en su valor medio y su fluctuación:

v = v + v′ w = w + w′ p = p + p′ (57)

Note que la raiz cuadrada de u′ 2 es una medida de la intensidad de la turbulencia. Tampoco es cero
los productos medios de las fluctuaciones como u′ v ′ , u′ w′ , etc.
De las ecuaciones de Navier-Stokes (ecuaciones 6, 7 y 8), la ecuación de cantidad de movimiento en x
promediada en el tiempo, la cual contiene el producto de las componentes de fluctuación, se expresa como:
     
∂ ∂ ∂u ′ 2 ∂ ∂u ′ ′
∂ ∂u ′ ′
∂u
− (p + γh) + µ − ρu + µ − ρu v + µ − ρu w = ρ (58)
∂x ∂x ∂x ∂y ∂y ∂z ∂z ∂t

Los términos −ρu′ 2 , −ρu′ v ′ y −ρu′ w′ , representan la aceleración convectiva y matemáticamente son
análogos al esfuerzo. Estos se denominan los esfuerzos de Reynolds. Note que los términos (e.g. µ ∂u ∂x ) son
los esfuerzos cortantes medios viscosos. Los esfuerzos de Reynolds son los responsables del intercambio
de cantidad de movimiento y la capacidad de mezcla en un flujo turbulento. En un flujo turbulento los
esfuerzos de Reynolds son superiores a los esfuerzos viscosos, mas aun, en la capa lı́mite. Los esfuerzos de
Reynolds se calculan empiricamente con ayuda de experimentos.
Para flujo en la dirección x, el esfuerzo turbulento mas importante es −ρu′ v ′ , por lo que la ecuación 58
se convierte:
∂ ∂τ ∂u
− (p + γh) + ≈ρ (59)
∂x ∂y ∂t
donde τ = τl + τt = µ ∂u ′ ′
∂y − ρu v es el esfuerzo de corte total, τl es el eufuerzo laminar y τt es el esfuerzo
cortante turbulento.

6.2 Longitud de mezcla de Prandtl

El esfuerzo cortante aparente para flujo turbulento, incluyendo los efectos viscosos, se expresa como:

du
τ = (µ + η) (60)
dy
Prandtl encontró expresiones para determinar u′ y v ′ en función de la longitud de mezcla l y el gradiente
de velocidad du
dy en donde u es la velocidad media temporal en un punto y y es la distancia normal a u
medida desde la frontera (ver figura 13)

Figure 13: Esquema de la teorı́a de longitud de mezcla (tomado de [5]).

De la figura 13, l y u′ se relacionan como:

14
 du
u′ ∼ l (61)
dy
Esto significa que los cambios de la velocidad dependen de los cambios en la velocidad media temporal
en dos punto separados por l, donde l es el tamaño promedio de los remolinos responsables de la mezcla. l
es tambien la distancia recorrida por una partı́cula antes de que su cantidad de movimiento sea modificada.
Teniendo en cuenta que u′ y v ′ están correlacionados:

du
v ′ ∼ u′ ∼ l (62)
dy
Reemplazando las expresiones de la ecuación 62 en τt = −ρu′ v ′ (esfuerzo turbulento), se obtiene una
ecuación que define la longitud de mezcla l:
 2
2 du
τt = ρl (63)
dl
Note que l absorbe el factor de proporcionalidad y el signo.
Para flujo turbulento, se puede encontrar una expresión similar a la ley de viscosidad de Newton como:

du
τt = η (64)
dy
donde η es la viscosidad de remolino que no es solo una propiedad del fluido ya que depende del
movimiento de este y de su densidad, y generalmente es mayor µ. η también se considera como un
coeficiente de transferencia de cantidad de movimiento. Igualando las ecuaciones 63 y 64, se tiene:

du
η = ρl2 (65)
dy
Note que l cambia con respecto a la distancia a la pared y y se hace cero en la frontera. Von Karman
a través del estudio de relaciones de similitud en flujo turbulento, encontró una relación para l como:

du/dy
l=κ (66)
d2 u/dy 2
donde κ es una constante universal (constante de Von Karman) en un flujo turbulento sin importar la
configuración de la frontera o el número de Reynolds. Se ha encontrado que κ ≈ 0.40.

6.3 Distribuciones de velocidad

En flujos de turbulentos, la distribución de velocidades en las zonas cercanas a la pared del conducto se
puede dividir en tres capas (ver Figura 14): Capa viscosa (o laminar), capa de traslape y zona turbulenta
exterior.

Capa laminar
En la capa viscosa mas cercana a la pared, el esfuerzo cortante es constante e igual al esfuerzo en la pared
τ0 . La distribución de velocidades en la capa viscosa se relaciona con la viscosidad absoluta y el esfuerzo
cortante a partir de la ley de viscosidad de Newton y se expresa como:
τ0 µu u
= =ν y ≤ δ′ (67)
ρ ρy y
donde δ ′ es el espesor de la capa viscosa de la pared.
Definiendo la velocidad de esfuerzo de corte (u∗ ) como u∗ = τ0 /ρ y reemplazando en la ecuación 67,
p

se tiene:

15
Figure 14: Para flujo turbulento en conductos a) distribución de esfuerzos y b) distribución de velocidades
(tomado de [5]).

u
u∗ 2 = ν (68)
y
dividiento por u e invirtiendo los términos en la ecuación 68, se tiene

u u∗ y
= y ≤ δ′ (69)
u∗ ν

La ecuación 69 muestra una relación lineal entre u y y en la capa viscosa (o laminar) y es conocida
como la ley de pared. A través de experimentos, se ha encontrado que δ ′ :

ν
δ′ = 5 (70)
u∗

En la practica, se ha encontrado que una superficie es hidráulicamente lisa si δ ′ > 6ε e hidráulicamente

rugosa si δ ′ < 3ε, donde ε es la rugosidad en unidades de L del material del conducto.

Capa de traslape
En la capa de traslape (ver figura 14), los esfuerzos turbulentos son dominantes mientras que el esfuerzo
laminar puede ser despreciable. Por lo tanto el esfuerzo cortante turbulento τt (Ecuación 63) es aproxi-
madamente igual al esfuerzo cortante en la pared τ0
 2
2 du
τ0 = ρl (71)
dy
Von Karman encontró que la longitud de mezcla l en p
cercanı́a a la pared es proporcional a la distancia
desde la pared y e igual l ≈ κy. Reemplazando l y u∗ = τ0 /ρ en la ecuación 71, se tiene que:

du 1 dy
= (72)
u∗ κ y
integrando la ecuación 72:
u 1
= ln y + A (73)
u∗ κ
donde A es una constante de integración. Si se reemplaza la expresión para u de la ecuación 73 en la
2
ecuación 66, se tiene que l es proporcional a y (note que ddyu2 es negativo por que el gradiente de velocidad
disminuye en la medida que y aumenta).

16
 Para y = δ ′ se tiene que u = uw que es la velocidad en donde el flujo cambia de laminar a turbulento.
Reemplazando en la ecuación 69, se tiene:

uw u∗ δ ′
= =N (74)
u∗ ν
donde N es un número de Reynolds crı́tico cuando el flujo cambia de laminar a turbulento. Sustituyendo
u = uw y y = δ ′ en la ecuación 73 para determinar el valor de A, se tiene que:
 
uw 1 ′ 1 Nν 1
A = ∗ − ln δ = N − ln ∗
= N − (ln N + ln ν − ln u∗ ) (75)
u κ κ u κ
Reemplazando la ecuación 75 en la ecuación 73, se tiene:

yu∗
 
u 1 1
∗
= ln +N − ln N (76)
u κ ν κ

Nikuradse mediante experimentos de flujo turbulento para conductos con paredes lisas y graficando
u/u∗ vs ln(yu∗ /ν) encontró que la constante universal adimensional κ = 0.40 y N − κ1 ln N = 5.5.

Zona turbulenta
En la zona turbulenta de la figura 14 para tuberı́as de radio r0 y reemplazando y = r0 , para el cual la
velocidad u = um (en el centro de la tuberı́a), en la ecuación 73, se tiene que la constante es:
um 1
A= ∗
− ln r0 (77)
u κ
Reemplazando en la ecuación 73:

um − u 1 r0
∗
= ln (78)
u κ y

La ecuación 78 es conocida como la ley del deficit de velocidad y se aplica a conductos lisos y rugosos.
Para tubos rugosos, yw = mϵ′ , donde ϵ′ es una altura tı́pica de la rugosidad y m es un coeficiente de
forma que depende de la naturaleza de la rugosidad, y la velocidad es uw (velocidad donde el flujo cambia
de laminar a turbulento), determinando el valor de la constante en la ecuación 73, se tiene:
uw 1
A= ∗
− (ln m + ln ϵ′ ) (79)
u κ
reemplazando la ecuación 79 en la ecuación 73, se tiene:

u 1 y uw 1
= ln ′ + ∗ − ln m (80)
u∗ κ ϵ u κ

donde el termino uuw∗ − κ1 ln m es una constante que depende del tipo de rugosidad. Nikuradse realizando
experimentos con tuberı́as rugosas (granos de arena aheridos a las paredes de la tuberı́a) y considerando
que la rugosidad de la tuberı́a ϵ′ es igual al diámetro del grano de arena, encontró que κ = 0.40 y
uw 1
u∗ − κ ln m = 8.48.
Note que la ley logarı́tmica de la ecuación 76 es la que tiene la aplicación mas amplia. Fı́jese que la
capa viscosa es muy pequeña en flujos turbulentos.
Prandtl desarrolló una formula sencilla para la distribución exponencial de velocidad para flujo turbu-
lento en tuberı́as:

17
  n
u y
= (81)
um r0

donde n varia con el número de Reynolds. Esta ecuación empı́rica es valida solo para ciertas distancias
de la pared. Para R < 100000, n = 1/7, para valores mayores de R, n decrece. Ambas ecuaciones 81 y
76 tienen la falla de predecir un valor du/dy ̸= 0 en el centro de la tuberı́a.

Example 3. Encuentrese una expresión aproximada para la distribución de longitud de mezclado en flujo
turbulento en un tubo a partir de la ley exponencial de velocidad de Prandtl para n = 1/7.

Example 4. En un tubo circular de 14 cm diámetro fluye aire a una temperatura de T = 20 o C. Si el

flujo esta completamente desarrollado, y la velocidad en el centro de la tuberı́a es 5 m/s, determinar a) la
velocidad de fricción u∗ y b) el esfuerzo de corte en la pared. Asumir flujo turbulento y la ocurrencia de la
ley logarı́tmica.

6.4 Pérdidas de energı́a

En flujos incompresibles y turbulentos a regimen permanente y uniforme en conductos cerrados de sección
transversal constante, el esfuerzo cortante en la pared τ0 se expresa como:

ρ
τ0 = λ V 2 (82)
2

donde λ es un coeficiente adimensional y V es la velocidad media. En conductos abiertos o cerrados no

circulares, el esfuerzo cortante no es constante sobre la superficie por lo que τ0 se calcula como el promedio
de los esfuerzos cortantes sobre la pared.
Si se analizan las fuerzas actuantes sobre un volumen de control de flujo en un conducto abierto o
cerrado (ver figura 15), la energı́a del flujo podrı́a proporcionarse por la caı́da de energı́a potencial , ası́
como por una caı́da en la presión p1 − p2 .

Figure 15: Fuerzas axiales sobre un volumen de control en un conducto (tomado de [5]).

Las pérdidas de energı́a entre dos secciones 1 y 2 en conductos se puede expresar usando la ecuación
de energı́a a partir de la ecuación de la ecuación de Bernoulli (ecuación 2):

p1 V12 p2 V22
+ + z1 = + + z2 + h1−2 (83)
γ 2g γ 2g

18
 donde h1−2 son las pérdidas de energı́a entre 1 y 2. Como la sección transversal en el conducto es
2
constante, el termino V2g se elimina a ambos lados, por lo que las perdidas se expresan como:
p1 − p2
h1−2 = + z 1 − z2 (84)
γ
Note que las pérdidas son proporcionales a los cambios en la cabeza de presión y en la energı́a potencial.
De acuerdo con la ecuación de conservación de la cantidad de movimiento lineal para el volumen de
control de la figura 15 en la dirección del flujo l, tenemos:
X
Fl = 0 = (p1 − p2 )A + γAL sin θ − τ0 LP (85)
donde L es la distancia entre 1 y 2, A es el área de la sección transversal en 1 y en 2, P es el perı́metro
mojado de la sección, el cual es la porción de la sección transversal en contacto con el fluido excluyendo la
superficie libre. Teniendo en cuenta que L sin θ = z1 − z2 , se tiene:

p1 − p2 τ0 LP
+ z1 − z2 = (86)
γ γA
Igualando las ecuaciones 84 y 86 y usando la ecuación 82, se tiene:

τ0 LP ρ LP LV2
h1−2 = =λ V2 =λ (87)
γA 2 γA R 2g
donde R es el radio hidráulico del conducto R = PA , para una tuberı́a R = D 4 donde D es el diámetro
del tubo. Note que el las perdidas de energı́a h1−2 tiene unidades de Newton-metro por newton o libras-pie
por libra. Teniendo en cuenta que las pérdidas de cabeza de energı́a son debido a la fricción, h1−2 = hf .
Si se expresan hf en términos de longitud, se tiene que la pendiente de la lı́nea de energı́a S se expresa:

hf λ V2
S= = (88)
L R 2g
Despejando la velocidad V de la ecuación 88:

2g √ √
r
V = RS = C RS (89)
λ

donde C es un coeficiente de fricción que depende de la rugosidad del material y del tamaño del conducto
y se encuentra de manera experimental. La ecuación 89 es conocida como la ecuación de Chezy y C es
el coeficiente de Chezy. Existen diferentes formulas para calcular C, para tuberı́as se tiene que λ = f /4,
reemplazando en la ecuación 87, se tiene:

L V2
hf = f (90)
D 2g

donde f es el factor de fricción que se ha determinado experimentalmente para tuberı́as. Esta ecuación
es aplicable a conductos abiertos (canales) de la siguiente forma:

8g √
r
V = RS (91)
f

19
6.5 Cálculo del factor de fricción
Al hacer un analisis del flujo turbulento, las variables que influyen en su comportamiento son: la velocidad
media del flujo V [LT −1 ], la viscosidad dinámica µ [M L−1 T −1 ] o cinemática ν [L2 T −1 ], la densidad del
fluido ρ [M L−3 ], el esfuerzo de corte en la pared τ0 [M L−1 T −2 ], una longitud caracterı́stica Lc que para el
caso especı́fico de flujo en tuberı́as de sección transversal circular es igual a D [L] y la rugosidad absoluta
del tubo ε [L]. De acuerdo con esto, podemos decir que en el flujo turbulento intervienen 6 variables y 3
dimensiones fundamentales: masa M , longitud L y tiempo T . Utilizando el teorema Pi de Buckingham y el
analisis dimensional se puede obtener una ecuación que relacione estas 6 variables: f (τ0 , ρ, µ, ν, D, V, ε) = 0.
De acuerdo con el teorema, si el numero de variables n = 6 y el numero de dimensiones m = 3, el número
de parametros adimensionales es n − m = 3. Aplicando el teorema y tomando como variables repetitı́vas
ε, ρ y V , los parámetros adimensionales son:

ρV D ρV 2 ε
Π1 = = Re Π2 = =E Π3 = (92)
µ τ0 D
2
donde E es el número de Euler. De acuerdo con lo anterior, f (Re, E, ε/D) = 0 o ρV τ0 = f (Re, ε/D).
Despejando τ0 , se tiene:
h  ε i
τ0 = ρV 2 f Re, (93)
D

−1
Note que al despejar f Re, Dε es equivalente a f Re, Dε ya que f es función de variables adimen-


h γR h ρgD
sionales. De la ecuación 87, τ0 = fL = f4L igualando esta expresión a la ecuación 93 y simplificando,
se tiene:

L V2 h  ε i
hf = 8f Re, (94)
D 2g D
Note que el término entre paréntesis en la ecuación 94 es una función de Re y ε/D, es el factor de
fricción (ver ecuación 90), por lo que la ecuación 94 se convierte en:

L V2
hf = f (95)
D 2g

la cual es conocida como la ecuación de Darcy-Weisbach. Esta ecuación se puede aplicar al flujo laminar
y al flujo turbulento, sin embargo, el factor de fricción en flujo turbulento depende no solo de Re si no de
la rugosidad ε del contorno.
Existen diferentes ecuaciones para definir el valor del coeficiente de fricción dependiendo de si la su-
perficie del conducto es lisa o rugosa. Estas ecuaciones han sido derivadas de análisis experimentales por
diferentes investigadores como Blassius, Prandtl, Von Karman, Nikuradse, Jain, etc.
Blasius encontró que, para superficies hidráulicamente lisas y flujo turbulento en el rango 4000 < Re <
100000, f se calcula como:

0.316
f= (96)
Re0.25

Nikuradse realizó experimentos usando tres tuberı́as con diferente diámetro a las cuales les adherı́a
uniformemente granos de arena de diferente diámetro (ε). Nikuradse encontró (ver figura 16) que cuando
la rugosidad relativa ε/D variaba desde 0.001 hasta 0.033, f presentaba un comportamiento complejo que
dependı́a del número de Reynolds y de la condición hidrodinámica del contorno. De los experimentos de
Nikuradse, se encontró:

20
Figure 16: Experimentos de Nikuradse con tubos con rugosidades formadas con diferentes diámetros de
arena (tomado de [5]).

ˆ Para flujo laminare (Re < 2100) los datos experimentales definen una lı́nea que satisface la ecuación 95,
donde:

64
f= (97)
Re

ˆ Para flujo turbulento (Re > 4000), los datos experimentales muestran una serie de curvas (ver
figura 16) definidas por ε/D, por lo que f es una función de Re y de la rugosidad del tubo. En esas
curvas se pueden definir tres zonas:

1. Una primera zona definida por una curva envolvente para superficies lisas (δ > ε), por lo que f
es independiente de la rugosidad.
2. Una segunda zona en donde las curvas de f vs Re son paralelas para diferentes valores de ε/D,
lo que evidencia que f es independiente de Re; si Re aumenta el valor de f permanece constante.
Esta zona define las superficies hidráulicamente rugosas y es conocida como la zona totalmente
rugosa. El espesor de la capa viscosa δ ′ disminuye a medida que incrementa Re.
3. Una tercera zona que define las superficies lisas (envolvente inferior). Se conoce como zona
rugosa de transición.

21
 Teniendo en cuenta el comportamiento de f , Prandtl y Von Karman presentaron las siguientes ecua-
ciones las cuales se ajustan muy bien a los datos experimentales obtenidos por Nikuradse:

ˆ Para Re > 4000 y tuberı́as lisas

1  p 
√ = 0.869 ln Re f − 0.8 (98)
f

ˆ Para Re > 4000 y tuberı́as rugosas

1 ε
√ = 1.14 − 0.869 ln (99)
f D

Sin embargo, teniendo en cuenta que los experimentos de Nikuradse se realizaron para tuberı́as con
granos de arena adheridos uniformemente a las paredes de la tuberı́a, dicha condición dista mucho de
la realidad en donde ni la rugosidad ni el diámetro son uniformes. Para estudiar el comportamiento en
tuberı́as comerciales en donde la rugosidad es irregular, C.F. Colebrook(1939), realizó experimentos con
el propósito de mostrar la aplicabilidad de las ecuaciones 98 y 99. Colebrook obtuvo que en la zona de
transición, existe un efecto de la no uniformidad de la rugosidad ε. Para tuberı́as comerciales de sección
circular, Colebrook obtuvo una ecuación implicita para f para la zona de transición:
 
1 ε 2.52
√ = −2 log + √ (100)
f 3.7D Re f

QD 4Q
En términos de caudal Q, Re = Aν = πDν , la ecuación 100 se expresa como:
  −2
ε 1.979νD
f = −2 log + √ (101)
3.7D Q f

L.F. Moody presentó un diagrama (ver figura 18) para determinar el factor de fricción f en tuberı́as
comerciales para flujos laminares y turbulento con base en las ecuaciones 96, 98, 99 y 100. El diagrama
de Moody sirve para determinar f en función de Re y ε/D. EL diagrama presenta una serie de curvas que
definen el comportamiento de flujo: laminar (Re < 2100, lı́nea recta con pendiente -1) de transición y turbu-
lento (Re > 4000). La zona de flujo turbulento esta compuesta por tres zonas: superficies hidráulicamente
lisas en donde f depende solamente de Re, superficies en transición de lisas a rugosas donde f depende de
Re y de ε/D y superficie hidráulicamente rugosa donde f depende de ε/D solamente.
El diagrama de Moody támbien puede ser utilizado en flujo a superficie libre (e.g. canales) utilizando
el radio hidráulico R en lugar del diámetro D.
En 1976, Swanee y Jain presentaron una ecuación explicita para el factor f como:
  −2
ε 5.74
f = 0.25 log + (102)
3.7D Re0.9

Esta ecuación es valida para los siguientes rangos de número de Reynolds y de rugosidad relativa ε/D:
4000 < Re < 1x108 y 0.01 < ε/D < 1x10−4 .

22
Figure 17: Diagrama de Moody para tuberı́as comerciales con flujo laminar y turbulento (tomado de [?]).

Example 5. Fluye petroleo (S = 0.8, µ = 0.0052 Pa.s) a través de una tuberı́a de 10 cm de diámetro
y rugosidad ε = 0.85 mm a razón de 40 l/s. Determinar la disipación de energı́a por unidad de longitud
en la tuberı́a y el esfuerzo de corte en la pared del conducto. Encontrar la magnitud de la velocidad a
una distancia radial de 2 cm desde el centro de la tuberı́a. Definir si la superficie se comporta como
hidraulicamente lisa, en transición o rugosa.

Example 6. Agua a 60 o F (ρ = 62.36 lbm/ft3 y µ = 7.536x10−4 lbm/ft.s ) fluye a través de una tuberı́a de
acero (ε = 7x10−6 ft) horizontal de 2 in de diametro y longitud 200 ft con un caudal 0.2 ft3 /s. Determinar
la caida de la cabeza de presión, la perdida de energı́a y la potencia requerida para bombear agua a través
de la tuberı́a.

23
 Figure 18: Rugosidades (ε) para tuberı́as comerciales (tomado de [?]).

6.6 Otras ecuaciones para calcular las pérdidas de energı́a

Ecuación de Hazen-Williams
Hazen y Williams presentaron la siguiente expresión para la velocidad media del agua en régimen turbu-
lento:

V = 0.85CH R0.63 Sf0.54 (103)

donde V es la velocidad media en m/s, CH es un coeficiente de rugosidad que depende del material de
la tuberı́a (ver tabla 1), R es el radio hidráulico en m y Sf es la pendiente de la lı́nea de energı́a Sf = hf /L,
donde L es la longitud en m.
Teniendo en cuenta que para tuberı́as de sección circular donde R = D/4, la ecuación 103 se puede
expresar en términos de hf como:
 1.85
Q
hf = 10.654L (104)
CH D2.63
donde Q es el caudal a través del conducto en m3 /s.

Ecuación de Manning
Robert Manning en 1880 obtuvo una expresión en el sistema Inglés de unidades para determinar la velocidad
media del agua en tuberı́as:

24
 Material del conducto CH
PVC 150
Fundición asfaltada 140
Eternit - Asbesto - Cemento 140
Acero 130−140
Hierro Forjado 130−140
Fundición 130
Hormigón 120
Acero liso 120
Madera 120
Fibra de vidrio 110

Table 1: Coeficientes de rugosidad de Hazen-Williams, CH (tomado de [2]).

1.49 2 12
V = R 3 Sf (105)
η

donde η es un coeficiente de rugosidad que depende del material del conducto y de las propiedades
hidráulicas del flujo. Despejando de la ecuación 105 las pérdida de energı́a hf y expresandola en sistema
internacional, se tiene:

10.29Lη 2 Q2
hf = (106)
D16/3
Al observar las ecuacione 104 y 106, se puede afirmar que la disipación o pérdida de energı́a se puede
expresar de manera general como:

hf = CQn (107)
donde C es una constante que depende de un sistema de unidades, longitud, diámetro y rugosidad del
conducto, y el exponente n varia entre 1.0 y 2.0 depende del régimen de flujo. La tabla 2 define el valor
de C y de n en la ecuación 107 para el sistema Internacional y para el sistema Ingés de unidades.

Sistema Internacional Sistema Ingés

L(m), D(m), Q(m3 /s) L(pie), D(pie), Q(pie3 /s)
Ecuación de Darcy-Weisbach (n=2) C = 0.0827f DL5 C = 0.0252f DL5
h i1.85 h i1.85
Ecuación de Hazen-Williams (n=1.85) C = 10.654L CH D1 2.63 C = 4.72L CH D1 2.63

Table 2: Ecuaciones para la disipación de energı́a hf = CQn (tomado de [2]).

7 Pérdidas menores
En un sistema de tuberı́as, el flujo pasa a través de múltiples entradas, salidas, uniones, válvulas, codos,
bifurcaciones, expansiones, contracciones, etc (ver figura 19) en adición de los tramos rectos de las tu-
berı́as. Estos componentes del sistema causan separación y mezcla del flujo induciendo pérdidas de energı́a
adicional. Comparadas con las perdidas por fricción, hf , las pérdidas debido a estos accesorios, he , son
menores. Sin embargo, en algunos casos en donde existen muchos cambios de dirección y válvulas en un
tramo corto de tuberı́a, he llega a ser mayor que hf . Cuando una válvula esta totalmente abierta la pérdidas
de energı́a a través de esta son despreciable. Sin embargo cuando la válvula esta parcialmente abierta,

25
existen pérdidas de energı́a debido a la disminución de caudal, como consequencia de la disminución de
presión, a través de la válvula.

Figure 19: Accesorios de uso común en tuberı́as (tomado de [2]).

Hay muchas clases de válvulas usadas comercialmente (ver figura 20): la válvula de compuerta la cual
se desliza hacia abajo en la sección, la válvula de globo la cual obtura un agujero usando un mecanismo
especial, la válvula de angulo la cual es similar a la de globo pero sobre un angulo de 90o , la válvula de cheque
la cual permite flujo en una sola dirección y la válvula de disco la disco la cual cierra la sección transversal
con una compuerta circular. La válvula de globo debido a su configuración interna es usualmente la que
genera mayores pérdidas cuando esta totalmente abierta. Existen también las válvula de mariposa la cual
consiste en un disco que gita y obtura la sección transversal dependiendo del angulo de rotación (ver
figuraacce01). Las pérdidas en válvulas están datas por el fabricante y dependen de su grado abertura. La
figura 22 muestra una curva del coeficiente de pérdidas K en función de la relación abertura (h) diametro
(D).
La cuantificación de las pérdidas menores es complejo desde el punto de vista teórico. Por esto, estas
pérdidas son determinadas experimentalmente por los fabricantes de los accesorios. Las pérdidas menores,
he , son expresadas usualmente como:

V2
he = K (108)
2g

donde K es el coeficiente de pérdida o de resistencia que depende del número de Reynolds, del material
del cual esta hecho el accesorio y de la forma como se acopla el accesorio. Sin embargo, para Re > 105 se
ha demostrado que K es independiente del número de Reynolds.
Vennard-Street explican fı́sicamente el efecto de un accesorio sobre la energı́a disponible en el sistema
y sobre el flujo. Dichos efectos se analizando mediante la figura 23:

ˆ Zona BC: Aguas arriba del accesorio, la vena lı́quida se contrae y el flujo se acelera. La máxima
contracción de la vena liquida se presenta justo en la posición del accesorio o un poco aguas abajo
de este. Debido a la contracción, la presión disminuye.

ˆ Zona DE: En esta zona la presión del flujo aumenta ya que la velocidad del disminuye. El flujo se
desacelera generando la creación de remolinos que ocacionan una turbulencia de gran escala. Por lo
tanto una parte de la energı́a se pierde debido a la creación de remolinos.

26
Figure 20: Válvulas comerciales tı́picas: a) válvula de compuerta, b) válvula de globo, c) válvula de angulo,
d) válvula de cheque y e) válvula de disco (tomado de [6]).

ˆ Zona EF: En esta zona se disipan (desaparecen) los remolinos. Aguas abajo de esta zona se restablece
la condición de flujo desarrollado.

Note que a lo largo de la longitud AG (ver figura 23) támbien actúan las fuerzas de fricción. Dichas
pérdidas debido a la fricción se calculan asumiendo flujo desarrollado a lo largo de AG. Las pérdidas
globales a lo largo de AG son entonces la suma de las pérdidas debido al accesorio (pérdidas locales entre C
y D) y debido a la fricción (a lo largo de AG). Las pérdida total de energı́a en un sistema puede expresarse
como:

X X X Li Vi2 X Vj 2
hT = hf i + hej = fi + Kj (109)
Di 2g 2g
i j i j

donde i representa cada tuberı́a con diámetro constante y j representa cada componente que causa una
pérdida menor. Si el sistema analizado tiene diámetro constante y del mismo material, la ecuación 109, se
convierte en:
 
L X V2
hT = f + Kj  (110)
D 2g
j

donde V es la velocidad media en todo el sistema.

Los sistemas de tuberı́as comunmente contienen contracciones o expansiones subitas o graduales de la
sección de la tuberı́a con el fin de acomodar cambios de caudal o velocidad o de propiedades del fluido (e.g.
densidad). Las pérdidas son usualmente mayores cuando los cambios son súbitos o con un ángulo grande
debido a la separación del flujo (ver figura 24).

27
 Figure 21: Válvula de mariposa (tomado de [6]).

Aplicando la ley de conservación de cantidad de movimiento para el volumen de control en la figura 24

entre las sección 1 y 2, tenemos:

p1 A2 − p2 A2 = ρV2 (V2 A2 ) + ρV1 (−V1 A1 ) (111)

p1 −p2
teniendo en cuenta que de acuerdo con la ley de continuidad V1 A1 = V2 A2 , y despejando γ de la
ecuación 111, se tiene que:

p1 − p2 V2 2 − V2 V1
= (112)
γ g
Aplicando la ley de conservación de la energı́a entre las secciones 1 y 2 en la figura 24, se tiene:

V1 2 p1 V2 2 p2
+ = + + he (113)
2g γ 2g γ
p1 −p2
despejando γ de la ecuación 113 e igualando a la ecuación 112, se tiene:

V2 2 − V2 V1 V2 2 − V1 2
= + he (114)
g 2g

28
 Figure 22: Coeficientes K promedio para válvulas parcialmente abiertas. (tomado de [6]).

Figure 23: Comportamiento del flujo en una tuberı́a con accesorio (tomado de [2]).

despejando he de la ecuación 114 y usando la ecuación de continuidad, se tiene:

(V1 − V2 )2 V1 2 A1 2
 
he = = 1− (115)
2g 2g A2

La ecuación 115 determina la pérdidas locales en una expansión brusca, por lo que:
 2
A1
K= 1− (116)
A2

Note que la ecuación 115 representa la deducción de una expresión para calcular las pérdidas menores
en términos generales teniendo en cuenta que afirma que las pérdidas menores varian con el cuadrado de
la velocidad afectada por un coeficiente K, que como lo habı́amos afirmado, esta dado generalmente por
el fabricante.

29
 Figure 24: Representación de la pérdida de energı́a en una expansión subita. (tomado de [5]).

Analizando la ecuación 116, en el caso en que la tuberı́a descargue a un deposito (A1 << A2 ), K ≈ 1
lo que significa que la energı́a cinética del flujo se convierte en energı́a términa (calor).
Para el caso de la contracción brusca que se muestra en la figura 25, se hace un análisis similar al
realizado para el caso de la expansión brusca. Partiendo de la ecuación 115, se tiene:

Figure 25: Contracción brusca en una tuberı́a (tomado de [5]).

(V0 − V2 )2
he = (117)
2g
Teniendo en cuenta la vena lı́quida contraı́da en la sección 0, la ecuación de continuidad entre 0 y 2
es V0 Cc A2 = V2 A2 en donde Cc es un coeficiente de contracción que fue determinado inicialmente por
Weisbach con respecto a la relación A2 /A1 . Reemplazando en la ecuación 117, se tiene:
2
V2 2 V2 2

1
he = −1 =K (118)
Cc 2g 2g

En las figuras 26 y 27 se presentan valore de K para diferentes tipos de accesorios.

7.1 Longitud equivalente

Las pérdidas menores se pueden expresar en términos de la longitud equivalente Le de tubo con la misma
pérdida de cabeza para el mismo caudal ası́:

Le V 2 V2
f =K (119)
D 2g 2g
en donde K puede referirse a una pérdida de carga menor o a la suma de varias pérdidas. Al despejar
Le , se tiene:

30
 KD
Le = (120)
f

Por ejemplo, si las pérdidas menores en una tuberı́a de 12 pulg = 1 f t se suman y da K = 20, y si
f = 0.020 para la tuberı́a, a la longitud real de la tuberı́a se puede sumar 20(1/0.020) = 1000f t y esta
longitud adicional o equivalente causa la misma resistencia al flujo que las pérdidas menores.

Example 7. Encuéntrese el caudal que fluye por la tuberı́a en la figura con H = 10 m, y las pérdidas de

energı́a para un caudal de 60 L/s.

Example 8. Por un sistema de conducción de un modelo de laboratorio, fluye agua con un caudal de
0.0014 m3 /s. El modelo consiste en una tuberı́a de 25 mm de diámetro que se acopla, a través de una
reducción brusca, a una tuberı́a de 19 mm de diámetro. Para el caudal ya definido, la lı́nea resultante de
gradiente hidráulico presenta los valores ilustrados en la figura. Determinar el coeficiente K del accesorio

reductor.

8 Flujo normal a superficie libre

El flujo uniforme ocurre en largos y rectos canales con una pendiente y un ancho constante (ver figure 28). Si
consideramos una sección 1 aguas arriba y otra sección 2 aguas abajo en un tramo de canal, la profundidad
es constante (y1 = y2 ) y se conoce como la rofundidad normal (y1 = y2 = yn ). Esto implica que la velocidad
V1 = V2 = Vo es también constante. La pendiente del canal se puede calcular como So = tan θ, donde θ es
el angulo de inclinación del canal, considerado positivo para un flujo descendente. De acuerdo con esto, de
la ecuación de Bernoulli, tenemos:

hf = z1 − z2 = So L (121)
donde L es la longitud entre 1 y 2. La ecuación 121 indica que la pendiente de la linea de energı́a
Sf = hf /L es la misma pendiente del canal So . Esto quiere decir ademas que la perdida de energı́a es
únicamente potencial.

31
 Teniendo en cuenta que el flujo esta totalmente desarrollado, la ecuación de Darcy-Weisbach para el
calculo de pérdidas por fricción en tuberı́as:

L V2
hf = f (122)
D 2g
Para el caso de canales se debe tener en cuenta que para tuberı́as Rh = D/4, D = 4Rh = 4A/P , donde
Rh es el radio hidráulico, A es el area de la sección transversal del canal y P es el perimetro mojado del
canal (ver figura 29).
Combinando las ecuacione 121 y 122, se puede obtener una expresión para la velocidad en un canal
con flujo uniforme:
 1/2
8g 1/2
Vo = Rh So1/2 (123)
f
 1/2
8g
De la ecuación anterior, el termino f es constante en canales y depende de la rugosidad de este y
es conocido como el coeficiente (C) de Chezy que tiene dimensiones L1/2 T −1 . La ecuación 123 se convierte
en la ecuación de Chezy:
1/2 1/2
Vo = CRh So1/2 Q = CARh So1/2 (124)
Diversos valore del C de Chezy se han obtenido a través de la experimentacion, esto valores van desde
60 pie1/2 s−1 para canales pequeños y rugosos hasta 160 pie1/2 s−1 para canales grandes y superficie suave.
De acuerdo con la ecuación 87, el esfuerzo cortante de pared promedio τ0 para un canal con flujo
uniforme, se expresa como:

τ0 = ρgRh So (125)
La ecuación 123 y quizas, desde el punto de vista teórico, la ecuacion mas apropiada para el analisis
de flujo uniforme, en donde el f se estima usando la ecuación de Colebrook-White. Sin embargo Robert
Manning en 1891 haciendo experimento en canales reales, encontró:
 1/2 1/6
8g Rh
C= ≈α (126)
f n
donde n es un parámetro de rugosidad y α es un factor de conversión: para Sistema Internacional α = 1
y para Sistema Ingles α = 1.486. Reemplazando la ecuación 126 en la ecuación 124, se tiene la ecuación
de Manning para flujo uniforme en canales:

1.0 2/3 1/2

Vo = R S Vo en m/s, Sistema Internacional
n h o (127)
1.486 2/3 1/2
Vo = Rh So Vo en pie/s, Sistema Inglés
n

Valores experimentales de n se encuentran en la figure 30 para varias superficies de canales. n puede

variar desde 0.01 para canales en vidrio con superficies lisas hasta valores de n = 0.15 para llanuras de
inundación en donde existen múltiples tipos de vegetación incluyendo arboles. Teniendo en cuenta que la
relación entre C y el n en la ecuación 126 no es del todo exacta, se ha encontrado que n cambia en función
del radio hidráulico; esto no solo incluye los efectos de cambios en la profundidad si no también de cambios
en la sección debido a la erosión y depositación de sedimentos.

32
Example 9. Unos ingenieros encontraron que el canal rectangular mas eficiente, el que transporta la
mayor cantidad de flujo uniforme para un area dada, es aquel que fluye con una profundidad igual a la
mitad del ancho del canal. Si consideramos un canal hecho en ladrillo (n = 0.015) sobre una pendiente
So = 0.006, determine el ancho óptimo si el canal transporta un caudal de 100 pie3 /s.

8.1 Estimación de la profundidad normal

En muchos problemas de ingenierı́a, dado el Q y las caracterı́sticas del canal y del fluido es necesario
estimar la profundidad normal yn . Dicho estimación se hace solucionando la ecuación 127 para yn . En la
tabla 3 se presenta la geometrı́a para las secciones transversales tı́picas de canales.

Tipo Perimetro(P) Area(A) Variable

Rectangular  b + 2yn   byn  yn
2
yn
Triangular yn sin1θ1 + sin1θ2 1
+ 1
yn
  h 2 tan θ1 tan θ2 i
yn
Trapezoidal yn sin1θ1 + sin1θ2 + b 1 1
2 yn tan θ1 + tan θ2 + 2b yn
yn = R Si θ = 90
sin 2θ
R2


Circular 2Rθ θ− 2 yn = R(1 − cos θ) Si 0 < θ < 90
yn = R(1 + sin θ) Si 90 < θ < 180

Table 3: Propiedades geométricas para diferentes secciones transversales de un canal. La notación de la

ecuaciones se puede observar en figura 31.

Reemplazando Rh = PA en la ecuación 127 para cualquiera de las secciones de la tabla 3, se obtiene

una ecuación implicita en yn , con exepción de la seccion circular en donde la función es implicita para
0o < θ < 180o . La ecuación a resolver es:

α 2/3
f (yn ) = Q − A Rh So1/2 Para todas las secciones exepto la circular
n (128)
α 2/3 1/2
f (θ) = Q − A Rh So Para la sección circular
n

Para resolver la ecuación 128, se presenta a continuación un procedimiento con base en el método de
bisección:

1. Leer la siguiente informacion:

(a) Para la todas las secciones exepto la circular leer: n, α, So , Q, b, θ1 y θ2 . Note que para la
sección rectangular θ1 = θ2 = 90o , y para la sección triangular b = 0.
(b) Para la sección circular leer: n, α, So , Q, R y θ.

2. Definir los lı́mites de iteración, a y b, para el método de bisección:

(a) Para todas las secciones exepto la circular leer: a = 0.01 y b = 1000 (a ≤ yn ≤ b). Note que
yn > 0 y que sign(f (a)) ̸= sign(f (b)). Note que la función sign extrae el signo de f .
(b) Para la sección circular leer: a = 0.01o y b = 179.9o (a ≤ θ ≤ b). Note que 0o < θ < 180o , y que
sign(f (a)) ̸= sign(f (b)).

3. Si N <= N max = 10000 ir a 4. Si N > N max ir a 10.

33
 a+b
4. Calcular C = 2 .
b−a
5. Calcular f (c) usando la ecuación 128. Si |f (c)| < Error o 2 < Error ir a 9. Si no se cumple
alguna de estas condiciones ir a 8.

6. Si sign(f (c)) = sign(f (a)), a = c, si no b = c. Incrementar N = N + 1 e ir 3.

7. Imprimir:

(a) para todas las secciones exepto la circular: yn = c.

(b) para la sección circular: θ = c. Usar las ecuaciones de la tabla 3 para convertir θ en yn .

8. Ha fallado la iteración. Revisar datos de entrada y lı́mites iniciales de iteración.

Example 10. Un canal trapezoidal en asfalto (n = 0.016) y con So = 0.0015, transporta agua con un
caudal de 300 f t3 s−1 . Si el ancho del canal es b = 6 f t y θ1 = θ2 = 50o , calcular la profundidad normal
(yn ).

Example 11. Determinar para que relaciones de yn /D es posible transportar un caudal máximo y obtener
una velocidad máxima en un canal de sección circular.

8.2 Sección óptima de un canal

El diseño de canales en ingenierı́a requiere determinar la sección óptima (aquella que garantice la menor
resistencia al flujo) en un canal para unas condiciones dadas. Usando la ecuación de Manning (ecuación 127)
es posible maximizar el radio hidráulico (Rh ) para un área (A) y un caudal (Q) dado. Teniendo en cuenta
que Rh = A/P , para maximizar Rh debemos minimizar el perı́metro (P ).
Considerando el canal trapezoidal de la figura 32, el area y el perimetro se pueden calcular como:

A = by + βy 2 (129)

1/2
P = b + 2W = b + 2y 1 + β 2 (130)
p
1/2
donde β = cot θ y y = yn . Note que W = y 2 β 2 + y 2 = y β 2 + 1 . Asumiendo que b es un dato dado,
se despeja de la ecuación 129:

A
b= − yβ (131)
y
y de la ecuación 130:

1/2
b = P − 2y 1 + β 2 (132)
igualando las ecuaciones 131 y 132 y despejando para P , se tiene:

A
1/2
P = − yβ + 2y 1 + β 2 (133)
y
Para minimizar P , se obtiene dP/dy de la ecuación 133 tomando como constantes A y β:

34
 dP A
1/2
= − 2 − β + 2 1 + β2 (134)
dy y
igualando dP/dy = 0 y despejando y se optiene el valor de y que minimiza P :
" #1/2
A
y= (135)
2 (1 + β 2 )1/2 − β
Despejando A de la ecuación 135 y reemplazando en la ecuación 133 para P , la geometrı́a optima de
un canal trapezoidal para un angulo θ dado es:
h
1/2 i
1/2 y
A = y2 2 1 + β 2 −β P = 4y 1 + β 2 − 2βy Rh = (136)
2

La ecuación 136 indica que para un angulo θ dado, la sección transversal mas eficiente de un canal en
flujo uniforme ocurre cuando el radio hidráulico es igual a la mitade la profundidad normal. Igualando las
ecuaciónes 130 y 136 para P , se tiene una expresión para la base inferior del canal b:
p 
b = 2y 1 + β2 − β (137)

y la base superior bs p 
bs = b + 2βy = 2y 1 + β2 (138)

Teniendo en cuenta que cuando θ = 90o β = 0 por lo que tenemos un canal rectangular, reemplazando
esto en las ecuaciones 136 se tiene la sección más eficiente para un canal rectangular:

y
A = 2y 2 P = 4y Rh = b = 2y (139)
2

Para encontrar la profundidad correcta, es necesario solucionar las ecuaciones anteriores en conjunto
con la ecuación 127 para un caudal dado.
Es posible encontrar el mejor angulo θ del canal trapezoidal para un área y una profundidad dada.
Evaluando dP/dβ para A y y constante se tiene:

dP
−1/2
= −y + 2βy 1 + β 2 (140)
dβ
igualando la ecuación 141 a cero se tiene:

1/2 1
2β = 1 + β 2 β = cot θ = √ θ = 60o (141)
3
Esto indica que para una profundidad dada, el angulo θ óptimo es 60o , o que la sección óptima (la que
maximiza el flujo) es la mitad de un hexágono.
Si se circunscribe una sección circular de radio R en la sección trapezoidal donde θ = 60o , se tiene
que R = 12 bs sin θ. Reemplazando bs en esta ecuación y sabiendo que β = √13 , se tiene que R = y. Esto
demuestra que una sección circular tiene su máxima eficiencia cuando y = R (semicirculo).
Para una sección triangular (b = 0 en la figura 32), el área y el perı́metro son:

A = βy 2 (142)

1/2
P = 2W = 2y 1 + β 2 (143)

35
  1/2
A
Despejando y de la ecuación 142 se tiene que y = β . Reemplazando en la ecuación 143, el
perı́metro es:
 1/2
A
1/2
P = 2W = 2 1 + β2 (144)
β
derivando P con respecto a β en la ecuación 144, se tiene:

2 1/2
"
#
dP −1/2 1 1 + β
= 2A1/2 β 1 + β 2 β −1/2 −


(145)
dβ 2 β 3/2
Haciendo dP/dβ = 0, se tiene:

1/2
β 1/2 1 + β2
= (146)
(1 + β 2 )1/2 2β 3/2
Solucionando la ecuación 146, se tiene que β = 1. Esto quiere decir que una sección triangular con
θ = 45o es la mas eficiente.

Example 12. Para un canal rectangular en ladrillo (n = 0.015) diseñado para transportar 5 m3 /s y
con una pendiente de So = 0.001, determinar las mejores dimensiones para y y b. Realizar los mismos
cálculos si la sección transversal del canal es la mitad de un hexágono y si es un semicirculo. Comparar
los resultados.

8.3 Cálculo de la profundidad normal en cauces naturales

En la naturaleza, los canales por los cuales transcurre flujo posseen secciones irregulares como en el caso
de rios y quebradas. Dicha sección cambia en el tiempo debido a procesos de erosion y sedimentación en
el fondo y en las orillas producto de cambios en la hidrodinámica del flujo. El cálculo de la profundidad
normal se hace con base en la solución de la ecuación de Manning implicita para yn siguiendo un método
numérico. Si tenemos una sección irregular de un canal como el mostrado en la figura 33, conformada por
una serie de m puntos p(x, y) con coordenadas espaciales x y y, y cuyos coeficientes de Manning estan
dados entre cada par de puntos pi y pi+1 (segmento) para un total de m − 1 coeficientes, el área mojada
de la sección para un valor de yn se puede calcular como:
m−1 m−1
X X xi+1 − xi
A= Ai = [(yn − yi ) + (yn − yi+1 )] (147)
2
i=1 i=1

, perimetro se puede calcular como:

m−1
X m−1
Xp
P = Pi = (xi − xi+1 )2 + (yi − yi+1 )2 (148)
i=1 i=1

y el coeficiente de rugosidad de Manning como:

Pm−1
n i Ai
n = Pi=1
m−1 (149)
i=1 Ai
Note que la ecuación 147 cálcula el área mediante la sumatoria de las areas individuales (trapecios)
Ai formadas entre cada par de puntos pi y pi+1 , el perimetro en la ecuación 148 es la sumatoria de las

36
distancias Euclidianas entre cada par de puntos pi y pi+1 y el coeficiente de rugosidad de Manning es el
promedio ponderado de los m − 1 coeficientes de rugosidad usando como pesos las areas Ai .
A continuación se muestra un procedimiento para el cálculo de yn en canales con sección irregular
basado en el método de bisección:

1. Leer la siguiente información: n, α, So , Q, xs, ys y ns. Note que xs y ys son dos vectores con las
coordenadas x y y respectivamente para un numero m de puntos que conforman la seción transversal.
ns es un vector de coeficientes de Manning para los m − 1 segmentos que conforman la seccion.

2. Definir los lı́mites de iteración, a y b, para el método de bisección: a = 0.01 y b = 1000 (a ≤ yn ≤ b).
Note que yn > 0 y que sign(f (a)) ̸= sign(f (b)). Note que la función sign extrae el signo de f .

3. Si N <= N max = 10000 ir a 4. Si N > N max ir a 10.

a+b
4. Calcular C = 2 .

5. Con base en las coordenadas datas (xs y ys), interpolar un nuevo conjunto de puntos para las
profundidades normales a, b y c.

6. Calcular f (a), f (b) y f (c) usando la ecuación ecuación 128 usando la nuevas secciones interpoladas
para a, b y c. Note que para evaluar la función es necesario, para cada nuevo conjunto de puntos,
calcular el área (ecuación 147), el perimetro (ecuación 148) y la rugosidad de Manning (ecuación 149.
b−a
7. Si |f (c)| < Error o 2 < Error ir a 9. Si no se cumple alguna de estas condiciones ir a 8.

8. Si sign(f (c)) = sign(f (a)), a = c, si no b = c. Incrementar N = N + 1 e ir 3.

9. Imprimir: yn = c.

10. Ha fallado la iteración. Revisar datos de entrada y lı́mites iniciales de iteración.

References
[1] Yunus Cengel and John Cimbala. Ebook: Fluid mechanics fundamentals and applications (si units).
McGraw Hill, 2013.

[2] CA Duarte Agudelo. Mecánica de fluidos e hidráulica. Universidad Nacional de Colombia. Facultad
de Ingenierı́a., 2011.

[3] I Shames. Fluid mechanics and application. editorial graw hill, 2010.

[4] RL Street, GZ Watters, and JK Vennard. Elementary fluid mechanics, (1995). John Wiley and Sons,
New York, ISBN 0, 471(01310):3.

[5] Victor L. Streeter and E. Benjamin. Wylie. Fluid mechanics /. Tata McGraw-Hill,, New Delhi,, c1983.

[6] Frank M White. Fluid mechanics. New York, 1990.

37
Figure 26: Valores de K para diferentes tipos de accesorios para flujo turbulento (tomado de [1]).

38
 Figure 27: Valores de K para diferentes tipos de accesorios para flujo turbulento (tomado de [1]).

Figure 28: Tipos de flujo en un flujo a superficie libre en un canal (tomado de [6]).

Figure 29: Geometrı́a y notación en un canal inclinado, a) vista lateral y b) sección transversal (tomado
de [6]).

39
Figure 30: Valores de n de Manning para diferentes tipos de materiales en canales (tomado de [6]).

Figure 31: Notación en varios tipos de canales.

40
Figure 32: Geometria de un canal trapezoidal (tomado de [6]).

Figure 33: Sección transversal de un cauce irregular (e.g. rio).

41