Ingeniería hidráulica en México, vol. XXI, núm. 1, pp.

29-42, enero-marzo de 2006

Perfil del flujo sobre un vertedor

tipo cimacio y su comprobación
Gilberto Sotelo-Ávila

Universidad Nacional Autónoma de México

El ingeniero que se enfrenta al diseño de un vertedor tipo cimacio utiliza la información experimental
de diversos autores y dependencias para determinar su capacidad de descarga, el perfil del fondo,
el perfil del flujo del agua sobre el cuerpo del cimacio y la distribución de la presión sobre la espal-
da del mismo. La información utilizada es sencilla, confiable, fácil de obtener y muy conocida por los
especialistas en este tipo de obras.
Por su parte, varios autores han publicado diversas aproximaciones para predecir el perfil de la
superficie del agua y la distribución de la presión sobre el vertedor, utilizando para ello redes de flujo,
hipótesis simplificadoras en la teoría del potencial, la teoría del flujo curvilíneo poco profundo de
Dressler y la teoría del vórtice libre. El autor de este trabajo muestra que la aplicación de la teoría del
flujo sobre fondos curvos permite una predicción segura, tanto de la posición de la sección crítica
en la proximidad de la cresta del cimacio, como del perfil del flujo y la distribución de la presión en
el fondo, cuando se conoce su perfil geométrico y las condiciones en que opera.

Palabras clave: cimacios vertedores, vertido sobre cimacios, flujo curvilíneo sobre cimacios.

curva teórica para el coeficiente de descarga. Comparó

sus resultados con los experimentales obtenidos por
Antecedentes

Fawer (1937) inicia la investigación del flujo sobre un varios autores en vertedores con distintos radios de
vertedor de cresta redondeada mediante la elaboración curvatura en la cresta y encontró algunas coincidencias.
de una red de flujo según los lineamientos muy conoci- Knapp (1960) aceptó que las líneas de corriente que
dos de Prasil (1926), comprobando la teoría con resulta- cruzan una sección ortogonal a ellas en un flujo curvilí-
dos experimentales. Concluyó que su solución gráfica neo poseen el mismo centro de curvatura del fondo y
proporciona el caudal máximo para una carga sobre la que la velocidad en cada línea sigue la ley del vórtice
cresta dada de acuerdo con el principio de Belanger libre. Con ello estableció la influencia de la curvatura del
(1849), expresamente limitado a vertedores de cresta fondo en curvas verticales cóncavas o convexas, deter-
ancha donde las líneas de corriente son rectas. minó la distribución de la presión en la sección y el valor
Jaeger (1956) presentó las que él llamó ecuaciones de la energía total del movimiento en la misma. Knapp
básicas del flujo curvilíneo gradualmente variado en ca- aplicó su teoría al caso de un vertedor tipo cimacio y
nales, en las que incluyó un término que toma en cuen- encontró algunas coincidencias, con el apoyo de los
ta, de modo aproximado, la influencia de la curvatura del resultados experimentales del US Bureau of Reclamation
fondo en un flujo bidimensional a superficie libre. Jaeger (1963).
admitió que era suficiente considerar un valor medio del Sivakumaran et al. (1981) presentaron un compendio
radio de curvatura para fines prácticos. Con base en las de las ecuaciones no lineales de Dressler (1978) para
ecuaciones antes mencionadas, Jaeger analizó el flujo flujos de poca profundidad con curvatura y las aplicaron
sobre un vertedor de cresta redondeada y determinó la sobre la región de la cresta de un cimacio tipo WES para
distribución de la velocidad en la sección vertical sobre obtener la relación carga-gasto y la distribución de la
la cresta y la expresión para el caudal vertido. Con ésta presión en dicha región, sin llegar a obtener el perfil del
y el principio de energía específica mínima obtuvo una flujo. Encontraron buena coincidencia con los experi-
expresión para el tirante crítico en dicha sección y una mentos realizados con el mismo tipo de cimacios en la

29
 Sotelo-Ávila, G., Perfil del flujo sobre un vertedor tipo cimacio y su comprobación

Estación Experimental de Vías Fluviales WES (1978) del

Cuerpo de Ingenieros de los Estados Unidos de
Ilustración 1. Condiciones de operación de un cimacio en
términos adimensionales y pérdida de energía por entrada
Norteamérica.
al mismo.

ζ
Objetivo Hd Nivel del energía
he sobre el cimacio
Hd
El vertido del agua sobre un cimacio equivale al que
ocurre en un canal rectangular con fondo de curvatura Sección
crítica
convexa variable conocida, donde es importante deter- H0 Ec H
Hd dc Hd
minar el coeficiente de descarga, la posición de la sec- Cresta Hd
Hd
x
ción crítica en la proximidad de la cresta, el perfil de la Hd
0 A
superficie libre del agua y la distribución de la presión en
ζ
Hd d
el fondo, para cualquier condición de operación, es de- B Hd
θ q
cir, para cargas o gastos de vertido iguales o distintos a
ζ = - 1 x 1.85
los que se utilizan para el diseño geométrico del perfil. Hd 2 ( Hd )
Para su análisis se emplea la información experimental
publicada por distintos autores o dependencias.
El propósito de este trabajo consiste en aplicar el mo-
delo de flujo curvilíneo desarrollado por Sotelo y Esca-
lante (2001 a y b) a las diversas condiciones de opera- el cual debe ser constante para cada caudal unitario
ción del vertedor tipo cimacio, conocido como tipo WES, vertido, q = Q/b, donde Q es el caudal total vertido y b
para tener resultados que puedan compararse con los la longitud de cresta del cimacio. El valor de q (por ende
experimentales publicados por la estación experimental de F ) varía con la carga H0 de vertido o de operación del
de vías fluviales, Water Experiment Station (1978) para cimacio, medida también desde la cresta hasta la super-
dicho vertedor y comprobar o no la bondad del análisis. ficie libre aguas arriba. El coeficiente de descarga (o de
vertido), C = q/H0 , en m /2 /s, se relaciona entonces con
3/2 1

Diseño geométrico e hidráulico F mediante la expresión:

La estación experimental de vías fluviales, Waterways  H 3 2

C
Experiment Station (WES), U. S. Army Corps of Engineers F=  o (3)
(1978), realizó pruebas experimentales exhaustivas en 2 g  Hd 
modelo hidráulico, con el fin de establecer criterios para
determinar el perfil del flujo vertido sobre cimacios y la donde g es la aceleración de la gravedad.
distribución de la presión sobre ellos para distintas La ilustración 2 muestra los valores experimentales
condiciones de operación. La ecuación adimensional de C para el perfil tipo WES, y la ilustración 3, los valo-
que define el perfil geométrico del cimacio, con para- res de F calculados de la ecuación (3). Por ejemplo, en
mento vertical aguas arriba y velocidad de aproximación el cuadro 1 se presentan los valores calculados de F y
insignificante, es: los experimentales de C para tres cargas distintas de
operación.
1.85
1 x  La curva de los valores experimentales de C se ajus-
(1)
ζ
Hd 2  Hd 
=−   ta bien a la ecuación:

con el origen del sistema coordenado ( x, ζ ) en la cresta, 0.1316

como se muestra en la ilustración 1, siendo Hd la carga C = 0.49867 2 g H0 Hd ( ) (4)
adoptada para el diseño geométrico del cimacio, medi-
da como el desnivel entre dicha cresta y la superficie li-
bre aguas arriba. Cuadro 1. Valores experimentales de C y los calculados de F
Un parámetro adimensional relacionado con el flujo y
para tres cargas de operación.

el perfil geométrico del cimacio es:

H0/Hd 0.5 1.0 1.33
q
F= (2)
C (m1/2/s) 2.004 2.222 2.28
2 g Hd3 2 F 0.16 0.5016 0.7895

30 ingeniería hidráulica en méxico /enero-marzo de 2006

 Sotelo-Ávila, G., Perfil del flujo sobre un vertedor tipo cimacio y su comprobación

Ilustración 2. Coeficiente de descarga C en un cimacio tipo Ilustración 3. Valores calculados de F a partir de los
WES, con velocidad de llegada despreciable y distintas cargas experimentales de C. Valores calculados del coeficiente
de operación. de pérdida ε según la ecuación (23).

1.4
0 0.01 0.02 0.03 0.04
3
H0 1.4
Hd H0
3
1.2 Hd

ε = 0.03053
ε = 0.01314
1.2

1.0 1.0
F

0.8
0.8
Experimental
0.6

0.6
0.4

0.4

ε = 0.00974
0.2

0
0.2 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8
F

0
1.7 1.8 1.9 2.0 2.1 2.2 2.3 2 2
 q (8)

C
H = ζ + d cos θ + 
κ
 1− κ d ln 1− κ d
( ) ( )  2g
 

con errores menores del 0.6% en el intervalo: 0.3 ≤ expresión que permite determinar el perfil de la superfi-
H0 / Hd ≤ 1.33. Por tanto, de la ecuación (3) resulta: cie libre de la lámina vertiente sobre el cimacio, cuya
geometría se obtiene de la carga Hd de diseño.
1.6316 La carga Hd sirve para adimensionalizar la ecuación
F = 0.49867 H0 Hd (5)
( ) (8), donde se sustituye también la ecuación (2) y se ex-
presa en cualquiera de las formas siguientes:
2
Ecuaciones por aplicar
H cos θ F κ Hd

(9a)
ζ
La energía del flujo (por unidad de peso y tiempo) en una
 
Hd Hd κ Hd  1− χ ln 1− χ 
− = χ+
sección ortogonal al fondo curvo de un canal rectangu-   ( ) ( )
lar se obtiene de la ecuación (64) del artículo de Sotelo y
2
Escalante (2001a), que se expresa en la forma: κ Hd  H κ Hd F κ Hd

ζ    (9b)
cos θ  Hd Hd  cos θ  1 − χ ln 1 − χ 
 − = χ+
1 Q2 ( ) ( )
H = ζ + d cos θ + (6)
 

IQ2 1− κ d
2 2g donde χ = κd. La suma de términos antes del signo igual
( ) representa la energía específica (adimensional) en cada
donde IQ es una integral que para la sección rectangular sección ortogonal al fondo, siendo H la altura de la línea
quedó valuada mediante la ecuación (26) del artículo de de energía respecto de la cresta, la cual es constante
Sotelo y Escalante (2001b), y resultó: para cualquier punto del campo de flujo.
Como ocurre en los canales de fondo plano, la ecua-
IQ = − b κ ln 1− κ d
( ) ( ) (7) ción (9) tiene dos raíces χ (alternas) para cada valor de
H/Hd y de x/Hd, una en régimen subcrítico y otra en su-
donde κ es la curvatura del fondo. Por tanto, la energía percrítico, la primera mayor y la segunda menor que el
total del flujo con respecto del plano horizontal que pasa crítico en cada sección. Para la solución puede adop-
por la cresta es: tarse el método de Newton-Raphson.

ingeniería hidráulica en méxico /enero-marzo de 2006 31

 Sotelo-Ávila, G., Perfil del flujo sobre un vertedor tipo cimacio y su comprobación

La ecuación (66) del artículo de Sotelo y Escalante La sección crítica se localiza donde F = Fc , es decir,
(2001a) proporciona la carga de presión en el fondo cur- el número de Froude se iguala con el valor crítico expre-
vo del canal, dada por: sado por la ecuación (14). Por tanto, se debe cumplir
que:
p0 −2  u2
= d cos θ +  1− κ d − 1 0 (10) 2 κ H  H

gρ
( ) 2g 2 1− χ c 
( ) d
 − 1
ζ  
 χ c cosθ  Hd Hd  
  −

(17a)
donde u0 es la velocidad local en el fondo de la sección 3
y está relacionada con IQ, de la ecuación (7), mediante el (1− χ c ) ln (1− χ c )
cociente u0 = Q / IQ. Al sustituir u0 e IQ y adimensionalizar χ c [ ln (1− χ c ) + 1]
=−

el resultado, se obtiene:
o bien:
2
p0 cos θ F κ Hd (11)

−2
 1− χ − 1 κ Hd  H 1 1− χ c ln 1− χ c
 
g ρ Hd
=
κ Hd
χ+
 ln 1 − χ
 ( ) ζ  ( ) ( ) (17b)
 ( ) 

 
cosθ  Hd
 −
Hd 
 = χc −
2 ln 1− χ c + 1 ( )
Al eliminar el término (cos θ/κ Hd ) χ mediante la Esta ecuación permite calcular χc para cada valor del
ecuación (9a) , resulta que: término a la izquierda del signo igual, el cual depende a
su vez de H/Hd y de ζ /Hd .
2 La ecuación (17b) no ubica a la sección crítica y para
p0 H ζ   F κ Hd

 (12) localizarla se debe satisfacer una expresión cuyo desa-
g ρ Hd  Hd Hd   ln 1 − χ
= − −
( ) rrollo se expone en el Anexo al final del trabajo y que es:


 x H  dζ
donde p0 es la presión en cualquier punto del perfil ,  = 1− χ c ( )
sólido del cimacio (de ordenada x y abscisa ζ ), ρ es la  Hd Hd  dx
ƒ
(18)
densidad del líquido y F está definido por la ecuación (2).
Por otra parte, las ecuaciones (5c) y (23) del artículo   H cosθ  Hd2 (dκ / dx) = 0
+2  − 3
ζ 
de Sotelo (2004) definen el número de Froude y su valor H
−
H κ Hd 
χc 
κ Hd
crítico en un canal rectangular, a saber:
  d d

La curvatura κ en cada punto y el término d κ/dx se

u0 obtienen de la ecuación (1), de modo que al adimen-
(13)
sionalizar resultan las expresiones siguientes:
F =
g dcosθ

dζ  x 0.85
3 1/ 2 tanθ = = − 0.925  (19)
 1− κ d ln 1− κ d

( ) ( ) dx  Hd 

(14)

Fc = −
 κ d ln 1− κ d + 1

[ ( ) ] 
0.78625
κ Hd = −
Se substituye u0 de la ecuación (13) en la (6) y se obtiene: 0.15  1.7 
3 /2 (20)
( x / Hd )  1+ 0.855625 x / Hd( ) 
 
 
F2 
H = ζ + 1+ dcosθ (15)

1.7
Hd2 dκ / dx 0.15  1+ 15.40125 x / Hd
 
2
 2 1 − χ ( ) ( )  (21)
κ Hd x / Hd 1+ 0.855625 x / H
=−
1.7 
 d ( ) 
de la cual se despeja F y en términos adimensionales
resulta: La curvatura adopta dos valores sobre la cresta del
cimacio (x = ζ = 0): κ Hd = –2, como se observa de su
 κH  H 1/ 2 geometría y la que se obtiene de hacer x/Hd = 0 en la
F = 2 1− χ   − 1 (16) ecuación (19). Para x/Hd > 0, la solución de la ecuación
d ζ  
( )
 χ cos θ  Hd Hd  
 −
(9b) cambia para cada valor de H/Hd y de x/Hd ; además

32 ingeniería hidráulica en méxico /enero-marzo de 2006

 Sotelo-Ávila, G., Perfil del flujo sobre un vertedor tipo cimacio y su comprobación

de que hay que localizar la sección crítica. Se puede donde el primer término es constante y permite obtener
demostrar que cuando el perfil del fondo es convexo y el valor de χc que satisface a la ecuación (18), con el
asimétrico, como es el caso, la sección crítica se despla- fin de localizar la sección crítica para el valor de F con-
za hacia aguas abajo de la cresta y para localizarla se siderado.
utiliza la ecuación (18). Para el segundo modo, H/Hd es constante y es prefe-
Para comprobar la ecuación (9b) por la vía experi- rible resolver χc de la ecuación (17b) que satisfaga la
mental con un cimacio de geometría diseñada con la ecuación (18), para después calcular F de la ecuación
carga Hd , se pueden seguir dos caminos: (3), que corresponda a la energía H/Hd dada.
1. Se establece un caudal en la instalación con una Un resumen de los cálculos se presenta en los
carga adimensional de operación H/Hd y después se cuadros 2a, 2b y 2c destacando sólo los resultados fina-
mide el perfil de la superficie libre con diferentes absci- les para cada modo. Por ejemplo, para el primer modo,
sas x, se calcula F de la ecuación (9b), el cual se com- con H0 / Hd = 0.5 y F = 0.16, la sección crítica se
prueba con q y la definición dada por la ecuación (2). localiza en x/ Hd = 0.059302; para ella, χc = –0.41249,
2. También se puede medir q para conocer F, así H/ Hd = 0.495128, al satisfacer la ecuación (18).
como el perfil de la superficie libre para obtener χ; El valor F = 0.16 se obtuvo experimentalmente para
después con cos θ, κ Hd y la ecuación (9b), se calcula la H0 / Hd = 0.5, como puede verse en el cuadro 1 o en la
carga relativa de operación H/Hd , la cual se comprueba ilustración 3. La diferencia que se observa respecto del
con las mediciones. teórico (H/ Hd = 0.495128) antes mencionado, puede in-
Sivakumaran et al. (1981) analizaron parcialmente el terpretarse de tres maneras distintas: a) que exista un
caso utilizando el primer camino y encontraron algunas error como resultado de las imperfecciones en la teoría y
concordancias. El autor de este trabajo siguió el segun- el experimento (para los valores anteriores, del orden del
do camino. uno por ciento); b) que se produzca una pérdida de
energía por efecto de la aceleración del agua, la turbu-
lencia y el desarrollo de la capa límite, y que no fue
considerada en la teoría; y c) que ocurra un efecto
Procedimiento de cálculo

El primer paso consistió en encontrar la ubicación y ca- combinado de las dos anteriores.
racterísticas de la sección crítica que se forma en la pro- La segunda posibilidad es la única manera factible
ximidad de la cresta, para determinar la relación que hay de valuar numéricamente englobando a las otras dos; es
entre el caudal y la energía. Para ello se consideraron decir, se acepta que ocurra la pérdida de energía he / Hd
dos modos equivalentes de localizar la sección crítica en que se acota en la ilustración 1. Dicha pérdida se puede
el cimacio, a saber: a) obtener la energía específica mí- valuar como he = ε H0 , donde ε ≤ 1, de manera que
nima (crítica) necesaria para que fluya un caudal dado; H0 = H + he. Así, se tiene que: H0 = H + ε H0 ; ε = 1
b) determinar el caudal máximo que puede fluir con una – (H/H0); o bien, en forma adimensional se tiene:
energía específica dada.
Ambos modos se utilizaron para definir la relación
entre el gasto q y la carga H0 sobre la cresta, como se ε = 1− ( H / Hd ) /( H0 / Hd ) (23)
muestra en la ilustración 1, aceptando que H0 es diferen-
te de la energía que posee el flujo en la región del cima- donde ε es el coeficiente de pérdida en términos de la
cio y sólo para las tres condiciones de operación antes carga aguas arriba.
mencionadas. Para el primer modo se eligieron los cau- Para el segundo modo de abordar el problema se
dales (representados por F ) y se calcularon las alturas presentan los valores calculados de F para las cargas
de energía H/Hd correspondientes a la sección crítica. H/ Hd utilizadas, en los cuadros 2a, 2b y 2c; dichos
Para el segundo modo se eligieron las cargas adimen- valores son mayores que los que se muestran en el
sionales y se calcularon los caudales críticos (máximos) cuadro 1. Por ejemplo, cuando el nivel de la energía para
correspondientes, representados por F. el flujo sobre el cimacio es H/ Hd = 0.5, se obtiene F =
En ambos casos son válidas las ecuaciones (6) a 0.162573, que es mayor que el valor dado en el cuadro
(21), pero para el primer modo es conveniente combinar 1. Esto significa que dicho gasto es posible sólo si la
las ecuaciones (9b) y (17b). Así, resulta que: carga aguas arriba es mayor; es decir, si H0 / Hd > H/ Hd
3 = 0.5, donde H/ Hd = 0.504995, debiendo interpretarse
κ Hd 2 1 1− χ c ln 1− χ c
[( ) ( )] la diferencia entre ambas como la pérdida antes expues-
F κ Hd
( ) (22) ta, para después calcular ε de la misma manera que en
cosθ 2 ln 1− χ c + 1
=−
( ) el modo anterior.

ingeniería hidráulica en méxico /enero-marzo de 2006 33

 Sotelo-Ávila, G., Perfil del flujo sobre un vertedor tipo cimacio y su comprobación

Cuadro 2a. Posición y características del flujo en la sección crítica para distintos gastos y cargas de operación en el cimacio.

Q = const.; H = Hc H = const.; Q = Qmáx

F = 0.16 H / Hd = 0.5

dκ 2
Hd2 (F κ H )
χc
χc H κHd 2
dx
d
F κ Hd
ζ
x
ζ
Hd tan θ cos θ κ Hd κ Hd H / Hd Hd Hd
−
cos θ
( )
κ Hd ƒ en la
Hd cos θ κHd 2 ƒ en la κ Hd  H
F κ Hd ec. 18
ζ 
ec. 1 ec. 19 ec. 20 ec. 21 cos θ
( ) ec. 18 cos θ  Hd Hd 
 − 
(F κ H )
d
2

0.05 - 0.001959 - 0.072488 0.997383 -1.222647 -1.225855 -3.266579 0.038269 - 0.423443 0.501959 - 0.427837
- 0.046912 0.49497 0.615329 0.037336
- 0.043386

0.06 - 0.002745 - 0.084639 0.996437 -1.186282 -1.190523 -2.802294 0.036026 - 0.411739 0.502745 - 0.415869
- 0.04289 0.495127 0.59853 - 0.000856
0.002815

0.059302 - 0.002686 - 0.083801 0.996507 -1.188616 -1.192782 -2.829297 0.036168 - 0.41249

- 0.04314 0.495128
0.0

0.059748 - 0.002724 - 0.084337 0.996462 -1.187120 -1.191335 -2.811953 0.502724 - 0.416141 - 0.044373
- 0.598912 0.0 0.037246
F = 0.162573

Cuadro 2b. Posición y características del flujo en la sección crítica para distintos gastos y cargas de operación en el cimacio
(continuación).

Q = const.; H = Hc H = const.; Q = Qmáx

F = 0.5016 H / Hd = 1.0

dκ 2
Hd2 (F κ H )
χc
χc H κHd 2
dx
d
F κ Hd
ζ
x
ζ
Hd tan θ cos θ κ Hd κ Hd H / Hd Hd Hd
−
cos θ
( )
κ Hd ƒ en la
Hd cos θ κHd 2 ƒ en la κ Hd  H
F κ Hd ec. 18
ζ 
ec. 1 ec. 19 ec. 20 ec. 21 cos θ
( ) ec. 18 cos θ  Hd Hd 
 − 
(F κ H )
d
2

0.098 - 0.006804 - 0.128435 0.991853 -1.086972 - 1.095901 - 1.952869 0.297271 - 0.768045 1.006804 - 0.778452
- 0.325779 0.986855 - 1.103357 - 0.501545
- 0.002339

0.099 - 0.006933 - 0.129548 0.991713 - 1.084859 -1.093924 - 1.940299 1.006933 - 0.777116

- 1.101508 - 0.002840

0.097375 - 0.006724 - 0.127739 0.99194 - 1.088302 - 1.097145 - 1.960876 0.297998 - 0.768876

- 0.326947 0.986856
0.0

0.09824 - 0.006834 - 0.128702 0.991819 - 1.086464 - 1.095425 - 1.949829 1.006834 - 0.778031 - 0.340176
- 1.102912 0.0 0.310546
F = 0.512915

34 ingeniería hidráulica en méxico /enero-marzo de 2006

 Sotelo-Ávila, G., Perfil del flujo sobre un vertedor tipo cimacio y su comprobación

Cuadro 2c. Posición y características del flujo en la sección crítica para distintos gastos y cargas de operación en el cimacio
(continuación).

Q = const.; H = Hc H = const.; Q = Qmáx

F = 0.7895 H / Hd = 1.33

dκ 2
Hd2 (F κ H )
χc
χc H κHd 2
dx
d
F κ Hd
ζ
x
ζ
Hd tan θ cos θ κ Hd κ Hd H / Hd Hd Hd
−
cos θ
( )
κ Hd ƒ en la
Hd cos θ κHd 2 ƒ en la κ Hd  H
F κ Hd ec. 18
ζ 
ec. 1 ec. 19 ec. 20 ec. 21 cos θ
( ) ec. 18 cos θ  Hd Hd 
 − 
(F κ H )
d
2

0.11 - 0.008425 - 0.141686 0.990111 -1.062683 - 1.073297 - 1.819848 0.703901 - 0.988031 1.338425 - 1.019683
- 0.755495 1.289330 - 1.436526 - 0.031353
0.010679

0.12 - 0.009896 - 0.152562 0.988562 -1.043986 - 1.056066 - 1.733346 0.679351 - 0.972969 1.339896 - 1.004086
0.717439 1.289360 - 1.415019 0.007038
0.015046

0.115796 - 0.009264 - 0.148007 0.989224 -1.051695 - 1.063152 - 1.76744 0.689420 - 0.97918

- 0.732959 1.28939
0.0

0.118059 - 0.009602 - 0.150462 0.988869 -1.047520 - 1.059311 - 1.748709 1.339602 - 1.007013 - 0.805506
- 1.419055 0.0 0.760405
F = 0.832453

Cuadro 3. Valores críticos para tres condiciones de operación del cimacio.

a) Para Q = constante, H = Hc

H0 χc = κ dc H dc
F κ Hd C
χc
Hd Hd Hd Hd
ε

0.50 0.1600 0.059302 - 0.412490 0.495128 0.009743 -1.188615633 0.347034 2.0342

1.00 0.5016 0.097375 - 0.768876 0.986856 0.013144 -1.088301875 0.706490 2.2663
1.33 0.7895 0.115796 - 0.979180 1.289390 0.030534 -1.051695345 0.931049 2.3885

b) Para H = constante, Q = Qc
H Fc
H0
χc = κ dc
χc
Hd Hd Hd
ε

0.50 0.059748 - 0.416141 0.162573 0.504995 0.009891

1.00 0.098240 - 0.778131 0.512915 1.013400 0.013220
1.33 0.118059 - 1.007013 0.832453 1.379000 0.035500

En los cuadros 3a y 3b se presenta el resumen de los

cálculos antes descritos para las condiciones críticas y
Perfil del flujo y distribución de la presión

los valores del factor ε de pérdida. Puede observarse El cálculo de los perfiles del flujo debe entonces realizar-
que la energía del flujo en la zona del cimacio (H/Hd ) es, se para los valores mostrados en los cuadros 3a y 3b, ya
en cualquier caso, siempre menor que la que existe que corresponden a los obtenidos experimentalmente
como carga aguas arriba (H0 / Hd ). La ilustración 3 mues- para las tres cargas de operación. Dicho cálculo se pre-
tra la representación gráfica de los valores de ε según senta en los cuadros 4a, 4b y 4c, donde las característi-
los cuadros 3a y 3b, y se observa que existe congruen- cas geométricas del perfil del cimacio se obtuvieron de
cia entre ambos modos de solución. las ecuaciones (19) y (20), las raíces χ i de la ecuación
(22) y los valores de la presión de la ecuación (12).

ingeniería hidráulica en méxico /enero-marzo de 2006 35

 Sotelo-Ávila, G., Perfil del flujo sobre un vertedor tipo cimacio y su comprobación

Cuadro 4a. Cálculo del perfil del flujo y carga de presión sobre la espalda de un cimacio para tres cargas de operación.

F = 0.16 F = 0.5016 F = 0.7895

H/Hd = 0.495128 H/Hd = 0.986856 H/Hd = 1.28939

H H H
= – –
ζ χ ζ χ ζ χ
tan θ Hd Hd Hd Hd Hd Hd
x κ Hd
cos θ κ Hd d/Hd d/Hd d/Hd
ζ
Hd Hd d ζ /dx cos θ
(F κ Hd)2 p/gρHd (F κ Hd)2 p/gρHd (F κ Hd)2 p/gρHd

- 0.2 - 0.042978 0.539781 0.879986 - 5.0 - 5.68191 0.538106 - 2.931999 1.029834 - 5.623513 1.332368 - 7.281528
0.640000 0.5864 6.290064 1.124703 15.582756 1.456306
0.196694 - 0.729887 - 2.154400

- 0.15 - 0.023030 0.314485 0.953939 - 2.0 - 2.09657 0.518158 - 0.964279 1.009886 - 1.894269 1.312420 - 2.471585
0.1024 0.482139 1.006410 0.947135 2.493241 1.235792
0.293496 0.118785 - 0.297101

- 0.1 - 0.010102 0.204124 0.979796 - 2.0 - 2.041241 0.505230 - 0.884856 0.996958 - 1.781193 1.299492 - 2.338312
0.1024 0.442428 1.006410 0.890596 2.493241 1.169156
0.250357 0.035067 - 0.416261

- 0.05 - 0.002506 0.100504 0.994987 - 2.0 - 2.010076 0.497635 - 0.833849 0.989362 - 1.712621 1.291896 - 2.258586
0.1024 0.416924 1.006410 0.856310 2.493241 1.129293
0.219178 - 0.021258 - 0.494792

0.0 0.0 0.0 1.00 - 2.0 -2.0 0.495128 - 0.815595 0.986856 - 1.689368 1.28939 - 2.231832
0.1024 0.407797 1.006410 0.844684 2.493241 1.115916
0.207252 - 0.041429 - 0.522501

0.059302 - 0.002686 - 0.083801 0.996507 - 1.188616 - 1.192782 0.497815 - 0.412490 0.989542 - 0.875635 1.292076 - 1.167732
crítico 0.036167 0.347034 0.355466 0.736685 0.880617 0.982431
0.19457 0.090935 - 0.179093

0.097375 - 0.006723 - 0.127739 0.99194 - 1.088302 - 1.097145 0.501852 - 0.358121 0.993579 - 0.768876
crítico 0.030321 0.329064 0.296993 0.706491
0.178254 0.077486

Cuadro 4b. Cálculo del perfil del flujo y carga de presión sobre la espalda de un cimacio para tres cargas de operación (continuación).

0.115796 - 0.009264 - 0.148007 0.989224 - 1.051695 - 1.063152 1.298654 - 0.97918

crítico 0.689420 0.931049
- 0.180614

0.2 - 0.025461 - 0.235514 0.973369 - 0.92308 - 0.948334 0.520589 - 0.270981 1.012317 - 0.598773 1.314851 - 0.809320
0.021813 0.293562 0.214384 0.648668 0.531108 0.876761
0.141222 0.03865 - 0.195733

0.4 - 0.091787 - 0.424515 0.920492 - 0.703576 - 0.764348 0.586915 - 0.173367 1.078643 - 0.403248 1.381118 - 0.554756
0.012673 0.246408 0.124548 0.573141 0.308551 0.78848
0.091135 - 0.006475 - 0.203069

0.6 - 0.194334 - 0.599198 0.857796 - 0.535785 - 0.624606 0.689463 - 0.114149 1.181190 - 0.277725 1.483724 - 0.388431
0.007349 0.213049 0.072226 0.518351 0.178931 0.724975
0.060470 - 0.021288 - 0.177686

0.8 - 0.330892 - 0.765188 0.794173 - 0.407233 - 0.512776 0.826020 - 0.076339 1.317748 - 0.193003 1.620282 - 0.274013
0.004245 0.187457 0.041725 0.473939 0.103369 0.672867
0.041542 - 0.022047 - 0.142268

36 ingeniería hidráulica en méxico /enero-marzo de 2006

 Sotelo-Ávila, G., Perfil del flujo sobre un vertedor tipo cimacio y su comprobación

Cuadro 4c. Cálculo del perfil del flujo y carga de presión sobre la espalda de un cimacio para tres cargas de operación (continuación).

F = 0.16 F = 0.5016 F = 0.7895

H/Hd = 0.495128 H/Hd = 0.986856 H/Hd = 1.28939

H H H
= – –
ζ χ ζ χ ζ χ
tan θ Hd Hd Hd Hd Hd Hd
x κ Hd
cos θ κ Hd d/Hd d/Hd d/Hd
ζ
Hd Hd d ζ /dx cos θ
(F κ Hd)2 p/gρHd (F κ Hd)2 p/gρHd (F κ Hd)2 p/gρHd

1.0 - 0.5000 - 0.925 0.734099 - 0.311047 - 0.423712 0.995128 - 0.051953 1.486856 - 0.135634 1.7894 0.195120
0.002477 0.167025 0.024343 0.436056 0.060305 0.627302
0.029599 - 0.017856 - 0.108677

1.2 - 0.700576 - 1.080055 0.679389 - 0.239905 - 0.353119 1.195704 - 0.036063 1.687432 - 0.096681 1.989966 - 0.140672
0.001473 0.150323 0.014481 0.402997 0.035874 0.586366
0.02183 - 0.012766 - 0.080911

1.4 - 0.931766 - 1.231262 0.630440 - 0.187316 - 0.297118 1.426894 - 0.025558 1.918622 - 0.07003 2.221156 - 0.102884
0.000898 0.136445 0.008828 0.373861 0.021870 0.549254
0.016604 - 0.008267 - 0.059364

1.6 - 1.192867 - 1.379252 0.586984 - 0.148191 - 0.252461 1.687995 - 0.018489 2.179723 - 0.051579 2.482257 - 0.0764
0.000562 0.124765 0.005525 0.348059 0.013688 0.515551
0.012972 - 0.004731 - 0.043150

1.8 - 1.483283 - 1.524486 0.548486 - 0.118787 - 0.216572 1.978412 - 0.013640 2.470139 - 0.038622 2.772673 - 0.057606
0.000361 0.114828 0.003550 0.325137 0.008795 0.484950
0.010372 - 0.002113 - 0.031128

Ilustración 4. Valores adimensionales teóricos y experimentales para tres cargas de operación.

1.4 0.6
Puntos experimentales H0 Puntos experimentales
1.2 p =0.5
Hd H0 / Hd = 0.5
H0 H0 / Hd = 0.5 gρH0 0.4
=1.33 H0 / Hd = 1.0
ζ 1.0 Hd H0 / Hd = 1.0 H0 / Hd = 1.33
Hd 0.8 0.2
H0 / Hd = 1.33 Valor teórico de p
Críticos
0.6 Perfil teórico 0
de la superficie
0.4 del agua
H0
=0.5 -0.2
0.2 Hd H0
=1.0
0 -0.4 Hd
H0 H0
=1.0 =1.33
-0.2 Hd Hd
0.5 Hd -0.6
-0.4 0.2 Hd 0.2 0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8
-0.6 x /Hd

-0.8
-1.0 b) Distribución de la presión sobre la espalda del cimacio
-1.2
-1.4
-1.6
-1.0 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 2.2
x /Hd

a) Perfiles de la superficie libre sobre la espalda del cimacio

ingeniería hidráulica en méxico /enero-marzo de 2006 37

 Sotelo-Ávila, G., Perfil del flujo sobre un vertedor tipo cimacio y su comprobación

en ella. Es obvio que en ningún caso dicha sección críti-

ca coincide con la vertical que pasa por la cresta, como
Resultados obtenidos

Los perfiles del flujo calculados para las tres cargas de ha sido común considerarlo en la práctica, con la acep-
operación en los cuadros 4a a 4c se dibujan con línea tación además de un flujo rectilíneo formado sobre dicha
llena en la ilustración 4a, donde también aparecen los cresta, que dista mucho de ser correcto. La profundidad
puntos que fueron obtenidos experimentalmente. La ilus- del flujo medida experimentalmente en dicha sección
tración 4a contiene asimismo el perfil del cimacio ajusta- vertical es cercana al valor de 0.75 H, distinto de 2H/3
do al sistema de ejes coordenados representado en la que resulta de la condición crítica del flujo rectilíneo.
ilustración 1. Se observa una excelente concordancia en- El modelo de flujo muestra, una vez más, que es
tre los perfiles teóricos y los experimentales para la re- capaz de predecir de manera confiable la distribución de
gión aguas abajo de la cresta, no así para la región antes la presión sobre el cuerpo del cimacio, variable de mu-
de ella, donde ocurren pequeñas diferencias. A pesar de cha importancia en el problema. Por lo que se refiere a
que en la región antes de la cresta los valores de χi son los valores teóricos, se determinaron valores positivos
mayores de dos para x/ Hd = – 0.2, se observa una bue- de la presión, casi iguales a cero y negativos, para una
na concordancia con los valores experimentales al efec- carga de operación menor, igual y mayor que la de dise-
tuar el trazo dentro de la ilustración. ño. Estos valores coinciden con los experimentales en
La ilustración 4b contiene los valores de x/ Hd en el magnitud, pero en puntos con posiciones ligeramente
eje horizontal y los de p/ gρHd en el eje vertical. En ella distintas. Las mayores diferencias en posición ocurren
se muestra una comparación gráfica entre los valores prácticamente sobre la cresta y antes de ella, es decir, el
teóricos correspondientes a los cuadros 4a, 4b, 4c y los mismo comportamiento observado para el perfil del flujo.
experimentales de la presión relativa para las mismas
tres cargas de operación, donde nuevamente se obser-
van las mayores discrepancias antes de la cresta, ya que
Anexo

los valores de χ se acercan o rebasan el valor límite (–2) Demostración de la ecuación (18).
acordado para la validez de la hipótesis del flujo poco La ubicación de la sección crítica depende de ζ c,
profundo. es decir, de la coordenada xc de la sección donde
Las mayores discrepancias se presentan ligeramente dH/dx = 0 y F = Fc. Para encontrarla, se sustituye
detrás de la cresta, pero una vez alcanzada la sección χ = κd en la ecuación (8) y ésta se convierte en:
crítica existe una gran coincidencia entre los valores
2g
teóricos y los experimentales.
(H − ζ) − 2gκcos
−2
[(1− χ) ln(1− χ)] =
κ q
2 2
q
3 2
θ
χ (A.1)

la cual se escribe de modo simplificado en la forma:

Conclusiones

La aplicación del modelo de flujo curvilíneo a curvas

verticales convexas de geometría conocida, como en un γ χ = A + Bχ
() (A.2)
cimacio, conduce efectivamente a muy buenos resulta-
dos, ya que, a pesar de las pequeñas diferencias entre donde:
los valores calculados y los valores experimentales, ha-
2g
bría que pensar que ambos contienen imperfecciones A= ( H − ζ) (A.3)
propias de cada método. Las mayores discrepancias se κ2 q 2
observan en el perfil del flujo antes de la cresta, donde
resulta difícil definir la sección ortogonal por la geometría 2gcosθ
B=− (A.4)
que el fondo tiene en esta zona. Además, como ya se κ3 q 2
explicó, se viola la hipótesis de flujo poco profundo que −2
sustenta la teoría. No obstante, puede decirse que las γ χ = 1− χ ln 1− χ
( ) [( ) ( )] (A.5)
diferencias no son importantes, pues en esa región no
se requiere mayor precisión en la determinación del perfil Las soluciones de la ecuación A.2 se muestran en la
del flujo. ilustración A.1 para los valores: κ > 0 en el primer cua-
Por lo que se refiere a la sección crítica en la proximi- drante del sistema coordenado [ χ, γ ( χ)]; κ < 0 en el
dad de la cresta, la teoría permite detectar su posición segundo. Las rectas de pendiente ±B intersectan a la
para cada valor de la carga de operación, así como las curva de ecuación A.2 en dos puntos de cada cuadran-
características de tirante, velocidad y energía específica te: χ± en régimen subcrítico, χ±0 en supercrítico.

38 ingeniería hidráulica en méxico /enero-marzo de 2006

 Sotelo-Ávila, G., Perfil del flujo sobre un vertedor tipo cimacio y su comprobación

Sólo el punto χ±c corresponde a la condición crítica y Al elevar al cuadrado, sustituir d = χ /κ y obtener el recí-
ésta se obtiene de la derivada de la ecuación (A.5) de la proco de los términos, para κ > 0 ( χ > 0), resulta:
siguiente manera:

2gcosθ 2 1+ ln 1− χ
[ ( )] χ
dγ χ
( ) = −2 −3  1− χ
[(1− χ) ln(1− χ)] − ln 1− χ 

κ q
χ=−
( ) 3 2 3
dχ  1− χ
−
 [(1− χ) ln (1− χ)]
o bien:
O bien, de las ecuaciones (A.4) y (A.6), la anterior se
dγ χ 2 1+ ln 1− χ
( ) = [ ( )] convierte en:
(A.6)
dχ 3
dγ χ
[(1− χ) ln(1− χ)] −B =
()
dχ
Para valores positivos de κ el valor mínimo de γ( χ)
se alcanza cuando d γ( χ)/d χ = 0, lo que ocurre para Cuando κ < 0 ( χ < 0) se invierte el signo de B, de
χ = 0.6321, como se muestra en la ilustración A.1. manera que cuando F = Fc se cumple que:
El número de Froude para el canal rectangular está
dado por la ecuación (13), donde se sustituye u0 toma-
1 dγ χ ()
do de la ecuación (7) y resulta: =1
B dχ
qκ / ln 1− χ
( )
F =− cualquiera que sea el signo de κ.
gdcosθ Se sustituye ahora la ecuación (A.5) en la (8) y la
energía total se expresa en la forma:
Cuando ocurre el estado crítico F = Fc y de la ecua-
ción (14) se obtiene:
( ) q2gκ
2 2
H = ζ + dcosθ + γ χ (A.8)
1/ 2
3
κq / ln 1− χ
) = − (1− χ) ln(1− χ) 
 
( La curvatura κ en una determinada sección es cons-
 χ [ ln(1− χ) + 1]  tante e independiente de d, pero θ es función de x y la
−
gdcos θ
 
c derivada respecto de la coordenada x de la ecuación
anterior resulta:

dH dς dd dθ
χ) = A + B χ; χ± subcrítico,
Ilustración A.1. Solución de γ (χ

+ cosθ − d senθ
χ±c crítico, χ±0 supercrítico.
dx dx dx dx
=
(A.9)
Canal rectangular

γ (χ)
dγ χ
( ) qgκ ddκx + q2gκ
2 2 2
20 ()
dx
+γ χ
18
A + Bχ para κ<0
16
2 2
14 Debido a que: κ = 1/R; R dθ = (dx) + (dς) =
12
2
10 dx 1+ dς / dx ; se tiene:
( )
A + Bχ para κ<0 Pendiente B
8
6
dθ  sen θ 2
Pendiente B
4 d senθ = κ d sen θ 1+ 
dx  cos θ 

2

-1.0 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1.0 sen θ dζ
χ
χ- χc- χ0- = χtanθ = χ
cos θ dx
χ0 χ+c
+
χ+ =χ
0.6321

ingeniería hidráulica en méxico /enero-marzo de 2006 39

 Sotelo-Ávila, G., Perfil del flujo sobre un vertedor tipo cimacio y su comprobación

Por otra parte, vale también que: Simbología

dd d  χ  1 dχ χ dκ b longitud de cresta del cimacio, en m
dx dx  κ  κ dx κ2 dx d tirante medido en dirección perpendicular al
=  = −

fondo, en m
dγ χ ( ) = d γ(χ) dχ dc tirante en la sección crítica, en m
dx dχ dx g aceleración de la gravedad, en m/s2
he pérdida de energía entre la existente aguas arriba
Al sustituir las tres últimas expresiones en la ecuación
del cimacio y la energía crítica, en m
(A.9) y agrupar términos, se obtiene:
p0 presión en el fondo de cualquier sección perpen-
dicular al fondo, en N/m2
dH dζ  q 2 κ3 d γ χ ( ) 

= 1− χ
( + 1+ ) q gasto unitario sobre la cresta (Q/b), en m2/s
dx dx  2gcosθ dχ  u0 componente de la velocidad local en el fondo de
la sección, en dirección tangente a dicho fondo,


cosθ dχ q 2 κ  2gcosθ  dκ en m/s

2γ χ − () x coordenada horizontal, cuyo valor cero coincide
κ dx 2g  q 2 κ3  dx
+ χ
con la cresta del cimacio, en m
C coeficiente de descarga del cimacio, en m1/2/s
Pero de las ecuaciones (A.2) a (A.4) resulta:
F parámetro adimensional relacionado con el flujo y
la geometría del cimacio, según ecuación (3)
dH dζ  1 dγ χ ( )  cos θ dχ + q2 κ H energía total del flujo (carga) respecto de un plano

= 1− χ
( + 1−
)
dx dx  B dχ  κ dx 2g horizontal que pasa por la cresta, en m
Hd carga sobre la cresta (energía total del flujo) elegi-

dκ da para el diseño del perfil geométrico del cima-
[ 2 A + 2Bχ + Bχ ]
dx cio, en m
H0 carga sobre la cresta con que opera el cimacio
Cuando dH/dx = 0, F = Fc , es decir, vale la ecuación
para verter el gasto Q, medida como el desnivel
(A.7); por tanto, la expresión anterior se convierte en:
entre la superficie libre aguas arriba y la cresta del
mismo cimacio
(1− χ c ) ddζx + q2gκ [2 A + 3Bχ c ] ddκx = 0
2
(A.10) Q gasto de operación vertido sobre el cimacio, en
m3/s
coeficiente de pérdida de energía entre la que
Se sustituyen los valores de A y B tomados de las
ε
existe aguas arriba del cimacio y la cresta del
ecuaciones (A.3) y (A.4) y se obtiene:
mismo, sin dimensiones
curvatura del fondo de la sección, en m-1
dζ q 2 κ  4 g 6 g cos θ  dκ
κ
(1− χc ) H−ζ − ( ) =0 densidad del agua, en kg/m3
dx 2g  q 2 κ2  dx
χc  ρ
q 2 κ3
+ 
θ ángulo de inclinación de la tangente al perfil del
cimacio en una sección cualquiera perpendicular
que al simplificar resulta:
al fondo
χ variable adimensional equivalente a κd
dζ  2Hd2  H ζ  3cos θ 2  dκ
 
valor de χ en la sección crítica
(1− χc ) dx  κ Hd  Hd Hd  κ H 2
+  − − Hd χ c 
dx
=0 χc
ζ coordenada vertical cuyo valor cero coincide con
 d  ( ) la cresta del cimacio, en m
Finalmente, en términos adimensionales vale la ζc desnivel entre la cresta y el fondo de la sección
expresión: crítica, en m

Recibido: 25/06/2004
(1− χc ) ddζx + 2 HH ζ  3cos θ  Hd2 dκ
= 0 (A.11)
 
−
Hd  κ Hd
−
 κ Hd dx
χc  Aprobado:11/05/2005
d

que es la condición que debe cumplirse para encontrar

la posición de la sección crítica e igual a la ecuación (18)
del texto.

40 ingeniería hidráulica en méxico /enero-marzo de 2006

 Sotelo-Ávila, G., Perfil del flujo sobre un vertedor tipo cimacio y su comprobación

Referencias SIVAKUMARAN, S. S., T. TINGSANCHALI, R. J., HOSKING.

Steady shallow flow over a spillway, Journal of fluid mechan-
BELANGER, J. B. CH., Notes sur le cours d’hydraulique. ics, 1981, núm. 111, pp. 411-420.
Memories cole Nationale des Ponts et Chaussées, 1849- SOTELO, G. y C., ESCALANTE. Ecuaciones generales del flujo
1850, pp. 32-33. impermanente en canales de fondo curvo. Ingeniería
DRESSLER, Robert F., New nonlinear shallow-flow equations Hidráulica en México. Segunda Época, abril-junio, 2001 a,
with curvature, Journal of Hydraulic Research, IAHR, 1978, vol. XVI, núm. 2, pp. 27-42.
vol.16, num. 3, pp. 205-222. SOTELO, G. y C., ESCALANTE. Ecuaciones del flujo imperma-
FAWER, C. Etude de quelques écoulements permanent à filets nente en canales trapeciales de fondo curvo. Ingeniería
courbes, Tesis Université de Laussane, Laussane, 1937. hidráulica en México. Segunda Época, octubre-diciembre,
JAEGER, Charles. Engineering Fluid Mechanics, Traducción de 2001 b, vol. XVI, núm. 4, pp. 45-51.
la edición en alemán de 1949. Londres, Glasgow. Blackie & SOTELO, G. Flujo crítico en canales trapeciales de fondo curvo.
Son Limited, 1949, 529 pp. Ingeniería hidráulica en México. Segunda Época, enero-
KNAPP F. H. Überfall und Durchfluss im Wasserbau. Karlsruhe, marzo, 2004, vol. XIX, núm. 1, pp. 27-40.
Verlag G. Braun, 1960, 671 pp. US ARMY CORPS OF ENGINEERS, Corps of Engineers.
PRÁSIL, Franz. Technische Hydrodynamik, Segunda edición, Hydraulic Design Criteria, Waterways Experiment Station,
Berlín, Springer Verlag, 1926. Vicksburg, Miss., 1978.

ingeniería hidráulica en méxico /enero-marzo de 2006 41

 Sotelo-Ávila, G., Perfil del flujo sobre un vertedor tipo cimacio y su comprobación

Abstract

SOTELO-ÁVILA, G. Flow profile on an overflow spillway and its verification. Hydraulic engineering in Mexico (in
Spanish). Vol. XXI, no. 1, January-March, 2006, pp. 29-42.

Many authors have shown different approaches to predict the flow profile and the pressure distributions of water
flowing on an overflow spillway. The authors have used flow nets, simplified hypothesis in the potential theory, the
shallow-flow equations with curvature by Dressler and other procedures based on the theory of vortex-free flow.
Author shows the application of the theory of flow on a curved bed allows predicting the location of the critical
section around the crest, the flow profile and the pressure distribution when the shape of the spillway is known.

Keywords: overflow spillways, curvilinear flow over a spillway.

Dirección institucional del autor:

Dr. Gilberto Sotelo-Ávila

Doctor en Ingeniería.
Profesor Emérito.
Universidad Nacional Autónoma de México,
Facultad de Ingeniería,
División de Ingeniería Civil, Topográfica y Geodésica,
Departamento de Ingeniería Hidráulica,
Ciudad Universitaria, Circuito Exterior,
Coyoacán, 04510, México, D. F. México,
teléfono: + (52) (55) 56 22 80 11,
fax: + (52) (55) 55 50 77 09,
[email protected]

42 ingeniería hidráulica en méxico /enero-marzo de 2006