Hidráulica Marítima - Teoria Lineal de Olas - Ejercicios
Hidráulica Marítima - Teoria Lineal de Olas - Ejercicios
Hidráulica Marítima - Teoria Lineal de Olas - Ejercicios
2π 2π
Siendo: σ=
T
k=
L
(
η( x , t) = a sin k x - σ t )
2 π 2 π
Onda 1: σ1 := Onda 2: σ2 :=
8s 10s
2 π 2 π
k1 := k2 :=
100m 80m
(
ηa( x, t) := a cos k1 x - σ1 t ) (
ηb ( x , t) := a cos k2 x + σ2 t )
De izquierda a derecha De derecha a izquierda
η1( x , 0s) 1
η1( x , 1s)
m
η1( x , 2s) 0
η1( x , 3s)
-1
-2
m
Se observa en el grafico que la progresion en el tiempo del mismo genera una onda cuasi-estacionaria,y esto
se debe principalmente a que las longitudes de ondas y sus periodos son similares, y son de sentidos
opuestos
La altura de la superficie libre a lo largo del tiempo para x = 50m
1
m
η1( 50m , t) 0
-1
-2
0 20 40 60 80
t
s
Si repetimos los cálculos para igual sentido de propagación
Onda 1: Onda 2:
(
ηc( x, t) := a cos k1 x - σ1 t ) (
ηd ( x , t) := a cos k2 x - σ2 t )
De izquierda a derecha De izquierda a derecha
2
Se observa en el grafico que la
η2( x , 0s) progresion en el tiempo del mismo
1 genera una onda cuasi-progresiva,y
η2( x , 1s) esto se debe principalmente a que
las longitudes de ondas y sus
m
-2
m
La altura de la superficie libre a lo largo del tiempo para x = 50m
1
m
η2( 50m , t) 0
-1
-2
0 20 40 60 80
t
s
P2) Sea una ola cuya longitud de onda en aguas profundas es L0 = 156m, una altura
H0 = 2m y una velocidad C0 = 15.6m/s, moviéndose en aguas abiertas con su cresta
paralela a sus contornos profundos.
Despreciando efectos de reflexión desde la playa.
Hallar:
a) Relación entre la altura de la ola en cualquier profundidad de agua y la altura en agua
profunda.
b) La altura de la ola cuando la profundidad es de 3m.
c) La tasa de energía por unidad de ancho de cresta transportada hacia la costa y la
energía total por unidad de ancho llevada a la costa en 1 hora.
a)
Los parámetros de la ola en aguas profundas son:
m
L0 := 156m H0 := 2m c0 := 15.6
s
L( h ) := root tanh 2 π
h h h
- , ξ
L0 ξ ξ
2 π
k( h) :=
L( h )
2 k0 h 0
(
tanh k0 h0 1 +)
H1 a1
( )
sinh 2k0 h 0
= = = Ks
H0 a0 2 k1 h 1
tanh ( k1 h1 ) 1 +
sinh ( 2k1 h 1 )
2 k h
Sabiendo que para aguas profundas k h >> 1, entonces lim = 0 y tanh ( k h ) = 1 , lo que resulta
k h ∞ senh( 2 k h )
en:
1
Ks( h) := para cualquier profundidad (teniendo en cuenta que
2 k( h ) h
tanh ( k( h ) h ) 1 +
venimos de aguas profundas)
sinh ( 2 k( h) h)
H( h ) := H0 Ks( h )
Graficando la altura de ola en función de la profundidad:
3
m
H ( h) 2
0
0 50 100 150
h
b)
Cuando la profundidad es de 3m, la altura de la ola es: H( 3m) = 2.5 m
c)
1 kgf H0
Siendo: cg0 := c0 γ := 1000 a :=
2 3 2
m
8 J
Et := Pm 3600s Et = 1.377 10
m
ξ := 1m
P3) Sea el gradiente de profundidad mostrado en la figura del enunciado.
Para una ola incidente de altura H1 = 1m y período T = 10seg, la cual incide con ángulo de θ1 = 30º,
determinar:
a) Qué tipo de aproximación es aplicable en cada zona para el cálculo de la celeridad.
b) Cómo varía el ángulo de incidencia al cambiar de zona (mediante la ley de snell) ilustrando con
esquema.
c) La altura de la ola en la zona III
Primero definimos los parametros para aguas profundas, siendo T := 10s , θ1 := 30deg , H1 := 1m :
2
g T
L0 := , entonces L0 = 156m
2π
Zona 1:
h 1 := 100m
h1 h1 h1 2π 1
L1 := root - tanh 2 π , ξ = 155.979 m por lo que k1 := = 0.04
L0 ξ ξ
L1 m
h1
entonces = 0.641 > 0.5, por lo que estamos en aguas profundas
L1
Zona 2:
h 2 := 10m
h2 h2 h2 2π 1
L2 := root - tanh 2 π , ξ = 92.356 m por lo que k2 := = 0.068
L0 ξ ξ
L2 m
h2
entonces = 0.108 > 0.05, por lo que estamos en aguas intermedias
L2
Zona 3:
h 3 := 2m
h3 h3 h3 2π 1
L3 := root - tanh 2 π , ξ = 43.692 m por lo que k3 := = 0.144
L0 ξ ξ
L3 m
h3
entonces = 0.046 < 0.05, por lo que estamos en aguas poco profundas
L3
Por ley de Snell, sabemos que se cumple la relacion
L2 ( ) entre los angulos incidentes de ola, por lo que
sin θ2
=
L1 sin( θ1 )
L2
podemos despejar el angulo incidente para la proxima onda como: θ2 := asin ( ) = 0.300557
sin θ1
L1
L3
De manera analoga para la zona 3: θ3 := asin sin θ2 ( ) = 0.140519
L2
a)
La aproximación que es aplicable en cada zona para el cálculo de la celeridad es:
g 1
Zona 1: aguas profundas, entonces c1 = y n 1 :=
k 2
g
Zona 2: aguas intermedias, entonces c2 = tanh ( k h)
k
b)
c)
Se puede demostrar que, para distintos cambios de profundidad se cumple la relacion:
HN n 1 L1 ( )
cos θ1
=
H1 n N LN cos( θN)
n 1 L1 ( )
cos θ1
H3 := H = 1.249 m
n 3 L3 cos( θ3 ) 1
ξ := 1m
P4) Calcular la refracción y el bajío de una ola de H1 = 2m, T = 6seg y con un ángulo de incidencia de θ1
= 20º en las condiciones del croquis. Indicar los valores de KS2, KR2 y H2
h1 h1 h1 2 π 1
L1 := root - tanh 2 π , ξ = 55.032 m k1 := = 0.114
L0 ξ ξ
L1 m
h2 h2 h2 L = 2 π 1
L2 := root - tanh 2 π , ξ = 38.082
2 m k2 := = 0.165
L0 ξ ξ
L2 m
1 2 k1 h 1 1 2 k2 h 2
n 1 := 1 + = 0.547 n 2 := 1 + = 0.829
2 sinh( 2 k1 h1 ) 2 sinh( 2 k2 h2 )
Las celeridades:
L1 m L2 m
c1 := = 9.172 c2 := = 6.347
T s T s
n 1 c1
ks21 := ks21 = 0.977
n 2 c2
n1 sin( θ1 )
θ2 := asin θ2 = 13.053 deg
n2
El coeficiente de refracción es:
( )
cos θ1
kr21 := kr21 = 0.982
cos( θ2 )
Para hallar las coordenadas del punto P en metros, debemos conocer la longitud de la onda incidente.
En aguas profundas
g 2
L0 := Ti L0 = 99.89 m
2 π
hi hi hi
Li := root - tanh 2 π , ξ = 53.071 m
L0 ξ ξ
HP = Hi K = 0.256m
P6.1) Calcular la altura de ola y direccion de incidencia a una profundidad h = 4m en
forma grafica y analitica, conociendo H a una profundidad de h = 25m. Dibujar el
perfil transversal de la playa.
Las batimetricas son rectas, paralelas y separadas 150m entre si.
Primero definimos los parametros para aguas profundas, suponiendo T0 := 10s , θ1 := 30deg , H1 := 2m :
2
g T0
L0 := entonces L0 = 156m
2π
h a5 := 8m h a6 := 6m h a7 := 4m h a8 := 2m h a9 := 0m
ξ := 1m
Zona 1:
h a1 ha1 h a1 2π 1
La1 := root - tanh 2 π , ξ = 130.353 m por lo que ka1 := = 0.048
L0 ξ ξ
La1 m
ha1
entonces = 0.192 < 0.5, por lo que estamos en aguas intermedias
La1
Zona 2:
h a2 ha2 h a2
La2 := root - tanh 2 π , ξ = 121.21 m
L0 ξ ξ
Zona 3:
h a3 ha3 h a3
La3 := root - tanh 2 π , ξ = 109.027 m
L0 ξ ξ
Zona 4:
h a4 ha4 h a4
La4 := root - tanh 2 π , ξ = 92.356 m
L0 ξ ξ
Zona 5:
h a5 ha5 h a5
La5 := root - tanh 2 π , ξ = 83.801 m
L0 ξ ξ
Zona 6:
h a6 ha6 h a6
La6 := root - tanh 2 π , ξ = 73.609 m
L0 ξ ξ
Zona 7:
h a7 ha7 h a7 2π 1
La7 := root - tanh 2 π , ξ = 60.946 m por lo que ka7 := = 0.103
L0 ξ ξ
La7 m
ha7
entonces = 0.066 < 0.5, por lo que estamos en aguas intermedias
La7
La3
Zona 3: θa3 := asin ( ) = 0.43146
sin θa2
La2
La4
Zona 4: θa4 := asin ( ) = 0.362114
sin θa3
La3
La5
Zona 5: θa5 := asin ( ) = 0.327249
sin θa4
La4
La6
Zona 6: θa6 := asin ( ) = 0.286238
sin θa5
La5
La7
Zona 7: θa7 := asin
La6
( ) = 0.235957
sin θa6 θa7 = 13.5deg
1 2 ka1 h a1
n a1 := 1 + = 0.718
2 sinh ( 2 ka1 h a1)
1 2 ka7 h a7
n a7 := 1 + = 0.948
2 sinh ( 2 ka7 h a7)
Se puede demostrar que, para distintos cambios de profundidad se cumple la relacion:
HN
=
n 1 L1
( )
cos θ1
H1 n N LN cos( θN)
Ha7 :=
na1 La1
( ) H = 2.403 m
cos θ1
na7 La7 cos( θa7) 1
b) h b1 := 25m h b2 := 15m h b3 := 10m h b4 := 6m
h b5 := 10m h b6 := 8m h b7 := 6m h b8 := 4m h b9 := 0m
ξ := 1m
Zona 1:
h b1 h b1 hb1 2π 1
Lb1 := root - tanh 2 π , ξ = 130.353 m por lo que kb1 := = 0.048
L0 ξ ξ
Lb1 m
hb1
entonces = 0.192 < 0.5, por lo que estamos en aguas intermedias
Lb1
Zona 2:
h b2 h b2 hb2
Lb2 := root - tanh 2 π , ξ = 109.027 m
L0 ξ ξ
Zona 3:
h b3 h b3 hb3
Lb3 := root - tanh 2 π , ξ = 92.356 m
L0 ξ ξ
Zona 4:
h b4 h b4 hb4
Lb4 := root - tanh 2 π , ξ = 73.609 m
L0 ξ ξ
Zona 5:
h b5 h b5 hb5
Lb5 := root - tanh 2 π , ξ = 92.356 m
L0 ξ ξ
Zona 6:
h b6 h b6 hb6
Lb6 := root - tanh 2 π , ξ = 83.801 m
L0 ξ ξ
Zona 7:
h b7 h b7 hb7
Lb7 := root - tanh 2 π , ξ = 73.609 m
L0 ξ ξ
Zona 8:
h b8 h b8 hb8 2π 1
Lb8 := root - tanh 2 π , ξ = 60.946 m por lo que kb8 := = 0.103
L0 ξ ξ
Lb8 m
hb8
entonces = 0.066 < 0.5, por lo que estamos en aguas intermedias
Lb8
Por ley de Snell, sabemos que se cumple la relacion
L2 ( ) entre los angulos incidentes de ola, por lo que
sin θ2
=
L1 sin( θ1 )
Lb2
podemos despejar el angulo incidente para la proxima onda como: θb2 := asin ( ) = 0.43146
sin θ1
Lb1
De manera analoga:
Lb3
Zona 3: θb3 := asin ( ) = 0.362114
sin θb2
Lb2
Lb4
Zona 4: θb4 := asin ( ) = 0.286238
sin θb3
Lb3
Lb5
Zona 5: θb5 := asin ( ) = 0.362114
sin θb4
Lb4
Lb6
Zona 6: θb6 := asin ( ) = 0.327249
sin θb5
Lb5
Lb7
Zona 7: θb7 := asin ( ) = 0.286238
sin θb6
Lb6
Lb8
Zona 8: θb8 := asin ( ) = 0.235957
sin θb7 θb8 = 13.5deg
Lb7
Los factores n son:
1 2 kb1 h b1
n b1 := 1 + = 0.718
2 sinh ( 2 kb1 hb1)
1 2 kb8 h b8
n b8 := 1 + = 0.948
2 sinh ( 2 kb8 hb8)
Se puede demostrar que, para distintos cambios de profundidad se cumple la relacion:
HN n 1 L1 ( )
cos θ1
=
H1 n N LN cos( θN)
n b1 Lb1 ( )
cos θ1
Hb8 := H = 2.403 m
n b8 Lb8 cos( θb8) 1
c) h c1 := 25m h c2 := 15m h c3 := 12m h c4 := 11m
h c5 := 10m h c6 := 8m h c7 := 6m h c8 := 4m h c9 := 0m
ξ := 1m
Zona 1:
h c1 hc1 h c1 2π 1
Lc1 := root - tanh 2 π , ξ = 130.353 m por lo que kc1 := = 0.048
L0 ξ ξ
Lc1 m
hc1
entonces = 0.192 < 0.5, por lo que estamos en aguas intermedias
Lc1
Zona 2:
h c2 hc2 h c2
Lc2 := root - tanh 2 π , ξ = 109.027 m
L0 ξ ξ
Zona 3:
h c3 hc3 h c3
Lc3 := root - tanh 2 π , ξ = 99.707 m
L0 ξ ξ
Zona 4:
h c4 hc4 h c4
Lc4 := root - tanh 2 π , ξ = 96.163 m
L0 ξ ξ
Zona 5:
h c5 hc5 h c5
Lc5 := root - tanh 2 π , ξ = 92.356 m
L0 ξ ξ
Zona 6:
h c6 hc6 h c6
Lc6 := root - tanh 2 π , ξ = 83.801 m
L0 ξ ξ
Zona 7:
h c7 hc7 h c7
Lc7 := root - tanh 2 π , ξ = 73.609 m
L0 ξ ξ
Zona 8:
h c8 hc8 h c8 2π 1
Lc8 := root - tanh 2 π , ξ = 60.946 m por lo que kc8 := = 0.103
L0 ξ ξ
Lc8 m
hc8
entonces = 0.066 < 0.5, por lo que estamos en aguas intermedias
Lc8
Lc3
Zona 3: θc3 := asin ( ) = 0.392446
sin θc2
Lc2
Lc4
Zona 4: θc4 := asin ( ) = 0.377776
sin θc3
Lc3
Lc5
Zona 5: θc5 := asin ( ) = 0.362114
sin θc4
Lc4
Lc6
Zona 6: θc6 := asin ( ) = 0.327249
sin θc5
Lc5
Lc7
Zona 7: θc7 := asin ( ) = 0.286238
sin θc6
Lc6
Lc8
Zona 8: θc8 := asin ( ) = 0.235957
sin θc7 θc8 = 13.5deg
Lc7
Los factores n son:
1 2 kc1 h c1
n c1 := 1 + = 0.718
2 sinh ( 2 kc1 h c1)
1 2 kc8 h c8
n c8 := 1 + = 0.948
2 sinh ( 2 kc8 h c8)
Se puede demostrar que, para distintos cambios de profundidad se cumple la relacion:
HN
=
n 1 L1
( )
cos θ1
H1 n N LN cos( θN)
Hc8 :=
nc1 Lc1
( ) H = 2.403 m
cos θ1
nc8 Lc8 cos( θc8) 1
Los perfiles transversales son, respectivamente:
a)
b)
c)
A modo de conclusion podemos decir que para todos los casos, sin importar el perfil transversal, se obtuvo la
misma direccion y altura para el punto de profundidad de 4m
P7) Dadas olas con altura Hi = 3m y período T = 7 seg con incidencia normal hacia un
montículo de pedregullo (rubble-mound breakwater) con pendiente 1/2.
Hallar el coeficiente de reflexión utilizando el gráfico Kr vs x (número de Iribarren).
Hi := 3m T0 := 7s
g 2
L0 := T0 L0 = 76.478 m
2 π
Para una escollera de pedregullo con pendiente 1/2 el parámetro de Iribarren es:
θ := atan
1
θ = 26.6 º
2
tan( θ)
ξ := ξ = 2.5
Hi
L0
Kr := 0.3
kg
Datos: PA := 8300 h := 9m T := 15s z := -8 m B := 150m X := 450m Y := 150m
2
m
2
g T
1) Para aguas profundas sabemos que: L0 := = 351.175 m
2π
tanh 2 π
h h h
=
L0 L L
L := root tanh 2 π
h h h
- , λ L = 137.129 m
L0 λ λ
2π 1
k := k = 0.046
L m
Debemos entrar en los graficos para obtener el coeficiente de difraccion, para ello expresamos las
coordenadas en funcion de la longitud de onda:
B X Y
= 1.094 = 3.282 = 1.094
L L L
Calculamos la altura de ola en el punto A, con el dato de la presion dada y sabiendo que PA = γ η Kz - z ( )
Suponemos que la presion dada corresponde a la maxima presion en el punto A, por lo que η = Ha/2
Definimos: cosh[ k ( h + z) ] kg
Kz := = 0.922 γ := 1000
cosh( k h) 3
m
HA1 PA 1 PA 1
η= = + z HA1 := 2 + z = 0.651 m
2 γ Kz γ
Kz
Suponiendo que la ola incide de manera prependicular al plano de los rompeolas, ayudandonos con el
siguiente grafico obtenemos el coeficiente de difracción
Adoptamos Kd := 0.1
Finalmente, teniendo la altura de la ola difractada, despejamos la altura de la ola incidente mediante:
HA HA1
Kd = HI := = 6.511 m
HI Kd
2) Suponiendo que la direccion de incidencia es la indicada por el enunciado, que la altura de la ola
incidente es la misma del punto anterior, buscamos el coeficiente de difraccion con los siguientes datos
geometricos:
( )
0.5
θ := atan
2 2 Y r
r := Y + X = 18.435 deg = 3.459
X L
HA2 := Kd HI = 0.651 m
λ := 1m
P9) Sea una dársena rectangular de aproximadamente 5m de ancho por 400 de largo,
habiendo en la misma una profundidad uniforme de 10m.
a) Calcular el período para los primeros cuatro modos normales, hallar las líneas
nodales e interpretar gráficamente los resultados.
b) Si un barco que genera olas con un período de 1/2seg. Usando la aproximación de
dársena angosta, determinar cuál es el modo dominante al entrar el barco, justificando
el resultado.
Datos: a := 400m b := 5m h := 10m
a) Dado que uno de los lados es mucho mayor que el otro, aproximaremos el periodo mediante la formula:
2 a
T( n ) :=
n g h
Entonces, para los 4 modos normales:
Ademas, sabemos que por la formula de Morian, la superficie libre en este caso esta dada por la ecuacion:
n π x
cos
H
η= cos( σ t)
2 a
n π x n π x
la cual se anula para cualquier t si cos
π
= 0 , o lo que es mismo = ( 2l + 1 )
a a 2
Para n := 1
π x π
= ( 2l + 1)
a 2
Por lo que si l = 0, la solucion es que x = 200m y es unica, ya que con cualquier otro l obtendriamos un valor de x
fuera de la darsena
Para n := 2
2π x π
= ( 2l + 1 )
a 2
Para l = 0, la solucion es x = 100m, para l = 1 la solucion es x = 300m , y son las unicas soluciones posibles
Para n := 3
3π x π
= ( 2l + 1 )
a 2
Para n := 4
4π x π
= ( 2l + 1 )
a 2
Para l = 0, la solucion es x = 50m, para l = 1 la solucion es x = 150m, para l = 2 la solucion es x = 250m, para l = 3
la solucion es x = 350m, y son las unicas soluciones posibles