CLASE 4
Método de Muskingum
vwL/K velocidad de la onda
Si vw>Dt/(2 x) subdividir el tramo en n partes
2Kx Dt K
No puede aplicarse a lechos de gran pendiente pues no
se consideran efectos de energía
K, x características del tramo
Tiempo de viaje de la onda K tiempo de subida de
I(t) para que x=cte
OTROS MÉTODOS
K=L/vw x ½ v/vw x = 0,2 - 0,3
Vw/v 1,67 1,44 1,33
Cauce rectangular parabólico triangular
ancho ancho
K tiempo entre C.G. de los hidrogramas de entrada
y salida del tramo tiempo entre Q máximos
1 dQ
= vw
B dy
CANAL LINEAL
Canal ficticio en el cual el tiempo T requerido para
trasladar un caudal Q de cualquier magnitud, en una
distancia x, es constante. No hay cambio de forma del
hidrograma: I = f(t) Q= f(t-T)
Además, en cada sección, v=cte para todo Q
v sección 1 v sección 2
relación entre área mojada y Q es lineal A = cQ
f(T)
T: tiempo de traslación, constante para sección dada.
Si entrada de duración Dt y volumen S se propaga
Q = S d(t,Dt)
1/Dt si 0tDt y t=t+T
0 resto t
d(t) es la función delta de Dirac cuando Dt>0 = HUI para
canal lineal= función impulso.
Un área de drenaje puede ser considerada análoga a
un canal lineal con Q espacialmente variado
Desfase (Canal) y Traslación (Embalse Lineal)
I(t) I(t-T)
Q(t-T)
S(t)=KQ(t-T)
dQ(t - T )
Q(t ) k = I (t )
dt
Si entrada es función impulso, se obtiene el
hidrograma unitario instantáneo del tramo
t -T
1 -
e K
tT
K
u(0,t)= u(t)
0 t<T
1/k
0 T
El HU para una lluvia de duración D resulta:
D t -T
1 -
(e - 1)e
K K
tT
K
h(D,t)=
0 t<T
Dooge (1973)
Dx 4 F 2 h0 Dx 2 F2 h
K= (1 - )( ) T= (1 - (1 - )( 0 ))
u0 9 4 S0 1,5u0 3 4 S0 Dx
MUSKINGUM DE 3 PARÁMETROS (K, x, a)
Para considerar aportes laterales que puedan existir en
el tramo O´Donnell (1985) plantea que éstos son
proporcionales a la entrada al tramo
S=K[x(1+a)I + (1-x)Q]
La ecuación de rastreo:
Q(t+Dt)= d1 I(t) + d2 I(t+Dt) + d3 Q(t)
d1, d2, d3 pueden determinarse por mínimos cuadrados
o gráficamente (= Musk2)
d1 d 2 d 3 1 d 2 d1d 3
K= Dt x = (1 - )
(1 - d 3 )( d1 d 2 ) 2 d1 d 2 d 3
d1 d 2 d 3 - 1
a=
1 - d3
Modelo de Clark Histograma Tiempo-área
(usando curvas isocronas)
HEC usa Ac=atn para 0t0,5tc
Ac=1-a(1-t)n para 0,5tct1tc
Ac es el área contribuyente acumulada dividida por
el área total y t es el tiempo en unidades del
tiempo de concentración, mientras n es un
parámetro que varía según forma de cuenca (=1 si
rectangular; =2 si rombo; =1,5 elíptica).
Clark combina HTA con embalse lineal
Km2
2,78 A Permite transformar
q= áreas en caudales (m3/s)
Dt
hrs
Ejemplo En la cuenca vertiente al embalse de Alhama de
Granada, de 54,3 km2, se han trazado las líneas
isocronas cada media hora, obteniéndose la relación
área-tiempo de la Tabla. Calcular el hidrograma unitario
de t=0,5 hrs según Clark.
Método de Muskingum-Cunge
Se han propuesto muchas variaciones al
método de la onda cinemática; una de éstas
corresponde al método propuesto por Cunge en 1969
que se basa en el método de Muskingum. La ecuación
de dicho método puede escribirse como:
j 1 j+1 j j
Qi+1 = C1 Qi + C 2 Qi + C3 Qi+1 (*)
donde los parámetros Ci se han definido anteriormente en función de los
parámetros del método de Muskingum K y X. Cunge mostró que al
considerar K y Dt constantes, entonces la ecuación * es una solución
aproximada de las ecuaciones de la onda cinemática También mostró
que la ecuación * puede considerarse una solución aproximada de la
ecuación de difusión modificada si:
DxDx
K= = **
ck dQ/dA
1 Q
X = (1 - )
2 B ck S 0 Dx
donde ck es la celeridad correspondiente al caudal Q y a
un ancho superficial B. El lado derecho de la expresión
** representa el tiempo que demora en propagarse en el
tramo Δx una descarga dada. Además, Cunge encontró
que para que existiera estabilidad numérica debía
cumplirse que 0X1/2.
Este método ofrece 2 ventajas sobre los
métodos de la onda cinemática estándar. En primer
lugar se usa una expresión algebraica para obtener la
solución, en vez de una aproximación de diferencias
finitas, de manera que se puede obtener el hidrograma
completo en un punto dado sin recurrir a resolver las
ecuaciones en cada punto del tramo como en el caso
de la onda cinemática. En segundo lugar, la solución de
la ecuación para x tenderá a mostrar menos
atenuaciones de las ondas de manera que permite una
mayor flexibilidad en la elección de los incrementos
temporales y espaciales que el método de la onda
cinemática.
Igualmente Chow et al. (1988) indican que este
método es preferible al método de la onda difusiva
debido a su simplicidad y por cuanto su precisión es
similar. Las desventajas que se indican es que este
método no permite manejar las perturbaciones de
aguas abajo que se propagan hacia aguas arriba y
que tampoco predice adecuadamente el hidrograma
de descarga en un límite aguas abajo cuando hay
grandes variaciones de la velocidad de la onda
cinemática, como ocurre en el caso de inundación
de planicies durante crecidas.
Coeficiente de rugosidad de Manning según uso
de suelo
Escurrimiento en embalses
