Metodo Mostkow
Metodo Mostkow
Metodo Mostkow
gb-
(6)
La energa especfica E
o
al inicio de la reja,
que se considera constante en cualquier
seccin del canal es igual a:
E
o
= o
i
y
o
+
v
c
2
2g
= y +
`
2gb``
(7)
Figura 2. Flujo sobre una reja de fondo.
Tomado de Sotelo(2002).
, de donde se obtiene que el caudal es:
= by2g(E
o
-y) (8)
Sustituyendo la ecuacin (8) en la (6), se
obtiene la siguiente ecuacin para determinar
el perfil de flujo sobre la rejilla:
Jx =
g b(3-2L
c
)
2L
c
- [
-dQ
dx
Jy (9)
El caudal Q
o
que conduce el ro es igual al de
la seccin inicial de la reja, y por lo tanto, de
la ecuacin de energa especfica en dicha
seccin se tiene que:
o
= by
o
2g (E
o
- oy
o
) (10)
Para la seccin al final de la reja, de la
ecuacin (8), el caudal Q
L
es igual a:
L
= by
L
2g (E
o
-y
L
) (11)
Por lo tanto, el caudal captado o retirado Q
r
,
por una reja de longitud L es igual a:
=
o
-
L
(12a)
, o a partir de las expresiones (11) y (12):
=
o
_1 -
L
_1-
j
L
E
c
c
_1-
oj
c
E
c
_ (12b)
Segn el trabajo realizado por Mostkow [8],
en el caso de rejas construidas con barras
paralelas a la corriente, como en la Figura 2c,
el flujo a travs de las aberturas es
aproximadamente vertical, la prdida de
energa es despreciable y la carga efectiva
sobre ellas es prcticamente igual a la
energa especfica E
o
con que el flujo llega a
la reja.
Mientras que en el caso de rejillas formadas
a base de una placa perforada (Figura 2d), la
direccin del flujo a travs de los orificios
tiene un ngulo apreciable con la vertical y se
ve afectado por los sedimentos incrustados
en las perforaciones, que producen una
prdida de energa apreciable por el cambio
de direccin, de inclinada eventualmente a
vertical. Dicha prdida es casi igual a la
carga de velocidad del flujo sobre la reja, por
lo que la carga efectiva es igual a la carga
esttica o se puede confundir con el tirante.
Expuestas estas razones, se hace necesario
presentar por separado las expresiones
especficas para determinar la longitud de la
rejilla segn su tipo.
2.1.1 Flujo de barras paralelas al flujo. En
general, la tasa de cambio del caudal unitario
derivado a travs de la rejilla est dada por
una expresin que depende de la carga
hidrulica sobre sta, afectada por un
coeficiente de descarga, el cual depende de
las caractersticas hidrulicas del flujo de
aproximacin, la geometra de la rejilla
(longitud, pendiente y orientacin respecto al
flujo), y las barras que forman la rejilla
(forma, tamao y espaciamiento).
En el caso de barras paralelas, la carga
hidrulica efectiva sobre las aberturas de la
reja es la energa especfica del flujo sobre
ellas y la variacin del caudal desviado a
travs de la reja se expresa de la siguiente
forma:
-d
dx
= eC
d
b2gE
o
= mb2gE
o
(13)
, donde:
E
o
=Energa especfica del flujo a la entrada
de la reja.
b=Ancho total de la rejilla, medido en
direccin perpendicular al flujo. (Se est
haciendo b
rejilla
=b
canal
=b)
m=coeficiente global de descarga (m=C
d
)
C
d
=Coeficiente de descarga a travs del
espacio entre dos barras consecutivas.
A
n
=rea neta de paso a travs de las barras.
=(1 - )
A
n
A
t
, es igual al cociente del rea
de paso entre rejas y su rea total o
coeficiente de porosidad afectado por el
factor de obturacin.
A
t
=rea total de rejas.
f=Coeficiente de obstruccin, originado por
las arenas y gravas que se incrustan y que
se toma de 5 a 30%.
En la reja de barras paralelas se considera
que el coeficiente de porosidad p, es igual a:
p = [
u
u+c
(14)
, donde e y a corresponden a los tamaos del
espesor y la abertura entre barras,
respectivamente, como se indica en la
Figura 2c.
Al sustituir la ecuacin (13) en la ecuacin (9)
y reorganizar, se obtiene que:
Jx =
3
2
j
E
c
-1
m_1-[
j
E
c
Jy (15)
Al integrar la ecuacin (15), y determinando
la constante de integracin de la condicin de
frontera a la entrada de la rejilla donde: x=0,
y=y
o
, la ecuacin para el clculo del perfil
de flujo a lo largo de una rejilla de barras
paralelas es igual a:
x =
1
m
_y
o
_1 -o [
L
c
- y _1 -[
L
c
_ (16)
De la condicin de frontera a la salida de la
rejilla, donde x=L
m
, se obtiene la ecuacin
para determinar la longitud mnima necesaria
para desviar todo el caudal Q
o
del ro a travs
de la reja al alcanzar el tirante y
L
=0 (Ver
Figura 2b):
I
m
=
c
m
_1 -o [
c
L
c
(17)
De esta forma, para una reja de longitud x=L
dada y caudal Q
o
mayor en el ro, de la
ecuacin (16) se obtiene el tirante y=y
L
al
final de la reja y de la ecuacin (12b) se
obtiene el caudal realmente captado con la
longitud L de la reja.
Para derivar una expresin que permita
determinar la longitud de la rejilla en trminos
del caudal captado, puede considerarse que
en la seccin de entrada a la rejilla: y=y
o
, y
que de la ecuacin (10), con Q
o
=Q
r
, se tiene
que al despejar y
o
, tras sustituir en la
ecuacin (17), se obtiene que:
I
m
=
r
mb2gL
c
(18)
, donde Q
r
es el caudal captado por la reja y
Q
o
el de la seccin inicial, es decir, el caudal
del ro.
2.1.2 Placa perforada. En este caso, la
carga efectiva sobre la placa perforada est
dada por el tirante y la variacin del caudal
desviado a travs de la placa, y se puede
expresar como:
-d
dx
= eC
d
b2gy = mb2gy (19)
, donde los trminos tienen igual
interpretacin que la dada para la ecuacin
(13).
Sustituyendo la ecuacin (19) en la ecuacin
(9), tras simplificar se obtiene:
Jx =
1
m
l
l
l
l
3
2
_[
j
E
c
d
_1-[
j
E
c
-
d
_[
j
E
c
-[
j
E
c
2
1
1
1
1
(20)
Luego de la integracin de la ecuacin (20),
determinando la constante de integracin de
la misma forma que en el caso de la rejilla de
placas paralelas, se obtiene que:
x =
L
c
m
_
1
4
orcscn j[1 - 2
L
c
-orcscn [1 - 2
u
c
L
c
[ -
3
2
__
L
c
[1 -
L
c
- _
u
c
L
c
[1 -
u
L
c
__ (21)
Nuevamente, para y
L
=0, se obtiene que
x=L
m
, tal que la longitud mnima necesaria
para la captacin del total de caudal, es igual
a:
I
m
=
L
c
m
_
3
2
_ _
u
|j
c
L
c
[1 -
u
|
L
c
-
1
4
orcscn [1 -
2
ui
c
L
c
_ +
n
8
_ (22)
Para el vertedero tirols Mostkow[8], propone
que la relacin entre el caudal Q
o
y la energa
especfica E
o
del flujo que llega a la rejilla,
est dada por la ecuacin de descarga de un
vertedero de cresta ancha:
o
= C
bE
o
1.5
(23)
La Tabla 1 muestra tambin valores de la
relacin E
o
/E
L
de energa especfica al
principio y final de la reja en caso de
captacin parcial, as como los valores
tpicos de los dems parmetros que
intervienen en el clculo de los dos tipos de
rejas analizadas.
Tabla 1. Valores experimentales de C
V
, y
o
/E
o
,
, E
o
/ E
L
, y C
d
para el flujo sobre rejas,
segn Mostkow (1959)
Tipo de
reja
De barras
paralelas
De placa
perforada
Inclinacin
Horizontal
Pendiente 1:5
Horizontal
Pendiente 1:5
C
v
1.55
1.55
1.634
1.528
y
o
/E
o
0.509
0.449
0.594
0.496
'
0.850
0.615
0.970
0.750
E
o
/E
L
0.640
0.609
0.815
1.130
C
d
0.497
0.435
0.800
0.750
3. DIMENSIONAMIENTO Y PERFIL DE
FLUJO SOBRE REJILLAS
Las ecuaciones las ecuaciones (16) y (21)
corresponden a las ecuaciones para el
clculo de los perfiles de flujo a lo largo de
rejillas de barras paralelas y de placa
perforada respectivamente, las cuales al ser
evaluadas para las condiciones de frontera
que definen en inicio y el fin de la rejilla
toman la forma presentada por las
ecuaciones (17) y (22) que permiten el
clculo de la longitud de la rejilla.
Si bien el clculo de la longitud de la rejilla no
reviste mayor complejidad pues una vez
determinada la pendiente de la rejilla es
posible conocer los valores de los diferentes
parmetros necesarios a partir de la Tabla 2,
el procedimiento de determinar un valor
definitivo para la longitud de la rejilla es un
procedimiento de ensayo y error que puede
tornarse tedioso en la medida que se deben
realizar diferentes evaluaciones de longitud
de rejilla.
Para facilitar esta evaluacin, en este trabajo
se construyeron las relaciones grficas para
la determinacin de la longitud de la rejilla de
barras paralelas o de placa perforada (las
cuales se presentan en el Anexo 1) en
trminos de la inclinacin de la rejilla, el
caudal unitario que pasa sobre ella, as como
de diferentes valores del factor de obturacin
f, para diferentes configuraciones de rejillas.
En la construccin de estos grficos se ha
considerado como lmite mximo para el
caudal unitario una magnitud de 0.5 m/s, tal
que, la longitud asociada, se considera la
longitud mxima recomendada en trminos
constructivos.
3.1 CLCULO DEL PERFIL DE FLUJO
SOBRE LA REJILLA
Para mostrar la importancia de determinar el
perfil de flujo para diferentes caudales en el
dimensionamiento de una rejilla, se
considerar un canal rectangular en el cual
se dispone de una bocatoma de fondo, con el
objeto de captar, para propsitos de
suministro a un sistema de acueducto, un
caudal unitario de 0.15 m/s, en una fuente
superficial donde el caudal unitario sobre la
rejilla puede alcanzar durante una creciente
un valor mximo de 4.5 m/s.
Del anlisis de sedimentos de la fuente, se
asigna un factor de obturacin del 25%, para
una rejilla de barras paralelas con espesor
e= 1 cm y de espaciamiento entre las barras
a= 1.5 cm, con una inclinacin del 10%.
Con los datos as definidos, de la Figura 5 se
obtiene una longitud mxima de la rejilla de
45 cm; luego, el clculo del perfil de flujo se
realiza con la ecuacin (16) y de la Tabla 1
se tiene que: C
v
=1.55, C
d
=0.435, =0.615.
A partir de la longitud de rejilla seleccionada,
la energa y la profundidad de flujo al inicio
de la rejilla en trminos de cualquier caudal
sobre ella est dadas por las ecuacin (23) y
el valor de y
o
/E
o
=0.449.
E
o
= u.747q
o
23
(24)
y
o
= u.S71q
o
23
(25)
De esta forma, en la ecuacin (16), con
x=L=0.45, y=y
L
, se obtiene que:
y
L
_1 -1.SS9 y
L
q
o
-23
+u.Su9 q
o
23
-u.u882 = u (26)
As mismo, de la ecuacin (12b), se tiene
que:
q
= q
o
_1 - S.2S6Sq
o
-23
y
L_
1 -
1.339
L
c
23
(27)
Tabla 2. Caractersticas de la rejilla para
caudal variable
qo
(m/s)
Eo
(m)
yo
(m)
Vo
(m/s)
yL
(m)
qr
(m/s)
qr/qo
0.15 0.211 0.105 1.43 0.000 0.150 1.000
0.3 0.335 0.166 1.80 0.054 0.193 0.644
0.6 0.531 0.264 2.27 0.157 0.241 0.401
0.9 0.696 0.346 2.60 0.248 0.278 0.309
1.2 0.844 0.419 2.86 0.334 0.307 0.255
1.8 1.105 0.549 3.28 0.498 0.347 0.193
2.4 1.339 0.665 3.61 0.647 0.385 0.160
3 1.554 0.772 3.89 0.792 0.412 0.137
3.75 1.803 0.895 4.19 0.962 0.447 0.119
4.5 2.036 1.011 4.45 1.130 0.472 0.105
De esta forma, a partir de la solucin de las
ecuaciones (24), (25), (26), (27), se
construye la Tabla 2 para diferentes caudales
unitarios entre 0.15-4.2 m/s.
La Figura 3, muestra los perfiles de flujo
sobre la rejilla para q
o
=0.15, 0.6 y 1.2 m/s, a
partir de la ecuacin (17) con m=0.196:
x = S.1u2u.8SS6y
o
- y
_
1 -
L
c
(28)
La Tabla 2 presenta la variacin del caudal q
r
captado para diferentes caudales de
circulacin sobre la rejilla q
o
, mientras que la
Figura 3 indica los niveles del flujo sobre la
rejilla para los diferentes caudales asociados.
De su anlisis conjunto se observa como un
valor exagerado de la longitud de la rejilla
genera una captacin de un caudal mayor al
requerido, lo cual tiene implicaciones en la
dinmica fluvial de la obra aguas abajo de la
estructura de captacin y que permite
disponer de valores ms confiables del
caudal de excedencia bajo el cual se deben
disear las posteriores obras de alivio, as
como los vertederos de la obra.
Figura 3. Perfiles de flujo para diferentes
caudales sobre la rejilla.
4. CONCLUSIONES
La condicin de flujo espacialmente variado
que se presenta sobre la rejilla de una
captacin de fondo no cuenta con una
solucin analtica que permita determinar la
variacin de la superficie de libre sobre la
estructura de captacin.
Esta situacin ha generado, que el estudio
del flujo sobre las rejillas se realice a partir
del anlisis experimental, produciendo
expresiones de clculo que tienen algunas
hiptesis que pueden ser cuestionables en
algunos puntos, y que ha generado que la
longitud de las rejillas se asigne en la
prctica ms con criterios empricos o
constructivos, que con base en la aplicacin
de metodologas que consideren aun con
limitaciones una aproximacin fsica
coherente del fenmeno.
En esta medida, las relaciones grficas
presentadas en el Anexo 1 que corresponden
al mtodo de Mostkow, facilitan
enormemente la labor de la determinacin de
la longitud de la rejilla, la cual se determina a
partir de los grficos construidos, una vez
que se define el tipo de rejilla que se desea
instalar nicamente en trminos del caudal
unitario que ser captado. De esta forma,
rpidamente, se puede conocer para
diferentes caudales del rgimen fluvial la
longitud mnima necesaria de la rejilla, de
forma tal que la obra de captacin se disee
de manera eficiente para el caudal que
realmente se desea captar, considerando a
su vez en el dimensionamiento de los
vertederos el rgimen de variabilidad del
caudal superficial mediante la inclusin del
anlisis de los perfiles de flujo.
Bajo este mismo enfoque, se considera a
futuro necesario evaluar el dimensionamiento
de los drenajes de carreteras que en
ocasiones presentan insuficiencias en la
captacin de los excedentes de precipitacin,
favoreciendo la desestabilizacin de las vas
y llegando inclusive a ser un motivo de la
accidentalidad vial.
5. REFERENCIAS BIBLIOGRFICAS
[1] Mizuyama, T., and Mizuno, H. Behavior of
debris flow at control structures. Proc., IAHR
Int. Workshop on Floods and Inundations
related to Large Earth Movements, C2.1
C2.12. 1994.
[2] Viparelli, C. Dissipatori a griglia di fondo.
Energ. Elettr., 7, p.509519. 1963.
[3] Righetti, M., Lanzoni, S. Experimental
Study of the Flow Field over Bottom Intake
Racks, Journal of Hydraulic Engineering, Vol.
134, No. 1, January 1, p.15-22. 2008.
[4] Sotelo A., G. Hidrulica de canales
abiertos. 1
ra
edicin, UNAM. 836 p. 2002.
[5] Ackers, P.A. A theoretical consideration of
side weirs as stormwater overflows.
Proceedings of the institution of Civil
Engineers, V.6, pp.305, Londres, 1957.
[6] De Marchi G. Saggio di teoria del
funzionamiento degli stramassi laterali.
LEnergia Elletrica, Miln, V.II, No. 11,
pp.849, Noviembre. 1934.
[7] Frank, J., Von Obering, E. Hydraulische
Untersuchungen fr das Tiroler Wehr. Der
Bauingenieur, 31, Helf 3. p 96-101. 1956.
[8] Mostkow, M. Sur le calcul des grilles de
prise deau (Theorical study of bottom type
water intake). La Houille Blanche, V.4, pp.
569-576. 1957.
[9] Noseda, G. Operation and design of
bottom intake racks. Proc., VI General
Meeting IAHR, Vol. 3, The Hague, C17-1
C17-11. 1955.
[10] Zamarin, E.A., Popow, K.W., Fandejew
W.W. Wasserbau. VEB Verlag fr
Bauinbenieur, Heft 3, pp.96-101. 1956.
[11] Hager, W. Discussion of Experimental
Study of the Flow Field over Bottom Intake
Racks Righetti and Lanzoni,S. Journal of
Hydraulic Engineering, January 2008. Vol.
134, No. 1, pp. 15-22. JHE, ASCE/October
2009/ p. 865-868.
[12] Venkataraman, P., Nasser, M. S., and
Ramamurthy, A. S. Flow behavior in power
channels with bottom diversion works. Proc.,
XVIII IAHR Conf., Vol. 4, Cagliari Italy, p.
115122. 1979.
[13] Brunella, S., Hager, W., Minor, H.E.
Hydraulics of Bottom Rack Intake. Journal of
Hydraulic Engineering, Vol. 129, No. 1,
January 1, p-2. 2003.
[14] Castillo Elsitdi, L.G; Lima Guamn, P.
Anlisis del dimensionamiento de la longitud
de reja y del canal lateral en una captacin
de fondo - XXIV congreso latinoamericano de
hidrulica, Punta del este, Uruguay,
Noviembre 2010.
[15] Subramanaya, K. Flow in open
channels. 3
rd
ed. McGraw-Hill. 549 p. 2009.
[16] French, R. Hidrulica de canales
abiertos. Mc-Graw Hill. 794 p. 1998.
[17] Chow, V. T. Hidrulica de canales
abiertos. McGraw Hill. 406 p. 1994
[18] Naudascher E. Hidrulica de Canales.
Ed. Limusa. 381 p. 2001.
[19] Hager W.H. Open channel hydraulics
of flows with increasing discharge. Journal of
Hydraulic Research. IAHR, V.21, No.3,
pp.177-193, 1983.
[20] Matern, Hernn Obras hidrulicas
rurales Editorial Universidad del Valle
p.323.
5. ANEXO. Relaciones grficas para la determinacin de la longitud de la rejilla en una
bocatoma de fondo
Se presentan a continuacin las relaciones
grficas de las ecuaciones (17) y (22), que
permiten determinar la longitud de rejilla L de
barras paralelas o de placa perforada, en
trminos del caudal unitario q que pasa sobre
ella segn el trabajo de Mostkow[8].
Las figuras 4, 5 y 6 se aplican a rejillas de
barras paralelas para inclinaciones de 0, 10
y 20, donde e es el espesor de las barras, a
es la abertura entre las barras y f
corresponde al porcentaje de obstruccin de
la rejilla el cual vara entre el 5% y el 30%..
.
Figura 4. Longitud de rejilla de barras inclinadas para =0
Figura 5. Longitud de rejilla de barras inclinadas para =10
Figura 6. Longitud de rejilla de barras inclinadas para =20
Las figuras 7, 8 y 9 se aplican a rejillas de
placa perforada para inclinaciones de 0, 10
y 20, donde p es el coeficiente de porosidad
de la rejilla, que se calcula en el caso de una
placa perforada como la relacin entre del
rea de paso del flujo a travs de la placa y
su rea total .
Figura 7. Longitud de rejilla para placa perforada con =0.
Figura 8. Longitud de rejilla para placa perforada con =10.
Figura 9. Longitud de rejilla para placa perforada con =20.