CURSO HIDROLOGÍA E HIDRÁULICA APLICADAS

2.2. ENERGIA ESPECÍFICA

2.2.1. DEFINICION DE ENERGIA ESPECÍFICA
Como se vio en la sección 2.1 cuando las líneas de flujo son rectilíneas, la
distribución de presiones en la sección es hidrostática, la carga piezométrica h
es constante en la sección, y se cumple que
p
h  z   zf  y cos2 
 ,
con z la cota desde un origen arbitrario, p la presión manométrica en el punto de
cota z, γ=gρ el peso unitario, zf la cota del fondo del canal en la sección, y el
tirante en la sección y θ el ángulo del fondo del canal con la horizontal.

La carga hidráulica H en una sección se define como

U2 Q2
H  zf  y cos   
2
 zf  y cos   
2
,
2g 2gA2
con α el coeficiente de forma de Coriolis del perfil de velocidades (en general se
asume α=1) y U=Q/A la velocidad media en la sección, con Q el caudal y A el
área de la sección perpendicular al flujo.
La energía específica E se define como la carga hidráulica en la sección referida
al fondo del canal
Q2
E  y cos2   
2gA2
Para un canal de pendiente pequeña (tanθ=S0<0.1) y tomando α=1, que son las
condiciones usuales en los canales que trabajaremos en el curso
U2 Q2
Ey y .
2g 2gA2
Nótese que para una geometría de sección conocida la energía específica es
solo función del tirante en la sección, ya que el área A es función de y.

Para canales rectangulares se define el caudal por unidad de ancho q=Q/b con
b el ancho del canal, por lo que la energía específica se puede expresar como
q2
Ey
2gy 2

2.2.2. CURVA DE ENERGIA

En esta sección se analizará la dependencia funcional de la energía específica
respecto al tirante. La energía específica para un canal rectangular, de pendiente
suave y con distribución aproximadamente uniforme de velocidades está dada
por E=y+q2/(2gy2), que se puede reescribir como (E-y)y2=q2/(2g). Esta ecuación

UdelaR – FI – IMFIA 2.2-0

CURSO HIDROLOGÍA E HIDRÁULICA APLICADAS

es de tercer grado y tiene una raíz negativa y 2 raíces reales positivas que se
denominan tirantes alternos.
Esta deducción puede realizarse para un canal de sección cualquiera y la
existencia de los tirantes alternos puede generalizarse a cualquier tipo de
sección.

Al graficar el tirante contra la energía específica para un caudal dado, resulta una
curva con dos asíntotas y un mínimo. Se distinguen dos ramas que se acercan
a las asíntotas: la rama superior, con asíntota E=y, y la rama inferior con asíntota
horizontal y=0. En la rama superior de la curva la componente de velocidad es
más pequeña, predominando la componente debida al tirante, mientras que en
la rama inferior predomina la componente de velocidad de la energía específica.

Figura 2.2.1: Curva energía específica – tirante para distintos caudales.

Para encontrar el tirante yc para el que se da el mínimo de la curva de energía

específica E, basta con derivar la expresión de la energía respecto al tirante e
igualar a cero. Para un canal de pendiente y sección cualquiera
dE Q2 dA Q2 Q2 B U2
 cos   
2
 cos   
2
B  cos   
2
 cos   
2
,
dy gA3 dy gA3 gA2 A gD
definiendo la profundidad hidráulica D=A/B con B el ancho superficial. En la
derivación anterior se asumió que α no depende de y. Resulta entonces que la
U 2
energía mínima se da cuando  1 . Por tanto el número de Froude para
gD cos2 
una sección de forma y pendiente cualquiera se define como
 Q 2B U 2
Fr 2   ,
gA3 cos2  gD cos2 

UdelaR – FI – IMFIA 2.2-1

CURSO HIDROLOGÍA E HIDRÁULICA APLICADAS

y la energía mínima se da cuando Fr=1, que corresponde a la condición de flujo

crítico. Por tanto para un caudal dado el flujo crítico es la condición de flujo que
lleva ese caudal con la menor energía específica posible.

Para un canal de pendiente pequeña (S0<0.1) y tomando α=1 se obtiene

Q 2B U 2
Fr 
2
 .
gA3 gD

Obsérvese que la rama superior de la curva (E,y) corresponde a flujos con Fr<1
(flujos subcríticos) y la rama inferior a números de Fr>1 (flujos supercríticos).

Para el caso particular de canales rectangulares el tirante crítico yc cumple que

1
 q2  3
yc   
g 
y la energía especifica crítica cumple que
3
Ec  yc
2 .
Para verificar que este extremo relativo de la curva es efectivamente un mínimo
se debe cumplir que la derivada segunda sea positiva. Para el caso particular de
un canal rectangular d2E/dy=3q2/(gy4)=3/yc>0.

La energía específica mínima se da entonces para el tirante crítico yc. Para

energías específicas mayores a esta energía mínima existen dos tirantes que
tienen la misma energía (tirantes alternos). Ambos tirantes son capaces de llevar
el caudal dado con la energía específica disponible.

En el caso rectangular los tirantes alternos verifican que

y A -1+ 1 + 8FrA2
=
yB 2 .

Para niveles de energía específica menores a la energía mínima no existen

tirantes capaces de llevar el caudal dado con tan poca energía.

2.2.3. CURVA CAUDAL ESPECÍFICO – TIRANTE

En este caso se estudia como varía el caudal específico en función del tirante
para un nivel de energía dado. Para el caso de un canal de sección rectangular
se puede escribir
q2  2gEy 2  2gy 3

UdelaR – FI – IMFIA 2.2-2

CURSO HIDROLOGÍA E HIDRÁULICA APLICADAS

Esta curva cumple q=0 para dos tirantes diferentes (y=0 e y=E). Luego para cada
caudal específico, menor a un valor máximo de q, se tienen dos valores de tirante
para cada nivel de energía.

Derivando q2 respecto a y e igualando a 0, se encuentran los extremos relativos

de q2 esta curva resulta que
dq2 dq
 2q  2gy  2E  3y   0
dy dy
se anula para dos raíces y=0 e y=2/3E.

De y=0 resulta q=0 que es un mínimo relativo. La otra condición y=2/3E coincide
con la condición de flujo crítico, e implica
8
q2  g E3
27
Obsérvese además que en este último caso resulta ser un máximo ya que
d 2q 2
 2gE
dy 2 .
Por tanto el caudal que conduce un canal en una sección para una energía
especifica dada es máximo para la condición de flujo crítico.

Figura 2.2.2: Curva caudal por unidad de ancho – tirante para distintas energías
específicas.

Se destaca entonces que el flujo crítico verifica simultáneamente ser:

 La condición de flujo para la cual circula un caudal dado con el mínimo

nivel de energía específica.
 La condición de flujo para la cual con un nivel de energía específica dado
circula el máximo caudal.

UdelaR – FI – IMFIA 2.2-3

CURSO HIDROLOGÍA E HIDRÁULICA APLICADAS

2.2.4. CONTROL DEL FLUJO

La celeridad de una onda se define como la velocidad de la onda relativa al medio
en la cual viaja. Para un canal prismático de sección rectangular de ancho b, la
celeridad c de una onda de pequeña amplitud se puede deducir a partir de aplicar
las ecuaciones de continuidad y cantidad de movimiento a una región de fluido
que incluye la perturbación.

Figura 2.2.3a: Esquema de volumen de control de fluido en un referencial

absoluto.

aquí la onda se propaga con velocidad ci sobre un flujo de velocidad media Ui,
por lo que en dicho referencial se observa la onda moviéndose con velocidad
(U+c) i.

Luego es conveniente expresar el problema en un referencial que se mueve con

la velocidad de la onda (U+c) i de forma de ver la onda fija,

-c dU-c

Figura 2.2.3b: Esquema de un volumen de control de fluido en un referencial

relativo, que se mueve acompañando la onda de pequeña amplitud.

y escribir las ecuaciones en este nuevo referencial inercial. Obtenemos para la

continuidad
cyb + (dU-c) (y+dy) b = 0
entonces
(dU-c) (y+dy) = -yc (1)

UdelaR – FI – IMFIA 2.2-4

CURSO HIDROLOGÍA E HIDRÁULICA APLICADAS

y despreciando los términos de segundo orden se encuentra

dUy = cdy (2)

Y para la cantidad de movimiento

-ρ c2y b + ρ (dU-c)2 (y+dy)b = ρg y2/2 b - ρg (y+dy)2/2 b

usando la ecuación (1) se obtiene

-cy2 - yc (dU-c) = g y2/2 - g (y+dy)2/2

luego

- yc dU = - g ydy -dy2/2

Usando la ecuación (2) y despreciando los términos de segundo orden

c2=gy

resulta que la celeridad de la onda es

c  gy

Para el caso general de un canal cualquiera es posible demostrar que de forma

más general la celeridad de una onda de pequeña amplitud vale

c  gD

Como puede observarse la celeridad de la onda de pequeña amplitud resulta ser

el denominador de la expresión del número de Froude

U U
Fr   .
gD c

En el caso de flujo supercrítico se tiene U>c, por lo que la velocidad absoluta de

la onda hacia aguas arriba w   U  c es negativa y nunca alcanzará a remontar
el flujo.

En el caso de flujo subcrítico se tiene U<c, por lo que la velocidad absoluta de la

onda hacia aguas arriba w   U  c es positiva y en este caso la onda remonta
el flujo.

En la situación particular de flujo crítico U=c, y la celeridad de la onda de pequeña

amplitud coincide exactamente con la velocidad del flujo.

UdelaR – FI – IMFIA 2.2-5

CURSO HIDROLOGÍA E HIDRÁULICA APLICADAS

Se dice por tanto que el flujo subcrítico está controlado desde aguas abajo, ya
que las perturbaciones que introduzca al flujo aguas abajo de una cierta sección
terminarán finalmente incidiendo en la condición del flujo en esa sección.

Por el contrario el flujo supercrítico se dice está controlado desde aguas arriba,
ya que el flujo en una determinada sección nunca se enterará de las
perturbaciones que sufra el flujo en una sección ubicada aguas abajo de la
primera.

2.2.5. TRANSICIONES
Una transición en canal puede definirse como un cambio en la sección del canal;
por ejemplo un cambio en el ancho del canal o en la elevación del fondo. Estos
cambios geométricos pueden desarrollarse en largas distancias o pueden ser
súbitos.

Una transición usualmente se diseña de forma que las pérdidas de carga sean
pequeñas, de forma que en los desarrollos que siguen se adoptará la hipótesis
de conservación de carga en la transición.

TRANSICIÓN POR PRESENCIA DE UN ESCALÓN

Si se desprecian las pérdidas de carga entre dos secciones (A) y (B)
HA  HB
Y se puede escribir
EA  HA  HB  EB  zf
Entonces si se conoce la carga en (A) y el cambio en la elevación del fondo
Δzf=zfB - zfA, se puede calcular la energía específica en (B) EB como
EB  EA  zf

UdelaR – FI – IMFIA 2.2-6

CURSO HIDROLOGÍA E HIDRÁULICA APLICADAS

Figura 2.2.4: Ejemplo de un perfil de flujo subcrítico sobre un escalón suave

(izquierda) y la correspondencia de los perfiles A y B en la curva de energía
específica – tirante (derecha). Asimismo, se indica gráficamente la máxima altura
del escalón (∆𝑧𝑚𝑎𝑥 ) que no altera la condición del flujo aguas arriba.
Para la energía específica EB existen dos posibles tirantes, uno subcrítico (B) y
otro supercrítico (B’). Si el flujo en (A) era subcrítico, como en la figura, el flujo
pasará de (A) a (B). Pasa a (B) y no a (B’), porque para pasar a (B’) debería o
bien incrementar el caudal por unidad de ancho y luego disminuirlo o bien perder
carga y luego recuperarla para llegar a (B’), claramente ninguno de esos dos
procesos son físicamente posibles.

Analicemos que sucede con los tirantes y la cota de la superficie libre aguas
abajo de un escalón si se desprecian las pérdidas de carga. Derivando la carga
respecto a la progresiva x a lo largo del canal
dH d  Q2  dzf dy Q2 dA
0   f
z  y    
dx dx  2gA2  dx dx gA3 dx
.
Si el canal es prismático, o sea que la forma del canal no depende de x, la
variación de A con x será solo asociada con la variación de y
dA dA dy dy
 B
dx dy dx dx ,
con B el ancho superficial. Notando además que
Q2B Q2 U2
Fr 2    ,
gA3 gA2D gD
se encuentra que
dy 1 dzf

dx 
Fr  1 dx
2

UdelaR – FI – IMFIA 2.2-7
CURSO HIDROLOGÍA E HIDRÁULICA APLICADAS

y
d d dy
zw   zf  y   Fr 2
dx dx dx .
Por lo que los cambios en el tirante y y en la cota de la superficie libre zw son en
el mismo sentido.

Debe notarse que Fr2 debe ser evaluado en x para poder proyectar la evolución
de y a lo largo del canal.

Entonces:
 Flujo subcrítico Fr<1 y el fondo se eleva dzf/dx>0 entonces el tirante y la
cota de la superficie libre descienden.
 Flujo subcrítico Fr<1 y el fondo desciende dzf/dx<0 entonces el tirante y
la cota de la superficie libre aumentan.
 Flujo supercítico Fr>1 y el fondo se eleva dzf/dx>0 entonces el tirante y la
cota de la superficie libre aumentan.
 Flujo supercítico Fr>1 y el fondo desciende dzf/dx<0 entonces el tirante y
la cota de la superficie libre descienden.

Nótese que existe un límite superior para la altura que puede tener el escalón
dados un caudal y tirante aguas arriba yA. Esta altura máxima de escalón está
dada en la figura por Δzf,max=EA-Ec. Si el incremento en el nivel de fondo supera
ese valor límite resulta que la condición de energía con la que viene el flujo no
es suficiente para sortear la transición. En dicho caso el flujo se altera en la
sección aguas arriba de la sección de transición (A), y en la sección de aguas
abajo (B) se genera una condición de flujo de mínima energía específica (flujo
crítico). Esta alteración puede darse mediante la disminución de caudal
específico o por el aumento del nivel aguas arriba del escalón.

Ejemplo 1:
Un canal rectangular de ancho b=3 m transporta un caudal Q=4.5 m3/s. En una
sección determinada del mismo existe un escalón de altura zf = 0.25 m, para
luego volver a la cota de fondo original. El tirante aguas arriba del escalón es de
1.66 m. Se pide:
 Determinar el tirante arriba del escalón y aguas abajo del mismo.
 Determinar la altura máxima que puede tener el escalón para no afectar
el flujo aguas arriba de la misma.

UdelaR – FI – IMFIA 2.2-8

CURSO HIDROLOGÍA E HIDRÁULICA APLICADAS

Figura 2.2.5: Esquema del perfil de flujo para el ejemplo 1.

Solución: Según las ecuaciones presentadas anteriormente, la energía aguas

arriba del escalón, en la sección 1, puede calcularse como
Q2 4.52
E1  y1   1.66   1.70m
2gA2 2  9.8   3  1.66 
2

Como no existe perdidas de carga en el escalón se debe verificar, para la sección

2 encima del escalón, que H1=H2, entonces
E2  E1  zf  1.45m ,
Q2 4.52 0.115
E2  y 2   y2   y2   E1  zf  1.45
2  9.81  3  y 2 
2 2
2gA2 y 22
,
y entonces

y23  1.45y22  0.115  0 .

Esta última ecuación tiene 3 soluciones

1.39m

y 2  0.32m
 0

La solución negativa no tiene sentido físico. De las dos soluciones positivas debe
analizarse cuál de las dos corresponde a la solución del problema. Para ello se
calculará el tirante crítico para saber qué tipo de flujo era el original
Q2 b Q2 Q2
Fr 2    1  y c  3  0.61m
gA3 gb2 y 3 gb2
Por lo tanto y1>yc, entonces el flujo original es subcrítico. Como el escalón no
provoca cambio de flujo (hecho que se comprobara a continuación), entonces
y2=1,39 m.

Luego el fondo desciende nuevamente a su nivel original (Sección 3) y como no

hay cambio de flujo el tirante vuelve a su valor en la Sección 1 y se obtiene que
y3= y1=1.66 m.
UdelaR – FI – IMFIA 2.2-9
CURSO HIDROLOGÍA E HIDRÁULICA APLICADAS

Para la resolución de la segunda parte del problema, se volverá a recurrir a el

tirante crítico calculado anteriormente yc=0.61 m.

Sustituyendo en la ecuación de energía, se obtiene

Q2
Ec  yc   0.92m
2gAc2

Utilizando nuevamente la conservación de la carga entre dos secciones, aguas

arriba y aguas abajo del escalón
E1  zf ,max  Ec  zf ,max  0.79m
Como Δzf,max>Δzf, entonces la elevación del fondo no introduce cambio de flujo
como se había asumido en la primera parte del problema.

UdelaR – FI – IMFIA 2.2-10

CURSO HIDROLOGÍA E HIDRÁULICA APLICADAS

TRANSICIÓN POR CONTRACCIÓN.

Figura 2.2.6: Esquema de planta (izquierda) y corte (derecha) de flujo subcrítico

a través de una contracción en el canal.

Recordando del punto anterior que la carga respecto a la progresiva a lo largo

del canal x se conserva
dH d  Q2  dzf dy Q2 dA
0   f
z  y    
dx dx  2gA2  dx dx gA3 dx
.
Estudiando ahora el caso de un canal de sección rectangular (A=by) cuya ancho
b varía con x. Asumiendo que la cota del fondo zf no varía
dy Q2 d dy Q2  y db dy  dy  y db dy 
0  3  by    3 b      Fr 2   
dx gA dx dx gA  b dx dx  dx  b dx dx  ,
entonces
dy  Fr 2   y  db
 
dx  1  Fr 2   b  dx
.
2
Debe notarse que Fr , b e y deben ser evaluados en x para poder proyectar la
evolución de y a lo largo del canal.

Por lo tanto, para un canal de sección rectangular:

 Flujo subcrítico Fr<1 y el canal se angosta db/dx<0 entonces el tirante y
la cota de la superficie libre descienden.
 Flujo subcrítico Fr<1 y el canal se ensancha db/dx>0 entonces el tirante y
la cota de la superficie libre aumentan.
 Flujo supercrítico Fr>1 y el canal se angosta db/dx<0 entonces el tirante
y la cota de la superficie libre aumentan.
 Flujo supercrítico Fr>1 y el canal se ensancha db/dx>0 entonces el tirante
y la cota de la superficie libre decienden.

UdelaR – FI – IMFIA 2.2-11

CURSO HIDROLOGÍA E HIDRÁULICA APLICADAS

Observar que en este caso también existe un límite para el cual se puede
contraer la sección. Si el ancho se reduce más allá de ese valor la condición de
energía con la que viene el flujo no es suficiente para superar la contracción y se
perturbará el flujo aguas arriba de la transición, de forma similar a la vista
anteriormente para el caso de transición en el nivel de fondo.

Figura 2.2.7: Ejemplo de un perfil de flujo subcrítico afectado por una contracción
del canal (izquierda) y la correspondencia de los perfiles 1 y 2 en la curva de
energía específica – tirante (derecha). Se observa que al aumentar el caudal por
unidad de ancho (q), se mueve hacia la derecha la curva de energía específica.
Asimismo, se indica gráficamente la máxima contracción (𝑏𝑚𝑖𝑛 ) que no altera la
condición del flujo aguas arriba.

Ejemplo 2:
Un canal rectangular de ancho b1 = 3 m transporta un caudal de 4.5 m3/s. En una
sección determinada del mismo se produce una contracción puntual del ancho
pasando a b2 = 1,5 m, para luego volver al ancho original. El tirante aguas arriba
de la contracción es de 1.66 m. Se pide:
 Determinar el tirante en la contracción.
 Determinar el ancho mínimo que puede tener la contracción para no
afectar el flujo aguas arriba de la misma.

UdelaR – FI – IMFIA 2.2-12

CURSO HIDROLOGÍA E HIDRÁULICA APLICADAS

Figura 2.2.8: Esquema de planta (izquierda) y corte (derecha) para el perfil de

flujo del ejemplo 2.

Solución: Según las ecuaciones presentadas anteriormente en este capítulo, la

energía aguas arriba de la contracción puede calcularse de la siguiente manera
Q2 4.52
E1  y1   1.66   1.70 m
2gA12 2  9.81  3  1.66 
2
.
Como no existe perdidas de carga en la contracción y la cota del fondo no varía
se debe verificar que E1=E2, entonces
Q2 4.52 0.46
E2  y 2   y 2   y 2   1.70  E1
2gA22 2  9.81 1.5  y 2 
2
y 22
entonces
y23  1.704y22  0.46  0 .
Esta última ecuación tiene 3 soluciones
1.50m

y 2  0.66m
 0

La solución negativa no tiene sentido físico. De las dos soluciones positivas debe
analizarse cuál de las dos corresponde a la solución del problema. Para ello se
calculará el tirante crítico para saber qué tipo de flujo era el original.
Q2 b Q2 Q2
Fr 2    1  y c  3  0.61m
gA3 gb2 y 3 gb2
Por lo tanto y1>yc, entonces el flujo original es subcrítico. Como la contracción no
provoca cambio de flujo (hecho que se comprobara en la parte siguiente del
problema), entonces y2=1.5 m.

Luego el ancho del canal vuelve a su valor original (Sección 3) b3 = b1 = 3 m y

como no hay cambio de régimen el tirante vuelve a su valor en la Sección 1 y se
obtiene que y3= y1=1.66 m.

Para la resolución de la segunda parte del problema se volverá a recurrir a la

condición de Fr=1 para calcular el tirante critico en la contracción.
Recordando que como se conserva la carga y no hay cambios en la cota del
fondo
E1  E2 .
Imponiendo que el flujo sea critico en la sección contraída y recordando que para
un canal rectangular Ec=3/2yc
3
1.70m  E1  E2  Ec  yc
2 .
Entonces
yc  1.14m

UdelaR – FI – IMFIA 2.2-13

CURSO HIDROLOGÍA E HIDRÁULICA APLICADAS

Sustituyendo en la ecuación de Fr=1, se puede calcular el ancho mínimo que

puede tener la contracción sin afectar el flujo aguas arriba
Q2 b Q2
Fr 2   1
gA3 gb2 y 3
Entonces
Q2
bmin   1.19m
gyc3
Como bmin<b2, entonces la contracción no introduce cambio de flujo como se
había asumido en la primer parte del problema.

2.2.6. SECCION DE CONTROL

Una sección de control se define como una sección en la que existe una
condición que determina una relación única entre el caudal circulante y el tirante
en la misma. Por lo tanto si se conoce el tirante en dicha sección es posible
calcular el gasto a través de ella y viceversa. Las secciones en que se da el flujo
crítico son secciones de control.

Existen varios dispositivos de medición del caudal circulante en un canal que

fuerzan la ocurrencia del mismo ya sea por una contracción, una elevación del
lecho o ambas. En estos dispositivos, a partir de la medida del tirante se puede
determinar el caudal que circula. Ejemplo de estos dispositivos son los
medidores Parshall y los vertederos de cresta ancha.

En el apartado anterior se mostró que la condición de flujo crítico (mínima

energía específica) se puede producir en un canal ya sea aumentando su nivel
de fondo, reduciendo su ancho, o mediante una combinación de ambas
acciones.

Para la situación en la que existe un escalón de fondo, idealmente la condición

de flujo crítico resulta de imponer que el tirante sea mínimo en la sección
(dy/dx=0, d2y/dx2<0). Notando que
dy dzf  1 
 0
dx dx  Fr 2  1  .
Entonces será necesario que dzf/dx=0 y la sección de control, para el caso ideal
en que no existen perdidas de carga ni curvatura de las líneas de flujo, será en
el punto más alto del escalón, donde dzf/dx=0.

Para la situación en la que existe una contracción en el ancho de la sección,

idealmente la condición de flujo crítico resulta de imponer (dy/dx=0, d2y/dx2<0).
Notando que

UdelaR – FI – IMFIA 2.2-14

CURSO HIDROLOGÍA E HIDRÁULICA APLICADAS

dy  Fr 2   y  db
  0
dx  1  Fr 2   b  dx
.
Entonces serán necesario que db/dx=0 y la sección de control, para el caso ideal
en que no existen perdidas de carga ni curvatura de las líneas de flujo, será en
el punto de ancho mínimo, donde db/dx=0.

Ejemplo 3:
En el canal del Ejemplo 1: canal rectangular de ancho b=3 m, que transporta un
caudal Q=4.5 m3/s, originalmente con un tirante de 1.66 m. Encuentra en una
sección determinada que el fondo asciende una altura zf = 1 m para luego volver
a la cota de fondo original. Se pide:
 Determinar el tirante encima del escalón, aguas arriba y aguas abajo del
mismo.

Figura 2.2.9: Perfil de flujo para el ejemplo 3.

Solución: De acuerdo con el Ejemplo 1 el escalón de zf = 1 m es mayor que el

escalón máximo determinado para un tirante de 1.66 m, Δzf,max= 0.79 m. Por
tanto el flujo aguas arriba del escalón será afectado de forma de ganar la
suficiente energía específica para poder sobrepasar el escalón. El flujo ganará
energía específica en la Sección 1 hasta lograr que el caudal pueda pasar por
encima del escalón con la energía mínima y por tanto sobre el escalón se
establecerá el flujo crítico. Entonces y2=yc=0.61 m, y E2=Ec=0.92 m.

Como se desprecian las perdidas de carga en el escalón se debe verificar que

H1=H2
E1  E2  zf  Ec  zf  1.92m ,
Q2 4.52 0.115
E1  y1   y1   y1   1.92
2  9.81  3  y1 
2 2
2gA1 y12
,
y entonces

UdelaR – FI – IMFIA 2.2-15

CURSO HIDROLOGÍA E HIDRÁULICA APLICADAS

y13  1.92y12  0.115  0 .

Esta última ecuación tiene 3 soluciones

1.89m

y1  0.26m
 0

Sabemos que el flujo original era subcritico en la Sección 1 por lo que se toma la
solución subcritica para el tirante en 1, y1=1.89 m.

Luego el fondo desciende nuevamente a su nivel original (Sección 3). En este

caso no se puede asegurar que no hay cambio de régimen ya que sobre el
escalón se pasa por flujo critico y aguas abajo el flujo puede tanto volver a la
condición de aguas arriba y se obtiene y3= y1=1.89 m o ir al tirante alterno, que
tiene la misma energía específica y3= y1alt=0.26 m. Cuál de estas dos soluciones
se de en la realidad dependerá de las condiciones del flujo aguas abajo de la
Sección 3, que en este problema no se especifican.

UdelaR – FI – IMFIA 2.2-16

CURSO HIDROLOGÍA E HIDRÁULICA APLICADAS

2.2.7. Apéndice

Introducción

Esta Sección se enfoca en la programación de herramientas de cálculo para el

tema abordado en el presente capítulo, energía específica; utilizando el lenguaje
de Octave (o Matlab).

En el capítulo anterior se describió la definición más general de una función en

los lenguajes de programación mencionados, definida como un archivo (.m) que
acepta argumentos de entrada y devuelve argumentos de salida. Sin embargo,
existen distintos tipos o clases de funciones. A continuación se presentan las
definiciones de función local y de funciones de funciones.

Función local: es una función definida dentro del archivo (.m) de otra función o
de un script. Solamente son visibles para el código de dicho archivo, y no pueden
ser llamadas desde la ventana de comandos ni desde otros archivos. En el
ejemplo abajo, la primera función del archivo (líneas de código 1 a 5) es
denominada función principal. Ésta es visible para otros archivos y también
puede ser ejecutada desde la ventana de comandos. Las funciones adicionales
en el archivo (líneas de código 6 a 9 y 10 a 18) son denominadas funciones
locales o sub-funciones, y pueden ser definidas en cualquier orden luego de la
función principal. Las funciones locales solamente son visibles para las otras
funciones en el mismo archivo. Aunque los comandos end no son necesarios
para las funciones, por buenas prácticas se pueden incluir (líneas 5, 9 y 18),
indicando el final de las respectivas funciones. En cambio el comando end de la
línea 17 es necesario para cerrar el bucle if.

1 function [avg, med] = mystats(x)

2 n = length(x);
3 avg = mymean(x,n);
4 med = mymedian(x,n);
5 % end

6 function a = mymean(v,n)
7 % MYMEAN: ejemplo de una función local.
8 a = sum(v)/n;
9 % end

10 function m = mymedian(v,n)
11 % MYMEDIAN: otro ejemplo de una función local.
12 w = sort(v);
13 if rem(n,2) == 1
14 m = w((n + 1)/2);
15 else

UdelaR – FI – IMFIA 2.2-17

CURSO HIDROLOGÍA E HIDRÁULICA APLICADAS

16 m = (w(n/2) + w(n/2 + 1))/2;

17 end
18 % end

Funciones de funciones: aunque suene redundante, se denomina así a una clase

de funciones que utilizan otra función como uno de sus argumentos de entrada.
Típicamente la optimización, búsqueda de raíces y resolución de ecuaciones
diferenciales ordinarias se encuentran dentro de esta clase de funciones. Para
diferenciar una función como argumento de entrada de una variable cualquiera,
se escribe el nombre de la función precedido por el símbolo @. Un ejemplo de
esta clase de funciones es la función predefinida fsolve, que resuelve un
sistema no lineal de ecuaciones de la forma: 𝐹(𝑥) = 0. Los argumentos de
entrada obligatorios son la función que se quiere minimizar (𝐹(𝑥)), y el valor
inicial de la variable (𝑥0 ) por el que comienza la búsqueda de la raíz del problema.
Más abajo se muestra un ejemplo de utilización de fsolve.

Es útil definir una función local cuando necesitamos utilizarla más de una vez
dentro del archivo que estamos escribiendo, pero a su vez es poco probable que
la utilicemos en otros archivos. A continuación se presentan dos ejemplos.

Cálculo de energía específica y tirante alterno

El siguiente código es una función que calcula la energía específica y el tirante

alterno en una sección de un canal rectangular de ancho y tirante conocidos,
para un caudal dado. Se puede observar que la función Eesp_rect es la función
principal (líneas de código 1 a 22), mientras que alt_rect (líneas 23 a 37 del
código) es una función local o sub-función auxiliar.

La función principal Eesp_rect tiene como argumentos de entrada y, b y Q,

definidos en las líneas 4 a 6 del código. Los argumentos de salida son E e yalt,
definidos en las líneas 8 y 9. Para calcular la energía específica E en la sección
del canal es necesario calcular algunos de sus parámetros geométricos, como
se muestra en la línea 12 del código. Por ser un canal rectangular, se utiliza la
función rect_geom descrita en la Sección 2.1.5 del capítulo anterior. A partir de
la velocidad media en la sección (línea 13), es posible calcular E (línea 14). Para
un canal de sección rectangular, el tirante alterno yalt se puede calcular
directamente a partir del número de Froude (líneas 15 a 18). Para el caso
rectangular las restantes líneas de código (19 a 37) no serían necesarias, pero
se incluyen como ayuda para generalizar el código para geometrías más
complejas, como la trapezoidal, donde el tirante alterno no puede ser despejado
directamente. Entonces, se utiliza la función fsolve para encontrar el tirante
alterno (línea 22 del código).

Como se mencionó arriba, fsolve requiere dos argumentos de entrada

obligatorios, una función para la cual se quiere hallar el cero, y el valor inicial de

UdelaR – FI – IMFIA 2.2-18

CURSO HIDROLOGÍA E HIDRÁULICA APLICADAS

la variable de interés por el que comienza la búsqueda. En el código abajo la

variable de interés es el tirante alterno. Como valor inicial para comenzar la
búsqueda se utiliza el tirante alterno calculado para el caso rectangular. La
función que se usa como argumento de entrada está definida como una función
local en el mismo archivo (líneas 23 a 37 del código). Esta función local
(alt_rect) calcula el error relativo entre la energía específica dada (E) y la
calculada con su variable de entrada y, teniendo en cuenta así la definición de
tirantes alternos, es decir, que poseen la misma energía específica. La energía
específica E, así como el caudal Q y el ancho de fondo del canal b son parámetros
dados a la función dentro de la variable de parámetros de entrada denominada
par (línea 20). El tirante y se define como una variable de entrada independiente
a par en la función alt_rect para que fsolve pueda iterar en y. Como se
observa en la línea de código 22, la sintaxis de fsolve es: yalt=fsolve(@(y)
alt_rect(y,par), yalt). El símbolo @ le indica a Octave (o Matlab) que dicha
variable de entrada es una función: entre paréntesis se indica la variable en la
que debe iterar, y a continuación el nombre de la función para la que se quiere
hallar el cero (si la función tuviese una única variable de entrada, no sería
necesario indicarle en qué variable iterar, solamente se pondría el símbolo @
precedido por el nombre de la función).

En resumen, fsolve va a iterar en el tirante y hasta que la función alt_rect

valga cero (numéricamente un valor muy bajo). El tirante de la iteración final es
el que devuelve como resultado en su variable de salida yalt. Es muy importante
el valor inicial con el que comienza la iteración, ya que si es muy similar al tirante
dado como variable de entrada a la función principal, fsolve va a converger a
éste y no a su alterno, ya que ambos poseen la misma energía específica. Para
evitar este inconveniente, en el código abajo se utiliza como valor inicial para la
iteración la solución del tirante alterno para el caso de canal rectangular.

Notar que para generalizar el código para un canal trapezoidal se debe agregar
como variable de entrada el factor horizontal m de los taludes de las paredes
laterales (1V:mH), y modificar algunas líneas del código para incorporar este
parámetro en los cálculos (por ejemplo la línea 12).

1 function [E,yalt]=Eesp_rect(y,b,Q)
2 % Función que calcula la energía específica y el tirante
%alterno para un canal rectangular para un tirante, ancho
%de canal y caudal dados
3 % INPUTS
4 % y tirante (m)
5 % b ancho (m)
6 % Q caudal (m^3/s)
7 % OUTPUTS
8 % E Energía específica (m)
9 % yalt tirante alterno de y (m)
10 %

UdelaR – FI – IMFIA 2.2-19

CURSO HIDROLOGÍA E HIDRÁULICA APLICADAS

11 %% Cálculo de la energía específica:

12 [B,A,P,R,yG,D]=rect_geom(y,b); % función que calcula
%parámetros geométricos de la sección (algunos de los
%cuales no son usados en esta función), debe estar en el
%mismo directorio.
13 U=Q./A; % velocidad media
14 E=y+U.^2/(2*9.8); % energía especifica
15 Fr2=(Q^2)*B./(9.8.*A.^3); % número de Froude al cuadrado
16 %
17 %% Búsqueda del tirante alterno:
18 yalt=y./(-1 + sqrt(1 + 8./Fr2))*2; % Solución para caso
%rectangular que se puede usar para el caso trapezoidal %como
primera estimación del tirante alterno para iniciar
%iteración en fsolve

19 % NOTA: para el caso rectangular las siguientes líneas de

%código no son necesarias. Se incluyen como ayuda para la
%tarea.
20 par=[Q,b,E]; % para el caso trapezoidal: armado de vector
%de parámetros que permite pasar los valores de Q, b y E
%dentro de la función que busca el tirante alterno

21 % Búsqueda del cero de la función:

22 yalt=fsolve(@(y) alt_rect(y,par), yalt); % para el caso
%trapezoidal: función de Octave que busca el cero de la
%función alt_rect. Ver help de Octave para fsolve.

23 function eE=alt_rect(y,par)
24 % Función auxiliar que calcula el error relativo entre la
%energía específica dada y la energía específica calculada
%usando el tirante alterno estimado, que se buscará
%minimizar.
25 % INPUTS
26 % par: vector de parámetros de entrada
27 % y: variable de entrada a la función cuyo valor se
%seleccionará para minimizar el error eE
28 % OUTPUTS
29 % eE error relativo
30 %
31 % descomposición del vector de parámetros en las variables
originales
32 Q = par(1);
33 b = par(2);
34 E = par(3);
35 [B,A,P,R,yG,D]=rect_geom(y,b); % función que calcula
%parámetros geométricos de la sección (algunos de los

UdelaR – FI – IMFIA 2.2-20

CURSO HIDROLOGÍA E HIDRÁULICA APLICADAS

%cuales no son usados en esta sub-función), debe estar en

%el mismo directorio.
36 U=Q/A; % velocidad media
37 eE=E/(y+U^2/(2*9.8))-1; % error relativo entre el valor de
%la energía específica dada y la calculada con el valor de
%la variable de entrada y.

Cálculo del tirante crítico

El siguiente código calcula el tirante crítico para un canal rectangular de ancho
de fondo conocido, y un caudal dado. Utiliza la función fsolve, de manera muy
similar al ejemplo anterior. En este caso se utiliza la definición de tirante crítico,
es decir, aquel tirante para el cual el número de Froude vale 1. Las únicas
variables de entrada necesarias son el caudal Q y el ancho de fondo b (líneas 3
y 4 del script).

Para un canal de sección rectangular, el tirante crítico yc se puede calcular

directamente (líneas 6 del script). Por lo tanto, para el caso rectangular las líneas
de código 8 y 9 no serían necesarias, pero se incluyen como ayuda para
generalizar el código para geometrías más complejas, como la trapezoidal,
donde el tirante crítico no puede ser despejado directamente. Entonces, se utiliza
la función fsolve para encontrarlo (línea 9 del código). Los dos argumentos de
entrada obligatorios de fsolve son la función para la cual se quiere hallar el cero
(froude_rect en este caso, cuyo código también se incluye más abajo), y el
valor inicial de la variable de interés por el cual comienza a iterar. Como valor
inicial para comenzar la búsqueda se utiliza el tirante crítico calculado para el
caso rectangular. La función auxiliar froude_rect calcula el error relativo entre
el número de Froude calculado con su variable de entrada y y el valor 1 (por la
definición de tirante crítico). El caudal Q y el ancho de fondo del canal b son
parámetros dados a la función froude_rect dentro de la variable de parámetros
de entrada denominada par (línea 8 del script y líneas 10 y 11 de la función). En
las últimas líneas del script (10 a 12) se calcula la velocidad y energía específica
correspondientes al tirante crítico.

Al igual que en el ejemplo anterior, notar que para generalizar el código para un
canal trapezoidal se debe agregar como variable de entrada el factor horizontal
m de los taludes de las paredes laterales (1V:mH), y modificar algunas líneas del
código para incorporar este parámetro en los cálculos.

Script para calcular el tirante crítico:

1 %% Datos de entrada
2 clear all% borra todas las variables

3 Q = 0.03;% caudal (m3/s)

4 b = 0.2;% ancho de fondo (m)

UdelaR – FI – IMFIA 2.2-21

CURSO HIDROLOGÍA E HIDRÁULICA APLICADAS

5 %% Cálculo tirante crítico,

6 yc=(Q^2/(9.8*b^2))^(1/3);% estimación inicial del tirante
%crítico suponiendo canal rectangular (m)
7 % En el caso rectangular esta estimación inical es la
%solución. Las siguientes líneas de código se incluye aquí
%como ayuda para el caso trapezoidal.
8 par=[Q b];% vector con parámetros
9 yc=fsolve(@(y) froude_rect(y,par), yc);% tirante critico
%(m)

10 [Bc,Ac,Pc,Rc,yGc,Dc]=rect_geom(yc,b);% cálculo de
%parámetros geométricos asociados a las condiciones de
%flujo crítico
11 Uc=Q/Ac;% velocidad en flujo crítico (m/s)
12 Ec=yc+Uc^2/(2*9.8);% energía específica en flujo crítico
%(m)

Función froude_rect:

1 function ec = froude_rect(y,par)
2 % Función auxiliar que calcula la diferencia entre 1 y el
número %de Froude
3 % calculado con el tirante.
4 % INPUTS
5 % par vector de parámetros de entrada
6 % y tirante de entrada con el que calculará el Froude
7 % OUTPUTS
8 % ec apartamiento del Froude respecto a 1

9 %% Descomposición del vector de parámetros de entrada

10 Q = par(1);% caudal (m3/s)
11 b = par(2);% ancho de fondo (m)

12 %% Función
13 [B,A,P,R,yG,D]=rect_geom(y,b);% función que calcula
%parámetros %geométricos de la sección
14 ec=(Q^2)*B./(9.8*A.^3)-1;% diferencia entre el Froude
%estimado %con y que tiende al yc con 1

UdelaR – FI – IMFIA 2.2-22