TEMA5 - Grado ICT Refracción, Difracción y Rotura 2012

TEMA V
REFRACCIN, DIFRACCIN Y ROTURA

Tema 5. Los niveles del mar y las variaciones del nivel medio. Convergencia con el oleaje
5.1
Niveles del mar y nivel medio del mar
5.2
Marea astronmica y otras fluctuaciones peridicas
5.3
Marea meteorolgica y otras variaciones de origen climtico
5.4
Propagacin y transformacin del oleaje: refraccin, difraccin y rotura
5.5
Efectos de la rotura del oleaje sobre las costas: erosin y transporte de

sedimentos
Vicente Negro Valdecantos

Doctor Ingeniero de Caminos, Canales y Puertos
Profesor Titular de Universidad
Curso Acadmico 2012 2013
REFRACCIN. Razonamiento terico del problema

La determinacin del oleaje es absolutamente imprescindible en la Ingeniera Martima,
Portuaria y Costera. Los datos de partida son generalmente obtenidos en aguas profundas,
d/L > 1/2, y es necesaria la transformacin del fenmeno ondulatorio hasta el rea de inters,
en zonas de relativa y escasa lmina de agua, zonas de transicin, profundidades reducidas e
incluso en procesos de descomposicin ondulatoria con fenmenos de rotura.
La modificacin de la onda con los procesos de refraccin o expansin frontal de la energa

del oleaje, difraccin o expansin lateral, shoaling o asomeramiento, disipacin, friccin,
reflexin y otros fenmenos que pueden afectar a la onda son elementos bsicos en el
anlisis. La onda es un ente matemtico, la ola es un ente fsico.
Los modelos de propagacin se basan en mltiples teoras.
Matemtica aplicada a la onda con su ecuacin de Laplace basada en las hiptesis de

incompresibilidad del fluido ( = constante), e irrotacionalidad del flujo (rotacional del
campo = 0), en la ecuacin de continuidad como operador diferencial divergencia y
la velocidad como funcin gradiente del potencial.
div u = 0, la divergencia de un campo vectorial es un escalar

u = - grad , el gradiente de un campo escalar es vectorial
div (grad ) = = 0 Ecuacin de Laplace, por operador de campo compuesto
Soluciones basadas en la ecuacin de pendiente suave de Berkoff, 1972 o "mild slope

equation".
Resoluciones de manera elptica
Resoluciones en aproximaciones parablicas
Modelos basados en las ecuaciones de conservacin
Refraccin, difraccin y rotura

Grado en Ingeniera Civil y Territorial
Obras Martimas
Conservacin de la accin de onda
Irrotacionalidad del nmero de onda
Estos elementos o herramientas tan sofisticadas han servido para, en cierta medida,
complementar a los modelos clsicos grficos basados en dibujos de evolucin mediante
semiavances, Iribarren, o la teora de rayos, Snell.
Iribarren
H1
=
H0
B0 2
; H c g = cte ; H2 c g cos = cte
B1
Altura de ola, m
Ancho del cuadriltero de avance, m
Subndice de profundidad indefinida, -
Subndice de profundidad considerada, -
Snell
sen 1
c1
sen 2
c2
; H12 cos 1 = H2 2 cos 2
Oblicuidad, ngulo que forma el frente con la batimetra,
Celeridad del movimiento, m/s
Hi
Altura de ola, m
1, 2
Indicadores de referencia, punto 1 y punto2
sen
= cte , H2 c g = cte
c
Los subndices representan la misma significacin que la teora clsica del Profesor Iribarren.

Obras Martimas
Han sido mltiples los intentos de obtener bacos de determinacin del coeficiente de
refraccin (Kr) y asomeramiento (Ks) y combinado (Kr x Ks) en funcin de criterios de batimetra
rectilnea, sensiblemente paralela a la lnea de costa, ms o menos homognea que permitan
definir,
K=f(
d
, 0 )
g . T2
obteniendo resultados como el proporcionado por el U.S. Army Coastal Engineering Research
Center 1.977, o baco del Shore Protection Manual.
Con hiptesis semejantes y en funcin de parmetros espectrales direccionales (Smax),

dependientes de la pulsacin angular y la velocidad del viento, y para oleajes tipo Sea y Swell
con pequeo y gran "decay", decaimiento, como representantes del peraltamiento de la ola,
Goda 1.985, en fondo paralelo - regular define el coeficiente de refraccin para un campo
irregular direccional de oleaje tomando en abscisas h/L0 y en ordenadas 0.
Mitsoyashu con hiptesis semejantes y en funcin de d/L > 0.50 y o define la variacin
angular para obtener el ngulo p y con ello, el coeficiente de refraccin correspondiente.
Kr =
c0 =

Obras Martimas
cos 0
cos i
gT
; celeridad indefinida
2
cp =
cp
gT
2d
th
2
Lp
celeridad a profundidad "p", m/s
sen p
cp
sen 0
c0
= cte
El hecho cierto es que la matemtica aplicada a la teora de ondas y los avances en las
tcnicas de resolucin de ecuaciones diferenciales, bien empleando elementos finitos o
diferencias finitas, sirvi para desarrollar el campo de la propagacin de los fenmenos
ondulatorios.
La clsica ecuacin de Laplace,

2
2
2

+
+
= 0
2
x2
z2
y
= 0 = diver ( grad ) ; u =
con las hiptesis de incomprensibilidad e irrotacionalidad y discretizacin,
(x , y , z , t ) = A ( x , y ) w ( z ) e ( - i * * t )
(z )=
ch k ( d + z )
2
;k =
ch k d
L
y la ecuacin de Helmholtz,
2
2
A
A
+
+ k2 A = 0
2
2
x
y

Obras Martimas
Esta ecuacin fue empleado en expansin lateral y fue bsico para estudiar la propagacin de
ondas sobre fondos suaves - variables, definiendo un potencial.
(x, y, z, t)
*
(x, y)
* (x, y) x w(z) x exp(- iwt)
A(x, y) x exp {is (x, y)}
A(x, y)
Amplitud
s(x, y)
Fase
Transformando la ecuacin de Laplace de 3D a 2D, integrando,

amplitud
-h
( z ) dz
La descripcin completa de los fenmenos de propagacin de ondas monocromticas lineales

sobre un fondo cualquiera, batimetra arbitraria fue descrita y formulada por Berkoff, 1972, con
la nica restriccin de taludes suaves, expresando su ecuacin de "pendiente suave", "mild
slope equation" mediante un potencial y una expresin matemtica de naturaleza elptica.
c c g + 2
cg
=0
c
siendo:
Celeridad de la onda, m/s
cg
Celeridad del grupo, m/s
Frecuencia, s-1

Obras Martimas

2
= g k th k d
-1
Nmero de onda, m
k=
2
L
Longitud de onda en profundidad reducida, m
L0
Longitud de onda en profundidad indefinida, m
Profundidad, m
Funcin potencial, -
La solucin numrica emple dos tendencias claramente diferenciadas, el mtodo de los

elementos finitos y los esquemas en diferencias finitas.
MEF
MDF
Berkoff,
1.972
Houston,
1.981
William, 1.980
Holmes, 1.980
Ambos modelos presentan un largo tiempo de clculo computacional con los inconvenientes
que esto conlleva.
La resolucin del problema elptico requera adems de un tiempo de computacin muy

elevado, condiciones de contorno en todo el dominio, as como reflexiones no deseadas en las
condiciones de absorcin, ante incidencia oblicua. Las grandes mallas conducen a enormes
matrices y como consecuencia, procesos laboriosos de clculo.

Obras Martimas
El conocimiento de la condicin de contorno en la costa casi siempre desconocida, obliga a

pensar en otros planteamientos de clculo que obviarn la misma y simplificarn notablemente
la resolucin numrica del problema.
Para simplificar el modelo elptico y descomponiendo la funcin potencial en una componente

de reflexin que se dispersa y en una difraccin, donde en direccin de avance es ms
reducido que en expansin lateral, se obtiene una ecuacin representativa del fenmeno de
propagacin, transformando el problema de elptico a parablico, empleando conceptos de
funcin potencial complejo.
Est simplificacin origina los modelos llamados "Parablicos" que funcionan en esquemas en
diferencias finitas; destacando:
Berkoff, Booij y Radder, 1.982

Radder, 1.982 (Ecuacin parablica de Radder)
Mei, 1.983
Kirby, 1.986
Equipo de la Universidad de Delaware (Mei, Kirby y Dalrymple) 1.988
REF - DIF
Vincent, 1.988
J.M. Grassa, REF - DIF, 1.988 90. Win Waves, 2.000 y siguientes
SMC, Sistema de Modelado Costero, Universidad de Cantabria, 2.001
La principal desventaja del modelo parablico se centra en la malla o celda de integracin,

sensiblemente paralela a la direccin de propagacin, lo que proporciona soluciones errneas
para batimetras irregulares, complejas, con alternancia de bajos, y giros pronunciados. La
ecuacin representativa de los fenmenos de propagacin resulta:
= i K
( C Cg K ) +
( C Cg
)
x
2 K C Cg x
2
y
C
g

Obras Martimas
donde "x" define ahora la direccin principal de propagacin y haciendo la hiptesis de

descomposicin del potencial en una componente dispersiva hacia adelante y una reflejada,
despreciando la componente de reflexin y suponiendo que los efectos de difraccin en la
direccin de propagacin son mucho menores que los producidos en la dimensin
perpendicular a la de avance, tal como se expona brevemente con anterioridad.
Tal como queda demostrado en la figura adjunta, para onda plana y pequeo valor del ngulo
, el error cometido es despreciable, para ngulos mayores de 45 los errores son ms
grandes, debiendo delimitar el campo de aplicacin de la aproximacin parablica.
Los modelos tipo Booij y Radder presentan una malla en la direccin de propagacin con eje X
en la direccin del oleaje y el eje Y en perpendicular antihorario. Los casos de oblicuidad
requieren la reorientacin de la malla de forma que el modelo Radder garantiza propagaciones
de onda 45; Booij maneja rangos de 60; debido al tratamiento en particiones de la matriz
elptica y el nmero de trminos en la aproximacin de las derivadas transversales. Kirby con
una aproximacin en minimax alcanz rangos de validez de 70.
As mismo, no es igual emplear estos esquemas en puertos que en playas, donde la

componente reflexin es de gran importancia en el primer caso y despreciable en el segundo,
de manera que puedan despreciarse frente al oleaje incidente.
Los modelos numricos de propagacin para ondas monocromticas lineales funcionan fuera
de rotura entre gamas d/L = 0.20 a d/L < 1/25, debiendo cambiar la funcin potencial para
caso de H/d H/L, fenmeno de rotura por fondo o peralte.
Los otros mtodos de clculo responden:
1.
Clsicos. Teora de Rayos

Obras Martimas
2.
Flujo de energa.
Teoras de Boussinesq, empleando la integracin de las ecuaciones de cantidad de

movimiento en el eje x, eje y, conservacin de masa, con una aproximacin en la componente
vertical dada por las ecuaciones de Boussinesq de aproximacin cuasilineal.
El modelo, aqu empleado, est basado en el estudio de Kirby y Dalrymple, Universidad de

Delaware, con aproximacin no lineal que incluye la celeridad de una onda en cada uno de los
dominios de existencia en funcin del modelo terico.
El clsico modelo REF - DIF incluye la teora de Dally y Dean para el proceso de rotura y post rotura asumiendo que tras sta existe una altura de ola estable proporcional a la profundidad y
que la tasa de disipacin de la energa en la zona de rompientes es proporcional a la
diferencia entre el flujo de energa real y el flujo de energa estable.
REF - DIF incluye tambin efectos de capas lmites laminares de superficie y fondo, capa
lmite de fondo turbulenta y propagacin sobre lechos porosos, preferentemente arena.
La tcnica de resolucin mediante el mtodo de Crank - Nicholson mediante esquema en

diferencias finitas implcito en dobles pasadas o barridos permiten progresar en la
determinacin de la respuesta numrica al problema.
El modelo parablico, de fondo de pendiente suave, restringido en el parmetro de Ursell a

valores inferiores a 40 y con direccin predominante del oleaje, permite obtener resultados
bastante aproximados al proceso refraccin - difraccin y a la combinacin del fenmeno
oleaje - corriente.
En dominios relativamente reducidos y efecto profundidad bastante importante, donde la

solucin terica es la onda solitaria y la onda cnoidal, las ecuaciones que modelan la
propagacin son las de Boussinesq, con soluciones numricas en 2D y fondo irregular.
Obras Martimas
10
Modelos actuales tipo Liu o Kirby permiten incluir efectos rotura, friccin e incidencia oblicua.
Los modelos de teora de rayos estn basados en las sencillas ecuaciones descritas
anteriormente y se basan en hiptesis muy sencillas de fsica general. stas son:
Primera Hiptesis: T, el perodo, es constante en la propagacin
sen 0
c T L0
c
= 0= 0
=
seni
ci T
Li
ci
sabiendo las celeridades se puede obtener,
1
2d
c0
=
= coth
= Ki
2d
L
ci
th
L
Por consecuencia,
seni =
sen 0
; Kr =
Ki
cos 0
sen 0
2
Ki
2
1-
Segunda Hiptesis: Flujo de energa entre dos ortogonales es constante
E0 =
1
g H0 2 cos 0
8
Ei =
1
g Hi2 cosi
8
Por tanto,

Obras Martimas
11
2
2
E0 = Ei ; H0 cos 0 = Hi cos i
Hi
=
H0
cos 0
= kr
cos i
Tercera Hiptesis: La celeridad solo depende de la profundidad, c = f(h)

Cuarta Hiptesis: El avance es ortogonal en cualquier profundidad
Quinta Hiptesis: La batimetra es rectilnea y paralela
Con estos criterios, se obtiene el coeficiente de refraccin por teora clsica de ortogonales.
Esta situacin puede complementarse con el concepto de asomeramiento o shoaling,
obtenido mediante el modelo clsico de Blunt siguiente:
Ks =
CG0
=
CG
1
4 h
2 h
L
th
1 +
4 h
L
sh
L
obteniendo el tradicional y clsico Kr x Ks mediante esta teora. Estos modelos siempre se

fundamentan en los principios siguientes:
El perodo es constante en la propagacin
El flujo de energa entre dos ortogonales se mantiene constante
La direccin de avance de la onda es perpendicular en cada momento de avance
La celeridad solamente depende de la profundidad

Obras Martimas
12
Se desprecian efectos corrientes, reflexiones, vientos,...
La batimetra es rectilnea, paralela, en la misma configuracin que la lnea de costa
Los cambios en fondo se consideran graduales, no existiendo discontinuidades
La refraccin, como consecuencia, es la expansin frontal de la energa del oleaje por efecto
fondo (kR), pudiendo presentar una componente de rayo (kr) y un efecto grupo (ks). Si
suponemos que el fondo es homogneo, todo el tren de ondas avanza con la misma celeridad
y se modifica en la propagacin de manera semejante, concepto que se define como
asomeramiento o shoaling (ks).
La existencia de irregularidades, variaciones de profundidad, entre otros, condiciona la

refraccin (kr).
El concepto global kR resulta como kr x ks. La refraccin o expansin frontal es un fenmeno

de muchas longitudes de onda; al contrario, que la difraccin o expansin lateral, que es un
fenmeno de muy pocas longitudes de onda.
Merece destacar tambin en los modelos denominados clsicos, los conceptos de Iribarren
fundamentados en los siguientes aspectos:
1.-
Definicin del proceso de refraccin, d/L < 1/2.
2.-
Definicin del cuadriltero de avance, nL0/2, nL0/2.
3.-
Seleccin de "n", recomendando valores entre 2 y 5
4.-
Definicin del avance en el proceso desde profundidad indefinida a reducida.

Obras Martimas
13
S0 =
5.-
n L0
2
Conociendo la ecuacin de onda y el avance, se define:
L =
2 S
2 d
g T2
; L = L0 th
; L0 =
2
n
L
S =
n L0
n d
th
2
S
n L0
2 n d
S
=
th
S
2
4
2
6.-
Esta situacin a perodo constante permite la propagacin y obtencin de la relacin

siguiente:
Hi
=
H0
B0
Bi
Estos esquemas se encuentran basados en la constancia del perodo y en la conservacin del

flujo de energa entre dos ortogonales.
Esta clsica descripcin del mtodo de Iribarren se observa en las figuras adjuntas, donde se
determina la curva de avances, en nuestra caso de semiavances, y en segundo lugar como se
aplica en tramos sucesivos de viaje desde profundidades indefinidas a reducidas pasando por
transicin del frente.
Debe recordarse que la teora de ondas responde a entes matemticos que transportan
energa, mientras que la solas son entes fsicos que transportan masa, diferencia fundamental
en los dos conceptos que se manejan en este captulo de modificaciones de las ondas. Los
grficos siguientes son suficientemente explicativos, adjuntando el clsico de Goda de Kr y
Shutto de Ks.
Obras Martimas
14
Grfico de Goda de refraccin (1985). Grfico de Shuto de shoaling - asomeramiento (1973)

Obras Martimas
15
Esquema general de Iribarren de expansin frontal de energa del oleaje

Obras Martimas
16
Esquema general de la propagacin de oleaje. Teora de ortogonales
Diferencias entre los modelos elpticos y parablicos de propagacin de oleaje

Obras Martimas
17
DIFRACCIN
Aunque resulta complejo separarlo del fenmeno anterior (refraccin y asomeramiento), la
difraccin es la expansin lateral de la energa del oleaje cuando un tren de ondas incide
sobre un obstculo emergido o sumergido a una determinada profundidad.
Al encontrarse con el dique se modifica la configuracin del tren de ondas y se crea una zona I
de Incidencia y reflexin. Esta zona es la de mxima agitacin. La onda contina sin estar
perturbada en su zona II o de Incidencia. Cuando la onda disminuye su energa por el efecto
obstculo, estamos en la zona III o de Alimentacin. La zona IV de expansin est formada
por la aparicin verdadera de la expansin lateral de energa cedida por la zona de
alimentacin. Finalmente, aparece una zona abrigada o V de calma. La planta de la difraccin
queda planteada en la figura adjunta.
Esquema de difraccin o expansin lateral de energa. Iribarren, 1.954

Obras Martimas
18
Para su comprensin terica, el fluido se considera perfecto y cumple la ecuacin de

continuidad. No hay cesin de energa entre la zona I y la II. El morro del dique es el emisor de
ondas, cediendo la energa a lo largo de las lneas de fase. La celeridad transversal coincide
con la de las ondas longitudinales, debido a la constancia del perodo y al anlisis de la
profundidad en cada celda o paraleleppedo de avance. Tanto el mtodo del Profesor Iribarren
como de Wiegell son de indudable valor prctico, aplicable a un tren monocromtico que en su
avance se encuentra un dique.
Forma en alzado de la expansin lateral de Iribarren
Esquema de difraccin de Wiegell

Obras Martimas
19
ROTURA O DESCOMPOSICIN DE LA ONDA. EFECTO FONDO

1.-
Munk - Mac Cowan
H
= 0.78
d
2.-
Gumbak
H
= 0.80 ; < 0.20
d
H
= 0.87 + 0.63 ; 0.20 < < 0.66
d
H
= 1.20 ; > 0.66
d
3.-
Weggel
H
=
d
b
; a = 43.70 ( 1 - exp ( 19 i )
ad
1+
g T2
b=
4.-
1.56
1 + exp ( 19.5 i )
Battjes
H
2 Hs
= 0.50 + 0.40 th ( 33 s0 ) ; s0 =
d
g T z2
5.-
Kimshi - Saeki
H
= 5.618 m 0.40
d
Obras Martimas
20
6.-
Collins
H
= 0.72 + 5.60 m
d
7.-
Miche (1949), Latoine (1962), Longuet - Higgins y Fenton (1974)
H
= 0.88 ; Miche
d
H
= 0.73 ; Latoine
d
H
= 0.827 ; Longuet - Higgins
d
Hasta un nmero de treinta modelos de rotura por fondo.
EFECTO FORMA.
1.-
Miche en profundidades indefinidas
H
1
=
L
7
2.-
Miche en profundidades reducidas
H
2
= 0.142 th ( k h ) ; K =
L
L
L0 =

Obras Martimas
g T2
2d
; L = L0 th
2
L
21
tag
2 Hs
tag
, s om =
, =
2
g Tz
H
som
L
siendo,
H1/3
Altura de ola significante o promedio del tercio de olas ms altas, m
H1/10 Altura de ola promedio del dcimo de olas ms altas, m

H1/50 Altura de ola promedio de 1/50 x N olas ms altas, m
Hmax Altura mxima del registro de N olas, m
Hd
Altura de ola de diseo, m
Hb
Altura de ola en rotura, m
Nmero de Iribarren, -
som
Peralte adimensional, -
tag Talud de la plataforma,

i
Pendiente del terreno natural, -
Pendiente del terreno natural, - (la notacin vara segn formulaciones)
Profundidad de la lmina de agua, m
Profundidad de la lmina de agua, m (la notacin vara segn formulaciones)
L0
Longitud de onda en profundidades indefinidas, m
Longitud de onda en profundidades reducidas, m
Aceleracin de la gravedad, m/s
Perodo ondulatorio, s
Nmero de onda, -
a, b Constantes de ajuste
Nmero pi
Nmero de olas activas de un registro de oleaje, -
Los Ejercicios presentados a continuacin representan casos prcticos y reales de estudio.
Se adjuntan igualmente los grficos de Goda y Weggel explicativos de la rotura.

Obras Martimas
22
Esquema de Goda para la determinacin de la ola rota en funcin de la pendiente de la

plataforma, el perodo ondulatorio y la altura de ola H0

Obras Martimas
23

Obras Martimas
24
Ejercicio 1
El morro de una estructura de abrigo convencional se encuentra situado a la profundidad de
diez metros en una zona de 4.50 metros de carrera de marea. El perodo del oleaje es de 15
segundos, y la altura de ola significante en profundidades indefinidas es de 14.25 metros. El
viento generador del temporal ms desfavorable es 26 m/s. El temporal tiene direccin NW
teniendo una orientacin angular del frente con relacin a la batimetra en aguas profundas de
22.50, 0 = 22.50.
Aplicando el baco de Goda - Mitsoyasu para un oleaje swell totalmente desarrollado,

calcular el coeficiente de refraccin - shoaling.
Aplicando la frmula de Blunt determinar el coeficiente de shoaling.
Determinar la altura de ola de diseo en funcin de la tipologa estructural y del posible anlisis
de rotura, tanto por peralte (H/L) o forma como por fondo (H/h).

Obras Martimas
25
FORMULACIN DE REFRACCIN Y SHOALING:

Coeficiente de refraccin segn la ley de Snell
seni =
sen 0
; Kr =
Ki
cos 0
sen 0
2
Ki
2
1-
Coeficiente de asomeramiento segn el modelo clsico de Blunt
4d
2d
th
L
L
4 d 4 d
sh
+
L
L
sh
Ks =
Modelo de Iribarren de refraccin
Hi
=
H0
B0
Bi
donde,
Kr
Coeficiente de refraccin
B0
Ancho del cuadriltero de avance en profundidades indefinidas, m
Bi
Ancho del cuadriltero de avance en profundidades reducidas, m
Angulo del frente con la batimetra en profundidades indefinidas,
Angulo del frente con la batimetra en profundidades reducidas,
Subndice de profundidades reducidas
Subndice de profundidades indefinidas
Ks
Coeficiente de shoaling, -
Profundidad de la lmina de agua, m
L0
Longitud de onda en profundidades indefinidas, m
Longitud de onda en profundidades reducidas, m
Perodo ondulatorio, s

Obras Martimas
26

2
Aceleracin de la gravedad, m/s
sh
Seno hiperblico, -
th
Tangente hiperblica, -
cth
Cotangente hiperblica, -
L0 =
g T2
2 d
; L = L0 th
2
L
CUANDO HAY REFRACCIN
d/L0 > 0.50
No hay refraccin
d/L0 < 0.50
Hay refraccin
d/L0
Entre 0.50 a 0.20 la refraccin es muy poco pronunciada
d/L0
Desde 0.20 a 0.050, la refraccin es muy acusada
SOLUCIN EJERCICIO 1
T = 15 segundos, L0 = 351 metros, Lpleamar = 173 metros, Lbajamar = 143 metros
Ki,p = 2.01, Ki,b = 2.43; Kr,pleamar = 0.97; Kr,bajamar = 0.97
Para la Bajamar Mxima Viva Equinoccial, se obtiene:
4
4d
2d
= 0.878 ; sh
= 1 ; th
= 0.413 ; k s = 1.12
L
L
L
Para la Pleamar Mxima Viva Equinoccial, se obtiene:
4
4d
2d
= 1.053 ; sh
= 1.259 ; th
= 0.483 ; k s = 1.061
L
L
L
Empleando Snell y Blunt, kr x ks = 1.04 en P.M.V.E. y 1.08 en B.M.V.E., adoptando 1.05
Obras Martimas
27
BACO DE GODA
d/L0 = 0.042 en pleamar y 0.028 en bajamar. Smax = 25, swell con bajo decay. Kr > 0.96 en
ambos casos
El coeficiente de shoaling calculado en pleamar resulta 0.84.
La altura de ola resultante ser:
Hs, indefinida x Ks x Kr = Hs, pie de dique

H0 = 14.25 x 1.05 = 14.96 m; Hs = 15.00 metros
Figura 1

Obras Martimas
baco de Goda de refraccin
28
baco 2
Diagrama de refraccin - shoaling del Shore Protection Manual
Empleando este modelo, resulta:
h
= 0.00657 ; kr k s = 1.05
g T2
h
= 0.00453 ; k r k s = 1.08
g T2
Adoptando un valor medio de 1.05, se obtiene una ola de clculo a pie de dique de 15.00
metros.

Obras Martimas
29
Ejercicio 2
Calcular con el procedimiento de Goda la altura de ola en rotura a una profundidad de 10
metros, en una zona de playa con pendiente 1/30 cuando sobre ella inciden oleajes de 10
segundos.
Ejercicio 3
Calcular lo mismo que el concepto anterior pero con la aproximacin de Osterdoff y Madsen
Ejercicio 4
Utilizando los bacos de Goda estimar la altura de ola significante en rotura a 8.00 metros de
profundidad cuando un swell H0'= 6.00 metros, T1/3 = 15 segundos, cuando incide sobre una
playa de pendiente 1/30.

Obras Martimas
30
Ejercicio 2
Empleando la frmula de Goda, previo al empleo de bacos, se obtiene:
4
hb
( 1 + 15 (tag ) 3 )
Hb = 0.17 L0 1 exp ( 1.50
L0
Para el caso que nos ocupa,
L0 =
g T2
10
h
; L0 = 156 m ; b =
= 0.064
2
156
L0
Entrando en la frmula, se obtiene, Hb = 7.80 metros
baco de Goda, lmite de rotura de oleaje regular

Obras Martimas
31
Ejercicio 3
La frmula de Osterdoff y Madsen, 1.979, resulta:
2 hb
Hb
= 0.14 th ( 0.80 + 5 tag )
Lb
Lb
Hb
= 0.14 th
Lb
2 hb
1.30
Lb
; tag < 0.10
; tag > 0.10
La longitud de onda en profundidad indefinidas es de 156 metros. La longitud de onda en

profundidades reducidas es 92.50 metros. La pendiente es 1/30, se toma la primera parte de
la frmula de Osterdoff y Madsen,
Hb = 92.50
( 0.14 )
= 7.50 metros
Ejercicio 4
Utilizando los bacos de Goda, la longitud de onda en profundidades indefinidas es 351
metros. El peralte equivalente y la profundidad relativa se obtienen:
*
6
h
8
H0
=
= 0.0171 ;
=
= 1.33
*
351
6
L0
H0
Entrando en el baco de pendiente 1/30, se obtiene la relacin:
H1
3
*
0
= 0.87 ; H 1 = 0.87 6 = 5.20 metros

3
Se adjuntan los bacos representativos.

Obras Martimas
32
Diagramas para la estima de la altura de ola en la zona surf o de rompiente de Goda

Obras Martimas
33
Diagramas para la estima de la altura de ola en la zona de surf o de rompiente de Goda

Obras Martimas
34
Ejercicio 5
Sea un oleaje de perodo 20 segundos y altura de ola de 1.00 metro propagndose sobre un
fondo en pendiente y formando un ngulo de 30 grados con la batimtrica de - 20.00 metros.
Se pide, definir la oblicuidad y la altura de ola a 10 metros de calado.
Ejercicio 6
Calclese la altura de ola de una onda unitaria a una profundidad de 20 metros cuando
alcance la profundidad de 10 metros. Se supone que la oblicuidad inicial respecto al fondo es
de 0 grados. Es de aplicacin la teora lineal.
Ejercicio 7
Mismo ejemplo anterior pero con oblicuidad de 30 grados.
NOTA: Aplquese la teora lineal de ondas para aguas poco profundas, es decir,
c =

Obras Martimas
g d = cg
35
Ejercicio 5
La ley de Snell puede escribirse de la forma siguiente,
sen 1
c1
sen 2
c2
; c1 =
g d1 ; c 2 =
g d2
La expresin en alturas de ola resulta,

2
2
H1 cos 1 = H2 cos 2
El ngulo crtico se define,
c2
sen 1 < 1
c1
Por tanto, en este caso resulta, c1 = 14.00 m/s, c2 = 9.90 m/s, empleando la Ley de
Snell,
sen 2
sen 30
=
9.90
14
2 = 20.70 grados
La altura de ola resulta, H2 = 1.00 x (cos 30/cos 20.70)0.50 = 0.96 metros
Ejercicio 6
La celeridad a diez metros de profundidad resulta, 9.90 m/s. La celeridad a 20.00 metros de
profundidad es de 14.00 m/s.
2
2
H c g = cte ; H c g cos = cte
2
2
1 14 = H2 9.90 ; H2 = 1.19 m

Obras Martimas
36
Ejercicio 7
La relacin es:
2
2
1 14 0.866 = H2 9.90 0.9355 ; H2 = 1.14 m
Ejercicio 8
Empleando los diversos esquemas de rotura por fondo definir la altura de ola en rotura en un
punto con seis metros de profundidad, 3.40 metros de carrera de marea y la sobrelevacin
meteorolgica dada por una borrasca de gradiente de presin de 953 milibares. La pendiente
del lecho es 1/20. El perodo es de 12 segundos.
Solucin del Ejercicio 8

Empleando el criterio de sobrelevacin por gradiente clsico:
p = 0.010 ( 1013 P a ) ; Pa milibares

n
La sobrelevacin sale 1013 - 953 = 60 cm. Por tanto, la lmina de agua es la profundidad con
la carrera de marea y la sobrelevacin, es decir, 6.00 + 3.40 + 0.60 = 10.00 metros.
a.-
Modelo de Mac Cowan, H = 0.78 x 10 = 7.80 metros
b.-
Modelo de Miche, H = 0.88 x 10 = 8.80 metros
c.-
Modelo de Latoine, H = 0.73 x 10 = 7.30 metros
d.-
Modelo de Longuet - Higgins y Fenton, H = 0.827 x 10 = 8.27 metros
e.-
Modelo de Collins, H = 10 x (0.72 + 5.60/20) = 10.00 metros

Obras Martimas
37
f.-
Modelo de Weggel, a = 26.79, b = 1.13, H = 9.49 metros
g.-
Modelo de Gumbak, Nmero de Iribarren inferior a 0.20, H = 8.00 metros
Se observa la notable variacin entre los distintos modelos, si bien, se puede decir,
H
< ; 0.70 < < 1.00
h
Ejercicio 9
Determinar mediante el mtodo de las ortogonales, para un frente de ondas en profundidades
indefinidas de 10 segundos de perodo y 5 metros de amplitud de onda, la amplitud de la
misma a la sonda o batimtrica - 15.00 metros de profundidad, sabiendo que las isobatas son
lneas rectas, paralelas entre ellas, el perodo permanece constante y el frente forma en
indefinidas un ngulo de 45 grados.
Ejercicio 10
Calcular el ngulo del frente a profundidad de 15 metros, el coeficiente de refraccin individual
y el coeficiente de refraccin - shoaling mediante el baco del Shore Protection Manual

Obras Martimas
38
Ejercicio 9
FORMULACIN DE REFRACCIN Y SHOALING:
Coeficiente de refraccin segn la ley de Snell
seni =
sen 0
; Kr =
Ki
cos 0
sen 0
2
Ki
2
1-
Coeficiente de asomeramiento segn el modelo clsico de Blunt
4d
2d
th
L
L
4d 4d
sh
+
L
L
sh
Ks =
Modelo de Iribarren de refraccin
Hi
=
H0
B0
Bi
Kr
Coeficiente de refraccin, -
B0
Ancho del cuadriltero de avance en profundidades indefinidas, m
Bi
Ancho del cuadriltero de avance en profundidades reducidas, m
Angulo del frente con la batimetra en profundidades indefinidas,
Angulo del frente con la batimetra en profundidades reducidas,
Subndice de profundidades reducidas, -
Subndice de profundidades indefinidas, -
Definimos los parmetros de la onda, es decir, la longitud de onda en profundidades

indefinidas y a 15 metros, resultando:

Obras Martimas
39

2
L0 = 1.56 x T = 156 metros
L, tanteando en la ecuacin trascendente, 109 metros
La constante Ki sale 1.43, que al elevarla al cuadrado es 2.048.
El ngulo del frente con la batimetra en indefinidas es 45 grados. Aplicando la ley de Snell, se
obtiene un coeficiente de refraccin de 0.90177, por lo que la altura de ola es de 4.51 metros.
Ejercicio 10
Empleando el SPM
El baco del Shore Protection Manual tiene en abscisas d/gT2 y 0 valiendo respectivamente,
0.015 y 45 grados.
El ngulo obtenido en profundidad de 15 metros es 30 grados, y el cociente de las funciones

coseno a la potencia 0.50, resulta:
Kr =
cos 0
; Kr =
cos i
cos 45
=
cos 30
0.707
= 0.9035
0.866
En este caso, Hi = 0.9035 x 5 = 4.51 metros
El coeficiente de refraccin - shoaling resulta 0.84, y por tanto, la ola es de 4.20 metros.

Obras Martimas
40
4797000
4796800
4796600
4796400
4796200
4796000
4795800
4795600
4795400
4795200
4795000
4794800
4794600
550600
550800
551000
551200
551400
551600
551800
552000
552200
552400
552600
552800
553000
553200
553400
Batimetra real reproduciendo un can submarino

Obras Martimas
41

4796500
4796400
4796300
4796200
4796100
4796000
4795900
4795800
4795700
4795600
4795500
4795400
551100
551200
551300
551400
551500
551600
551700
551800
551900
552000
552100
552200
552300
552400
552500
552600
Ejemplo de propagacin mediante un modelo numrico de refraccin - difraccin

Febrero 2013

Obras Martimas
42

TEMA5 - Grado ICT Refracción, Difracción y Rotura 2012

Cargado por

Copyright:

Formatos disponibles

TEMA5 - Grado ICT Refracción, Difracción y Rotura 2012

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

TEMA5 - Grado ICT Refracción, Difracción y Rotura 2012

Cargado por

Copyright:

Formatos disponibles

TEMA V

REFRACCIN, DIFRACCIN Y ROTURA

Niveles del mar y nivel medio del mar

Marea astronmica y otras fluctuaciones peridicas

Marea meteorolgica y otras variaciones de origen climtico

Propagacin y transformacin del oleaje: refraccin, difraccin y rotura

Efectos de la rotura del oleaje sobre las costas: erosin y transporte de

Vicente Negro Valdecantos

Curso Acadmico 2012 2013

REFRACCIN. Razonamiento terico del problema

La modificacin de la onda con los procesos de refraccin o expansin frontal de la energa

Los modelos de propagacin se basan en mltiples teoras.

Matemtica aplicada a la onda con su ecuacin de Laplace basada en las hiptesis de

div u = 0, la divergencia de un campo vectorial es un escalar

Soluciones basadas en la ecuacin de pendiente suave de Berkoff, 1972 o "mild slope

Resoluciones de manera elptica

Resoluciones en aproximaciones parablicas

Modelos basados en las ecuaciones de conservacin

Refraccin, difraccin y rotura

Curso Acadmico 2012 2013

Conservacin de la accin de onda

Irrotacionalidad del nmero de onda

Ancho del cuadriltero de avance, m

Subndice de profundidad indefinida, -

Subndice de profundidad considerada, -

; H12 cos 1 = H2 2 cos 2

Oblicuidad, ngulo que forma el frente con la batimetra,

Celeridad del movimiento, m/s

Indicadores de referencia, punto 1 y punto2

Refraccin, difraccin y rotura

Curso Acadmico 2012 2013

Con hiptesis semejantes y en funcin de parmetros espectrales direccionales (Smax),

Refraccin, difraccin y rotura

Curso Acadmico 2012 2013

celeridad a profundidad "p", m/s

La clsica ecuacin de Laplace,

Refraccin, difraccin y rotura

Curso Acadmico 2012 2013

* (x, y) x w(z) x exp(- iwt)

A(x, y) x exp {is (x, y)}

Transformando la ecuacin de Laplace de 3D a 2D, integrando,

La descripcin completa de los fenmenos de propagacin de ondas monocromticas lineales

Celeridad de la onda, m/s

Celeridad del grupo, m/s

Refraccin, difraccin y rotura

Curso Acadmico 2012 2013

Longitud de onda en profundidad reducida, m

Longitud de onda en profundidad indefinida, m

La solucin numrica emple dos tendencias claramente diferenciadas, el mtodo de los

La resolucin del problema elptico requera adems de un tiempo de computacin muy

Refraccin, difraccin y rotura

Curso Acadmico 2012 2013

El conocimiento de la condicin de contorno en la costa casi siempre desconocida, obliga a

Para simplificar el modelo elptico y descomponiendo la funcin potencial en una componente

Berkoff, Booij y Radder, 1.982

La principal desventaja del modelo parablico se centra en la malla o celda de integracin,

Refraccin, difraccin y rotura

Curso Acadmico 2012 2013