JRM
JRM
JRM
DE CARTAGENA
ESCUELA UNIVERSITARIA DE
INGENIERÍA TÉCNICA CIVIL
CARACTERIZACIÓN DE RESALTOS
RESALTOS
HIDRÁULICOS LIBRES Y SUMERGIDOS A
PARTIR DE MEDICIONES DE VELOCIDADES
VELOCIDADES
INSTANTÁNEAS CON EQUIPO
EQUIPO DOPPLER AGUAS
ABAJO DE UNA COMPUERTA. ANÁLISIS
TEMPORAL Y FRECUENCIAL.
Dirigido por:
D. LUIS G. CASTILLO ELSITDIÉ
Doctor Ingeniero de Caminos, C. y P.
D. FCO. JAVIER PÉREZ DE LA CRUZ
Ingeniero de Caminos, C. y P.
2
INDICE
3
5 CONCLUSIONES ................................................................................................... 151
ANEXO I. RESULTADOS DEL FILTRADO A PARTIR DE LAS CAMPAÑAS
DE MUESTREO ..................................................................................................... 154
ANEXO II. GRÁFICAS DE VELOCIDAD FRENTE CALADO ......................... 184
ANEXO III. CARACTERISTICAS ASOCIADAS A LOS RESALTOS
HIDRÁULICOS ANALIZADOS .......................................................................... 197
ANEXO IV. DESCOMPOSICIONES FRECUENCIALES POR MEDIO DEL
ALGORITMO NUMÉRICO DE LA FFT ........................................................... 202
4
1 INTRODUCCIÓN Y OBJETIVOS
Por último se extraen las conclusiones más importantes y se planifican y realizan las
campañas de medición.
Realizamos un estudio del comportamiento del agua en una canal rectangular en el que
se han formado diferentes resaltos hidráulicos libres y sumergidos.
A continuación se muestra una descripción del canal y de los elementos utilizados para
la adquisición de datos con un equipo ADV (Acoustic Doppler Velocimeter).
5
1.2.1 Descripción de los elementos del canal
La fotografía 1.1 muestra los principales elementos que componen el canal sobre el que
realizamos el estudio:
6
Figura 1.3. Depósito regulador del calado del agua con compuerta reguladora del calado del agua
en el canal y depósito de almacenamiento (parte izquierda de la fotografía 1.1).
7
Figura 1.5. Caudalímetro
Donde:
1. Tubería de aspiración.
2. Bomba de alimentación.
4. Tubería de impulsión.
5. Interruptor.
6. Depósito regulador del calado del agua con compuerta reguladora del calado del
agua en el canal.
7. Depósito de almacenamiento.
12. Caudalímetro.
8
1.2.2 Descripción de los elementos para la adquisición de datos con ADV
9
Descripción de cada uno de ellos:
2. Limnímetro.
3. Compuerta.
8. Cinta métrica
9. Medidor de nivel
10
La fotografía 1.7 a muestra los principales elementos de la sonda 3D de 16 MHz de la
marca SONTEK. Se ha de comentar que el resto de sondas utilizadas comentadas
anteriormente tienen los elementos parecidos:
Transmisor acústico
La figura 1.8 a y 1.8 b muestra como captan las mediciones las sondas que toman
medidas hacia abajo y medidas hacia arriba:
11
Figura 1.8 b. Sonda tomando medidas hacia arriba
12
1.3 Objetivos :
Adquisición de datos por medio del equipo Doppler, validación de datos, filtrado
numérico y la comprobación, contrastación teórico-experimental de las distribuciones
de velocidad e índice turbulentos.
13
2 RESUMEN
Una vez hecho el análisis crítico, se extraen las conclusiones más importantes y se
planifican y realizan las campañas de medición.
Los puntos principales que se han tratado en la realización de este proyecto son los
siguientes:
Configuración de la sonda
14
Q caudalímetro Q real (m3/h) Q real (m3 /s)
(m3 /h)
4,50 5,3790 0,0015
5,50 6,5150 0,0018
7,00 8,2190 0,0023
8,50 9,9230 0,0028
9,50 11,0590 0,0031
10,50 12,1950 0,0034
11,00 12,7630 0,0035
12,00 13,8990 0,0039
15
Criterio establecido para obtener la longitud de resalto
L j / y1 = 9,75 ( F1 − 1) 1, 01
(2.1)
16
Figura 2.3 b. Esquema de la situación de las secciones en un resalto hidráulico sumergido.
Filtrado de datos.
1. Detectar las series que presentan puntos anómalos (también llamados “spikes”)
para tratarlas y descartarlos.
17
Relación entre profundidades inicial y final.
Hemos utilizado para el presente proyecto fin de carrera la ecuación de Bélanguer para
los calados conjugados, ya que nuestros resaltos pertenecen al tipo de flujo “no
desarrollado”:
y2
= ( 8 F12 + 1 − 1) / 2
y1 (2.2)
= +
(2.3)
Lj puede ser interpretado como la longitud de la zona necesaria para que la disipación de
energía en el resalto se complete, y la pérdida de energía HL entre el inicio (x=0) y el
final (x=Lj) del resalto puede ser expresada mediante la siguiente ecuación
unidimensional:
H L = H 1 − H 2 = (v1 / 2 g + y1 ) − (v 2 / 2 g + y 2 )
2 2
(2.4)
Siendo:
1 Es la energía cinética al principio del resalto, donde v1 es la
18
H2 = Energía al final del resalto sumergido en función de:
1 Es la energía cinética al final del resalto, donde v2 es la
velocidad al final del resalto libre y g la aceleración de la gravedad 9.81 m/s2
y2 1
2(1 − ) + (1 − ) F12
HL y1 ( y 2 / y1 ) 2
= (2.5)
H1 2 + F12
0,70
Ley teórica de
0,60 HL/H1(con y1/y2
teórico)
0,50 Compuerta (2009)
0,40
HL/H1
Vicente (2008)
0,30
0,20
Márquez (2006)
0,10
0,00
1,0 2,0 3,0 4,0 5,0 6,0 7,0
F1
19
resultado de la comparación entre los valores teóricos y los experimentales es bueno,
pero exceptuando que para algunos valores, los resultados experimentales se desvían de
la ley teórica.
H 2 (8F12 + 1) 3 / 2 − 4F12 + 1
= (2.6)
H1 8F12 (2 + F12 )
H2
% de disipación de energía= (1 − ) × 100 (2.7)
H1
H 1 = (v1 / 2 g + y1 ) y H 2 = (v 2 / 2 g + y 2 )
2 2
(2.8)
En la figura 2.5 se presentan los valores medidos de disipación de energía, junto con los
de Vicente (2008) y Márquez (2006). Como podemos observar, los resultados
experimentales, salvo alguna excepción, son mayores que los teóricos.
20
60
50
% de disipación de energía
40 Ley teorica
30 Compuerta
(2009)
Márquez
20 (2006)
Vicente (2008)
10
0
1 2 3 4 5 6
F1
v2
% reducción de velocidad máxima experimental= (1 − ) × 100 (2.9)
v1
Los resultados obtenidos se representan junto con los de Márquez (2006) y Vicente
(2008) en la figura 2.6. y se puede observar que los resultados experimentales se
agrupan de una forma satisfactoria, y mantienen una buena relación, obteniéndose el
siguiente ajuste:
21
100,00
90,00
80,00 Ley teórica de
% Disipación de energía
disipación de
70,00 energía 2009
Compuerta (2009)
60,00
50,00 Márquez (2006)
40,00
Vicente (2008)
30,00
20,00
10,00
0,00
2,00 3,00 4,00 5,00 6,00 7,00
F1
Resaltos sumergidos.
567 se interpreta como la longitud de la zona requerida para disipar el total de la energía
perdida en el resalto sumergido entre el inicio (x=0) y el final (x=567 ) del resalto
expresada mediante la siguiente ecuación unidimensional:
89 = + :; − + :<
(2.11)
Donde:
Energía cinética al principio del resalto, siendo v0 la velocidad en
Energía cinetica al final del resalto, donde v4 es la velocidad al
22
y4 Energía potencial al final del resalto.
, ? 1 A, ⁄ 1 B
=
=>
= , ⁄ 1?
(2.12)
2.6.
ecuaciones 2.12 y 2.3 se muestra en la figura 2.7 a con líneas continuas para cada valor
de considerado.
1.5<y3/y0<2.5
2.5<y3/y0<3.5
3.5<y3/y0<4.5
0,50
4.5<y3/y0<6
6<y3/y0<8
2,4 2,9 3,4 3,9 4,4 4,9 Ley teórica para y3/y0=8
F0
Figura 2.7a. Pérdida relativa de energía en resaltos sumergidos 85 ⁄80 = E A00 , :3B
:
0
23
Se puede observar que los datos medidos nos salen por encima de las leyes teóricas, es
decir, nos dan una mayor pérdida de energía, Pero esto es solo aparente, puesto que los
resaltos sumergidos se alejan de la configuración teórica.
Para ver la diferencia que hay entre la configuración de flujo obtenida y la teórica,
ajustamos el valor del calado conjugado y3, despejando de la ecuación 2.12, para
ponerlo en su posición correcta en la gráfica.
1.5<y3/y0<2.5
2.5<y3/y0<3.5
3.5<y3/y0<4.5
0,50 4.5<y3/y0<6
6<y3/y0<8
0,40
HL/H0
2,4 2,9 3,4 3,9 4,4 4,9 Ley teórica para y3/y0=8
F0
24
Los resultados experimentales obtenidos los hemos representado junto con los de
menores, para los que no se han podido conseguir grandes valores de sumergencia,
considerando la limitación de la dimensión del canal de nuestro laboratorio.
0,6 F0=2.45
F0=2.93
0,5 F0=3.33
F0=4.307
F0=4.318
0,4
F0=2,87 Márquez 2006
HL/H0
0,2
F0=4,57 Márquez 2006
y3/y0
Figura 2.7 c. Pérdida de energía relativa para resaltos sumergidos: GH ⁄GI = J KI ,
LM
LI
Como hemos comentado para los resaltos libres, tenemos otra forma de evaluar la
disipación de energía en los resaltos, mediante la ecuación 2.6.
Para calcular la disipación de energía en los resaltos sumergidos hemos utilizado las
siguientes expresiones:
25
Siendo:
En la figura 2.8 se han representado los datos obtenidos en el laboratorio junto con los
de Márquez (2006) y la ley teórica. Podemos observar que los datos experimentales
tienen mayor disipación de energía que la obtenida con la ley teórica.
También, si comparamos ésta gráfica con la de disipación de energía para los resaltos
libres (ver figura 2.5) se comprueba que, para los resaltos sumergidos tenemos mayor
disipación de energía que para los valores de los resaltos libres.
60,00
50,00
% Disipación de energía
40,00 Disipación de
energía para
velocidades medias
Compuerta (2009)
30,00
Márquez (2006)
20,00
10,00
0,00
1,00 2,00 3,00 4,00 5,00 6,00
F0
26
Comparando las figuras 2.6 y 2.9, veremos que la reducción de energía por velocidades
máximas es similar para un mismo número de Froude, tanto en el caso de los resaltos
libres como para los resaltos sumergidos, ya que las velocidades máximas medidas en el
perfil de la sección 4 son similares en ambos casos. Este hecho es contradictorio con la
teoría, ya que la reducción de las velocidades máximas debería de ser mayor en los
resaltos sumergidos, ya que el calado será superior (siendo superior, a su vez, el área de
la sección) y a mayor área e igualdad de caudal, la velocidad obtenida será inferior.
90,00
80,00
70,00
% Disipación de energía
60,00
50,00
20,00
10,00
0,00
2,00 2,50 3,00 3,50 4,00 4,50 5,00 5,50
F0
Resaltos libres.
En la figura 2.10 de acuerdo con los datos experimentales obtenidos por nuestras
mediciones junto con los de Márquez (2006) y Vicente (2008), se aprecia una clara
correlación, donde también se representa la ley de Ohtsu et al. (1990) que sigue la
ecuación 2.17:
27
Lj HL
log10 = −1,71 + 1,58 ; en (2,3≤ F1 ≤9,5) y (0,14≤ HL/H1 ≤0,71) (2.17)
HL H1
RhL C
100
Vicente (2008)
Márquez (2006)
Lj/HL
Ajuste propuesto
Lj/HL Márquez
(2006)
Ajuste propuesto
Lj/HL Vicente
(2008)
Ajuste propuesto
1 Lj/HL Compuerta
(2009)
0,0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1,0
HL/H1
Comparando la ley teórica propuesta por Ohtsu et al. (1990) con los resultados
experimentales se observa que los datos experimentales para los resaltos siguen la
misma tendencia y pendiente que las teóricas, lo que nos indica que el criterio de
longitud Lj cumple el criterio establecido.
Márquez (2006):
Lj HL
log10 = −1,71 + 1,49 ; para (2,87≤ F1 ≤4,96) (2.18)
HL H1
Vicente (2008):
Lj HL
log10 = −1,71 + 1,53 ; para (1,96≤ F1 ≤5,45) (2.19)
HL H1
28
Ajuste propuesto por compuerta (2009).
Nuestro ajuste es igual al obtenido por Ohtsu et al. (1990). Esto puede ser debido a que
hemos empleado la ecuación de Silvester (1964) como criterio a la hora de establecer la
longitud del resalto y este es un valor teórico.
Analizando los datos experimentales, se aprecia una buena correlación que se muestra
en la figura 2.11, con la siguiente ecuación:
!OD =VR = −1,71 => + 1,58; para ,2,3 ≤ 0D ≤ 101 : ,1 ≤ :; ⁄:D ≤ 201 (2.21)
9 =
>
100,00
Compuerta (2009)
Márquez (2006)
Lsj/HL
Ajuste propuesto
Lsj/HL Márquez
(2006)
Ajuste propuesto
1,00 Lsj/HL
Compuerta (2009)
0,00 0,20 0,40 0,60 0,80 1,00
HL/H0
Analizando las leyes teóricas junto con los resultados experimentales obtenidos por
Márquez (2006) y en nuestro proyecto se observa que la nube de puntos experimentales
sigue la misma pendiente que las leyes teóricas verificándose que el criterio de longitud
del resalto, Lsj seguido por este proyecto es el correcto. A continuación se muestra la ley
29
de ajuste para los valores de Márquez (2006), la cual es similar a la obtenida por
nuestros datos experimentales, Compuerta (2009).
Márquez (2006)
Compuerta (2009)
30
u 1 y y
= { ( )}1 / n , 0 ≤ ≤ k (2.24)
um k Y Y
u 1 1,177 y y
= exp[− { ( − k )}2 ], k ≤ ≤ 1,5 (2.25)
um 2 1− k Y Y
donde k = y máx / Y
Distribución de Velocidad K n
Ya que en los resaltos libres no se han podido medir muchos puntos cercanos a la
superficie debido a que no había calado suficiente, añadimos unos puntos teóricos
aplicando la formulación de Hagen & Vischer (1995) de velocidades máximas negativas
para los puntos ubicados en la superficie del resalto por medio de las ecuaciones dadas:
ecuación 2.26 y la ecuación 2.27:
N = \6 ] (2.26)
\6 = −N! D.;`
^?D,
(2.27)
31
Siendo:
X = la posición dentro del resalto
V2 =la velocidad al final del resalto.
32
ohtsu a
ohtsu b
7,5 Ley propuesta UPCT
X/Lj=0.5 (R.L.0) A
7,0 X/Lj=0.75 (R.L.0)A
X/Lj=0.5 (R.L.1) A
X/Lj=0.75 (R.L.1)A
6,5 X/Lj=0.5 (R.L.2)A
X/Lj=0.75 (R.L.2)A
X/Lj=0.25( R.S.1) (4.37<y4/y0<13.10) A
6,0 X/Lj=0.5 (R.S.1) (4.47<y4/y0<13.10) A
X/Lj=0.75(R.S.1) (4.37<y4y0<13.10) A
X/Lj=0.25 (R.S.2)( 5.90<y4/y0<13.195) A
5,5 X/Lj=0.5 (R.S.2) (5.90<y5/y0<13.195) A
X/Lj=0.75 (R.S.2) ( 5.90< y5/y0<13.195) A
5,0 X/Lj=0.25 (R.S.3)(y4/y0= 6.06) A
X/Lj=0.5 (R.S.3)(y4/y0= 6.06) A
X/Lj=0.75 (R.S.3) (y4/y0= 6.06) A
4,5 X/Lj=0.5 (R.L.0) C
X/Lj=0.25 (R.L.0)C
X/Lj=0.75 (R.L.0)C
4,0 X/Lj=0.25 (R.L.1) C
y/Y
X/Lj=0.5 (R.L.1)C
X/Lj=0.75 (R.L.1) C
3,5 X/Lj=0.5 (R.L.2)C
X/Lj=0.75(R.L.2)C
X/Lj=0.5 (R.L.3)C
3,0 X/Lj=0.75(R.L.3)C
X/Lj=0.25 (R.S.1)( 2.62<y4/y0< 4.48)
2,5 X/Lj=0.5(R.S.1) (2.62<y4/y0< 4.48) C
X/Lj=0.75 (R.S.1) ( 2.62<y4/y0<4.48) C
X/Lj=0.25 (R.S.2)( 2.76<y4/y0< 4.92 C
2,0 X/Lj=0.5 (R.S.2) (2.76< y4/y0<4.92) C
X/Lj=0.75 (R.S.2) (2.76<y4/y0< 4.92) C
X/Lj=0.25(R.S.3) (2.92<y4/y0<5.52) C
1,5 X/Lj=0.5(R.S.3) (2.92<y4/y0<5.52) C
X/Lj=0.75 (R.S.3) (2.92<y4/y0<5.52) C
X/Lj=0.25 (R.S.4) (3.05<y4/y0<6.2) C
1,0 X/Lj=0.5 (R.S.4) (3.05<y4/y0<6.2) C
X/Lj=0.75 (R.S.4)( 3.05<y4/y0<6.2) C
0,5 X/Lj=0.25(R.S.5) (3.02<y4/y0<4.12) C
X/Lj=0.5 (R.S.5) (3.02<y4/y0< 4.12) C
X/Lj=0.75(R.S.5) ( 3.02<y4/y0<4.12) C
0,0 X/Lj=0.25 (R.S.6) (3.45<y4/y0<5.06) C
X/Lj=0.5(R.S.6) (3.45<y4/y0<5.06) C
-0,8 -0,7 -0,6 -0,5 -0,4 -0,3 -0,2 -0,1 0,0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1,0 1,1 1,2 X/Lj=0.75 (R.S.6) (3.45<y4/y0<5.06) C
X/Lj=0.25(R.S.7) (3.88<y4/y0<6) C
X/Lj=0.5(R.S.7) (3.88<y4/y0<6) C
u/um X/Lj=0.75 (R.S.7) (3.88<y4/y0<6) C
Hagen & Vischer
Figura 2.13. Distribución de velocidades inicial para los resaltos hidráulicos libres y sumergidos en compuerta y aliviadero.
33
Como se puede observar en la figura 2.13 hay registros que se desvían de la
distribución. Estos datos han sido analizados y no se dan por válidos ya que se deben a
malas mediciones (ver figura 2.14). Por lo complejo de este fenómeno, al representar la
ley para las velocidades negativas los puntos siguen una tendencia, pero se encuentran
dispersos.
Hemos realizado un ajuste para resalto hidráulico libre para compuerta mediante la
ecuación 2.28, la cual se indica en la figura 2.15.
[
(2.28)
34
ohtsu a
ohtsu b
Ley propuesta UPCT
7,5 X/Lj=0.5 (R.L.0) A
X/Lj=0.75 (R.L.0)A
7,0 X/Lj=0.5 (R.L.1) A
X/Lj=0.75 (R.L.1)A
X/Lj=0.5 (R.L.2)A
6,5 X/Lj=0.75 (R.L.2)A
X/Lj=0.25( R.S.1) (4.37<y4/y0<13.10) A
X/Lj=0.5 (R.S.1) (4.47<y4/y0<13.10) A
6,0 X/Lj=0.75(R.S.1) (4.37<y4y0<13.10) A
X/Lj=0.25 (R.S.2)( 5.90<y4/y0<13.195) A
X/Lj=0.5 (R.S.2) (5.90<y5/y0<13.195) A
5,5 X/Lj=0.75 (R.S.2) ( 5.90< y5/y0<13.195) A
X/Lj=0.25 (R.S.3)(y4/y0= 6.06) A
X/Lj=0.5 (R.S.3)(y4/y0= 6.06) A
5,0 X/Lj=0.75 (R.S.3) (y4/y0= 6.06) A
X/Lj=0.5 (R.L.0) C
X/Lj=0.25 (R.L.0)C
4,5 X/Lj=0.75 (R.L.0)C
X/Lj=0.25 (R.L.1) C
X/Lj=0.5 (R.L.1)C
4,0 X/Lj=0.75 (R.L.1) C
y/Y
X/Lj=0.5 (R.L.2)C
3,5 X/Lj=0.75(R.L.2)C
X/Lj=0.5 (R.L.3)C
X/Lj=0.75(R.L.3)C
3,0 X/Lj=0.25 (R.S.1)( 2.62<y4/y0< 4.48)
X/Lj=0.5(R.S.1) (2.62<y4/y0< 4.48) C
X/Lj=0.75 (R.S.1) ( 2.62<y4/y0<4.48) C
2,5 X/Lj=0.25 (R.S.2)( 2.76<y4/y0< 4.92 C
X/Lj=0.5 (R.S.2) (2.76< y4/y0<4.92) C
X/Lj=0.75 (R.S.2) (2.76<y4/y0< 4.92) C
2,0 X/Lj=0.25(R.S.3) (2.92<y4/y0<5.52) C
X/Lj=0.5(R.S.3) (2.92<y4/y0<5.52) C
X/Lj=0.75 (R.S.3) (2.92<y4/y0<5.52) C
1,5 X/Lj=0.25 (R.S.4) (3.05<y4/y0<6.2) C
X/Lj=0.5 (R.S.4) (3.05<y4/y0<6.2) C
X/Lj=0.75 (R.S.4)( 3.05<y4/y0<6.2) C
1,0 X/Lj=0.25(R.S.5) (3.02<y4/y0<4.12) C
X/Lj=0.5 (R.S.5) (3.02<y4/y0< 4.12) C
X/Lj=0.75(R.S.5) ( 3.02<y4/y0<4.12) C
0,5 X/Lj=0.25 (R.S.6) (3.45<y4/y0<5.06) C
X/Lj=0.5(R.S.6) (3.45<y4/y0<5.06) C
0,0 X/Lj=0.75 (R.S.6) (3.45<y4/y0<5.06) C
X/Lj=0.25(R.S.7) (3.88<y4/y0<6) C
-0,8 -0,7 -0,6 -0,5 -0,4 -0,3 -0,2 -0,1 0,0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1,0 1,1 1,2 X/Lj=0.5(R.S.7) (3.88<y4/y0<6) C
X/Lj=0.75 (R.S.7) (3.88<y4/y0<6) C
u/um Hagen & Vischer
Figura 2.14. Distribución de velocidades analizada para los resaltos hidráulicos libres y sumergidos en compuerta y aliviadero excluyendo los puntos erroneos.
35
6,0
5,5
X/Lj=0.5 (R.L.0)
:
= 1,5 −0,69 # + 0,024 −7,5 #
X/Lj=0.75 (R.L.0)
4,5
f
X/Lj=0.25 (R.L.1)
4,0
X/Lj=0.5 (R.L.1)
X/Lj=0.5 (R.L.2)
3,0
y/Y
X/Lj=0.75(R.L.2)
2,5
X/Lj=0.5 (R.L.3)
2,0 X/Lj=0.75(R.L.3)
1,0 ohtsu b
0,0
-0,8 -0,7 -0,6 -0,5 -0,4 -0,3 -0,2 -0,1
0,1 0,0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1,0 1,1
u/um
36
Resaltos hidráulicos sumergidos.
Los datos experimentales se explican por la ecuación de resalto libre si y4/y0 → y2/ y0.
YZ
(2.30)
donde:
h = [ , siendo Y la distancia perpendicular a la solera, en la que =
YZ
Esta función la hemos acotado para x=0,68η para su límite inferior y 0 para su límite
superior.
37
En relación a los datos utilizados tenemos algunos en los que la relación y4/y0 es grande
para aplicar la ecuación anterior; y otros en los que la relación anterior se asemeja a
y2/y0, pudiendo aplicar la ley del resaltos sumergidos junto a resaltos libres.
38
X/Lj=0.25( R.S.1) (4.37<y4/y0<13.10) A
7,5 X/Lj=0.5 (R.S.1) (4.47<y4/y0<13.10) A
X/Lj=0.75(R.S.1) (4.37<y4y0<13.10) A
7,0 X/Lj=0.25 (R.S.2)( 5.90<y4/y0<13.195) A
X/Lj=0.5 (R.S.2) (5.90<y5/y0<13.195) A
6,5
X/Lj=0.75 (R.S.2) ( 5.90< y5/y0<13.195) A
1
= 1,48 h 7 [1 − erf,0.68 h1]
#
X/Lj=0.25 (R.S.3)(y4/y0= 6.06) A
6,0
X/Lj=0.5 (R.S.3)(y4/y0= 6.06) A
X/Lj=0.75 (R.S.3) (y4/y0= 6.06) A
5,5
X/Lj=0.25 (R.S.1)( 2.62<y4/y0< 4.48)
X/Lj=0.5(R.S.1) (2.62<y4/y0< 4.48) C
5,0
X/Lj=0.75 (R.S.1) ( 2.62<y4/y0<4.48) C
X/Lj=0.25 (R.S.2)( 2.76<y4/y0< 4.92 C
4,5
X/Lj=0.5 (R.S.2) (2.76< y4/y0<4.92) C
X/Lj=0.25(R.S.3) (2.92<y4/y0<5.52) C
39
Distribución de velocidad en dirección negativa (hacia atrás).
= 1,2 a,tcd + 1
c
[
(2.32)
40
7,5
X/Lj=0.25 (R.S.1)( 2.62<y4/y0< 4.48)
:
X/Lj=0.75 (R.S.1) ( 2.62<y4/y0<4.48)
= 1,2 −6.8# + 1
6,5
f
X/Lj=0.25 (R.S.2)( 2.76<y4/y0< 4.92 )
X/Lj=0.5(R.S.3) (2.92<y4/y0<5.52)
5,0
X/Lj=0.75 (R.S.3) (2.92<y4/y0<5.52)
X/Lj=0.5(R.S.6) (3.45<y4/y0<5.06)
2,0 X/Lj=0.75 (R.S.6) (3.45<y4/y0<5.06)
X/Lj=0.25(R.S.7) (3.88<y4/y0<6)
1,5
X/Lj=0.5(R.S.7) (3.88<y4/y0<6)
1,0 X/Lj=0.75 (R.S.7) (3.88<y4/y0<6)
ohtsu b
0,5
Ley propuesta UPCT
u/um
Figura 2.18. Distribució
Distribución de velocidades para resalto hidraulico
draulico sumergido en compuerta.
41
Los ajustes obtenidos para compuerta se muestran a continuación, en las tablas 2.3 y
2.4, diferenciando entre ajustes para resaltos libres (ecuación 2.33) y para resaltos
sumergidos (ecuación 2.34).
Resaltos libres:
Distribución de
velocidades negativas a b c d
Resaltos sumergidos:
Distribución de
velocidades negativas f g h
Cuando tomamos medidas con los aparatos de alta resolución de medición “acoustic
doppler velocimeters (ADV)” en el seno de un fluido, se observa que, sufren un ruido
que distorsiona la turbulencia.
El concepto de Kolmogoroff (1941) afirma que para una señal turbulenta con poca
contaminación por el ruido blanco sigue una pendiente de -5/3 en un rango de altas
frecuencias, 1Hz<f <10Hz y se encuentra horizontal a bajas frecuencias, f <1Hz. Tal y
como se muestra en la siguiente figura:
42
Figura 2.19. Señal turbulenta sin contaminación acústica que cumple el concepto de Kolmogoroff
(1941)
En este estudio veremos cómo afecta este ruido a las mediciones tomadas,
comparándolas unos puntos con otros, tomadas a distintas frecuencias de medición,
5Hz, 25Hz y 50Hz, tomadas en la sección 3 y sección 5, en dos configuraciones de
resalto, es decir, resalto libre y sumergido. Y en la normal “X”, “Y”, y la cortante “XY”,
cuyos espectros de energía se nombraran de la siguiente manera:
Para ello hacemos un análisis frecuencial de las mediciones por medio del algoritmo
numérico de la Transformada Rápida de Fourier (FFT), y las analizamos en gráficas
log-log.
43
Comenzando nuestro análisis vemos que a mayor frecuencia a de medición tenemos
el ruido a mayor descomposición frecuencial, tal y como se puede observar en la
figura 2.20, que en la cual hemos representado la sección 5 del resalto sumergido
medido a distintas frecuencias, en la Sx’.
0,1
0,01
S5RS 5Hz
Sx'
1E-06
0,001 0,01 0,1 1 10
frecuencia Hz
También comprobamos que hay variación de la energía con respecto al eje que se
analiza, tal y como se muestra en la figura. 2.21, (sección 3 del resalto sumergido,
medido a una frecuencia de 5Hz).
44
1
0,1
0,01 Sx'
Sij
Sx'Sy'
0,001
Sy'
0,0001 p -5/3
P = -5/3
0,00001
0,000001
0,001 0,01 0,1 1 10
F(Hz)
Figura 2.21. Descomposición frecuencial de la sección 3 sumergida para xsy y xy la xsy xy , a una
En este caso observamos que la energía de la z{y es superior a la z{y zLy , siendo
frecuencia a 5Hz.
0,1
0,01
Sx'
S3RL 5Hz
0,001 S3RL 25Hz
S3RL 50Hz
0,0001 p=-5/3
P = -5/3
1E-05
1E-06
0,001 0,01 0,1 1 10
frecuencia Hz
Figura 2.22. Espectro de energía en la sección 3 RL medidas a una frecuencia de 5Hz, 25Hz y 50Hz
45
En la figura 2.22 (espectro de energía en la sección 3 de un resalto libre medido a las
frecuencias de 5, 25, 50 Hz, en la Sx’) podemos comprobar que el ruido es bastante
elevado en la sección 3, y que distorsiona la señal turbulenta, pero tienden al
concepto de Kolmogoroff (1941) en el último tramo, (este hecho se puede verificar
sumergido, que la señal tiene bastante ruido, siendo este similar en z{y y zLy la
también en la figura 2.21, en la que se ve en distintos ejes, de la sección 3 del resalto
0,1
S5RL 5Hz
0,01 S5RL 25Hz
S5RL 50Hz
Sx'
0,001 p=-5/3
Distribución de
0,0001 tramos por
pendientes
P = -5/3
1E-05
1E-06
0,001 0,01 0,1 1 10
frecuencia Hz
Figura 2.23. Descomposición frecuencial en la sección 5 medida a distintas frecuencias.
Para una misma sección de resalto libre y sumergido, vemos que la energía es superior
en el resalto libre que en el resalto sumergido en la sección 3 igual pasa con el
ruido, tal y como se muestra en la figura 2.24, en la que se muestra la sección 3 de un
resalto libre y sumergido medida a una frecuencia de 5Hz, en la Sx’
46
1
0,1
0,01
Sx'
S3 LIBRE
0,001
S3 SUMERGIDO
p=-5/3
0,0001
1E-05
P = -5/3
1E-06
0,001 0,01 0,1 1 10
F(Hz)
0,1
0,01
Sx'
0,001
S5 LIBRE
S5 SUMERGIDO
0,0001 p=-5/3
1E-05 P = -5/3
1E-06
0,001 0,01 0,1 1 10
F(Hz)
47
3 LABOR EXPERIMENTAL Y FILTRADO DE DATOS
1. Exactitud.
2. Rapidez.
3. Medida de velocidades en 2D y 3D.
4. Calibración de fábrica invariable (por lo que no se requiere ninguna calibración
periódica).
5. Operación simple.
6. Cálculo directo de parámetros turbulentos, tales como tensión de Reynolds.
7. Rendimiento excelente en flujos lentos.
48
Los principales parámetros de calibración de la sonda son los siguientes:
Rango de velocidad:
Frecuencia de muestreo:
Se establece con el fin de obtener unas series de velocidades más homogéneas y unos
valores de correlación y de SNR mejores. El valor de frecuencia utilizado en la
adquisición de datos es de 5 Hz, y para la realización del análisis frecuencial, utilizamos
también frecuencias de 25Hz y 50Hz.
Salinidad:
Intensidad de señal:
49
Coeficiente de Correlación:
El ADV calcula valores de correlación para cada receptor acústico, con cada muestra de
velocidad. La correlación es expresada en porcentaje siendo la correlación perfecta de
un 100%, indicando la máxima fiabilidad en mediciones de velocidad con bajo ruido; la
correlación 0% indica que el valor de velocidad de salida está dominado por el ruido
(por lo que no hay señal coherente).
Los valores de correlación deben estar entre 70 y 100%. Valores por debajo del 70%
indican que el ADV esté operando en un régimen de medida difícil, es decir, la sonda
está fuera del agua, el SNR es demasiado bajo, o que puede pasarle algo al ADV. En
ciertos entornos (flujos altamente turbulentos), puede no ser posible lograr altos valores
de correlación.
Los valores de correlación bajos afectarán la variabilidad a corto plazo en los datos de
velocidad (por ejemplo, aumento del ruido), pero no predisponen las mediciones de
velocidad media.
Para mediciones de velocidad media, pueden usarse los valores de correlación entre el
30%-100%.
Los caudales que marca el caudalímetro cuando se pone en marcha la bomba no son los
correctos, por tanto, es necesario realizar un ajuste con los caudales reales.
Para ello, se realizó un aforo volumétrico a distintos caudales para obtener un gráfico
Qreal – Qcaudalímetro que despejara todas las dudas respecto al caudal que realmente circula
por el canal. Para ello se procedió al llenado de un recipiente de 22 litros, tomando el
tiempo de la operación, para diversos caudales marcados con el caudalímetro.
50
De esta forma, se calculó, a partir del volumen de agua recogido en un determinado
intervalo de tiempo, el caudal realmente trasegado en cada momento. Para realizar este
aforo se contó con la ayuda de una probeta de 1000 mililitros, tal como se muestra en
las fotografías.
Para cada uno de los caudales marcados por la bomba se realizaron tres mediciones y se
obtuvo el valor promedio de las mismas.
Una vez que tenemos los caudales reales aforados y los que marca el caudalímetro,
podemos conocer la relación que guardan mediante una recta de regresión.
Para realizar el aforo nos hemos basado en los datos obtenidos en los proyectos
anteriores, muy similares a los nuestros, añadiéndole también el caudal 2.5 m3/h, tal y
como hicieron en el año 2008.
En el gráfico también hemos incluido los datos de años anteriores con el fin de realizar
una comparación con los obtenidos este año. Los datos obtenidos por Vicente (2008) y
García (2008) están en color azul, y los obtenidos por Márquez (2006), en color rojo.
Sin embargo, los datos de otros años no se han considerado para obtener el ajuste por
regresión lineal, representado en color verde, que se realizó solo para nuestros datos.
51
Los caudales empleados para realizar el ajuste fueron 2,5, 4, 5, 6, 8, 10, 10,5, 11,5 y 12
m3/h.
16
14
12 Vicente
y = 1,136x + 0,267
(2008), Garcia
R² = 0,997 (2008)
10
Márquez (2006)
8
Q real(m3/h)
6
Corrección de
caudales 2009
4
2 Lineal
(Corrección de
0 caudales 2009)
0 5 10 15
Q caudalimetro(m3/h)
52
3.1.3 Descripción del canal:
El canal tiene una longitud de 5000 mm, y una sección transversal de 81x250mm. y las
partes fundamentales son: dos soportes situados a ambos lados del canal, los cuales
permiten la nivelación del canal, un depósito de almacenamiento de agua, una bomba
centrífuga, una compuerta aguas abajo y unos tranquilizadores de flujo.
• Soportes del canal: Los soportes del canal cumplen con una doble función, la
primera es apoyarlo, y la segunda es permitir su nivelación.
• Depósito de almacenamiento: Contiene el agua que recircula por el canal.
• Bomba de agua: La bomba de agua permite recircular el agua a distintos
caudales gracias a una válvula de regulación, el máximo caudal bombeado para
nuestro proyecto fue de 13,89 m3/h.
• Tranquilizadores de flujo: Se encuentran en ambos extremos del canal, la función de
los tranquilizadores que se encuentran aguas arriba de la compuerta es facilitar que
el flujo llegue tranquilizado a la compuerta y evitar así las oscilaciones que se
puedan producir en ese punto. Y el del aguas abajo tiene la finalidad de que no se
produzcan ondas de retorno que influyan a los resaltos.
• Compuerta (Aguas abajo): Cumple una función reguladora ya que permite
controlar el caudal de aguas abajo.
• Situación de la compuerta: Es el elemento con el cual modificamos el flujo del
caudal, ya que este pasará de supercrítico a subcrítico por medio de las distintas
configuraciones de resaltos hidráulicos. La distancia desde el depósito
tranquilizador de salida del agua hasta la compuerta es de 73 cm. La apertura
máxima utilizada es de 5,6 cm por las limitaciones del canal.
Con ayuda de un medidor de nivel se ha comprobado que la primera mitad del canal
situada en la parte izquierda de la siguiente fotografía está perfectamente nivelada y la
pendiente del canal es nula en ese tramo.
53
Figura 3.4. Situación del canal.
Por lo contrario, la otra mitad del canal, situada en la parte derecha de la fotografía, no
está totalmente nivelada porque la estructura del canal no lo permite al no disponer en
este extremo de un tornillo regulador de altura.
Dado el resultado de la nivelación del canal la parte útil considerada para nuestros
ensayos es la primera mitad situada a la izquierda de la fotografía.
54
Figura 3.5. Vista transversal del canal (Vicente 2008)
Para los resaltos sumergidos, hemos utilizado 5 caudales (1,5, 1,8, 2,8, 3,4, 3,9 l/s),
midiendo para cada caudal 3 sumergencias, exceptuando los caudales más altos (3,4 y
3,9 l/s), en los que solo hemos podido medir 2 sumergencias, ya que era imposible
medir la tercera, por desbordamiento del agua, aguas arriba de la compuerta.
El calado de las 3 sumergencias varía con el caudal, ya que para obtenerlas se toman,
para el resalto libre de referencia, los valores de la carga de agua (Htotal) aguas arriba de
la compuerta y el calado (y5) en una sección intermedia entre la compuerta y la mitad
del canal (ver figura 3.6 a). La diferencia obtenida entre los dos valores nos da el valor
de la tercera sumergencia, obteniéndose las otras dos a 2/3 (segunda sumergencia) y 1/3
(primera sumergencia) de la inicial (ver figura 3.6 b).
55
Figura 3.6 a. Obtención de Htotal - y5 a partir de un resalto libre.
En las siguientes imágenes podemos apreciar la formación de los resaltos, tanto libres,
como sumergidos.
56
2. Tres resaltos sumergidos estabilizados a pie de compuerta, primera
sumergencia, R.S.1, segunda sumergencia, R.S.2 y tercera sumergencia,
R.S.3.
57
• Sección 4: localizada a una distancia igual a la longitud del resalto.
58
3.3 Equipos de medida. Principales especificaciones.
En la toma de datos se emplearon velocímetros acústicos de alta frecuencia, también
llamados ADV (Acoustic Doppler Velocimeter). Estos instrumentos sirven para medir
velocidades en un medio acuoso y su funcionamiento está basado en el efecto Doppler.
Para mantener la sonda estática dentro del flujo de agua, nos ayudamos de un brazo
construido para ese fin.
Las patillas de la sonda deben encontrarse a más de 5 cm de la solera del canal, debido a
que al emitir el sonido, se pueden producir reflexiones en la solera que, si la distancia es
inferior a 5 cm, pueden interferir el sonido reflejado con el emitido, obteniéndose
medidas erróneas.
En condiciones de aireación, cuando una burbuja de aire de los remolinos del resalto,
atraviesa la trayectoria de muestreo, el aparato deja de medir obteniéndose un valor de
velocidad que no es el real.
59
Por esta razón es necesario realizar un filtrado de los datos, con el fin de desechar o
paliar las anomalías debidas a dicho fenómeno de aireación.
Para las mediciones realizadas en este proyecto fin de carrera se han utilizado sondas de
dos marcas diferentes: sondas ADV Sontek y sondas ADV Nortek.
60
3. Sonda 2D de 16 MHz MicroADV, que efectúa mediciones tomando como
referencia la solera del canal (medidas hacia abajo).
Los parámetros y equipamientos de esta sonda son prácticamente iguales a los de los
instrumentos Sontek, diferenciándose únicamente en el software de adquisición de
datos, ya que este funciona con el sistema operativo MS-DOS, y solo permite medir
hasta una frecuencia de 25Hz.
Para tratar los datos obtenidos con esta sonda, es decir, para poder pasarlos a formato
compatible con Excel, lo haremos mediante el programa “WinADV”, utilizando la
opción “process unfiltered data” y marcando la casilla “produce time series data file”,
61
dentro del menú “process”. Esto generará un archivo de salida (extensión “.Vu”) que se
ubicará en la misma carpeta en la que se encuentre el archivo ADV.
62
3.4 Filtrado digital de datos.
Como hemos comentado anteriormente, el equipo de medida de velocidades
instantáneas (ADV, Acoustic Doppler Velocimeter) necesita estar completamente
sumergido y presenta un buen rendimiento en el muestreo de flujos laminares y
turbulentos sin la presencia de aire.
Estos errores tenemos que eliminarlos y/o corregirlos verificando los registros, para lo
cual se realiza una filtración digital de la información.
1. Detectar las series que presentan puntos anómalos (también llamados spikes)
para tratarlas y descartarlos.
63
3.4.1 Filtrado por umbral de aceleración (Castillo, 2009).
En primer lugar escogemos la serie de datos de velocidad muestreados que vayan a ser
nuestro objeto de estudio.
u j F j gy j
aj = = = λaj g (3.1)
∆t ∆t
Donde:
Fj y j
λaj = ≥ 0 .5 (3.2)
∆t g
Para el final del resalto, la aceleración máxima de corte hallada no deberá ser menor a
media gravedad, pues de lo contrario podríamos alterar la serie tomando puntos válidos
por spikes.
Variables utilizadas:
del flujo es hacia aguas abajo; y toma el valor y uj cuando la dirección es hacia
64
El parámetro λ aj , como hemos dicho anteriormente, se establece en función de las
Figura 3.12. Principales parámetros del resalto hidráulico. Variación del número de Froude.
a i = (u i − u i −1 ) / ∆t , (3.3)
Donde
u i = velocidad de la serie registrada en el tiempo i .
∆t = intervalo de muestreo.
65
Siendo:
| a | = El valor absoluto de la aceleración de las velocidades discretas en
el intervalo de tiempo muestreado. La consideramos en valor absoluto
para actuar tanto en las aceleraciones positivas como en las negativas.
Este paso debería repetirse hasta que no se detectaran nuevos spikes, pero en nuestro
trabajo hemos aplicado una sola iteración.
Obtenida la serie filtrada por el umbral de aceleración, aplicamos el descarte del corte
progresivo de los límites inferior y superior en función de 5% y del 95% estadísticos
(Castillo 2008) si es necesario:
La diferencia entre la media Xmedia y A1 nos dará el valor Xmín = Xmedia - A1 que
utilizaremos para el cálculo de la amplitud general A = Xmáx – Xmín.
Así adquiriremos los valores de corte superior Xmáx,c e inferior Xmín,c que
utilizaremos para filtrar la serie finalmente.
66
En definitiva, los valores de velocidad que sean mayores que Xmáx,c y menores que
Xmín,c se considerarán fuera de lo aceptable por nosotros y serán sustituidos por ambos
respectivamente, obteniendo una serie casi limpia. Este método se aplica reiteradamente
hasta que no detectemos más spikes. Aunque se recomienda no hacer más de dos
filtrados para no alterar la serie original, establecemos la posibilidad de la aplicación
hasta cuatro filtrados para casos específicos.
67
Figura 3.14. Resumen de la hoja de filtrado.
y2
= ( 8 F12 + 1 − 1) / 2
y1
(3.5)
68
5. La velocidad teórica en la sección 1.
6. El número de Froude.
7. Número de Reynolds.
y2 calado
Q
Apertura y0 calado conjugado V1 teórica Froude Reynolds
Caudalímetro contraído Bélanguer
m3/h a(cm) cm cm m/s
4,5 2,10 1,23 6,92 1,50 4,32 14107
5,5 2,15 1,40 7,85 1,60 4,31 16553
7 2,90 2,05 7,92 1,37 3,07 18658
8,5 3,35 2,20 9,32 1,55 3,33 21985
9,5 3,90 2,70 9,16 1,40 2,73 22687
10,5 3,90 3,10 9,29 1,35 2,45 23618
11 4,10 2,80 10,49 1,56 2,98 25801
12 4,30 3,00 11,02 1,59 2,93 27300
Tabla 3.2. Características principales de los resaltos libres
La siguiente Tabla 3.3 incluye las principales características medidas para los resaltos
libres:
69
Q Q V1 V2 F1 Y1 Y2 Htotal Lj Reynolds
Caudalimetro Real m/s m/s cm cm cm cm
m3/h m3/s
4,5 0,00149 1,50 0,28 4,32 1,23 6,70 12,60 40,27 14107
5,5 0,00181 1,60 0,32 4,31 1,40 7,05 12,36 45,64 16553
7 0,00228 1,37 0,35 3,07 2,05 7,95 13,14 41,60 18658
8,5 0,00276 1,55 0,39 3,33 2,20 8,65 13,00 50,40 21985
9,5 0,00307 1,40 0,43 2,73 2,70 8,90 14,00 45,78 22687
10,5 0,00339 1,35 0,41 2,45 3,10 10,20 17,40 43,88 23618
11 0,00355 1,56 0,45 2,98 2,80 9,80 15,30 54,50 25801
12 0,00386 1,59 0,43 2,93 3,00 11,00 17,90 56,79 27300
Tabla 3.3. Características principales para los resaltos al pie de la compuerta
Siendo:
Primera sumergencia, h5=1/3(Htotal- y5)+ y5
Segunda sumergencia, h5=2/3(Htotal- y5)+ y5
Tercera sumergencia, h5= (Htotal- y5)+ y5
En la siguiente tabla se presentan los valores del caudalímetro, el caudal real, la y5, la
Htotal, y la h para los resaltos libres, así como el valor del calado de cada sumergencia en
la sección 5, siendo R.S.1 la primera sumergencia, R.S.2 la segunda sumergencia y R.S.3
la tercera sumergencia para los resaltos sumergidos.
70
Q Q Resalto libre resaltos sumergidos
Caudalímetro real R.L.0 R.S.1 R.S.2 R.S.3
Htotal y5 h y5+1/3*h y5+2/3*h y5+h
3 3
m /h m /s cm cm cm cm cm cm
4,5 0,0015 12,60 6,7 5,90 14,57 16,53 18,50
5,5 0,0018 12,36 7,1 5,26 14,11 15,87 17,62
8,5 0,0028 13,00 8,7 4,30 14,43 15,87 17,30
10,5 0,0034 17,40 10,35 7,05 19,75 22,10 24,45
12 0,0039 17,90 11,17 6,73 20,14 22,39 24,63
72
4 FUNDAMENTOS Y ANÁLISIS TEÓRICO Y EXPERIMENTAL.
4.1 Introducción.
El resalto hidráulico es un fenómeno conocido como método útil para disipar el exceso
de energía de flujos de alta velocidad.
Los saltos hidráulicos ocurren cuando hay un conflicto entre los controles que se
encuentran aguas arriba y aguas abajo, los cuales influyen en la misma extensión del
canal.
El salto hidráulico puede tener lugar sobre la superficie libre de un flujo homogéneo o
en una interfase de densidad de un flujo estratificado.
Este proyecto es una continuación del iniciado por Márquez (2006) y seguido por
Vicente (2008) y García (2008), en los cuales se realizaron estudios de resaltos
sumergidos y libres.
4. Dependiendo del desplazamiento del inicio del resalto, así podremos obtener
flujos desarrollados y no desarrollados.
73
profundidades inicial y final. Se comprueban y se proponen nuevas formulaciones para
los resaltos libres y sumergidos sobre la energía disipada y las condiciones para la
formación de los resaltos tras la compuerta, las características de la caída de velocidad y
la distribución de velocidades de dichos resaltos. También hacemos un análisis
frecuencial por medio de la Transformada Rápida de Fourier (FFT) de algunas de las
mediciones tomadas, tanto para resaltos libres como para los resaltos sumergidos
1. Disipar la energía del agua que fluye sobre presas, vertederos y otras estructuras
hidráulicas, con el fin de prevenir la socavación aguas abajo de las estructuras.
2. Recuperar altura o aumentar el nivel del agua en el lado de aguas abajo de una
canaleta de medición y mantener un nivel alto del agua en el canal de irrigación
o de cualquier estructura para distribución de aguas.
74
8. La identificación de condiciones especiales de flujo, como la existencia del flujo
supercrítico o la presencia de una sección de control para la medición de la
relación efectividad-costo del flujo mediante una estación de aforo.
3. Para F = 1,7 a 2,5, se desarrolla una serie de remolinos sobre la superficie del
resalto, pero la superficie del agua hacia aguas abajo permanece uniforme. La
velocidad a través de la sección es razonablemente uniforme y la pérdida de
energía es baja. Se presenta entonces el resalto débil.
4. Para F = 2,5 a 4,5, existe un chorro oscilante que entra desde el fondo del resalto
hasta la superficie y se devuelve sin ninguna periodicidad. Cada oscilación
produce una onda grande con periodo irregular, muy común en canales, que
puede viajar a lo largo de varios kilómetros causando daños a bancales de tierra y
escolleras de protección. Se produce entonces el resalto oscilante.
6. Para F = 9,0 y mayores, el chorro de alta velocidad choca con paquetes de agua
intermitentes que discurren hacia abajo a lo largo de la cara frontal del resalto,
75
generando ondas hacia aguas abajo, y puede prevalecer una superficie rugosa. La
acción del resalto es brusca pero efectiva debido a que la disipación de energía
puede alcanzar un 85%. Se produce entonces el
e resalto fuerte.
76
4.2 Conclusiones de estudios anteriores.
f (L / H L , H L / H* ) = 0 (4.1)
77
De acuerdo con lo dicho en el párrafo anterior, Lj es la longitud del resalto, siendo y1 la
profundidad de agua al inicio del resalto (x=0) e y2 la profundidad al final del resalto
(x= Lj).
En la sección x=Lrj se definió el final del rulo del resalto. También indica el final de la
zona de separación entre Lrj y Lj .
La curva de ajuste que obtuvo con sus datos Vicente (2008) es la siguiente:
L jc / y1 = 9,51( F1 − 1) 0 ,879
(4.2)
Se comprueba que ambas curvas de ajuste son similares a la obtenida por Silvester
(1964):
L j / y1 = 9,75 ( F1 − 1) 1, 01
(4.4)
78
Resaltos hidráulicos sumergidos.
Márquez (2006), tomó la longitud de los resaltos hidráulicos sumergidos tal y como se
ha descrito para los resaltos libres, es decir, el final del resalto sumergido (x=L ) se
sj
determinó como la primera sección donde casi no existe una variación de la velocidad
máxima del flujo. Así, L es la longitud del resalto, siendo y la profundidad del agua
sj 4
Además, la sección donde x=L es la primera en la que la superficie del agua llega a
sj
incrementa.
79
F=2,87
10
F=3,56
F=3,66
F=4,19
F=4,57
F=4,96
5
Lsj/y2
Línea de ajuste
experimental UPCT (R.S)
Línea de ajuste
experimental UPCT (R.L)
[2,5<F<5]
0 Línea teórica (Bureau of
Reclamation) para R.L
1,00 1,25 y4/y2 1,50 1,75 (6<F<11)
Figura 4.4. Características de longitud de los resaltos libres y sumergidos (Márquez 2006).
Según la ecuación 4.6, se puede obtener la ecuación 4.7 para el caso y /y =1:
4 2
De esta manera, Ohtsu hace coincidir en su ajuste este valor para el caso de resalto libre,
aunque el coeficiente 5,9 es válido para un rango de números de Froude de [6≤ F ≤11],
1
siendo este coeficiente más bajo, del orden 5,4 aproximadamente para los números de
Froude [2,5≤ F ≤5].
1
80
La curva de ajuste expresada mediante la ecuación 4.8 constituye el mejor ajuste a los
resultados experimentales obtenidos por Márquez (2006):
A partir de la ecuación 4.8 se puede establecer para el caso del resalto libre (y /y =1).
4 2
Comentando la figura 4.4 se observa que hasta una sumergencia del orden de y /y =1,20
4 2
las longitudes de resalto sumergido, aunque van aumentando, podrían ser inferiores a la
longitud de resalto libre definido por el Bureau of Reclamation [6< F <11]; pero a
1
partir de esta sumergencia, las longitudes de resalto sumergido superan los valores de
definición de L del Bureau of Reclamation y tienden a coincidir con la definición de
sj
81
4.2.2 Relación entre profundidades inicial y final.
∂V
∫ V (n .V )dS = ∫ n pdS − [ ∫ µ dS + ρ ∫ V ' (n0 .V ' )dS ] (4.10)
∂n I
0 I
S S S S
Donde
S: superficie de control
Lj
∂u
y2
donde:
A: área de la sección
82
y : coordenada vertical dirigida hacia arriba empezando desde el fondo del canal
q: caudal unitario ( q = Q / B )
β = ∫ (u / v) 2 dA / A
A
donde:
S f = P f /(γ y12 / 2 )
Lj
∂u
τ 0 = [µ + ρ ( −u ' v ')] y = 0 )
∂y
F1 = v1 / gy 1
y2
γ 2 = ∫ u' 2 dy / v12 y1
0
donde:
83
β1 = 1 + ∆β1
β 2 = 1 + ∆β 2
Harleman (1958) encontró que los errores debidos a la uniformidad asumida del flujo
( β1 = β 2 = 1) y a no tener en cuenta la turbulencia al comienzo y al final del resalto
( γ 1 = γ 2 = 0 ), son muy pequeños y tienden a cancelarse unos a otros, y que la tensión
de corte del fondo integrada Sf es el único término adicional importante a considerar.
Rajaratnam (1965), demostró que los datos obtenidos para diversos experimentos y de
diversa procedencia, mostraban mejor concordancia con la ecuación 4.13 que con la
ecuación 4.5. De acuerdo con esto, en un supuesto Caso “a” referido al inicio de un
resalto libre bajo las condiciones del flujo “no desarrollado”, justo debajo de la
compuerta, el valor experimental de y2/y1 para un F1 dado era algo más pequeño que el
valor de y2/y1 dado por la ecuación 4.5; pero en un supuesto Caso “b” referido al inicio
de un resalto libre bajo las condiciones del flujo totalmente desarrollado, se asumía que
el valor experimental podía coincidir con el valor de y2/y1 derivado de la ecuación 4.5.
Esto se explicaba por la suposición de que en el Caso “a”, los efectos de Sf podrían ser
grandes dado que la capa límite no se separa del fondo del canal, mientras que en el
Caso “b”, Sf ≈0, porque la capa límite podría estar separada del fondo del canal (flujo
totalmente desarrollado).
Ohtsu et al. (1990) considera que las diferencias comentadas en el anterior párrafo
demostradas por Rajaratnan (1965) son pequeñas y sugieren que para propósito de
diseño práctico, se aplique la ecuación de Bélanguer, a los dos resaltos tipo “a” y “b”.
84
Leutheusser y Kartha (1972) llevaron a cabo un estudio semi-empírico de resaltos
hidráulicos rechazados aguas abajo de una compuerta, estos corresponden a flujos
totalmente “desarrollados”.
Vicente (2008), comparando sus resaltos experimentales con respecto a las leyes
teóricas de Rajaratnan, Leutheusser y Kartha, y la ley teórica de los calados conjugados
de Bélanguer, comprobó que los datos se ajustan más a la ley teórica de Rajaratnan,
que a la ley teórica de los calados conjugados de Bélanguer.
Para asegurarse que los resultados que obtenían al pie de la compuerta eran correctos
Vicente (2008) realizó mediciones de calados contraídos aguas abajo de la compuerta
para todos los resaltos al pie, utilizando las gráficas de los coeficientes de contracción
(Rouse, Engineering Hydraulics, 1950) para obtener estos a partir de la apertura de
compuerta y la carga aguas arriba de ésta. Multiplicando estos coeficientes por la
apertura de compuerta se obtiene el calado contraído teórico e incluso la distancia a la
que debe ser medido y1 .
Se observó claramente que sus resultados experimentales caen por debajo de los
teóricos.
85
8 Experimental
7 Valor(1) según
Henry linea
continua
6
Valor(2) según
Henry linea
5 discontinua
Y2/Y1
Valor(3)
Condiciones
4 ideales
Belanguer
3
Rajarantnan
2
1 Leutheusser y
Kartha
1 2 3 4 5
F1
Figura 4.5. Comparación valores teóricos y reales y2 y F1 (Vicente, 2008).
y1
= +
(4.15)
Márquez (2006), demuestra que para y3/y0 en función de y4/y0 hay una buena
concordancia entre los valores teóricos calculados con la ecuación 4.15 y los resultados
experimentales.
86
4.3 Pérdida de energía en los resaltos hidráulicos.
Aplicando la ecuación de la energía del flujo medio para la zona de resalto (0≤x≤Lj), se
obtiene la ecuación 4.16. En este caso, los términos de tensión viscosa media y los
efectos de la tensión de Reynolds en la superficie de control no se tienen en cuenta:
∂u ∂u ∂ v
2 2
v1 v 1
α1 + λ1 y1 = α 2 2 + λ2 y 2 + [− ∫ {ρ (u ' 2 − v' 2 ) + ρ u ' v'( + )}dw] (4.16)
2g 2g γ q ( D ) ∂x ∂y ∂x
Donde:
De acuerdo con esto, la ecuación 4.16 puede ser expresada como la ecuación 4.17:
∂u ∂u ∂v
2 2 Lj y( x)
v1 v 1
+ y1 = 2 + y 2 + [ − ∫ ∫ { ρ (u ' − v' 2 ) + ρ u ' v '( + )}dydx ] (4.17)
2
2g 2g γ q 0 0
∂x ∂y ∂x
Prod.Turb.=Disip.Vis.+Convecc.Turb. (4.18)
donde:
87
Disip.Vis.= Disipación viscosa
y2
1 V '2
Convecc.Turb.=
γ q ∫ρ
0
2
udy
H1=H2+(Disip.Vis.+Convecc.Turb.) (4.19)
Donde:
2
v
H1 = 1 + y1
2g
2
v
H 2 = 2 + y2
2g
Siendo:
1 Es la energía cinética al principio del resalto, donde v1 es la
88
1 Es la energía cinética al final del resalto, donde v2 es la velocidad
Sin embargo, la ecuación 4.19 puede ser rescrita en la forma de la ecuación 4.20:
De acuerdo con esto, Lj puede ser interpretado como la longitud de la zona necesaria
para que la disipación de energía en el resalto se complete, y la pérdida de energía HL
entre el inicio (x=0) y el final (x=Lj) del resalto puede ser expresada mediante la
siguiente ecuación unidimensional:
H L = H 1 − H 2 = (v1 / 2 g + y1 ) − (v 2 / 2 g + y 2 )
2 2
(4.21)
y2 1
2(1 − ) + (1 − 2
) F12
HL y1 ( y 2 / y1 )
= (4.22)
H1 2 + F12
89
0,70
Ley teórica de
0,60 HL/H1(con y1/y2
teórico)
0,50 Compuerta (2009)
0,40
HL/H1
Vicente (2008)
0,30
0,20
Márquez (2006)
0,10
0,00
1,0 2,0 3,0 4,0 5,0 6,0 7,0
F1
H 2 (8F12 + 1) 3 / 2 − 4F12 + 1
= (4.23)
H1 8F12 (2 + F12 )
90
Resultados teóricos
Porcentaje de disipación de
F1 H2/H1 energía en velocidades
teóricas medias
2,50 0,82 17,52
3,00 0,74 25,67
3,50 0,67 32,90
4,00 0,61 39,14
4,50 0,56 44,48
5,00 0,51 49,06
5,50 0,47 53,00
H2
% de disipación de energía= (1 − ) × 100 (4.24)
H1
H 1 = (v1 / 2 g + y1 ) y H 2 = (v 2 / 2 g + y 2 )
2 2
(4.25)
En la figura 4.7 se presentan los valores medidos de disipación de energía, junto con los
de Vicente (2008) y Márquez (2006). Como podemos observar, los resultados
experimentales, salvo alguna excepción, son mayores que los teóricos.
91
60
50
% de disipación de energía
40 Ley teorica
30 Compuerta
(2009)
Márquez
20 (2006)
Vicente (2008)
10
0
1 2 3 4 5 6
F1
v2
% reducción de velocidad máxima experimental= (1 − ) × 100 (4.26)
v1
Los resultados obtenidos se representan junto con los de Márquez (2006) y Vicente
(2008) en la figura 4.8. y se puede observar que los resultados experimentales se
agrupan de una forma satisfactoria, y mantienen una buena relación, obteniéndose el
siguiente ajuste:
92
100,00
90,00
80,00 Ley teórica de
% Disipación de energía
disipación de
70,00 energía 2009
Compuerta (2009)
60,00
50,00 Márquez (2006)
40,00
Vicente (2008)
30,00
20,00
10,00
0,00
2,00 3,00 4,00 5,00 6,00 7,00
F1
Resaltos sumergidos.
100% = + 100%
}6}.}6. .Yu.
= ?= ,= ?= 1
(4.29)
8D = ; :; + D D /2O
Donde:
Energía cinética al principio del resalto, siendo v0 la velocidad en
93
8< = < :< + < < /2O
Donde:
Energía cinetica al final del resalto, donde v4 es la velocidad al
; = < = 1
D = < = 1
y2
1 V '2
Ind. Turb.=
γ q ∫ρ
0
2
udy
y4
1 £ '2
V
,convecc.Turb.1x = ρ ug dy§
γq 2
0 s¨s
94
8D ,8< y ,!. 1s¨967 en la ecuación 4.25.
También, el valor de Disip. Vis. se ha determinado por la sustitución de los valores
Al igual que para el caso de los resaltos libres los valores de N&. ]N.⁄,8D − 8< 1
resultan muy superiores a los de Convecc. Turb.⁄,8D − 8< 1 y por este motivo se puede
concluir que Ind. Turb.⁄,8D − 8< 1 es despreciable. Entonces,
De acuerdo con esto, 567 se interpreta como la longitud de la zona requerida para disipar
el total de la energía perdida en el resalto sumergido, y la pérdida de energía entre el
inicio (x=0) y el final (x=567 ) del resalto sumergido puede ser expresada mediante la
siguiente ecuación unidimensional:
89 = + :; − + :<
(4.32)
, ? 1 A, ⁄ 1 B
=
=>
= , ⁄ 1?
(4.33)
4.23.
ecuaciones 4.33 y 4.15 se muestra en la figura 4.9 a con líneas continuas para cada
valor de considerado.
95
1.5<y3/y0<2.5
2.5<y3/y0<3.5
3.5<y3/y0<4.5
0,50
4.5<y3/y0<6
6<y3/y0<8
2,4 2,9 3,4 3,9 4,4 4,9 Ley teórica para y3/y0=8
F0
Figura 4.9 a. Pérdida relativa de energía en resaltos sumergidos 85 ⁄80 = E A00 , :3B
:
0
Se puede observar que los datos medidos se nos salen por encima de las leyes teóricas,
es decir, nos dan una mayor pérdida de energía, Pero esto es solo aparente, puesto que
los resaltos sumergidos se alejan de la configuración teórica.
Para ver la diferencia que hay entre la configuración de flujo obtenida y la teórica,
ajustamos el valor del calado conjugado y3, despejando de la ecuación 4.33, para
ponerlo en su posición correcta en la gráfica.
96
1.5<y3/y0<2.5
2.5<y3/y0<3.5
3.5<y3/y0<4.5
0,50 4.5<y3/y0<6
6<y3/y0<8
Los resultados experimentales obtenidos los hemos representado junto con los de
Márquez (2006) y se comprueba que, aun siendo modificados, ajustan bien a las leyes
teóricas representadas para cada valor de 0D .
menores, para los que no se han podido conseguir grandes valores de sumergencia,
considerando la limitación de la dimensión del canal de nuestro laboratorio.
97
0,6 F0=2.45
F0=2.93
0,5 F0=3.33
F0=4.307
F0=4.318
0,4
F0=2,87 Márquez 2006
HL/H0
0,2
F0=4,57 Márquez 2006
y3/y0
Figura 4.9 c. Pérdida de energía relativa para resaltos sumergidos: GH ⁄GI = J KI ,
LM
LI
Como hemos comentado para los resaltos libres, tenemos otra forma de evaluar la
disipación de energía en los resaltos, mediante la ecuación 4.23.
Para calcular la disipación de energía en los resaltos sumergidos hemos utilizado las
siguientes expresiones:
Siendo:
En la figura 4.10 se han representado los datos obtenidos en el laboratorio junto con los
de Márquez (2006) y la ley teórica. Podemos observar que los datos experimentales
tienen mayor disipación de energía que la obtenida con la ley teórica.
También, si comparamos ésta gráfica con la de disipación de energía para los resaltos
libres (ver Figura.4.7) se comprueba que, para los resaltos sumergidos tenemos mayor
disipación de energía que para los valores de los resaltos libres.
98
60,00
50,00
% Disipación de energía
40,00 Disipación de
energía para
velocidades medias
Compuerta (2009)
30,00
Márquez (2006)
20,00
10,00
0,00
1,00 2,00 3,00 4,00 5,00 6,00
F0
Comparando las figuras 4.8 y 4.11, veremos que la reducción de energía por
velocidades máximas es similar para un mismo número de Froude, tanto en el caso de
los resaltos libres como para los resaltos sumergidos, ya que las velocidades máximas
medidas en el perfil de la sección 4 son similares en ambos casos. Este hecho es
contradictorio con la teoría, ya que la reducción de las velocidades máximas debería de
ser mayor en los resaltos sumergidos, ya que el calado será superior (siendo superior, a
su vez, el área de la sección) y a mayor área e igualdad de caudal, la velocidad obtenida
será inferior.
99
90,00
80,00
60,00
50,00
20,00
10,00
0,00
2,00 2,50 3,00 3,50 4,00 4,50 5,00 5,50
F0
100
4.4 Longitud del resalto hidráulico.
f (L / H L , H L / H* ) = 0 (4.37)
Resaltos libres.
En el caso del resalto libre, la ecuación 4.37 puede ser expresada como:
f (L j / H L , H L / H1 ) = 0 (4.38)
f ( L j / H L , H L / H1 ) = 0 (4.39)
En la figura 4.13 de acuerdo con los datos experimentales obtenidos por nuestras
mediciones junto con los de Márquez (2006) y Vicente (2008), y con la ecuación 4.39
se aprecia una clara correlación, donde también se representa la ley de Ohtsu et
al.(1990) que sigue la ecuación 4.40.
101
Lj HL
log10 = −1,71 + 1,58 ; en (2,3≤ F1 ≤9,5) y (0,14≤ HL/H1 ≤0,71) (4.40)
HL H1
RhL C
100
Vicente (2008)
Márquez (2006)
Lj/HL
Ajuste propuesto
Lj/HL Márquez
(2006)
Ajuste propuesto
Lj/HL Vicente
(2008)
Ajuste propuesto
1 Lj/HL Compuerta
(2009)
0,0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1,0
HL/H1
De acuerdo con Márquez (2006) y Vicente (2008) se puede interpretar que la pérdida de
energía relativa más grande HL/H1 se da cuando el remolino superficial y las
fluctuaciones de la velocidad turbulenta son mayores. Así, cuando el gradiente de la
línea de energía HL/H1 es mayor, entonces Lj/HL será menor.
Comparando la ley teórica propuesta por Ohtsu et al. (1990) con los resultados
experimentales se observa que los datos experimentales para los resaltos siguen la
misma tendencia y pendiente que las teóricas, lo que nos indica que el criterio de
longitud Lj cumple el criterio establecido.
102
Márquez (2006):
Lj HL
log10 = −1,71 + 1,49 ; para (2,87≤ F1 ≤4,96) (4.41)
HL H1
Vicente (2008):
Lj HL
log10 = −1,71 + 1,53 ; para (1,96≤ F1 ≤5,45) (4.42)
HL H1
Nuestro ajuste es igual al obtenido por Ohtsu et al (1990). Esto puede ser debido a que
hemos empleado la ecuación de Silvester (1964) como criterio a la hora de establecer la
longitud del resalto y este es un valor teórico.
Los ajustes propuestos por Márquez (2006) y Vicente (2008) son prácticamente
similares.
L j / H L = f ( F1 ) (4.44)
103
Resaltos hidráulicos sumergidos.
En el caso del resalto sumergido, la ecuación 4.38 puede ser expresada mediante la
siguiente forma:
Analizando los datos experimentales en base a la ecuación 4.45, se aprecia una buena
correlación que se muestra en la figura.4.14, con la siguiente ecuación:
!OD = = −1,71 => + 1,58; para ,2,3 ≤ 0D ≤ 101 : ,1 ≤ :; ⁄:D ≤ 201 (4.46)
9VR =
>
una expresión general que incluye el caso del resalto libre ,:; ⁄:D = 11. El hecho de
La ecuación 4.46 es idéntica a la ecuación 4.40 para el resalto libre, la cual constituye
que las ecuaciones 4.40 y 4.46 sean idénticas tiene una explicación coherente dado que,
como se ha demostrado, la pérdida de energía HL en resaltos sumergidos es mayor que
en resaltos libres y, en consecuencia, la longitud de resalto en los resaltos sumergidos
será algo mayor que en los resaltos libres.
100,00
Compuerta (2009)
Márquez (2006)
Lsj/HL
Ajuste propuesto
Lsj/HL Márquez
(2006)
Ajuste propuesto
1,00 Lsj/HL
Compuerta (2009)
0,00 0,20 0,40 0,60 0,80 1,00
HL/H0
104
Analizando las leyes teóricas junto con los resultados experimentales obtenidos por
Márquez (2006) y en nuestro proyecto se observa que la nube de puntos experimentales
sigue la misma pendiente que las leyes teóricas verificándose que el criterio de longitud
del resalto, Lsj seguido por este proyecto es el correcto. A continuación se muestra la ley
de ajuste para los valores de Márquez (2006), la cual es similar a la de Compuerta
(2009).
Márquez (2006)
Compuerta (2009)
105
4.5 Campos de velocidad media en los resaltos hidráulicos.
Resaltos hidráulicos libres.
• Distribución de velocidad en dirección negativa (hacia atrás)
106
En la figura 4.15 se representan los valores adimensionales de los registros obtenidos,
Y
YZ [
representado y , con una ley de distribución de velocidad, ecuación. 4.49 y una
u 1 y y
= { ( )}1 / n , 0 ≤ ≤ k (4.49)
um k Y Y
u 1 1,177 y y
= exp[− { ( − k )}2 ], k ≤ ≤ 1,5 (4.50)
um 2 1− k Y Y
donde k = y máx / Y
Distribución de Velocidad K n
Ya que en los resaltos libres no se han podido medir muchos puntos cercanos a la
superficie debido a que no había calado suficiente, añadimos unos puntos teóricos
aplicando la formulación de Hagen & Vischer de velocidades máximas negativas para
107
los puntos ubicados en la superficie del resalto por medio de la ecuación (4.50) y la
ecuación 4.51:
N = \6 ] (4.50)
\6 = −N!
^?D,
D,;`
(4.51)
Siendo:
X = la posición dentro del resalto
V2 =la velocidad al final del resalto.
Pero, ya que la ecuación dada por Hagen & Vischer tenía un pequeño error numérico, se
tuvo que realizar una pequeña modificación de la ecuación anterior, ecuación 4.51, a la
ecuación corregida, ecuación 4.52:
\6 = −N!
^?D,
,
(4.52)
108
ohtsu a
ohtsu b
7,5 Ley propuesta UPCT
X/Lj=0.5 (R.L.0) A
7,0 X/Lj=0.75 (R.L.0)A
X/Lj=0.5 (R.L.1) A
X/Lj=0.75 (R.L.1)A
6,5 X/Lj=0.5 (R.L.2)A
X/Lj=0.75 (R.L.2)A
X/Lj=0.25( R.S.1) (4.37<y4/y0<13.10) A
6,0 X/Lj=0.5 (R.S.1) (4.47<y4/y0<13.10) A
X/Lj=0.75(R.S.1) (4.37<y4y0<13.10) A
X/Lj=0.25 (R.S.2)( 5.90<y4/y0<13.195) A
5,5 X/Lj=0.5 (R.S.2) (5.90<y5/y0<13.195) A
X/Lj=0.75 (R.S.2) ( 5.90< y5/y0<13.195) A
5,0 X/Lj=0.25 (R.S.3)(y4/y0= 6.06) A
X/Lj=0.5 (R.S.3)(y4/y0= 6.06) A
X/Lj=0.75 (R.S.3) (y4/y0= 6.06) A
4,5 X/Lj=0.5 (R.L.0) C
X/Lj=0.25 (R.L.0)C
X/Lj=0.75 (R.L.0)C
4,0 X/Lj=0.25 (R.L.1) C
y/Y
X/Lj=0.5 (R.L.1)C
X/Lj=0.75 (R.L.1) C
3,5 X/Lj=0.5 (R.L.2)C
X/Lj=0.75(R.L.2)C
X/Lj=0.5 (R.L.3)C
3,0 X/Lj=0.75(R.L.3)C
X/Lj=0.25 (R.S.1)( 2.62<y4/y0< 4.48)
2,5 X/Lj=0.5(R.S.1) (2.62<y4/y0< 4.48) C
X/Lj=0.75 (R.S.1) ( 2.62<y4/y0<4.48) C
X/Lj=0.25 (R.S.2)( 2.76<y4/y0< 4.92 C
2,0 X/Lj=0.5 (R.S.2) (2.76< y4/y0<4.92) C
X/Lj=0.75 (R.S.2) (2.76<y4/y0< 4.92) C
X/Lj=0.25(R.S.3) (2.92<y4/y0<5.52) C
1,5 X/Lj=0.5(R.S.3) (2.92<y4/y0<5.52) C
X/Lj=0.75 (R.S.3) (2.92<y4/y0<5.52) C
X/Lj=0.25 (R.S.4) (3.05<y4/y0<6.2) C
1,0 X/Lj=0.5 (R.S.4) (3.05<y4/y0<6.2) C
X/Lj=0.75 (R.S.4)( 3.05<y4/y0<6.2) C
0,5 X/Lj=0.25(R.S.5) (3.02<y4/y0<4.12) C
X/Lj=0.5 (R.S.5) (3.02<y4/y0< 4.12) C
X/Lj=0.75(R.S.5) ( 3.02<y4/y0<4.12) C
0,0 X/Lj=0.25 (R.S.6) (3.45<y4/y0<5.06) C
X/Lj=0.5(R.S.6) (3.45<y4/y0<5.06) C
-0,8 -0,7 -0,6 -0,5 -0,4 -0,3 -0,2 -0,1 0,0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1,0 1,1 1,2 X/Lj=0.75 (R.S.6) (3.45<y4/y0<5.06) C
X/Lj=0.25(R.S.7) (3.88<y4/y0<6) C
X/Lj=0.5(R.S.7) (3.88<y4/y0<6) C
u/um X/Lj=0.75 (R.S.7) (3.88<y4/y0<6) C
Hagen & Vischer
Figura 4.16. Distribución de velocidades inicial para los resaltos hidráulicos libres y sumergidos en compuerta y aliviadero.
109
Como se puede observar en la figura 4.16 hay registros que se desvían de la
distribución. Estos datos han sido analizados y no se dan por válidos ya que se deben a
malas mediciones (ver figura 4.17)
Por lo complejo de este fenómeno, al representar la ley para las velocidades negativas
los puntos siguen una tendencia, pero se encuentran dispersos.
Hemos realizado un ajuste para resalto hidráulico libre para compuerta mediante la
ecuación 4.53, la cual se indica en la figura 4.18.
[
(4.53)
110
ohtsu a
ohtsu b
Ley propuesta UPCT
7,5 X/Lj=0.5 (R.L.0) A
X/Lj=0.75 (R.L.0)A
7,0 X/Lj=0.5 (R.L.1) A
X/Lj=0.75 (R.L.1)A
X/Lj=0.5 (R.L.2)A
6,5 X/Lj=0.75 (R.L.2)A
X/Lj=0.25( R.S.1) (4.37<y4/y0<13.10) A
X/Lj=0.5 (R.S.1) (4.47<y4/y0<13.10) A
6,0 X/Lj=0.75(R.S.1) (4.37<y4y0<13.10) A
X/Lj=0.25 (R.S.2)( 5.90<y4/y0<13.195) A
X/Lj=0.5 (R.S.2) (5.90<y5/y0<13.195) A
5,5 X/Lj=0.75 (R.S.2) ( 5.90< y5/y0<13.195) A
X/Lj=0.25 (R.S.3)(y4/y0= 6.06) A
X/Lj=0.5 (R.S.3)(y4/y0= 6.06) A
5,0 X/Lj=0.75 (R.S.3) (y4/y0= 6.06) A
X/Lj=0.5 (R.L.0) C
X/Lj=0.25 (R.L.0)C
4,5 X/Lj=0.75 (R.L.0)C
X/Lj=0.25 (R.L.1) C
X/Lj=0.5 (R.L.1)C
4,0 X/Lj=0.75 (R.L.1) C
y/Y
X/Lj=0.5 (R.L.2)C
3,5 X/Lj=0.75(R.L.2)C
X/Lj=0.5 (R.L.3)C
X/Lj=0.75(R.L.3)C
3,0 X/Lj=0.25 (R.S.1)( 2.62<y4/y0< 4.48)
X/Lj=0.5(R.S.1) (2.62<y4/y0< 4.48) C
X/Lj=0.75 (R.S.1) ( 2.62<y4/y0<4.48) C
2,5 X/Lj=0.25 (R.S.2)( 2.76<y4/y0< 4.92 C
X/Lj=0.5 (R.S.2) (2.76< y4/y0<4.92) C
X/Lj=0.75 (R.S.2) (2.76<y4/y0< 4.92) C
2,0 X/Lj=0.25(R.S.3) (2.92<y4/y0<5.52) C
X/Lj=0.5(R.S.3) (2.92<y4/y0<5.52) C
X/Lj=0.75 (R.S.3) (2.92<y4/y0<5.52) C
1,5 X/Lj=0.25 (R.S.4) (3.05<y4/y0<6.2) C
X/Lj=0.5 (R.S.4) (3.05<y4/y0<6.2) C
X/Lj=0.75 (R.S.4)( 3.05<y4/y0<6.2) C
1,0 X/Lj=0.25(R.S.5) (3.02<y4/y0<4.12) C
X/Lj=0.5 (R.S.5) (3.02<y4/y0< 4.12) C
X/Lj=0.75(R.S.5) ( 3.02<y4/y0<4.12) C
0,5 X/Lj=0.25 (R.S.6) (3.45<y4/y0<5.06) C
X/Lj=0.5(R.S.6) (3.45<y4/y0<5.06) C
0,0 X/Lj=0.75 (R.S.6) (3.45<y4/y0<5.06) C
X/Lj=0.25(R.S.7) (3.88<y4/y0<6) C
-0,8 -0,7 -0,6 -0,5 -0,4 -0,3 -0,2 -0,1 0,0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1,0 1,1 1,2 X/Lj=0.5(R.S.7) (3.88<y4/y0<6) C
X/Lj=0.75 (R.S.7) (3.88<y4/y0<6) C
u/um Hagen & Vischer
Figura 4.17. Distribución de velocidades analizada para los resaltos hidráulicos libres y sumergidos en compuerta y aliviadero excluyendo los puntos erroneos.
111
6,0
5,5
X/Lj=0.5 (R.L.0)
:
= 1,5 −0,69 # + 0,024 −7,5 #
X/Lj=0.75 (R.L.0)
4,5
f
X/Lj=0.25 (R.L.1)
4,0
X/Lj=0.5 (R.L.1)
X/Lj=0.5 (R.L.2)
3,0
y/Y
X/Lj=0.75(R.L.2)
2,5
X/Lj=0.5 (R.L.3)
2,0 X/Lj=0.75(R.L.3)
1,0 ohtsu b
0,0
-0,8 -0,7 -0,6 -0,5 -0,4 -0,3 -0,2 -0,1
0,1 0,0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1,0 1,1
u/um
112
Resaltos hidráulicos sumergidos.
• Distribución de velocidad en dirección positiva (hacia delante)
g⁄Z = E,:⁄f1
del rango 0,2 < x/Lsj < Lrsj/Lsj, como se puede ver en la figura 4.19:
(4.54)
Los datos experimentales se explican por la ecuación de resalto libre si y4/y0 → y2/ y0.
YZ
(4.55)
donde:
h = [ , siendo Y la distancia perpendicular a la solera, en la que =
YZ
113
La función error es la función primitiva de la distribución de Gauss o Normal:
erf,%1 = q r '
s
√p D
(4.56)
Esta función la hemos acotado para x=0,68η para su límite inferior y 0 para su límite
superior.
En relación a los datos utilizados tenemos algunos en los que la relación y4/y0 es grande
para aplicar la ecuación anterior; y otros en los que la relación anterior se asemeja a
y2/y0 , pudiendo aplicar la ley del resaltos sumergidos junto a resaltos libres. En nuestro
caso, hemos representado en la figura 4.20, la ecuación 4.55, ya que tenemos resaltos
donde la relación y4/y0 es grande. Por este motivo hemos representado los resaltos
sumergidos separados de los resaltos libres.
114
X/Lj=0.25( R.S.1) (4.37<y4/y0<13.10) A
7,5 X/Lj=0.5 (R.S.1) (4.47<y4/y0<13.10) A
X/Lj=0.75(R.S.1) (4.37<y4y0<13.10) A
7,0 X/Lj=0.25 (R.S.2)( 5.90<y4/y0<13.195) A
X/Lj=0.5 (R.S.2) (5.90<y5/y0<13.195) A
6,5
X/Lj=0.75 (R.S.2) ( 5.90< y5/y0<13.195) A
1
= 1,48 h 7 [1 − erf,0.68 h1]
#
X/Lj=0.25 (R.S.3)(y4/y0= 6.06) A
6,0
X/Lj=0.5 (R.S.3)(y4/y0= 6.06) A
X/Lj=0.75 (R.S.3) (y4/y0= 6.06) A
5,5
X/Lj=0.25 (R.S.1)( 2.62<y4/y0< 4.48)
X/Lj=0.5(R.S.1) (2.62<y4/y0< 4.48) C
5,0
X/Lj=0.75 (R.S.1) ( 2.62<y4/y0<4.48) C
X/Lj=0.25 (R.S.2)( 2.76<y4/y0< 4.92 C
4,5
X/Lj=0.5 (R.S.2) (2.76< y4/y0<4.92) C
X/Lj=0.25(R.S.3) (2.92<y4/y0<5.52) C
115
• Distribución de velocidad en dirección negativa (hacia atrás)
= 1,2 a,tcd + 1
c
[
(4.57)
116
7,5
X/Lj=0.25 (R.S.1)( 2.62<y4/y0< 4.48)
:
X/Lj=0.75 (R.S.1) ( 2.62<y4/y0<4.48)
= 1,2 −6.8# + 1
6,5
f
X/Lj=0.25 (R.S.2)( 2.76<y4/y0< 4.92 )
X/Lj=0.5(R.S.3) (2.92<y4/y0<5.52)
5,0
X/Lj=0.75 (R.S.3) (2.92<y4/y0<5.52)
X/Lj=0.5(R.S.6) (3.45<y4/y0<5.06)
2,0 X/Lj=0.75 (R.S.6) (3.45<y4/y0<5.06)
X/Lj=0.25(R.S.7) (3.88<y4/y0<6)
1,5
X/Lj=0.5(R.S.7) (3.88<y4/y0<6)
1,0 X/Lj=0.75 (R.S.7) (3.88<y4/y0<6)
ohtsu b
0,5
Ley propuesta UPCT
u/um
Figura 4.21. Distribucion de velocidades para resalto hi
hidraulico
draulico sumergido en compuerta.
117
Los ajustes obtenidos para compuerta se muestran a continuación, en las tablas 4.3 y
4.4, diferenciando entre ajustes para resaltos libres (ecuación 4.58) y para resaltos
sumergidos (ecuación 4.59).
Resaltos libres:
Distribución de
velocidades negativas a b c d
Resaltos sumergidos:
Distribución de
velocidades negativas f g h
118
4.6 Caída de la velocidad máxima
La ecuación es aplicable en la zona del resalto, dentro del rango 0,2 ≤ x/Lj ≤ 0,7.
En la figura 4.22 se representan las leyes teóricas para cada número de Froude según la
ecuación 4.60, junto con nuestros resultados experimentales, los de Vicente (2008) y los
de Márquez (2006) y se observa que los puntos medidos quedan por debajo de las líneas
teóricas correspondientes, por lo tanto no se puede considerar el ajuste
teórico-experimental como bueno. Se ha de considerar igualmente que la mejor
correlación entre los valores teóricos y los experimentales se obtiene para el primer
Froude medido (F1=1,73), muy similar a lo que ya le ocurrió a Márquez (2006) y a
Vicente (2008), en su estudio en el mismo canal.
Se ha representado también el ajuste medio de las leyes teóricas en el rango del número
de Froude para el rango de trabajo (1≤ F1 ≤ 5,5) ecuación 4.61 y la ley de ajuste que
obtuvo Márquez (2006), ecuación 4.62, Vicente (2008), ecuación 4.63 y con nuestros
datos experimentales ecuación 4.64:
Ajuste propuesto por Ohtsu et al. (1990) independiente del número de Froude :
119
Ajuste propuesto por nuestros resultados experimentales Compuerta 2009:
En las figuras 4.23 y 4.24, se representa la ecuación 4.65 para el caso “b” (flujo no
desarrollado) de Ohtsu et al. (1990), la cual se puede aplicar dentro del rango 0,2 ≤ x/Lj
≤ 0,7:
En la figura 4.23 se han representado las leyes teóricas para los resaltos libres según la
ecuación 4.65 y para todos los números de Froude que se encuentran en el rango de
trabajo.
Se ha representado también la ley teórica correspondiente a los resaltos sumergidos
según la ecuación 4.77, la cual es independiente del número de Froude. Junto a estas
leyes teóricas se han representado los resultados obtenidos experimentalmente en
nuestro proyecto, los de Márquez (2006), para los resaltos libres y los resaltos
sumergidos, y los de Vicente (2008) para resaltos libres.
Según las bases teóricas, los resultados experimentales de nuestros resaltos sumergidos
parámetros medidos :< ⁄:D en los resaltos sumergidos, como hemos dicho
deberían de ajustarse más a la ley teórica correspondiente a resaltos libres, porque los
anteriormente, no son suficientemente grandes, ya que cumple (:< ⁄:D < 20) .Sin
embargo, observando la representación se puede comprobar que no ocurre así, sino que
los valores experimentales son muy dispersos y se ubican entre la ley teórica para el
resalto libre y la ley teórica para el resalto sumergido.
120
Ley teórica propuesta
independiente de F
Ajuste Compuerta 2009
Ajuste Márquez 2006
Ajuste Vicente 2008
1,0 Ley teórica F1=5.5
Ley teórica F1=5
Ley teórica F1=4.5
0,9
Ley teórica F1=4
Ley teórica F1=3.5
0,1 F1=4.32
F1=4.3
F1=3.3
0,0 F1=3
0 5 10 15 20 25 30 F1=2.9
F1=2.8
x/y1 F1=2.6
F1=2.3
Figura 4.23. Relación °⁄L± = J,{⁄L± 1 resalto libre ó °⁄LI = J,{⁄LI 1 resaltos sumergidos.
F0=2.45 (R.S.)
121
En la figura 4.24 se representa el caso de la ymax, como hemos comentado en el apartado
anterior, estando dentro del rango 0,2 ≤ x/Lj ≤ 0,7, pudiéndose aplicar la ecuación 4.66
para el Caso “b”.
Se han representado las leyes teóricas según la ecuación 4.66 para todos los números de
Froude que se encuentran en nuestro rango de trabajo y junto a éstas se ha representado
la ley teórica correspondiente a los resaltos sumergidos según la ecuación (4.78)
(independiente del número de Froude). Comparándolas con los resultados obtenidos
experimentalmente por Márquez (2006), Vicente (2008) y los de este proyecto se
observa que los valores medidos se encuentran entre el caso teórico para resaltos libres
y el caso teórico para resaltos sumergidos. Según las bases teóricas, los resultados para
los resaltos sumergidos deberían ajustarse más a la ley teórica correspondiente a resaltos
medidos :< ⁄:D en los resaltos sumergidos no son suficientemente grandes, ya que
libres, por el mismo motivo que ha sido explicado para la figura 4.29: los parámetros
Figura 4.24. Relación °⁄L± = J,{⁄L± 1 resalto libre ó °⁄LI = J,{⁄LI 1 resaltos sumergidos.
F0=2.45 (R.S.)
122
Si se compara la figura 4.23 con la figura 4.24, observando los datos experimentales
obtenidos por Márquez (2006) y Vicente (2008) y el actual proyecto, se aprecia que
tienen un mejor ajuste con los teóricos esperados en el caso de la figura 4.24. Esto se
explica porque la distancia y máx , en la que u = u m , se obtiene directamente de los
um ∝ xα (4.67)
Y∝x (4.68)
dy u m 2
/( ) = 0 ó valor constate (4.69)
dx gY
Cuando las ecuaciones 4.67, 4.68 y 4.69 se satisfacen aproximadamente dentro del
rango 0,2 ≤ x/Lj ≤ 0,7, se ha obtenido una ley de similitud para la distribución de
velocidades.
u / v 2 = f ( y, y 2 ) (4.70)
123
Márquez F1=4.96
100,00 Márquez F1=4.57
Márquez F1=4.19
90,00
Márquez F1=3.66
Márquez F1=3.56
80,00
Márquez F1=2.87
Vicente F1=5.45
70,00
Vicente F1=4.34
Vicente F1=4.05
60,00
Vicente F1=3.89
y/y2 (%)
Vicente F1=3.77
50,00
Vicente F1=3.49
40,00 Vicente F1=2.97
Vicente F1=2.75
30,00 Vicente F1=1.73
F1=4.318
20,00 F1=4.307
F1=3.33
10,00 F1=3.066
F1=2.983
0,00 F1=2.929
F1=2.729
0,00 0,50 1,00 1,50 2,00
F1=2.446
u/v2
Figura 4.25. Distribución de velocidades para todos los resaltos libres a x=Lj
La caída de Z para el resalto sumergido está entre el caso de chorro de pared y el caso
de resalto libre. En otras palabras, si la profundidad aguas abajo :< ⁄:D
suficientemente grande, tanto que :< ⁄:D ≥ 20 (en este caso, la superficie del agua es
es lo
casi horizontal y :; ⁄:D = :< ⁄:D ), la caída Z coincide con el chorro de pared, como
se indica la ecuación 4.72:
Si la profundidad :< ⁄:D disminuye y cumple con :< ⁄:D < 20, la caída de Z descrita
por la ecuación 4.72 se transforma en la ecuación 4.73:
124
Además, si :< ⁄:D → : ⁄:D, la caída de Z coincide con la del resalto libre (ecuación
4.64). En los casos experimentales en los que se ha trabajado se cumple que :< ⁄:D <
20, por lo que en la figura 4.26 los resaltos sumergidos deberían ajustarse a la ley
teórica representada para el caso de los resaltos libres.
F0=4.318
1,0 F0=4.307
F0=3.33
0,9
F0=2.929
F0=2.446
0,8
Ley teórica independiente de F0
(y4/y0>20)
0,7 Ley teórica (R.L) F0=5
resalto libre dado ya que nuestros valores :< ⁄:D son inferiores a 20. Esta observación
resultados experimentales se ajustan mejor a la caída de velocidad correspondiente a un
coincide con las bases teóricas explicadas anteriormente. Aún se observa que la caída de
velocidad máxima obtenida experimentalmente tiende a ser inferior a la teórica para
todos los Froude ensayados.
resaltos sumergidos en un rango aproximado de 4 ≤ :< ⁄:D ≤ 10. La ley propuesta por
Para obtener un buen ajuste de los datos experimentales se ha propuesto una ley para
125
Z ⁄D = 4,%⁄:D 1D,et (4.74)
caída de velocidad máxima en los resaltos sumergidos en función de su valor :< ⁄:D,
También se ha representado la ecuación 4.76 clasificando los valores experimentales de
Márquez (R.S)
7<y4/y0<8
0,6 Márquez (R.S)
8<y4/y0<9
0,5
Márquez (R.S)
9<y4/y0<10
0,4
2<y4/y0<3
0,3
3<y4/y0<4
0,2
0,1 4<y4/y0<5
0,0 5<y4/y0<6
0 10 20 30 40 50 60
x/y0
Figura 4.27. Caída de velocidad máxima en resaltos sumergidos: ³´ ⁄µI = J,{⁄LI , L¶ ⁄LI 1
126
experimentales en los que se ha trabajado. Los resultados se muestran en las figuras
4.23 y 4.24 anteriormente representadas junto a los casos de resaltos libres.
0,2 ≤ x/Lj ≤ 0,7 y donde las ecuaciones 4.67, 4.68 y 4.69 son casi satisfechas, y se
obtiene la ley de similitud para la distribución de velocidades.
127
4.6 Analisis frecuencial de algunas de las mediciones por la Transformada
Rápida de Fourier.
4.7.1.1 Introducción.
Para una correcta adquisición, procesamiento e interpretación, de los datos de
laboratorio necesitamos técnicas especiales, las cuales se denominan “Técnicas de
Tratamiento de Señales”.
Hay diversas técnicas de procesamiento de señales que se han desarrollado para la
detección y la estimación.
La detección es la técnica que determina si una señal específica está presente en una
observación, mientras que con la estimación se obtienen los valores de los parámetros
que describen la señal.
Para facilitar la detección y la estimación de sistemas de señales, acudimos al análisis
armónico de datos de extensión finita. La transformada de Fourier es el mecanismo por
el cual somos capaces de realizar la descomposición frecuencial.
A pesar de las ventajas del uso de la transformada discreta de Fourier (DFT), sobre
todo con la implementación del algoritmo de la transformada rápida de Fourier (FFT),
existen limitaciones en su funcionamiento (ventanas y suavizado).
En el intento de resolver estas limitaciones, se han desarrollado algunos procedimientos
de estimación espectral alternativos, como los modelos autoregresivos (AR), medias
móviles (MA), autorregresivos de medias móviles (ARMA), descomposición espectral
de Pisarenko (PHD), método de máxima verosimilitud (MLM)…
Sin embargo, para bajas relaciones señal-ruido (SNR), los estimadores espectrales
modernos no son mejores que el procesamiento convencional FFT.
En este caso, hemos decidido usar como herramienta de análisis el método convencional
FFT, ya que la información que nos proporciona es suficiente.
A continuación describiremos los elementos básicos de una transformada de Fourier,
tanto para parámetro continuo como discreto, por medio de la cual vamos a ser capaces
de realizar una correcta adquisición de los datos del laboratorio.
128
4.7.1.2 Transformada de Fourier.
El análisis de las frecuencias de los procesos estocásticos, se fundamentan en el
concepto de la transformación de Fourier. Por lo tanto, a continuación presentamos los
elementos básicos de las transformadas de Fourier.
|%,'1| '
∞
∞
º = 2¾E tenemos:
El par de transformadas de Fourier en función de la frecuencia en Hertz (Hz)
129
su transformada en el dominio de la frecuencia. Así |¹,º1| representa la distribución
Nos define que la energía de una señal en el dominio del tiempo, es igual a la energía de
como:
Entonces:
… , %D = %,01,% = %, 1, … … % = %, 1, …
¹Ã,º1 = ∑¼
¨¼ %
}»
(4.86)
Donde: ºÅ[−º , º ]; º =
p
¹Ã ,º1 }» º
Y su transformada inversa:
% = q
»W
p » W
(4.87)
130
El par de transformadas de Fourier, en función de la frecuencia en (Hz), vienen
definidas por:
¹Ã,E1 = ∑¼
¨¼ %
} pÆ
(4.88)
Y su transformada inversa:
% = q» ¹g ,E1 } E
»W pÆ
W
(4.89)
Aunque se especifica que ¹g está definida en [−º , º ], la ecuación (4.86) define una
función periódica de periodo ,2¾⁄ 1, hecho que será utilizado frecuentemente, para el
cálculo, que se realiza en un intervalo [0,2º ].
La propiedad de la convolución es análoga a la de la transformación de Fourier. Si
ℎ = ∑¼
7¨¼ % O7 (4.90)
ℎ = % O
Entonces:
(4.93)
131
No obstante, la transformada discreta de Fourier, ecuación 4.86 es impracticable
Ç% È.
generalmente ya que solo se conocerán un número N, limitado, de puntos de sucesión
: = º % (4.94)
º = 0 N < 0 ó ≥ Ë (4.95)
¹Ê ,ÌΩ1 = ∑Ê
¨D %
}ΩÎ
; Ì = 0,1, … … … , Ë − 1 (4.96)
%Î = ∑Ê
¨D ¹Ê ,Ω1 ; Ì = 0,1, … … … , Ë − 1
}ΩÎ
Ê
(4.97)
Debe observarse que ¹Ê ,ÌΩ1 = ¹ Ï,Ì + Ë7 1ΩÐ, siendo j un número entero cualquiera.
Por ello, cuando ÌΩ supera a º , es decir, k es mayor a N/2, los valores de ¹
corresponden a las frecuencias negativas dadas por ,Ì − Ë1Ω.
En lo que sigue, adoptamos la notación de ecuación 4.96 y ecuación 4.97 sin mencionar
este hecho.
La DFFT ecuación 4.96 tiene la propiedad de convolución:
∑Ê
ΨD % OÎ = Ê ∑7¨D ¹ ,ÑΩ1Á,ÑΩ1
Ê }Ω7
(4.98)
132
Es muy importante observar que el cálculo de la DFT del producto de Ç% È, equivale a
la DFT del producto de Ç% È por la ventana Ǻ È por los que las propiedades ecuación
4.92 y 4.93 nos indican que, la DFFT es la convolución circular de la DFT de Ǻ È con
la ventana, circunstancia que es aprovechada en el cálculo para mejorar la precisión de
estimas.
El teorema de Parversal aplicado a la transformada finita discreta de Fourier (DFFT),
nos expresa en forma análoga al caso continuo, que:
∑Ê
¨D |% | = Ê ∑¨D |¹Ê ,Ω1|
Ê
(4.95)
En donde
N = Número de datos muestreados
Ω = Intervalo de frecuencia
función de la frecuencia.
133
4.7.2. Análisis del espectro de las mediciones con ADV.
4.7.2.1. Introducción.
Cuando tomamos medidas con los aparatos de alta resolución de medición “acoustic
doppler velocimeters (ADV)” en el seno de un fluido, se observa que, sufren un ruido
que distorsiona la señal turbulenta.
El concepto de Kolmogoroff (1941) afirma que una señal turbulenta con baja distorsión
por el ruido sigue una pendiente de -5/3 en un rango de altas frecuencias, 1Hz<f <10Hz
y se encuentra horizontal a bajas frecuencias, f <1Hz, tal y como se muestra en la
siguiente figura:
Figura 4.28. Señal turbulenta sin contaminación acústica que cumple el concepto de Kolmogoroff
(1941).
El ruido es producido cuando el movimiento de los átomos de la materia interfiere en la
señal del aparato, es decir, cuando el sonido de la señal emitida por el doppler es
inferior al sonido (ruido) en el ambiente, produciéndose así una contaminación en la
medida de las fluctuaciones de la velocidad turbulenta.
134
4.7.2.2. Análisis de la turbulencia de las series de velocidades mediante FFT.
4.7.2.2.1 Introducción.
El objeto del estudio es hacer un análisis frecuencial de las principales mediciones por
medio del algoritmo numérico de la Transformada Rápida de Fourier (FFT). Para ello,
tomamos mediciones en dos puntos, tal y como se muestra en la figura 4.22 para
realizar diferentes comparaciones, a una altura aproximada de 0.5 cm de la solera del
canal, para que la situación de cada uno sea lo más parecida posible. Las comparaciones
que realizamos son:
135
Con ello queremos comprobar cómo afecta el ruido al espectro obtenido a partir de las
mediciones, y si se cumple el concepto de Kolmogoroff (1941), en las diferentes
frecuencias de medición, dentro y fuera del resalto hidráulico, en resalto hidráulico libre
y sumergido, y en los diferentes ejes.
Para situar las mediciones de nuestro análisis diremos que la sección 3 estará situada a
una distancia de 43 cm aguas abajo de la compuerta, y la sección 5 estará a 115,5 cm de
la compuerta, siendo el caudal de 8,5 m3 /s en el cuadalímetro.
frecuencias para xsy , xy , y xsy xy , de la velocidad media y de la raíz cuadrática media
Para comenzar, se presentan los datos de las mediciones tomadas a las distintas
Tabla 4.5 a. RMS y velocidad media para cada una de las mediciones a las frecuencias 5Hz, 25Hz y
50Hz .
136
Q8,5 m3/s Resalto sumergido
S3
5 Sx' 105,57 0,5 1 10,07
5 Sx'Sy' 6,59 0,5 - 4,37
5 Sy' 4,2 0,5 4 2,92
25 Sx' 108 0,48 3 10,49
25 Sx'Sy' 0,2 0,48 - 5,17
25 Sy' 6,99 0,48 2 4,3
50 Sx' 98,78 0,47 2 14,44
50 Sx'Sy' 28,24 0,47 - 9,31
50 Sy' 6,44 0,47 1 8,93
S5
5 Sx' 20,64 0,43 2 4,91
5 Sx'Sy' -0,34 0,43 - 2,32
5 Sy' 0,29 0,43 1 1,73
25 Sx' 20,81 0,45 2 5,3
25 Sx'Sy' -1,38 0,45 - 2,82
25 Sy' 0,48 0,45 1 2,43
50 Sx' 21,62 0,41 1 6,67
50 Sx'Sy' -2,26 0,41 - 3,78
50 Sy' 0,75 0,41 4 3,78
Tabla 4.5 b. RMS y velocidad media para cada una de las mediciones a las frecuencias 5Hz, 25Hz y
50Hz.
Se muestra que los valores de las desviaciones son mayores en z{′ que en zL′ y z{′ zL′
, esto es debido a que es el eje en el que más velocidad se mide, y por lo tanto las
desviaciones de la velocidad turbulenta en el eje x son superiores. También se aprecia
que cuanto mayor es la frecuencia, mayor es la desviación. Esto se debe a que, al ser
medida con mayor frecuencia, se registran mayores fluctuaciones debidas a esa alta
frecuencia.
137
Para la realización de nuestro análisis frecuencial de los espectros de energía obtenidos
por medio del algoritmo de la transformada rápida de Fourier (FFT), utilizaremos el
programa Microsoft Excel.
El paquete de análisis de Fourier, calcula los coeficientes complejos, Xn, desde las series
de datos, xk. ello requiere que el numero de medidas en la serie de datos sea potencia de
2, siendo N=2m.
Las funciones son calculadas con un intervalo de tiempo entre medida y medida en
función de la frecuencia en Hz, T=1/f, siendo la f la frecuencia de la medición en Hz.
Analizamos las fluctuaciones de la velocidad turbulenta en el eje V’x , en el eje V’y y en
Donde:
3. ]′s es la flutuación de la velocidad turbulenta (en este caso para la normal X).
138
Figura 4.30. Fluctuaciones de la velocidad turbulenta.
139
De ellos, solo analizamos la mitad, ya que el análisis para N datos solo se hace para
N/2+1, según el teorema de Nyquist, el cual dice que los valores para una frecuencia f/2
son simétricos.
A continuación pasamos a analizar las distintas descomposiciones frecuenciales.
comporta la energía y el ruido en xsy , xy y xsy xy en un mismo punto de medición, y a
sección 3, tanto para el resalto libre como para el sumergido, estudiaremos cómo se
140
1
0,1
Sx'
0,01
Sx'Sy'
Sij
0,001 Sy'
p -5/3
0,0001
0,00001
P= -5/3
0,000001
0,001 0,01 0,1 1 10
F(Hz)
Figura 4.31 a. Descomposición frecuencial de la sección 3 del resalto libre para xsy y xy la xsy xy ,
0,1
Sx'
0,01
Sx'Sy'
Sij
0,001 Sy'
p -5/3
0,0001
0,00001
P= -5/3
0,000001
0,001 0,01 0,1 1 10
F(Hz)
Figura 4.31b. Descomposición frecuencial de la sección 3 del resalto sumergido para xsy y xy la
141
Si nos centramos en el ruido, comprobamos que va aumentando conforme aumenta la
frecuencia de descomposición. Este ruido comienza a notarse a una frecuencia de
0.03Hz, donde empieza a distorsionarse la señal, siendo más acusado este aumento a
partir de una frecuencia de 0.1Hz para las mediciones tomadas a una frecuencia de 5Hz.
Este ruido, en algunos casos llega a ser tan elevado que puede distorsionar la señal
obtenida por la velocidad turbulenta. Para conocer como el ruido interfiere en la señal
obtenida en nuestras descomposiciones espectrales, la analizamos con el criterio de
Kolmogoroff (1941), que como se ha explicado anteriormente, afirma que la señal de la
velocidad turbulenta con poco nivel de ruido adquiere una pendiente prácticamente
horizontal en frecuencias bajas, hasta 1Hz y, en frecuencias altas, en un rango
1Hz<f<10Hz aproximadamente, adquiere una pendiente de -5/3, para realizar la
comparación hemos dibujado en todas las gráficas una recta con esta pendiente. La
142
1
0,1
S3RL 25Hz
0,001 S3RL 50Hz
p=-5/3
0,0001
p = -5/3
1E-05
1E-06
0,001 0,01 0,1 1 10
frecuencia Hz
0,1
S3RS 5Hz
0,01
S3RS 25Hz
Sx'
S3RL 50Hz
0,001
p=-5/3
0,0001
1E-05 p = -5/3
1E-06
0,001 0,01 0,1 1 10
frecuencia Hz
Figura 4.32b. Descomposición frecuencial de la sección 3 RS medidas a las frecuencias de
5Hz, 25Hz y 50Hz.
En este caso observamos que la energía varia si una medición se ha hecho a una
frecuencia mayor o menor. Es decir, a una frecuencia de 50Hz (verde) tiene mayor
energía que una de 5Hz (azul), aunque las diferencias entre cada una de ellas no son
muy grandes.
143
En cuanto al ruido, diremos que será parecido en todas las frecuencias de medición,
pero parece algo más elevado a una frecuencia de medición más baja.
Este comienza a ser más intenso a una frecuencia de descomposición de 0.012Hz
cuando la medición se ha hecho a un frecuencia de 5 Hz, a una frecuencia de 0.3Hz
cuando las medición se ha hecho a la frecuencia de 25Hz y a una frecuencia de 0.5Hz
cuando la medición se ha hecho a una frecuencia de 50Hz, por lo que, a frecuencias de
medición más elevadas, empieza a aparecer el ruido a mayor descomposición
frecuencial. La distorsión ocasionada por el ruido la analizamos con el concepto de
Kolmogoroff (1941), y aunque la señal cumple con el concepto, vemos que hay una alta
distorsión para todas las frecuencias medidas, por lo que la señal turbulenta estará
bastante perturbada.
Concluimos que, la sonda capta un nivel de ruido parecido en todas las frecuencias de
medición, pero es algo más elevado a frecuencias de medición más bajas, pero, hay que
tener en cuenta una limitación, ya que las series a frecuencias altas no son
representativas en su totalidad, ya que estas representan a una parte de la serie, es decir,
las series hechas a las frecuencias de medición de 25Hz y a 50Hz no tienen solo 4096
datos, sino son de 8192 y 32768 datos respectivamente, pero el programa que hemos
usado para hacer estas descomposiciones, solo deja hacer la transformada rápida de
Fourier a un valor máximo de 4096 datos, por lo que es un limitación a tener en cuenta.
144
1
0,1
0,01
0,001 Vx
Sij
VxVy
0,0001
Vy
0,00001 p= -5/3
0,000001
P = -5/3
0,000000
0,001 0,01 0,1 1 10
F(Hz)
Figura 4.33 a. Espectro de energía en la sección 5 del resalto libre medido a una frecuencia de 5Hz.
0,1
0,01
0,001 Sx'
Sij
Sx'Sy'
0,0001 Sy'
0,00001 p=-5/3
P = -5/3
0,000001
0,000000
0,001 0,01 0,1 1 10
F(Hz)
Figura 4.33 b. Espectro de energía en la sección 5 del resalto sumergido medido a una frecuencia de
5Hz.
energía de la xsy es superior a la de xsy xy , siendo xsy xy ligeramente superior a la
en la energía, que ocurre lo mismo que hemos comentado en la sección 3, es decir, la
normal xy .
145
ruido observamos que xsy xy y xy tiene una alta distorsión
Centrándonos ahora en el ruido,
0,1
0,01
S5RL 5Hz
Sx'
1E-05 P = -5/3
1E-06
0,001 0,01 0,1 1 10
frecuencia Hz
146
1
0,1
S5RS 5Hz
0,01
S5RS 25Hz
Sx'
S5RS 50Hz
0,001
p=-5/3
0,0001
P = -5/3
1E-05
1E-06
0,001 0,01 0,1 1 10
frecuencia Hz
En este caso observamos que, para la xsy , la medición hecha a una frecuencia de 25Hz
y 50Hz tiene mayor nivel de energía, que otra medida a una frecuencia de 5Hz.
En cuanto el ruido, cabe decir que, por lo general, las pendientes que adquieren los
espectros de energía tienden al concepto de Kolmogoroff (1941), es decir, tienen una
pendiente de -5/3, comenzando esta en la frecuencia de descomposición de 0.09Hz para
una medición hecha a una frecuencia de 5Hz, la medición a una frecuencia más elevada,
como es la de 25Hz, comienza a una frecuencia de descomposición más alta, 0.3 Hz, y
para una medición hecha a 50Hz, la pendiente empieza a 0.8 Hz, por lo que se acaba
concluyendo lo dicho para la sección 3. A mayor frecuencia de medición, mayor es la
frecuencia de descomposición a la que comienza a aparecer el ruido
147
y posteriormente haremos comparaciones de su comportamiento en las diferentes
siendo z{y zLy ligeramente superior a la normal zLy , y también varia con la frecuencia
de medición a la que hemos realizado, ya que, a mayor frecuencia de medición,
tenemos mayor cantidad de energía (aunque es un diferencia poco notoria), y el
ruido comienza a una frecuencia de descomposición más alta.
A continuación vamos a estudiar como varía la energía de una sección a otra, es decir, la
sección 3 y sección 5, y entre ambas configuraciones de resalto, libre y sumergido.
Comenzamos comparando ambas secciones en un resalto libre, que es donde más
diferencia encontramos, ya que la sección 3, como hemos comentado anteriormente, se
sección 5, se encuentra fuera del resalto (flujo turbulento, subcrítico), en la xsy (eje
encuentra dentro del resalto hidráulico (flujo altamente turbulento, inestable) y la
0,1
0,01 S3 LIBRE
S5 LIBRE
Sx'
0,001 p=-5/3
0,0001
1E-05
1E-06
0,001 0,01 0,1 1 10
F(Hz)
Figura 4.35. Descomposición frecuencial de la sección 3 y 5 del resalto libre, en la z{y a una
frecuencia de 5Hz.
148
agua tiene mayor movimiento comparado con la sección 5 debido a las características
del flujo en cada una de ellas.
A continuación pasamos a comparar el comportamiento de la energía en la sección 3 de
un resalto libre, y la misma en otro sumergido, ya que es la sección donde habrá más
diferencia en su comportamiento debido a que está en el seno del resalto hidráulico,
porque la sección 5 tiene características similares. Como ejemplo ponemos una
descomposición espectral de la Sx’ medida a una frecuencia de 5Hz:
0,1
0,01 S3 LIBRE
Sx'
S3 SUMERGIDO
0,001
p=-5/3
0,0001
1E-05
1E-06
0,001 0,01 0,1 1 10
F(Hz)
Por último vamos a centrarnos en el ruido. Para comenzar cabe decir unas cuantas
afirmaciones que hemos obtenido, tanto de la sección 3 como de la sección 5. La
primera es que, a mayor frecuencia de medición, en nuestro caso 50Hz, menor
distorsión producida por el ruido (aunque esa diferencia es pequeña), la distorsión
149
producida por el ruido será similar en Sx’, Sy’ y Sx’Sy’ de un mismo punto a una
misma frecuencia de medición, y por último es, analizando las series de acuerdo con
el concepto de Kolmogoroff (1941), vemos que la señal de los espectros obtenidos en la
sección 5 se encuentra menos distorsionada por el ruido, ya que la señal adquiere
una pendiente de -5/3 a frecuencias altas, cumpliéndose así el concepto de Kolmogoroff
(1941), pero por lo contrario, la sección 3, aunque la señal tiende a este concepto, está
más distorsionada. A continuación vamos a mostrar una descomposición frecuencial
que servirá como ejemplo y afirmará lo dicho, es en la Sx’, hecha a una frecuencia de
25Hz, en la sección 3 y 5 de un resalto libre:
0,1 -5/3
0,01 S3 LIBRE
Sx'
S5 LIBRE
0,001 p=-5/3
0,0001
1E-05
1E-06
0,001 0,01 0,1 1 10
F(Hz)
Figura 4.37. Descomposición frecuencial medida a 25Hz, de la sección 3 y 5 del resalto sumergido.
Esto es debido a la diferencia de las características del flujo en una zona y en otra,
evidentemente, la señal está más perturbada en la sección 3, donde es un flujo inestable,
que en la sección 5, donde el flujo es subcrítico
Otra comparación que vamos a realizar es en la misma sección, estando libre y
sumergida, para ello hemos escogido la sección 3, que es donde hay mayor diferencia en
cuanto al flujo. Como ejemplo hemos puesto una descomposición frecuencial medida a
25Hz, de la Sx’ :
150
1
0,1
0,01
S3 LIBRE
Sx'
0,001
S3 SUMERGIDO
P = -5/3 p=-5/3
0,0001
1E-05
1E-06
0,001 0,01 0,1 1 10
F(Hz)
Figura 4.38. Descomposición frecuencial medida a 25Hz, en la Sx’ , de la sección 3 de los dos
resaltos, libre y sumergido.
En este caso observamos que ambas cumple con el concepto de Kolmogoroff (1941),
pero se aprecia que tenemos mayor cantidad de ruido en la sección 3 del resalto
libre que en la sección 3 del sumergido, este hecho es debido a que en el libre hay
mayor movimiento, debido al menor calado con respecto al sumergido, ya que a
igualdad de caudal entre los dos resaltos, habrá mayor velocidad en el libre. Para ver
más gráficas, mirar el anexo IV.
151
5 CONCLUSIONES.
Configuración de la sonda
Una vez contrastados los datos obtenidos en este proyecto fin de carrera con trabajos
experimentales anteriores y hecho el análisis crítico correspondiente, se extraen las
siguientes conclusiones en los diversos campos de estudio que se han tratado.
Campos de estudio
L j / y1 = 9,75 ( F1 − 1) 1, 01
(5.1)
152
• Relación entre profundidades inicial y final
y2 (5.2)
= ( 8 F12 + 1 − 1) / 2
y1
= +
(5.3)
Respecto a la caída máxima de velocidad se aprecia una reducción de los valores experi,
mentales aproximadamente un 30-40% en comparación con la ley de Ohtsu.
Referente al análisis frecuencial de las principales mediciones por medio del algoritmo
numérico de la Transformada Rápida de Fourier (FFT) concluimos que:
153
1. A mayor frecuencia de medición, hay mayor desviación.
6. La energía de xsy es superior a la de xsy xy , siendo xsy xy ligeramente superior a
la normal xy . .
154
ANEXO I
RESULTADOS DEL
FILTRADO A PARTIR DE LAS
CAMPAÑAS DE MUESTREO
155
En este primer anexo se muestran las tablas con los resultados del filtrado a partir de las
campañas de muestreo. En la primera columna se muestra la denominación de cada
punto que se describe a continuación:
A. Las tres primeras cifras de cada archivo muestran el valor del caudal (m3/h),
procedente de la lectura del caudalímetro.
El siguiente cuadro refleja las características de cada uno de estos resaltos, con el caudal
correspondiente, la apertura de compuerta y el valor de cada sumergencia:
Para los caudales 7, 9,5 y 11 m3/h no se tomaron medidas en resaltos sumergidos. Para
los caudales más altos, es decir para los caudales 10,5 m3/h y 12 m3/h solo se pudieron
156
tomar medidas en dos sumergencias, ya que la tercera no se podía hacer por desborde
del agua, aguas arriba de la compuerta.
E. La letra última indica el punto de muestreo de cada sección, siendo la letra “a”
la correspondiente al punto de muestreo más cercano a la solera del canal para
cada una de las secciones. A continuación se muestra un ejemplo de la manera
de nombrar cada punto:
157
En la segunda columna se muestra el valor de velocidad media considerado tras el
proceso del filtrado de la muestra original se denomina ux.
158
RESALTOS LIBRES
Datos filtrados del resalto hidráulico libre a pie de compuerta, Q=4.5.
Archivo ux (cm/s) Distancia al % Calado Filtros Calado X (cm) Lj (cm)
límite (cm) superior
Sección 3
0452103a 36,20 0,55 8,80 1 6,25 33,80 14,49
Sección 4
0452104a 39,76 0,47 7,01 1 6,7 43,80 14,49
0452104b 43,29 0,72 10,75 1 6,7 43,80 14,49
Sección 5
0452105a 24,89 0,57 8,38 1 6,8 96,30 14,49
0452105b 27,04 0,95 13,97 1 6,8 96,30 14,49
159
Datos filtrados del resalto hidráulico libre a pie de compuerta, Q=7.
Archivo ux (cm/s) Distancia al % Calado Filtros Calado X (cm) Lj (cm)
límite (cm) superior
Sección 2
0072902a 65,47 0,49 7,10 1 6,90 24,90 16,40
0072902b 55,19 0,59 8,55 1 6,90 24,90 16,40
Sección 3
0072903a 36,20 0,42 5,45 3 7,70 35,20 16,40
0072903b 42,76 0,98 12,73 2 7,70 35,20 16,40
0072903c 49,05 1,35 17,53 1 7,70 35,20 16,40
0072903d 51,00 2,50 32,47 1 7,70 35,20 16,40
Sección 4
0072904a 25,80 0,34 4,28 1 7,95 45,70 16,40
0072904b 32,84 1,09 13,71 1 7,95 45,70 16,40
0072904c 36,38 1,41 17,74 1 7,95 45,70 16,40
0072904d 36,37 1,80 22,64 1 7,95 45,70 16,40
0072904e 43,29 2,29 28,81 1 7,95 45,70 16,40
Sección 5
0072905a 28,41 0,38 4,75 1 8 96,30 14,49
0072905b 30,07 1,05 13,13 1 8 96,30 14,49
0072905c 34,29 1,51 18,88 1 8 96,30 14,49
0072905d 31,61 1,94 24,25 1 8 96,30 14,49
0072905e 33,85 2,62 32,75 1 8 96,30 14,49
160
Datos filtrados del resalto hidráulico libre a pie de compuerta, Q=8.5.
Archivo ux (cm/s) Distancia al % Calado Filtros Calado X (cm) Lj (cm)
límite (cm) superior
Sección 2
0853302a 57,37 0,50 6,02 3 8,30 30,20 14,09
0853302b 38,17 1,35 16,27 2 8,30 30,20 14,09
Sección 3
0853303a 29,31 0,30 3,57 1 8,40 42,90 14,09
0853303b 35,14 0,32 3,81 4 8,40 42,90 14,09
0853303c 42,71 0,97 11,55 2 8,40 42,90 14,09
0853303d 49,54 2,87 34,17 1 8,40 42,90 14,09
Sección 4
0853304a 29,38 0,27 3,12 1 8,65 55,50 14,09
0853304b 35,43 0,89 10,29 1 8,65 55,50 14,09
0853304c 37,27 1,23 14,22 1 8,65 55,50 14,09
0853304d 38,69 2,31 26,71 1 8,65 55,50 14,09
0853304e 38,01 3,33 38,50 1 8,65 55,50 14,09
Sección 5
0853305a 34,77 0,53 6,09 1 8,7 98,00 14,09
0853305b 37,28 1,59 18,28 1 8,7 98,00 14,09
0853305c 38,05 2,6 29,89 1 8,7 98,00 14,09
161
Datos filtrados del resalto hidráulico libre a pie de compuerta, Q=9.5
Archivo ux (cm/s) Distancia al % Calado Filtros Calado X (cm) Lj
límite (cm) superior
Sección 1
0953901a 61,09 0,57 8,99 1 6,34 17,077 45,77
Sección 2
0953902a 61,16 0,59 7,05 1 8,37 28,52 45,77
0953902b 81,59 1,86 22,22 1 8,37 28,52 45,77
0953902c 37,88 2,39 28,55 1 8,37 28,52 45,77
0953902d 7,81 8,27 98,81 4 8,37 28,52 45,77
Sección 3
0953903a 49,47 0,54 6,27 1 8,61 39,97 45,77
0953903b 64,85 1,06 12,31 1 8,61 39,97 45,77
0953903c 63,71 1,71 19,86 1 8,61 39,97 45,77
0953903d 61,82 2,85 33,10 1 8,61 39,97 45,77
0953903e 44,18 3,53 41,00 1 8,61 39,97 45,77
0953903f 16,22 8,27 96,05 4 8,61 39,97 45,77
Sección 4
0953904a 38,98 0,33 3,71 1 8,90 51,41 45,77
0953904b 47,50 1,67 18,76 1 8,90 51,41 45,77
0953904c 38,63 3,11 34,94 1 8,90 51,41 45,77
0953904d 33,54 4,41 49,55 1 8,90 51,41 45,77
0953904e 25,73 8,27 92,92 3 8,90 51,41 45,77
Sección 5
0953905a 38,18 0,72 8,04 1 8,96 113,46 45,77
0953905b 38,77 1,09 12,17 1 8,96 113,46 45,77
0953905c 39,60 2,02 22,54 1 8,96 113,46 45,77
0953905d 39,88 3,13 34,93 1 8,96 113,46 45,77
0953905e 40,18 3,78 42,19 1 8,96 113,46 45,77
0953905f 35,40 8,27 92,30 1 8,96 113,46 45,77
162
Datos filtrados del resalto libre a pie de compuerta,Q=10.5
Archivo ux (cm/s) Distancia al % Calado Filtros Calado X (cm) Lj
límite (cm) superior
Sección 2
1053902a 85,89 0,37 4,40 1 8,40 26,85 43,88
1053902b 74,04 0,93 11,07 1 8,40 26,85 43,88
1053902c 53,79 1,72 20,48 1 8,40 26,85 43,88
1053902d -0,94 8,27 98,45 4 8,40 26,85 43,88
Sección 3
1053903a 57,57 0,43 4,57 1 9,40 37,82 43,88
1053903b 68,00 0,93 9,89 1 9,40 37,82 43,88
1053903c 67,50 1,91 20,32 1 9,40 37,82 43,88
1053903d 62,76 2,66 28,30 1 9,40 37,82 43,88
1053903e 57,92 2,98 31,70 1 9,40 37,82 43,88
1053903f 7,02 8,27 87,98 4 9,40 37,82 43,88
Sección 4
1053904a 36,70 0,40 3,98 1 10,05 48,79 43,88
1053904b 45,99 1,10 10,95 1 10,05 48,79 43,88
1053904c 49,06 1,54 15,32 1 10,05 48,79 43,88
1053904d 25,51 8,27 82,29 4 10,05 48,79 43,88
Sección 5
1053905a 35,94 0,39 3,77 1 10,35 112,14 43,88
1053905b 36,95 0,79 7,63 1 10,35 112,14 43,88
1053905c 37,70 1,31 12,66 1 10,35 112,14 43,88
1053905d 36,78 2,02 19,52 1 10,35 112,14 43,88
1053905e 36,77 2,47 23,86 1 10,35 112,14 43,88
1053905f 36,30 3,07 29,66 1 10,35 112,14 43,88
1053905g 35,75 3,62 34,98 1 10,35 112,14 43,88
1053905h 39,85 8,25 79,71 1 10,35 112,14 43,88
163
Datos filtrados del resalto libre a pie de compuerta,Q=11
Archivo ux (cm/s) Distancia al % Calado Filtros Calado X (cm) Lj
límite (cm) superior
Sección 1
0114101a 106,17 0,36 5,67 1 6,35 19,63 54,50
0114101b 65,90 0,79 12,44 1 6,35 19,63 54,50
Sección 2
0114102a 69,08 0,76 8,89 1 8,55 33,25 54,50
0114102b 85,13 1,03 12,05 1 8,55 33,25 54,50
0114102c 54,80 2,89 33,80 1 8,55 33,25 54,50
0114102d 3,33 8,27 96,73 4 8,55 33,25 54,50
Sección 3
0114103a 44,03 0,85 9,04 1 9,40 46,88 54,50
0114103b 52,73 1,11 11,81 3 9,40 46,88 54,50
0114103c 61,13 2,11 22,45 1 9,40 46,88 54,50
0114103d 54,51 3,36 35,74 1 9,40 46,88 54,50
0114103e 27,75 8,27 87,98 4 9,40 46,88 54,50
Sección 4
0114104a 39,95 0,42 4,29 1 9,80 60,50 54,50
0114104b 42,87 1,14 11,63 1 9,80 60,50 54,50
0114104c 44,35 1,77 18,06 4 9,80 60,50 54,50
0114104d 46,95 2,63 26,84 4 9,80 60,50 54,50
0114104e 46,83 3,43 35,00 2 9,80 60,50 54,50
0114104f 48,26 4,31 43,98 3 9,80 60,50 54,50
0114104g 38,46 8,27 84,39 4 9,80 60,50 54,50
Sección 5
0114105a 40,31 0,37 3,70 1 9,99 118,00 54,50
0114105b 41,54 1,05 10,51 1 9,99 118,00 54,50
0114105c 41,99 1,59 15,92 1 9,99 118,00 54,50
0114105d 42,50 2,37 23,72 1 9,99 118,00 54,50
0114105e 42,50 3,20 32,03 1 9,99 118,00 54,50
0114105f 43,06 4,18 41,84 1 9,99 118,00 54,50
0114105g 41,30 8,27 82,78 1 9,99 118,00 54,50
164
Datos filtrados del resalto libre a pie de compuerta, Q=12
Archivo ux (cm/s) Distancia al % Calado Filtros Calado X (cm) Lj
límite (cm) superior
Sección 1
0124301a 88,58 0,69 8,21 1 8,40 20,36 56,79
0124301b 68,60 1,08 12,86 1 8,40 20,36 56,79
0124301c -4,71 8,27 98,45 4 8,40 20,36 56,79
Sección 2
0124302a 81,11 0,51 5,41 1 9,43 34,56 56,79
0124302b 70,50 2,04 21,63 1 9,43 34,56 56,79
0124302c 62,73 2,78 29,48 1 9,43 34,56 56,79
0124302d 42,85 3,51 37,22 1 9,43 34,56 56,79
0124302e 35,92 4,12 43,69 1 9,43 34,56 56,79
0124302f 5,12 8,27 87,70 4 9,43 34,56 56,79
Sección 3
0124303a 51,12 0,54 5,14 1 10,50 48,76 56,79
0124303b 57,14 1,25 11,90 1 10,50 48,76 56,79
0124303c 69,11 2,51 23,90 1 10,50 48,76 56,79
0124303d 63,50 3,55 33,81 1 10,50 48,76 56,79
0124303e 46,37 4,68 44,57 1 10,50 48,76 56,79
0124303f 33,83 5,14 48,95 1 10,50 48,76 56,79
0124303g 33,38 8,25 78,57 4 10,50 48,76 56,79
Sección 4
0124304a 50,02 1,03 9,32 1 11,05 62,95 56,79
0124304b 51,72 4,01 36,29 1 11,05 62,95 56,79
0124304c 51,36 4,60 41,63 1 11,05 62,95 56,79
0124304d 51,98 5,29 47,87 1 11,05 62,95 56,79
0124304e 49,38 5,77 52,22 1 11,05 62,95 56,79
0124304f 43,75 5,93 53,67 1 11,05 62,95 56,79
0124304g 39,74 8,25 74,66 4 11,05 62,95 56,79
0124304h 38,44 8,83 79,91 4 11,05 62,95 56,79
Sección 5
0124305a 42,28 0,36 3,22 1 11,17 119,23 56,79
0124305b 43,41 0,98 8,77 1 11,17 119,23 56,79
0124305c 44,40 1,66 14,86 1 11,17 119,23 56,79
0124305d 43,69 2,39 21,40 1 11,17 119,23 56,79
0124305e 44,09 3,22 28,83 1 11,17 119,23 56,79
0124305f 44,46 3,78 33,84 1 11,17 119,23 56,79
0124305g 45,39 4,65 41,63 1 11,17 119,23 56,79
0124305h 45,23 5,34 47,81 1 11,17 119,23 56,79
0124305i 42,33 8,25 73,86 1 11,17 119,23 56,79
0124305j 42,42 9,61 86,03 1 11,17 119,23 56,79
165
RESALTOS SUMERGIDOS
Datos filtrados del resalto hidráulico sumergido a pie de compuerta (1ªSum.), Q=4.5.
Archivo ux (cm/s) Distancia al % Calado Filtros Calado X (cm) Lsj (cm)
límite (cm) superior
Sección 2
0452112a 41,39 0,44 5,87 2,00 7,50 13,60 14,49
Sección 3
0452113a 21,65 0,35 4,16 2 8,30 33,80 14,49
0452113b 33,25 1,10 13,08 3 8,30 33,80 14,49
0452113c 24,29 1,67 19,86 2 8,30 33,80 14,49
0452113d 0,88 8,27 98,34 4 8,30 33,80 14,49
Sección 4
0452114a 27,16 2,03 23,88 4 8,50 43,80 14,49
0452114b 12,42 8,27 97,29 2 8,50 43,80 14,49
Sección 5
0452115a 17,57 0,56 6,51 1 8,60 96,30 14,49
0452115b 18,56 0,95 11,05 1 8,60 96,30 14,49
0452115c 19,45 1,54 17,91 1 8,60 96,30 14,49
0452115d 20,15 2,15 25,00 3 8,60 96,30 14,49
0452115e 19,12 8,27 96,16 3 8,60 96,30 14,49
166
Datos filtrados del resalto hidráulico sumergido a pie de compuerta (2ªSum.), Q=4.5.
Archivo ux (cm/s) Distancia al % Calado Filtros Calado X (cm) Lsj (cm)
límite (cm) superior
Sección 1
0452121a 94,64 0,36 6,67 1 5,40 13,60 14,49
0452121b 95,18 0,66 12,22 1 5,40 13,60 14,49
Sección 2
0452122a 99,17 0,39 5,62 3 6,94 23,70 14,49
0452122b 95,70 0,67 9,65 1 6,94 23,70 14,49
Sección 3
0452123a 11,37 0,22 2,64 4 8,40 33,80 14,49
0452123b 42,03 1,81 21,55 3 8,40 33,80 14,49
0452123c -6,57 8,27 98,45 4 8,40 33,80 14,49
Sección 4
0452124a 31,06 1,89 20,77 3 9,10 43,80 14,49
0452124b 29,17 2,61 28,68 4 9,10 43,80 14,49
0452124c 2,24 8,27 90,88 4 9,10 43,80 14,49
Sección 5
0452125a 10,55 0,42 4,53 1 9,27 96,30 14,49
0452125b 16,02 1,24 13,38 1 9,27 96,30 14,49
0452125c 17,00 2,11 22,76 1 9,27 96,30 14,49
0452125d 18,45 2,61 28,16 4 9,27 96,30 14,49
0452125e 16,61 8,27 89,21 3 9,27 96,30 14,49
167
Datos filtrados del resalto hidráulico sumergido a pie de compuerta (3ªSum.), Q=4.5.
Archivo ux (cm/s) Distancia al % Calado Filtros Calado X (cm) Lsj (cm)
límite (cm) superior
Sección 1
0452131a 97,58 0,38 4,56 1 8,34 13,60 14,49
0452131b 72,98 0,86 10,31 1 8,34 13,60 14,49
0452131c 34,77 1,40 16,79 1 8,34 13,60 14,49
0452131d 0,15 2,34 28,06 3 8,34 13,60 14,49
0452131e -12,25 8,27 99,16 4 8,34 13,60 14,49
Sección 2
0452132a 15,35 0,34 4,02 3 8,45 23,70 14,49
0452132b 70,00 0,88 10,41 1 8,45 23,70 14,49
0452132c 37,58 1,48 17,51 1 8,45 23,70 14,49
0452132d 10,24 3,11 36,80 4 8,45 23,70 14,49
0452132e -17,80 8,27 97,87 4 8,45 23,70 14,49
Sección 3
0452133a 50,77 0,41 4,39 2 9,35 33,80 14,49
0452133c 40,41 1,64 17,54 1 9,35 33,80 14,49
0452133d 32,80 2,99 31,98 3 9,35 33,80 14,49
0452133e -12,77 8,27 88,45 4 9,35 33,80 14,49
0452133f -12,23 9,20 98,40 4 9,35 33,80 14,49
Sección 4
0452134a 33,24 1,82 17,17 3 10,60 43,80 14,49
0452134b 28,65 3,17 29,91 3 10,60 43,80 14,49
0452134c -2,63 8,27 78,02 4 10,60 43,80 14,49
0452134d -5,31 9,18 86,60 4 10,60 43,80 14,49
Sección 5
0452135a 6,54 0,41 3,81 1 10,76 96,30 14,49
0452135b 12,00 1,01 9,39 1 10,76 96,30 14,49
0452135c 12,66 1,44 13,38 1 10,76 96,30 14,49
0452135d 12,66 1,87 17,38 1 10,76 96,30 14,49
0452135e 13,07 2,36 21,93 1 10,76 96,30 14,49
0452135f 11,61 2,85 26,49 1 10,76 96,30 14,49
0452135g 14,50 4,74 44,05 3 10,76 96,30 14,49
0452135h 17,80 8,27 76,86 2 10,76 96,30 14,49
0452135i 18,16 9,15 85,04 2 10,76 96,30 14,49
168
Datos filtrados del resalto hidráulico sumergido a puie de compuerta (1ªSum.), Q=5.5.
Archivo ux (cm/s) Distancia al % Calado Filtros Calado X (cm) Lsj (cm)
límite (cm) superior
Sección 2
0552112a 61,74 0,38 4,551 1 8,35 26,14 45,68
0552112b 62,31 0,84 10,060 1 8,35 26,14 45,68
0552112c 56,96 1,07 12,814 1 8,35 26,14 45,68
0552112d -6,9 8,27 99,042 4 8,35 26,14 45,68
Sección 3
0552113a 32,86 0,43 5,119 1 8,40 37,56 45,68
0552113b 36,13 0,67 7,976 1 8,40 37,56 45,68
0552113c 37,04 1,06 12,619 1 8,40 37,56 45,68
0552113d 39,71 1,39 16,548 1 8,40 37,56 45,68
0552113e 39,67 1,66 19,762 1 8,40 37,56 45,68
0552113f 38,3 1,97 23,452 1 8,40 37,56 45,68
0552113g 38,02 2,15 25,595 1 8,40 37,56 45,68
0552113h 4,4 8,27 98,452 4 8,40 37,56 45,68
Sección 4
0552114a 20,18 0,38 4,368 1 8,70 48,97 45,68
0552114b 22,22 0,58 6,667 1 8,70 48,97 45,68
0552114c 24,14 1,08 12,414 1 8,70 48,97 45,68
0552114d 25,45 1,27 14,598 2 8,70 48,97 45,68
0552114e 27,36 1,62 18,621 1 8,70 48,97 45,68
0552114f 27,46 1,97 22,644 1 8,70 48,97 45,68
0552114g 27,9 2,27 26,092 1 8,70 48,97 45,68
0552114h 27,93 2,54 29,195 1 8,70 48,97 45,68
0552114i 12,4 8,27 95,057 4 8,70 48,97 45,68
Sección 5
0552115a 15,43 0,38 4,318 1 8,8 112,23 45,68
0552115b 15,05 0,66 7,500 1 8,8 112,23 45,68
0552115c 19,21 1,02 11,591 1 8,8 112,23 45,68
0552115d 21,17 1,25 14,205 1 8,8 112,23 45,68
0552115e 21,03 1,56 17,727 1 8,8 112,23 45,68
0552115f 21,38 1,88 21,364 1 8,8 112,23 45,68
0552115g 21,43 2,16 24,545 1 8,8 112,23 45,68
0552115h 21,54 2,57 29,205 1 8,8 112,23 45,68
0552115i 17,44 8,27 93,977 1 8,85 112,23 45,68
169
Datos filtrados del resalto hidráulico sumergido a pie de compuerta (2ªSum), Q=5.5
Archivo ux (cm/s) Distancia al % Calado Filtros Calado X (cm) Lsj (cm)
límite (cm) superior
Sección 0
0552120a 95,27 0,59 7,02 1 8,40 3,30 45,68
0552120b 75,00 0,86 10,24 1 8,40 3,30 45,68
0552120c 1,79 8,27 98,45 4 8,40 3,30 45,68
Sección 1
0552121a 73,45 0,61 7,30 1 8,36 14,72 45,68
0552121b 67,51 0,73 8,73 1 8,36 14,72 45,68
0552121c 66,96 0,83 9,93 1 8,36 14,72 45,68
0552121d 59,77 0,95 11,36 1 8,36 14,72 45,68
0552121e -8,85 8,27 98,92 4 8,36 14,72 45,68
Sección 2
0552122a 48,16 0,64 7,62 1 8,40 26,14 45,68
0552122b 68,45 0,80 9,52 1 8,40 26,14 45,68
0552122c 59,49 0,96 11,43 1 8,40 26,14 45,68
0552122d 20,88 1,38 16,43 1 8,40 26,14 45,68
0552122e -14,55 8,27 98,45 4 8,40 26,14 45,68
Sección 3
0552123a 36,40 0,68 7,56 1 9,00 37,56 45,68
0552123b 46,78 0,98 10,89 1 9,00 37,56 45,68
0552123c 44,59 1,20 13,33 1 9,00 37,56 45,68
0552123d 44,85 1,54 17,11 1 9,00 37,56 45,68
0552123e 41,88 1,81 20,11 1 9,00 37,56 45,68
0552123f 40,21 2,11 23,44 1 9,00 37,56 45,68
0552123g 35,81 2,45 27,22 1 9,00 37,56 45,68
0552123h -11,85 8,27 91,89 4 9,00 37,56 45,68
Sección 4
0552124a 12,76 0,35 3,50 1 10,00 48,98 45,68
0552124b 22,90 0,89 8,90 1 10,00 48,98 45,68
0552124c 25,17 1,33 13,30 1 10,00 48,98 45,68
0552124d 27,34 1,88 18,80 1 10,00 48,98 45,68
0552124e 29,45 2,39 23,90 1 10,00 48,98 45,68
0552124f 26,78 2,84 28,40 1 10,00 48,98 45,68
0552124g 22,35 3,33 33,30 1 10,00 48,98 45,68
0552124h 17,85 3,47 34,70 1 10,00 48,98 45,68
0552124i 20,07 8,27 82,70 1 10,00 48,98 45,68
0552124j 11,83 8,69 86,90 2 10,00 48,98 45,68
0552124k 11,30 8,88 88,80 3 10,00 48,98 45,68
170
Sección 5
0552125a 6,81 0,46 4,47 1 10,30 112,24 45,68
0552125b 13,36 0,90 8,74 1 10,30 112,24 45,68
0552125c 14,35 1,48 14,37 1 10,30 112,24 45,68
0552125d 15,53 1,99 19,32 1 10,30 112,24 45,68
0552125e 16,37 2,49 24,17 1 10,30 112,24 45,68
0552125f 12,10 2,97 28,83 1 10,30 112,24 45,68
0552125g 10,93 3,48 33,79 1 10,30 112,24 45,68
0552125h 10,15 3,73 36,21 1 10,30 112,24 45,68
0552125i 21,73 8,27 80,29 1 10,30 112,24 45,68
0552125j 21,73 8,32 80,78 1 10,30 112,24 45,68
171
Datos filtrados del resalto hidráulico sumergido a pie de compuerta (3ªSum), Q=5.5
Archivo ux (cm/s) Distancia al % Calado Filtros Calado X (cm) Lsj (cm)
límite (cm) superior
Sección 0
0552130a 133,68 0,40 4,49 1 8,90 3,30 45,68
0552130b 100,68 0,67 7,53 1 8,90 3,30 45,68
0552130c 92,02 0,96 10,79 1 8,90 3,30 45,68
0552130d 41,22 1,25 14,04 1 8,90 3,30 45,68
0552130e 12,78 1,59 17,87 1 8,90 3,30 45,68
0552130f 8,89 1,77 19,89 1 8,90 3,30 45,68
0552130g -0,81 2,10 23,60 1 8,90 3,30 45,68
0552130h -5,19 2,31 25,96 1 8,90 3,30 45,68
0552130i -1,68 8,27 92,92 4 8,90 3,30 45,68
0552130j -2,54 8,68 97,53 4 8,90 3,30 45,68
Sección 1
0552131b 22,63 0,55 6,40 1 8,60 14,72 45,68
0552131c 103,13 0,81 9,42 1 8,60 14,72 45,68
0552131d 52,09 1,25 14,53 1 8,60 14,72 45,68
0552131g 28,00 1,92 22,33 1 8,60 14,72 45,68
0552131h -19,12 8,27 96,16 4 8,60 14,72 45,68
Sección 2
0552132a 32,95 0,38 4,18 1 9,10 26,14 45,68
0552132b 82,64 0,65 7,14 1 9,10 26,14 45,68
0552132c 74,52 0,95 10,44 1 9,10 26,14 45,68
0552132d 40,27 1,31 14,40 1 9,10 26,14 45,68
0552132f 36,01 1,77 19,45 1 9,10 26,14 45,68
0552132g 36,37 2,08 22,86 1 9,10 26,14 45,68
0552132h -20,47 8,27 90,88 4 9,10 26,14 45,68
Sección 3
0552133a 21,62 0,40 3,77 1 10,60 37,56 45,68
0552133b 52,78 0,80 7,55 1 10,60 37,56 45,68
0552133c 57,06 1,14 10,75 1 10,60 37,56 45,68
0552133e 45,05 2,34 22,08 1 10,60 37,56 45,68
0552133f 44,70 2,79 26,32 1 10,60 37,56 45,68
0552133g -10,31 8,27 78,02 4 10,60 37,56 45,68
0552133h -12,79 8,69 81,98 4 10,60 37,56 45,68
0552133i -13,93 9,31 87,83 4 10,60 37,56 45,68
172
Sección 4
0552134a 16,69 0,38 3,17 1 12,10 48,98 45,68
0552134b 22,51 0,89 7,42 1 12,10 48,98 45,68
0552134c 27,00 1,35 11,25 1 12,10 48,98 45,68
0552134d 34,47 1,90 15,83 1 12,10 48,98 45,68
0552134e 35,52 2,66 22,17 1 12,10 48,98 45,68
0552134h 28,92 3,82 31,83 1 12,10 48,98 45,68
0552134i 6,44 8,27 68,92 3 12,10 48,98 45,68
0552134j -2,05 8,66 72,17 4 12,10 48,98 45,68
0552134k -2,59 8,96 74,67 3 12,10 48,98 45,68
Sección 5
0552135a 12,07 0,37 3,01 1 12,25 112,24 45,68
0552135b 13,74 1,02 8,29 1 12,25 112,24 45,68
0552135c 14,94 1,78 14,47 1 12,25 112,24 45,68
0552135d 14,75 2,39 19,43 1 12,25 112,24 45,68
0552135f 14,24 3,96 32,20 1 12,25 112,24 45,68
0552135g 15,39 4,60 37,40 1 12,25 112,24 45,68
0552135h 15,85 5,04 40,98 1 12,25 112,24 45,68
0552135i 14,44 8,27 67,24 3 12,25 112,24 45,68
0552135j 18,49 10,24 83,25 1 12,25 112,24 45,68
0552135k 17,55 10,86 88,29 1 12,25 112,24 45,68
0552135k 17,55 10,86 88,29 1 12,25 112,24 45,68
173
Datos filtrados del resalto hidráulico sumergido a pie de compuerta (1ªSum.), Q=8.5
Archivo ux (cm/s) Distancia al % Calado Filtros Calado X (cm) Lsj (cm)
límite (cm) superior
Sección 1
0853311a 44,52 0,63 8,58 1 7,34 17,6 50,4
0853311b 119,64 1,02 13,90 4 7,34 17,6 50,4
Sección 2
0853312a 41,05 0,55 6,17 3 8,92 30,3 50,4
0853312b 43,65 1,09 12,22 2 8,92 30,3 50,4
0853312c 41,00 1,63 18,27 2 8,92 30,3 50,4
0853312d 37,43 1,97 22,09 3 8,92 30,3 50,4
Sección 3
0853313a 28,06 0,72 7,44 4 9,68 42,9 50,4
0853313b 29,06 1,05 10,85 4 9,68 42,9 50,4
0853313c 32,45 1,78 18,39 3 9,68 42,9 50,4
0853313d 32,27 2,34 24,17 2 9,68 42,9 50,4
0853313e 29,63 2,87 29,65 3 9,68 42,9 50,4
0853313f 16,85 8,27 85,43 3 9,68 42,9 50,4
Sección 4
0853314a 21,09 0,53 5,31 2 9,95 55,49 50,4
0853314b 23,64 0,91 9,11 3 9,95 55,49 50,4
0853314c 24,39 1,40 14,01 3 9,95 55,49 50,4
0853314d 25,45 1,87 18,72 3 9,95 55,49 50,4
0853314e 25,46 2,39 23,92 3 9,95 55,49 50,4
0853314f 23,46 2,91 29,13 2 9,95 55,49 50,4
0853314g 26,67 8,27 82,78 3 9,95 55,49 50,4
Sección 5
0853315a 20,68 0,55 5,39 1 10,2 115,5 50,4
0853315b 20,82 1,06 10,39 1 10,2 115,5 50,4
0853315c 21,88 1,55 15,20 2 10,2 115,5 50,4
0853315d 21,77 2,03 19,90 2 10,2 115,5 50,4
0853315e 22,04 2,51 24,61 1 10,2 115,5 50,4
0853315f 19,90 2,96 29,02 3 10,2 115,5 50,4
0853315g 29,47 3,91 38,33 3 10,2 115,5 50,4
0853315h 29,82 8,27 81,08 2 10,2 115,5 50,4
174
Datos filtrados del resalto hidráulico sumergido a pie de compuerta(2ªSum.),Q=8.5
Archivo ux (cm/s) Distancia al % Calado Filtros Calado X (cm) Lsj (cm)
límite (cm) superior
Sección 0
0853320a 135,25 0,35 4,17 1 8,40 5,10 50,39
0853320b 135,72 1,06 12,62 2 8,40 5,10 50,39
0853320c 113,41 1,44 17,14 1 8,40 5,10 50,39
0853320d -6,21 8,27 98,45 2 8,40 5,10 50,39
Sección 1
0853321a 111,27 0,34 4,11 1 8,28 17,60 50,39
0853321b 100,16 0,99 10,67 3 8,28 17,60 50,39
0853321c 72,60 1,64 15,31 4 8,28 17,60 50,39
0853321d -19,21 8,27 71,54 4 8,28 17,60 50,39
Sección 2
0853322a 82,75 0,47 5,06 2 9,28 30,30 50,39
0853322b 79,23 1,16 12,50 4 9,28 30,30 50,39
0853322c 66,73 1,69 18,21 4 9,28 30,30 50,39
0853322d 69,89 2,85 30,71 3 9,28 30,30 50,39
0853322e -14,76 8,27 89,12 4 9,28 30,30 50,39
Sección 3
0853323a 47,40 0,42 3,92 4 10,71 42,90 50,39
0853323b 54,64 1,04 9,71 2 10,71 42,90 50,39
0853323c 54,13 1,57 14,66 3 10,71 42,90 50,39
0853323d 52,85 2,19 20,45 3 10,71 42,90 50,39
0853323e 40,90 2,82 26,33 3 10,71 42,90 50,39
0853323g 30,86 3,90 36,41 3 10,71 42,90 50,39
0853323h -0,11 8,27 77,22 4 10,71 42,90 50,39
Sección 4
0853324a 25,98 0,42 3,63 4 11,56 55,50 50,39
0853324b 30,77 0,96 8,30 4 11,56 55,50 50,39
0853324c 33,91 1,56 13,49 3 11,56 55,50 50,39
0853324d 37,47 2,20 19,03 4 11,56 55,50 50,39
0853324e 35,30 2,85 24,65 4 11,56 55,50 50,39
0853324g 36,80 4,27 36,94 4 11,56 55,50 50,39
0853324h 35,63 4,89 42,30 3 11,56 55,50 50,39
0853324i 15,96 8,27 71,54 3 11,56 55,50 50,39
0853324j 11,40 9,82 84,95 4 11,56 55,50 50,39
0853324k 8,91 10,38 89,79 4 11,56 55,50 50,39
175
Sección 5
0853325a 21,88 0,45 3,85 3 11,70 115,50 50,39
0853325b 22,62 1,03 8,80 3 11,70 115,50 50,39
0853325c 23,82 1,55 13,25 2 11,70 115,50 50,39
0853325d 24,22 2,18 18,63 3 11,70 115,50 50,39
0853325e 22,24 2,72 23,25 4 11,70 115,50 50,39
0853325g 25,12 4,11 35,13 3 11,70 115,50 50,39
0853325h 25,07 4,59 39,23 3 11,70 115,50 50,39
0853325i 25,38 5,43 46,41 3 11,70 115,50 50,39
0853325j 25,22 8,27 70,68 3 11,70 115,50 50,39
0853325k 25,44 9,00 76,92 4 11,70 115,50 50,39
0853325l 25,45 9,92 84,79 3 11,70 115,50 50,39
176
Datos filtrados del resalto hidráulico sumergido a pie de compuerta (3ªSum.), Q=8.5
Archivo ux (cm/s) Distancia al % Calado Filtros Calado X (cm) Lsj (cm)
límite (cm) superior
Sección 0
0853330a 134,71 0,91 9,47 3 9,61 5,10 50,40
0853330b 118,09 1,50 15,61 2 9,61 5,10 50,40
0853330c 35,85 2,12 22,06 3 9,61 5,10 50,40
0853330d 8,82 2,64 27,47 4 9,61 5,10 50,40
0853330e -4,73 3,22 33,51 4 9,61 5,10 50,40
0853330f -7,41 3,65 37,98 4 9,61 5,10 50,40
0853330g -9,19 8,27 86,06 4 9,61 5,10 50,40
0853330h -1,86 9,30 96,77 4 9,61 5,10 50,40
Sección 1
0853331a 124,14 0,40 4,17 1 9,59 17,60 50,40
0853331b 110,56 0,95 9,91 2 9,59 17,60 50,40
0853331c 42,97 2,27 23,67 1 9,59 17,60 50,40
0853331d 21,36 2,89 30,14 3 9,59 17,60 50,40
0853331e 7,29 3,44 35,87 4 9,59 17,60 50,40
0853331f -25,24 8,27 86,24 4 9,59 17,60 50,40
Sección 2
0853332a 100,54 0,41 3,97 3 10,34 30,30 50,40
0853332b 94,56 1,02 9,86 2 10,34 30,30 50,40
0853332c 87,93 2,25 21,76 2 10,34 30,30 50,40
0853332d 66,28 2,76 26,69 3 10,34 30,30 50,40
0853332e 8,18 3,33 32,21 4 10,34 30,30 50,40
0853332f 4,82 3,86 37,33 4 10,34 30,30 50,40
0853332g -23,46 8,27 79,98 4 10,34 30,30 50,40
0853332h -19,40 9,50 91,88 4 10,34 30,30 50,40
Sección 3
0853333a 42,56 0,37 3,19 4 11,10 42,89 50,40
0853333b 65,97 1,00 8,62 4 11,10 42,89 50,40
0853333c 80,80 2,20 18,97 4 11,10 42,89 50,40
0853333d 69,50 2,82 24,31 4 11,10 42,89 50,40
0853333e 33,97 4,67 40,26 3 11,10 42,89 50,40
0853333f -8,65 8,27 71,29 4 11,10 42,89 50,40
0853333g -11,78 9,42 81,21 4 11,10 42,89 50,40
177
Sección 4
0853334a 22,06 0,43 3,33 4 12,90 55,50 50,40
0853334b 38,65 0,96 7,44 3 12,90 55,50 50,40
0853334c 42,15 1,60 12,40 3 12,90 55,50 50,40
0853334d 32,11 3,06 34,26 4 12,90 55,50 50,40
0853334e 27,13 5,61 43,49 4 12,90 55,50 50,40
0853334f 10,05 8,27 64,11 4 12,90 55,50 50,40
0853334g 2,69 9,78 75,81 4 12,90 55,50 50,40
0853334h -3,61 10,71 83,02 4 12,90 55,50 50,40
Sección 5
0853335a 19,04 0,62 4,70 4 13,20 115,50 50,40
0853335b 20,07 1,56 11,82 4 13,20 115,50 50,40
0853335c 21,19 2,36 17,88 3 13,20 115,50 50,40
0853335d 15,69 3,21 24,32 3 13,20 115,50 50,40
0853335e 22,14 4,04 30,61 3 13,20 115,50 50,40
0853335f 21,92 4,97 37,65 3 13,20 115,50 50,40
0853335g 22,40 5,74 43,48 4 13,20 115,50 50,40
0853335h 22,35 6,85 51,89 3 13,20 115,50 50,40
0853335i 21,82 8,27 62,65 3 13,20 115,50 50,40
0853335j 22,05 9,71 73,56 4 13,20 115,50 50,40
0853335k 22,07 10,87 82,35 3 13,20 115,50 50,40
0853335l 22,03 11,44 86,67 4 13,20 115,50 50,40
178
Datos filtrados del resalto hidráulico sumergido a pie de compuerta (1ªSum.), Q=10.5.
Archivo ux (cm/s) Distancia al % Calado Filtros Calado X (cm) Lsj (cm)
límite (cm) superior
Sección 1
1053911a 150,34 0,10 1,19 1 8,40 10,43 22,11
1053911b 71,72 1,40 16,67 1 8,40 10,43 22,11
1053911c -26,42 8,27 98,45 4 8,40 10,43 22,11
Sección 2
1053912a 101,31 0,35 3,85 1 9,09 15,96 22,11
1053912c 123,98 1,33 14,63 1 9,09 15,96 22,11
1053912d 95,65 2,21 24,31 1 9,09 15,96 22,11
1053912e -27,39 8,27 90,98 3 9,09 15,96 22,11
Sección 3
1053913a 74,13 0,38 3,42 1 11,10 21,49 22,11
1053913b 95,22 0,91 8,20 1 11,10 21,49 22,11
1053913c 72,00 2,17 19,55 1 11,10 21,49 22,11
1053913d 70,14 2,19 19,73 1 11,10 21,49 22,11
1053913e 50,29 3,04 27,39 1 11,10 21,49 22,11
1053913f -8,16 8,27 74,50 4 11,10 21,49 22,11
Sección 4
1053914a 43,60 0,45 3,59 1 12,52 27,02 22,11
1053914b 55,95 1,32 10,54 1 12,52 27,02 22,11
1053914c 56,72 2,31 18,45 1 12,52 27,02 22,11
1053914e 55,86 4,13 32,99 1 12,52 27,02 22,11
1053914f 48,52 4,67 37,30 1 12,52 27,02 22,11
1053914g 2,19 8,27 66,05 3 12,52 27,02 22,11
1053914h -1,63 10,20 81,47 4 12,52 27,02 22,11
Sección 5
1053915a 25,26 0,32 2,52 1 12,7 101,26 22,11
1053915b 29,24 1,27 10,00 1 12,7 101,26 22,11
1053915c 30,18 2,35 18,50 1 12,7 101,26 22,11
1053915d 25,82 3,39 26,69 1 12,7 101,26 22,11
1053915e 31,2 4,42 34,80 1 12,7 101,26 22,11
1053915f 31,52 5,33 41,97 1 12,7 101,26 22,11
1053915g 31,32 6,15 48,43 1 12,7 101,26 22,11
1053915h 30,62 8,27 65,12 3 12,7 101,26 22,11
1053915i 30,02 9,54 75,12 1 12,7 101,26 22,11
1053915j 30,15 10,18 80,16 2 12,7 101,26 22,11
1053915k 30,03 11,04 86,93 2 12,7 101,26 22,11
179
Datos filtrados del resalto hidráulico sumergido a pie de compuerta (2ªSum.), Q=10.5.
Archivo ux (cm/s) Distancia al % Calado Filtros Calado X (cm) Lsj (cm)
límite (cm) superior
Sección 1
1053921a 96,89 1,35 16,05 1 8,41 10,43 22,11
1053921b 33,05 2,88 34,24 1 8,41 10,43 22,11
1053921c -31,47 8,27 98,34 4 8,41 10,43 22,11
Sección 2
1053922a 134,17 0,10 1,00 1 9,98 15,96 22,11
1053922b 137,95 0,35 3,51 1 9,98 15,96 22,11
1053922c 19,52 2,64 26,45 1 9,98 15,96 22,11
1053922e 1,83 4,52 45,29 1 9,98 15,96 22,11
1053922f -28,76 8,27 82,87 4 9,98 15,96 22,11
Sección 3
1053923a 116,90 0,02 0,17 1 12,00 21,49 22,11
1053923c 94,41 2,04 17,00 1 12,00 21,49 22,11
1053923d 23,30 3,57 29,75 1 12,00 21,49 22,11
1053923g 10,90 5,23 43,58 2 12,00 21,49 22,11
1053923h -13,13 8,50 86,67 4 12,00 21,49 22,11
Sección 4
1053924a 67,57 0,60 4,43 1 13,54 27,02 22,11
1053924b 75,63 1,56 11,52 1 13,54 27,02 22,11
1053924c 79,75 2,08 15,36 1 13,54 27,02 22,11
1053924d 64,07 3,22 23,78 1 13,54 27,02 22,11
1053924e 38,90 4,38 32,35 1 13,54 27,02 22,11
1053924g 32,59 4,88 36,04 1 13,54 27,02 22,11
1053924h 21,99 5,43 40,10 2 13,54 27,02 22,11
1053924i 0,72 10,50 77,55 2 13,54 27,02 22,11
Sección 5
1053925a 26,36 0,78 5,67 1 13,75 101,26 22,11
1053925b 27,2 1,57 11,42 1 13,75 101,26 22,11
1053925c 29,07 2,65 19,27 1 13,75 101,26 22,11
1053925d 29,94 3,68 26,76 1 13,75 101,26 22,11
1053925e 30,37 4,68 34,04 1 13,75 101,26 22,11
1053925g 30,81 5,71 41,53 1 13,75 101,26 22,11
1053925h 30,82 6,78 49,31 1 13,75 101,26 22,11
1053925i 30,52 7,82 56,87 1 13,75 101,26 22,11
1053925j 27,88 8,27 60,15 3 13,75 101,26 22,11
1053925j 27,76 9,45 68,73 2 13,75 101,26 22,11
1053925k 27,47 10,52 76,51 2 13,75 101,26 22,11
1053925l 27,2 11,49 83,56 2 13,75 101,26 22,11
180
Datos filtrados del resalto hidráulico sumergido a pie de compuerta (1ªSum.), Q=12.
Archivo ux (cm/s) Distancia al % Calado Filtros Calado X (cm) Lsj (cm)
límite (cm) superior
Sección 1
0124311a 127,69 0,32 3,68 1 8,70 20,40 20,60
0124311c 125,81 0,77 8,85 1 8,70 20,40 20,60
0124311d 118,16 0,94 10,80 1 8,70 20,40 20,60
0124311e 101,07 1,27 14,60 1 8,70 20,40 20,60
0124311h 85,31 2,11 24,25 1 8,70 20,40 20,60
0124311i -20,47 8,27 95,06 4 8,70 20,40 20,60
Sección 2
0124312a 86,78 0,28 2,80 1 10,00 34,60 20,60
0124312b 93,73 0,83 8,30 1 10,00 34,60 20,60
0124312e 83,39 2,36 23,60 1 10,00 34,60 20,60
0124312f 71,37 2,91 29,10 1 10,00 34,60 20,60
0124312g 59,87 3,44 34,40 1 10,00 34,60 20,60
0124312h 3,67 8,27 82,70 4 10,00 34,60 20,60
Sección 3
0124313a 51,91 0,45 3,78 1 11,90 48,70 20,60
0124313b 59,04 1,10 9,24 1 11,90 48,70 20,60
0124313c 61,50 1,95 16,39 1 11,90 48,70 20,60
0124313d 64,92 2,61 21,93 3 11,90 48,70 20,60
0124313e 64,48 3,25 27,31 1 11,90 48,70 20,60
0124313f 59,52 3,87 32,52 1 11,90 48,70 20,60
0124313g 56,54 4,55 38,24 1 11,90 48,70 20,60
0124313h 48,41 5,35 44,96 1 11,90 48,70 20,60
0124313i 16,67 8,27 69,50 1 11,90 48,70 20,60
Sección 4
0124314a 35,04 0,38 3,01 1 12,62 62,90 20,60
0124314b 39,20 1,17 9,27 1 12,62 62,90 20,60
0124314c 41,39 1,96 15,53 1 12,62 62,90 20,60
0124314d 44,16 2,63 20,84 1 12,62 62,90 20,60
0124314e 40,33 3,44 27,26 1 12,62 62,90 20,60
0124314f 46,44 4,23 33,52 1 12,62 62,90 20,60
0124314g 45,59 4,95 39,22 1 12,62 62,90 20,60
0124314h 43,08 5,73 45,40 1 12,62 62,90 20,60
0124314i 31,54 8,27 65,53 3 12,62 62,90 20,60
0124314j 21,25 11,10 87,96 3 12,62 62,90 20,60
0124314k 20,60 11,18 88,59 4 12,62 62,90 20,60
181
Sección 5
0124315a 31,81 0,43 3,36 1 12,8 101,13 20,60
0124315b 32,98 1,17 9,14 1 12,8 101,13 20,60
0124315c 33,62 2,11 16,48 1 12,8 101,13 20,60
0124315d 32,03 3,24 25,31 1 12,8 101,13 20,60
0124315e 34,62 4,1 32,03 1 12,8 101,13 20,60
0124315f 34,58 4,81 37,58 1 12,8 101,13 20,60
0124315g 35,09 5,58 43,59 1 12,8 101,13 20,60
0124315h 34,69 6,14 47,97 1 12,8 101,13 20,60
0124315i 32,16 8,27 64,61 1 12,8 101,13 20,60
0124315j 32,08 10,2 79,69 1 12,8 101,13 20,60
0124315k 32,22 11,42 89,22 1 12,8 101,13 20,60
182
Datos filtrados del resalto hidráulico sumergido a pie de compuerta (2ªSum.), Q=12.
Archivo ux (cm/s) Distancia al % Calado Filtros Calado X (cm) Lsj (cm)
límite (cm) superior
Sección 1
0124321a 132,25 0,50 5,62 1 6,70 20,40 20,60
0124321b 127,44 0,84 9,44 1 6,70 20,40 20,60
0124321c 122,33 0,98 11,01 1 6,70 20,40 20,60
0124321d 110,81 1,26 14,16 1 6,70 20,40 20,60
0124321g -30,21 8,27 92,92 4 6,70 20,40 20,60
Sección 2
0124322a 101,71 0,66 6,23 1 10,60 34,50 20,60
0124322d 91,46 2,48 23,40 1 10,60 34,50 20,60
0124322e 76,14 3,14 29,62 1 10,60 34,50 20,60
0124322f 59,38 3,84 36,23 1 10,60 34,50 20,60
0124322g 54,72 4,08 38,49 1 10,60 34,50 20,60
0124322h -8,50 8,90 78,02 4 10,60 34,50 20,60
Sección 3
0124323a 59,03 0,67 5,40 1 12,40 48,70 20,60
0124323b 64,67 1,50 12,10 1 12,40 48,70 20,60
0124323c 68,83 2,43 19,60 1 12,40 48,70 20,60
0124323d 67,15 3,28 26,45 3 12,40 48,70 20,60
0124323e 62,41 4,05 32,66 1 12,40 48,70 20,60
0124323f 54,07 4,88 39,35 1 12,40 48,70 20,60
0124323g 45,23 5,59 45,08 1 12,40 48,70 20,60
0124323h 15,51 10,49 84,60 2 12,40 48,70 20,60
0124323i 5,22 11,79 95,08 4 12,40 48,70 20,60
Sección4
0124324a 42,15 0,76 5,71 1 13,32 62,90 20,60
0124324b 46,61 1,94 14,56 1 13,32 62,90 20,60
0124324c 49,17 2,77 20,80 1 13,32 62,90 20,60
0124324d 49,72 3,89 29,20 1 13,32 62,90 20,60
0124324e 47,01 4,84 36,34 1 13,32 62,90 20,60
0124324f 43,47 5,95 44,67 1 13,32 62,90 20,60
0124324g 40,57 6,38 47,90 1 13,32 62,90 20,60
0124324h 29,43 9,75 73,20 1 13,32 62,90 20,60
0124324i 18,71 12,11 90,92 1 13,32 62,90 20,60
183
Sección 5
0124325a 33,29 0,58 4,28 1 13,55 101,13 20,60
0124325b 35,13 1,74 12,84 1 13,55 101,13 20,60
0124325c 35,23 2,75 20,30 1 13,55 101,13 20,60
0124325d 35,80 3,59 26,49 1 13,55 101,13 20,60
0124325e 36,18 4,59 33,87 1 13,55 101,13 20,60
0124325f 36,49 5,82 42,95 1 13,55 101,13 20,60
0124325g 36,26 6,74 49,74 1 13,55 101,13 20,60
0124325h 32,58 9,91 73,14 1 13,55 101,13 20,60
0124325i 32,62 12,49 92,18 1 13,55 101,13 20,60
184
ANEXO II
GRÁFICAS DE VELOCIDAD
FRENTE A CALADO
185
Una vez recogidos en el anexo anterior los valores de velocidad medidos en cada punto
de muestreo, se representan a continuación los perfiles de velocidad clasificados por
secciones muestreadas y caudal de trabajo para cada resalto hidráulico analizado.
Las velocidades más altas se localizan, lógicamente, al inicio del resalto, debido
fundamentalmente a la influencia del chorro principal saliente de la compuerta, y por
este mismo motivo se localizan a porcentajes de calado bajos. Se observa generalmente
un descenso rápido de las velocidades con el incremento del calado.
186
Perfiles de velocidad a pie de compuerta en resalto
libre, Q=4.5m³/s.
16
14
12 Sección 3
% Calado
10
Sección 4
8
6
Sección 5
4
0
0 10 20 30 40 50
Velocidad (cm/s)
25 Sección 3
20
Sección 4
15
10
Sección 5
5
0
0 10 20 30 40 50 60
Velocidad (cm/s)
187
Perfiles de velocidad a pie de compuerta en resalto
libre, Q=7m³/s.
35
30
Sección 2
25
% Calado
Sección 3
20
15 Sección 4
10
Sección 5
5
0
0 20 40 60
Velocidad (cm/s)
40
35
30 Sección 2
% Calado
25 Sección 3
20 Sección 4
15
Sección 5
10
0
0 20 40 60
Velocidad (cm/s)
188
Perfiles de velocidad a pie de compuerta en resalto
libre, Q=9.5m3/h
120
100
% Calado
80
Sección 1
60 Sección 2
Sección 3
40
Sección 4
20 Sección 5
0
0 20 40 60 80 100
Velocidad (cm/s)
100,00
% Calado
80,00
Sección 2
60,00 Sección 3
Sección 4
40,00 Sección 5
20,00
0,00
0 20 40 60 80 100
velocidad (cm/s)
189
Perfiles de velocidad a pie de compuerta en resalto
sumergido,Q=11m3/h
120
100
% Calado
80
Sección 1
60 Sección 2
Sección 3
40
Sección 4
20
Sección 5
0
0 20 40 60 80 100 120
Velocidad (cm/s)
100
% calado
80 Sección 1
Sección 2
60
Sección 3
40 Sección 4
Sección 5
20
0
0 20 40 60 80 100
Velocidad (cm/s)
190
Perfiles de velocidad a pie de compuerta en resalto
sumergido (1ªSum.),Q=4.5 m3/s.
120
100
% calado
80
60 Seccion 2
Sección 3
40 Sección 4
Sección 5
20
0
0 10 20 30 40 50
Velocidad (cm/s)
70 Sección 1
% Calado
Sección 2
60
Sección 3
50
Sección 4
40
Sección 5
30
20
10
0
-10 40 90 140
Velocidad (cm/s)
191
Perfiles de velocidad a pie de compuerta en resalto
sumergido (3ªSum), Q=4.5 m³/s.
120
100
Sección 0
80 Sección 1
% Calado
Sección 2
60 Sección 3
Sección 4
40 Sección 5
20
0
-25 -5 15 35 55 75 95 115 135
Velocidad (cm/s)
60 Sección 3
50 Sección 4
40
Sección 5
30
20
10
0
-10 10 30 50 70 90
Velocidad (cm/s)
192
Perfiles de velocidad a pie de compuerta en resalto
sumergido (2ªSum.), Q=5,5m3/h.
100
90
80 Sección 0
70
% Calado
Sección 1
60 Sección 2
50 Sección 3
40 Sección 4
Sección 5
30
20
10
0
-20 0 20 40 60 80 100
Velocidad ( cm/s)
70
60 Sección 2
50
Sección 3
40
Sección 4
30
20 Sección 5
10
0
-20 0 20 40 60 80 100 120 140
Velocidad (cm/s)
193
Perfiles de velocidad a pie de compuerta en resalto
sumergido (1ªSum.), Q=8,5m3/h.
100
90
80 Sección 1
70 Sección 2
% Calado
60
Sección 3
50
Sección 4
40
Sección 5
30
20
10
0
-20 0 20 40 60 80 100 120
Velocidades (cm/s)
70
60 Sección 2
50 Sección 3
40
Sección 4
30
Sección 5
20
10
0
-20 30 80 130
Velocidades (cm/s)
194
Perfiles de velocidad a pie de compuerta en resalto
sumergido (3ªSum.), Q=8,5m3/h.
100
90 Sección 0
80 Sección 1
% Calado
70
Sección 2
60
Sección 3
50
Sección 4
40
30 Sección 5
20
10
0
-30 20 70 120
Velocidades (cm/s)
100 Sección 0
Sección 1
80 Sección 2
% Calado
Sección 3
60 Sección 4
Sección 5
40
20
0
-30 20 70 120 170
Velocidades (cm/s)
195
Perfiles de velocidad a pie de compuerta en resalto
sumergido (2ªSum.), Q=10.5 m³/s.
120
100
Sección 0
80
Sección 1
% Calado
Sección 2
60
Sección 3
Sección 4
40
Sección 5
20
Velocidades (cm/s)
60
Sección 3
50
Sección 4
40
Sección 5
30
20
10
0
-50 0 50 100 150 200
Velocidades (cm/s)
196
Perfiles de velocidad a pie de compuerta en resalto
sumergido (2ªSum.), Q=12 m³/s.
100
90
80 Sección 1
70 Sección 2
60
% Calado
Sección 3
50
Sección 4
40
Sección 5
30
20
10
0
-50 0 50 100 150
Velocidad (cm/s)
197
ANEXO III
CARACTERÍSTICAS
ASOCIADAS A LOS
RESALTOS HIDRÁULICOS
ANALIZADOS
198
En este tercer anejo se recogen las características asociados a cada uno de los resaltos
hidráulicos analizados.
En la tabla 1 se presentan los principales valores de los resaltos libres en los que hemos
medido velocidades, se incluye el caudal marcado por el caudalímetro, la apertura de
compuerta para cada caudal, que determinará el nº de Froude 0 =
,y el tipo de
Rh=A/P, y ν= viscosidad cinética del agua, que para una Tª=20ºC, tiene un valor de
1,06 10-6 m2/s.
199
- En la 5ª, el número de Froude correspondiente al pie del resalto.
- En la 6ª, el calado y1 al principio del resalto.
- En la 7ª, el calado y2 al final del resalto.
- En la 8ª, la carga de agua, Htotal, aguas arriba de la compuerta.
- En la 9ª, la longitud del resalto Lj .
- En la 10ª, el Re correspondiente, los datos más importantes son los que
tengan un Re > 25.000.
Q Q V1 V2 F1 Y1 Y2 Htotal Lj Reynolds
Caudalímetro Real m/s m/s cm cm cm cm
m3/h m3/s
4,5 0,00149 1,50 0,28 4,32 1,23 6,70 12,60 40,27 14107
5,5 0,00181 1,60 0,32 4,31 1,40 7,05 12,36 45,64 16553
7 0,00228 1,37 0,35 3,07 2,05 7,95 13,14 41,60 18658
8,5 0,00276 1,55 0,39 3,33 2,20 8,65 13,00 50,40 21985
9,5 0,00307 1,40 0,43 2,73 2,70 8,90 14,00 45,78 22687
10,5 0,00339 1,35 0,41 2,45 3,10 10,20 17,40 43,88 23618
11 0,00355 1,56 0,45 2,98 2,80 9,80 15,30 54,50 25801
12 0,00386 1,59 0,43 2,93 3,00 11,00 17,90 56,79 27300
200
Primera sumergencia, h5=1/3(Htotal- y5)+ y5
Segunda sumergencia, h5=2/3(Htotal- y5)+ y5
Tercera sumergencia, h5= (Htotal- y5)+ y5
En la siguiente tabla se presentan los valores del caudal del caudalímetro, el caudal real,
la y5, la Htotal, y la h.
El calado de cada sumergencia en la sección 5, siendo R.S.1 la primera sumergencia,
R.S.2 la segunda sumergencia y R.S.3 la tercera sumergencia.
Los valores que aparecen tachados en la anterior tabla, son de sumergencias que no han
podido hacerse, ya que debido a las dimensiones del canal, no permitía por que el agua
se desborda, aguas arriba de la compuerta, por lo que se tomó como valor de la segunda
sumergencia el máximo para que el canal no se desbordara.
En la tabla que se muestra a continuación, se muestran esos valores de los calados
máximos de la segunda sumergencia del caudal 10,5 y 12.
Resalto
Q Q sumergido
Caudalímetro real R.S.2
m3/h m3/s h5 Máxima
10,5 0,0034 13,75
12 0,0039 13,55
201
En la siguiente tabla, se muestran características principales de los resaltos sumergidos.
En la 1ª columna, es el caudal que marca el caudalímetro.
En la 2ª, es el caudal real.
En la 3ª, el número de sumergencia para cada caudal.
En la 4ª, la velocidad en la apertura de compuerta. ]D =
202
ANEXO IV
DESCOMPOSICIONES
FRECUENCIALES POR
MEDIO DEL ALGORITMO
NUMÉRICO DE LA
TRANSFORMADA RÁPIDA
DE FOURIER FFT
203
En este último anexo, se muestran todas las descomposiciones frecuenciales hechas por
medio de la transformada rápida de Fourier, el objetivo es analizar como el ruido
perturba a la señal turbulenta en las diferentes mediciones, para ello comparamos la
descomposición frecuencial con el concepto de Kolmogoroff (1941) que afirma que
para una señal turbulenta sin ruido, el espectro tendría que seguir la siguiente estructura:
204
ESPECTROS DE ENERGÍA DE z{′ zL′ Y LA CORTANTE z{′ zL′ EN UN
MISMO PUNTO.
0,1
0,01
Sx'
Sij
0,001 Sx'Sy'
Sy'
0,0001
p -5/3
0,00001
0,000001
0,001 0,01 0,1 1 10
F(Hz)
0,1
0,01
Sx'
0,001
Sx'Sy'
Sij
0,0001
Sy'
0,00001 P -5/3
0,000001
0,000000
0,001 0,01 0,1 1 10
F(Hz)
205
Espectros de energía en la S3(libre) a una
frecuencia de 50Hz
10
0,1
Sij
0,01 Vx
VxVy
0,001
Vy
0,0001 p=-5/3
0,00001
0,01 0,1 1 10
F(Hz)
0,1
0,01
0,001
Sij
Vx
0,0001 VxVy
0,00001 Vy
0,000001 p= -5/3
0,000000
0,001 0,01 0,1 1 10
F(Hz)
206
Espectros de energía en la S5(libre) a una
frecuencia de 25Hz
1
0,1
0,01
0,001
0,0001 Sx'
Sij
0,00001
Sx'Sy'
0,000001
Sy'
0,000000
p=-5/3
1E-08
0,001 0,01 0,1 1 10
F(Hz)
0,1
0,01
Sx'
0,001
Sij
Sx'Sy'
0,0001
Sy'
0,00001
p=-5/3
0,000001
0,000000
0,01 0,1 1 F(Hz) 10
207
Espectro de energía en la S3( sum) a una
frecuencia de 5Hz
1
0,1
0,01
0,001
Sij
0,0001 Sx'
0,00001 Sx'Sy'
0,000001 Sy'
p=-5/3
0,000000
1E-08
0,001 0,01 0,1 1 10
F(Hz)
1
0,1
0,01
Sij
0,001 Sx'
Sx'Sy'
0,0001
Sy'
0,00001
p=-5/3
0,000001
0,001 0,01 0,1 1 10
F(Hz)
208
Espectro de energía en la sección 3(sum) a una
frecuencia de 50Hz
10
0,1
0,01 Sx'
Sij
Sx'Sy'
0,001
Sy'
0,0001 p=-5/3
0,00001
0,01 0,1 1 10
F(Hz)
0,0001 Sx'
Sx'Sy'
0,00001
Sy'
0,000001
p=-5/3
0,000000
0,001 0,01 0,1 1 10
F(Hz)
209
Espectro de energía en la S5(SUM) a una
frecuencia de 25Hz
1
0,1
0,01
0,001
0,0001 Sx'
Sij
0,00001 Sx'Sy'
Sy'
0,000001
p=-5/3
0,000000
0,001 0,01 0,1 1 10
F(Hz)
0,0001
Sx'Sy'
0,00001
Sy'
0,000001
p=-5/3
0,000000
0,01 0,1 1 10
F(Hz)
210
ESPECTROS DE ENERGÍA A DISTINTAS FRECUENCIAS EN UNA MISMA
SECCIÓN.
0,1
0,01
Sx'
S3RL 5Hz
0,001 S3RL 25Hz
1E-06
0,001 0,01 0,1 1 10
frecuencia Hz
0,1
0,01
Sx'
S3RS 5Hz
0,001 S3RS 25Hz
1E-06
0,001 0,01 0,1 1 10
frecuencia Hz
211
Espectro de energía en la S5(libre) medida a
f.5Hz, 25Hz y 50Hz
0,1
0,01
Sx'
1E-06
0,001 0,01 0,1 1 10
frecuencia Hz
0,1
Sx'
0,01
S5RS 5Hz
0,001 S5RS 25Hz
1E-06
0,001 0,01 0,1 1 10
frecuencia Hz
212
Espectro de energía en la S3(libre) medidas a las
f. 5Hz, 25Hz y 50Hz
1
0,1
0,01
Sx'Sy'
1E-06
0,001 0,01 0,1 1 10
frecuencia Hz
0,1
0,01
Sx'Sy'
1E-06
0,001 0,01 0,1 1 10
frecuencia Hz
213
Espectro de energía en la S5 (libre) medidas a
f.5Hz, 25Hz y 50Hz
0,1
0,01
Sx'Sy'
1E-06
0,001 0,01 0,1 1 10
frecuencia Hz
0,1
0,01
Sx'Sy'
1E-06
0,001 0,01 0,1 1 10
frecuencia Hz
214
Espectro de energía en la S3(libre) medida a f.
5Hz, 25Hz y 50Hz
1
0,1
0,01
Sy'
1E-06
0,001 0,01 0,1 1 10
frecuencia Hz
0,1
0,01
Sy'
1E-06
0,001 0,01 0,1 1 10
frecuencia Hz
215
Espectro de energía en la S5(libre) medidas a f.
5Hz, 25Hz y 50Hz
1
0,1
0,01
Sy'
1E-06
0,001 0,01 0,1 1 10
frecuencia Hz
0,1
0,01
Sy'
1E-06
0,001 0,01 0,1 1 10
frecuencia Hz
216
ESPECTROS DE ENERGÍA EN LA SECCIÓN 3 Y SECCIÓN 5 DE UNA
MISMA CONFIGURACIÓN DE RESALTO.
1
Sx'
0,1
0,01
0,001 S3 LIBRE
0,0001 S5 LIBRE
p=-5/3
1E-05
1E-06
0,001 0,01 0,1 1 10
F(Hz)
0,1
0,01
Sx'
0,001 S3 LIBRE
0,0001 S5 LIBRE
p=-5/3
1E-05
1E-06
0,001 0,01 0,1 1 10
F(Hz)
217
Espectro de energía en S3-S5 (libre) a 50 Hz
0,1
Sx'
0,01
S3 LIBRE
0,001
S5 LIBRE
0,0001 p=-5/3
1E-05
1E-06
0,01 0,1 1 10
F(Hz)
0,1
0,01
Sx'
0,001
S3 SUM
0,0001 S5 SUM
p=-5/3
1E-05
1E-06
0,001 0,01 0,1 1 10
F(Hz)
218
Espectro de energía S3-S5 (sum) a 25 Hz
0,1
0,01
Sx'
S3 SUM
0,001
S5 SUM
0,0001 p=-5/3
1E-05
1E-06
0,001 0,01 0,1 1 10
F(Hz)
1
Sx'
0,1
0,01
S3 SUM
0,001 S5 SUM
p=-5/3
0,0001
1E-05
1E-06
0,01 0,1 1 10
F(Hz)
219
Espectro de energía de S3-S5
S3 S5 ( libre) a 5 Hz
1
0,1
0,01
Sx'Sy'
0,001
S3 LIBRE
0,0001
S5 LIBRE
1E-05
1E-06
0,001 0,01 0,1 1 10
0,1
0,01
Sx'Sy'
0,001
S3 LIBRE
0,0001 S5 LIBRE
1E-05
1E-06
0,001 0,01 0,1 1 10
F(Hz)
220
Espectro de energía en S3-S5 ( libre) a 50 Hz
1
0,1
0,01
Sx'Sy'
0,001
S3 LIBRE
S5 LIBRE
0,0001
p=-5/3
1E-05
1E-06
0,001 0,01 0,1 1 10
F(Hz)
0,1
0,01
Sx'Sy'
0,001
S3 SUM
0,0001 S5 SUM
p=-5/3
1E-05
1E-06
0,001 0,01 0,1 1 10
F(Hz)
221
Espectro de energía en S3-S5 ( sum) a 25 Hz
1
0,1
0,01
Sx'Sy'
0,001
S3 SUM
0,0001 S5 SUM
p=-5/3
1E-05
1E-06
0,001 0,01 0,1 1 10
F(Hz)
0,1
0,01
Sx'Sy'
0,001
S3 SUM
S5 SUM
0,0001
p=-5/3
1E-05
1E-06
0,01 0,1 1 10
F(Hz)
222
Espectro de energía de S3-S5 (libre) a 5 Hz
1
0,1
0,01
0,001
S3 LIBRE
Sy'
0,0001 S5 LIBRE
p=-5/3
1E-05
1E-06
0,001 0,01 0,1 1 10
F(Hz)
0,1
0,01
Sx'
0,001
S3 LIBRE
0,0001 S5 LIBRE
p=-5/3
1E-05
1E-06
0,001 0,01 0,1 1 10
F(Hz)
223
Espectro de energía en S3-S5 (libre) a 50 Hz
1
0,1
0,01
Sy'
0,001
S3 LIBRE
S5 LIBRE
0,0001
p=-5/3
1E-05
1E-06
0,01 0,1 1 10
F(Hz)
0,1
0,01
Sy'
0,001
S3 SUM
0,0001 S5 SUM
p=-5/3
1E-05
1E-06
0,001 0,01 0,1 1 10
F(Hz)
224
Espectro de energía de S3-S5 (sum) a 25 Hz
1
0,1
0,01
Sy'
0,001
S3 SUM
S5 SUM
0,0001
p=-5/3
1E-05
1E-06
0,001 0,01 0,1 1 10
F(Hz)
1
Espectro de energía en S3-S5 (sum) a 50 Hz
0,1
0,01
0,001 S3 SUM
Sy'
S5 SUM
0,0001 p=-5/3
1E-05
1E-06
0,01 0,1 1 10
F(Hz)
225
ESPECTROS DE ENERGÍA EN UNA MISMA SECCIÓN, RL Y RS
0,1
0,01
Sx'
0,001 S3 LIBRE
S3 SUMERGIDO
0,0001 p=-5/3
1E-05
1E-06
0,001 0,01 0,1 1 10
F(Hz)
0,1
0,01
Sx'
0,001 S3 LIBRE
S3 SUMERGIDO
0,0001 p=-5/3
1E-05
1E-06
0,001 0,01 0,1 1 10
F(Hz)
226
Espectro de energía de las S3(libre)-S3( sum) a
50 Hz
0,1
0,01
Sy'
0,001 S3 LIBRE
S3 SUMERGIDO
0,0001 p=-5/3
1E-05
1E-06
0,01 0,1 1 10
F(Hz)
0,1
0,01
Sx'
0,001
S5 LIBRE
0,0001 S5 SUMERGIDO
p=-5/3
1E-05
1E-06
0,001 0,01 0,1 1 10
F(Hz)
227
Espectro de energía en las S5(libre)-S5( sum) a
25 Hz
1
0,1
0,01
Sx'
0,001
S5 LIBRE
0,0001 S5 SUMERGIDO
p=-5/3
1E-05
1E-06
0,001 0,01 0,1 1 10
F(Hz)
0,1
0,01
Sx'
S5 LIBRE
0,001 S5 SUMERGIDO
p=-5/3
0,0001
1E-05
1E-06
0,01 0,1 1 10
F(Hz)
228
Espectro de energía en las S3(libre)-S3( sum) a 5
Hz
1
0,1
Sx'Sy'
0,01
0,001
S3 LIBRE
0,0001 S3 SUMERGIDO
p=-5/3
1E-05
1E-06
0,001 0,01 0,1 1 10
F(Hz)
0,1
0,01 S3 LIBRE
Sx'Sy'
S3
0,001 SUMERGIDO
p=-5/3
0,0001
1E-05
1E-06
0,001 0,01 0,1 1 10
F(Hz)
229
Espectro de energía en S3(libre)-S3( sum) a 50
Hz
1
0,1
0,01
Sx'Sy'
0,001
S3 LIBRE
0,0001 S3 SUMERGIDO
p=-5/3
1E-05
1E-06
0,01 0,1 1 10
F(Hz)
0,1
0,01
Sx'Sy'
0,001
S3 LIBRE
S3 SUMERGIDO
0,0001
p=-5/3
1E-05
1E-06
0,001 0,01 0,1 1 10
F(Hz)
230
Espectro de energía en S5(libre)-S5( sum) a 25
Hz
1
0,1
0,01
Sx'Sy'
0,001
S3 LIBRE
0,0001 S3 SUMERGIDO
p=-5/3
1E-05
1E-06
0,001 0,01 0,1 1 10
F( Hz)
0,1
0,01
Sx'Sy'
0,001
S3 LIBRE
0,0001 S3 SUMERGIDO
p=-5/3
1E-05
1E-06
0,01 0,1 1 10
F( Hz)
231
Espectro de energía de S3(libre)-S3( sum) a 5 Hz
1
0,1
0,01
Sy'
0,001
S3 LIBRE
S3 SUMERGIDO
0,0001
p=-5/3
1E-05
1E-06
0,001 0,01 0,1 1 10
F(Hz)
0,1
0,01
Sy'
0,001
S3 LIBRE
0,0001 S3 SUMERGIDO
p=-5/3
1E-05
1E-06
0,001 0,01 0,1 1 10
F(Hz)
232
Espectro de energía en S3(libre)-S3( sum) a
50Hz
1
0,1
0,01
Sy'
0,001
S3 LIBRE
0,0001 S3 SUMERGIDO
p=-5/3
1E-05
1E-06
0,01 0,1 1 10
F(Hz)
0,1
0,01
Sy'
0,001
S5 LIBRE
S5 SUMERGIDO
0,0001
p=-5/3
1E-05
1E-06
0,001 0,01 0,1 1 10
F(Hz)
233
Espectro de energía de S5(libre)-S5( sum) a
25Hz
1
0,1
0,01
Sy'
0,001
S5 LIBRE
0,0001 S5 SUMERGIDO
p=-5/3
1E-05
1E-06
0,001 0,01 0,1 1 10
F(Hz)
0,1
0,01
Sy'
0,001
S5 LIBRE
0,0001 S5 SUMERGIDO
p=-5/3
1E-05
1E-06
0,01 0,1 1 10
F(Hz)
234
BIBLIOGRAFÍA
235
Doppler”. Proyecto de Fin de Carrera. Ingeniería Técnica de Obras Públicas-
Especialidad Hidrología. UPCT.
- NAUDASCHER, E. (2002). “Hidráulica de canales. Diseño de estructuras”. Ed.
Limusa.
- RAJARATNAM, N. (1965). “The Hydraulic Jump as a Wall Jet”. Proc. ASCE,
Hy., Sept.
- RANGA RAJU, K.G. (1998) “Flow Through Open Channels”. Ed. Tata
McGraw-Hill.
- SILVESTER, R. (1964). “Hydraulic Jump in All Shapes on Horizontal Channels”.
ASCE. Journal of the Hydraulic Division, vol. 90, no. Hy. 1, January, pp. 23-55.
- SONTEK ADEV FIELD (2001). “Acoustic Doppler Velocimeter”. Technical
Documentation.
- VEN TE CHOW (1998) “Hidráulica de canales abiertos” Ed. McGraw-Hill.
- VICENTE, J. (2008). “Caracterización de resaltos hidráulicos libres a partir de
velocidades instantáneas con equipo Doppler. Flujo aguas abajo de una
compuerta”. Proyecto Fin de Carrera. Ingeniería Técnica de Obras Públicas-
Especialidad Hidrología. UPCT.
- VISCHER, D.L. AND HAGER, W.H. (1995) “ Energy Dissipators” Hydraulic
Desig Considerations. Hydraulic Structures Design Normal Nº 9 IAHR.
236