Mecánica de Olas

TRANSPORTE LITORAL COSTERO
CAPITULO 3. MECANICA DE OLAS
3.1 Introducción
Las olas generadas por el viento son la principal fuente de energía a la zona litoral, son las
responsables de los procesos de erosión y deposición de sedimentos en las playas, determinando de
esta manera su geometría y composición. Las olas generalmente derivan su energía de las
tormentas en el mar adentro. Una cantidad significante de esta energía puede disiparse antes de que
las olas alcancen el lugar de diseño. La energía en la region cercana a la costa ocurre en un rango
amplio de frecuencias. Una distribución aproximada de esta energía se muestra en la Figura 3-1.
Las olas de interés caen dentro de la banda de olas de gravedad (períodos entre 3 y 25 s). La
energía de estas olas forma playas; sortea sedimentos de fondo en la cara de la playa; transporta
materiales de fondo hacia la costa, mar adentro, y paralelo a la costa; y ejerce fuerzas sobre
estructuras costeras. Otros estudios han sugerido que las olas de la banda de infragravedad también
juegan un papel importante en los procesos litorales (e.g., Komar & Holman 1986). El
entendimiento de la acción de las olas es fundamental para la comprensión de los procesos que
tienen lugar en la zona costera. Por ejemplo, la determinación de las condiciones apropiadas de
oleaje es necesaria antes de proceder a la estimación de corrientes y transporte de sedimentos y a la
planeación y diseño de estructuras costeras.
Las olas oceánicas son complejas e irregulares, pero haciendo algunas consideraciones teóricas se
idealizan. Aqui examinaremos brevemente algunas teorías de Oleaje Regular que nos
proporcionarán fórmulas matemáticas simples las cuales pueden aplicarse a las olas del mar.
Veremos las suposiciones envueltas en su formulación y las limitaciones resultantes de su
aplicación.
Las teorías suponen trenes de olas largos, suaves y uniformes de tal forma que pueden generarse
aproximadamente en un laboratorio. Las olas de Swell, distintas a las del área de generación, a
menudo se aproximan a esta condición de tal forma que uno podría aplicar las teorías para describir
los movimientos de las olas y la transferencia de energía. Otras veces las olas reales son de crestas
cortas y más aleatorias e irregulares que las tratadas por las teorías. Posteriormente veremos como
este Oleaje Irregular se mide y analiza.
La sección de Oleaje Regular delinea los principios fundamentales que rigen la mecánica del
movimiento del oleaje, esencial en la planeación y diseño de trabajos en la costa. La sección de
Oleaje Irregular discute las teorías estadísticas y probabilísticas que se aplican.
La descripción detallada de las ecuaciones básicas de la mecánica de olas puede consultarse en

varios libros de texto (por ej, Kinsman 1965; Stoker 1957; Ippen 1966; Le Méhauté 1976; Phillips
1977; Crapper 1984; Mei 1991; Dean & Dalrymple 1991). La sección de Oleaje Regular solamente
ofrece una introducción a la mecánica de las olas, y se enfoca en teorías sencillas de oleaje para
ingenieros costeros. Se discuten métodos para estimar los perfiles de olas superficiales, movimiento
carlos r torres iio/uabc 4/18/2008 1

de las partículas, energía de la ola, y transformaciones del oleaje debidas a la interacción con el
fondo y con estructuras.
Figura 3-1. Distribución aproximada de la energía de las ondas oceánicas superficiales (de Kinsman
1965).
Teorías más Importantes
i) Teoría Lineal de Airy: olas de amplitud pequeña en aguas profundas

ii) Stokes: olas de amplitud finita en aguas profundas, intermedias y someras.
iii) Olas Cnoidales: olas de amplitud finita en aguas intermedias y someras.
iv) Ola Solitaria: crestas aisladas o solitarias de amplitud finita que se mueven en aguas someras.
La teoría de oleaje más sencilla es de primer orden, de amplitud finita o la teoría de Airy, que
llamaremos aqui Teoría Lineal. Muchos problemas ingenieriles pueden manejarse sin dificultad
con esta teoría. Por convenciencia, los métodos de predicción en ingeniería costera generalmente se
basan en olas sencillas. En algunas situaciones, las teorías sencillas ofrecen estimaciones aceptables
de las condiciones de oleaje.
Cuando las olas llegan a ser grandes o viajan hacia la costa en la región de aguas someras, las
teorías de oleaje de orden superior a menudo se requieren para describir el fenómenos del oleaje.
Estas teorías reprsentan olas alineales. La teoría lineal que es válida cuando las olas son
infinitesimalmente pequeñas y su movimiento es pequeño también nos ofrece alguna información
para ondas periódicas de amplitud finita (alineal). Sin embargo, la teoría lineal no puede tomar en

cuenta el hecho de que arriba del nivel medio del mar las olas son mayores que los valles. Los
resultados obtenidos de las diferentes teorías deben de interpretarse cuidadosamente para su uso en
el diseño de proyectos costeros o para la descripción del ambiente costero.
En la práctica es importante diferenciar entre dos tipos de olas superficiales. Mar y mar de leva
(swells). Mar se refiere a olas de período corto que se están generando por los vientos. Mar de leva
se refiere a olas que se han movido fuera de la zona de generación. En general, el mar de leva
consta de olas más regulares, con crestas bien definidas y períodos relativamente largos.
El crecimiento de las olas oceánicas no es indefinido. El punto en que dejan de crecer se llama una
Condición de mar completamente desarrollado (a fully developed sea condition). La energía del
viento se imparte al agua dando como resultado el crecimiento de las olas; Sin embargo, después de
un punto, la energía impartida al agua se disipa por las olas que rompen. Las olas de Mar son de
crestas cortas e irregulares y sus períodos caen en el rango entre 3 a 25 seg. Estas olas generalmente
tienen períodos y longitudes cortos y su superficie aparece más perturbada que la del mar de leva.
Las olas aparecen más ordenadas cuando ya no están bajo la influencia del viento.
3.2 Oleaje Regular

3.2.1 Introducción
Las teorías de olas son aproximaciones a la realidad. Estas pueden describir ciertos fenómenos bajo
ciertas condiciones que satisfacen las suposiciones hechas en su derivación. También pueden fallar
en describir otros fenómenos que violan esas suposiciones. Dado que el diseño y protección costera
depende de la habilidad para predecir los perfiles superficiales de las olas y su movimiento, así
como de la exactitud de las predicciones, debe tenerse especial cuidado de que la teoría adoptada
asegure que el fenómeno de oleaje de interés sea descrito razonablemente por ella.
3.2.2 Definición de parámetros de la ola
Una onda progresiva puede ser representada por una variable espacial x, una temporal t o por su
combinación (la fase, θ = κx − ωt). Los valores de θ varían entre 0 y 2π. La Figura 3-2 ilustra los
parámetros que definen una onda progresiva que pasa por un punto fijo en el oceáno. Una ola de
forma permanente que se propaga sobre un fondo horizontal puede caracterizarse completamente
por su altura, Η, su longitu L y la profundidad a la que se encuentra h.
Al intervalo de tiempo entre el paso de dos crestas (o valles) sucesivas por el mismo punto se le
denomina período de la ola, Τ. L es la distancia horizontal entre dos puntos idénticos sobre dos
crestas o valles sucesivos.
Otros parámetros incluyen la frecuencia angular (ω = 2π/Τ), el número de onda κ = 2π/L, la

velocidad de fase o celeridad de la ola C = L/Τ = ω/κ, la esbeltéz de la ola δ = Η/L, la profundidad
relativa h/L, y la altura relativa Η/h. Los parámetros adimensionales κH and κh, preferidos en
trabajos de investigación, pueden sustituirse fácilmente por Η/L y h/L, respectivamente, ya que solo
difieren por un factor de 2π de aquellos utilizados por los ingenieros.

Figura 3-2. Definición de términos – onda elemental senoidal progresiva
3.2.3 Teoría Lineal o de Airy
En 1776, Laplace presentó el primer tratamiento satisfactorio para olas de amplitud pequeña en
aguas profundas arbitrarias. Airy, en 1845 desarrolló una teoría de olas no-rotacionales que viajan
sobre un fondo horizontal en aguas de cualquier profundidad. Para su obtención, las ecuaciones de
movimiento se linealizan y, por esta razón, generalmente esta teoría se denomina Teoría Lineal del
Oleaje. Una descripción teórica más completa de las olas puede obtenerse como la suma de muchas
aproximaciones sucesivas, donde cada término adicional de la serie es una corrección del
precedente.
Si no tomamos en cuenta la viscosidad del agua, y si consideramos olas fuera del área de
generación, entonces no habrá pérdida ni ganancia de energía. Más aun, si el movimiento es
irrotacional, entonces la velocidad horizontal u y la vertical w de las partículas, pueden obtenerse
de un potencial de velocidad, ϕ(x, z, t), tal que;
u = ∂ϕ/∂x; w = ∂ϕ/∂z (2.1)
Si la densidad del agua de mar, ρ, es constante, entonces;
∂u/∂x + ∂w/∂z = 0 (2.2)
o, en términos de ϕ(x, z, t), queda como;
∇2ϕ = 0 (2.3)

la cual se conoce como Ecuación de Laplace. Las condiciones de frontera son las siguientes:
a) no puede pasar fluído a través del piso oceánico:
w(x, −z, t) = 0
b) la presión en la superficie tiene que ser igual a la atmosférica. Esto se aplica mediante la
ecuación de Bernoulli;
p/ρ = ∂ϕ/∂t – (u2 + w2)/2 − gz = 0 (2.4)
Suponiendo que la pendiente de la ola, Η << L, es pequeña y que h >> Η, se obtiene una solución
en la cual la elevación de la superficie, η, está dada por;
η( x, t) = Η/2 cos(κx − ωt) (2.5)
O sea, el perfil resultante es senoidal, con período Τ, longitud, L, y altura Η . Una cantidad asociada
es la frecuencia f (= 1 / Τ)
Un resultado fundamental obtenido en la solución, lo constituye la relación de dispersión;
ω2 = gκ tanh(κh)
donde g es la aceleración debida a la gravedad.
Sustituyendo los equivalentes de κ y ω, queda;
L = (gΤ2/2π) tanh (2πh/L ) = gΤ/ω tanh(κh) (II-1-10)
que es la expresión general para la relación de dispersión de ondas de gravedad o "wind waves". La
utilización de esta ecuación involucra alguna dificultad pues la desconocida L aparece en ambos
lados de la ecuación. Pueden utilizarse los valores tabulados de h/L (SPM 1984) para simplificar su
solución. Eckart (1952) da una expresión de aproximación correcta dentro del 10%. Esta expresión
está dada por
L ≈ (gΤ2/2π)(tanh (4π2 h /gΤ2))1/2
El error máximo del 10% ocurre cuando h/L ≈ 1/2. Existen otras aproximaciones para resolver la
ecuación II-1-10 (Hunt 1979; Venezian and Demirbilek 1979; Wu and Thornton 1986; Fenton and
McKee 1990).
Si κh →∞, que es el límite para aguas profundas, entonces tanh(κh) se aproxima a la unidad
quedando;
L = (gΤ2/2π)
y su velocidad de fase, será ;

Cο = Lο/Τ = (gΤ/2π)
Es decir, en aguas profundas la longitud y la velocidad de fase dependen del período, ambos siendo
mayores para las olas más largas. El límite para aguas profundas es; h/L > 1/2. Cuando usamos
estas expresiones y h/L = 0.25 el error resultante será del 9%.
Aunque las aguas profundas ocurren a una profundidad infinita, tanh(κh), para la mayoría de
propósitos prácticos se aproxima a la unidad a un h/L mucho menor que infinito. Para un h/L = 1/2
(i.e., cuando la profundidad es la mitad de la longitud de onda), tanh(κh) = 0.9964.
El otro límite es cuando κh → 0 (i.e. h << L), por lo que la función tanh(κh) → κh. Por tanto, la
longitud de onda se reduce a;
L = Τ(gh)1/2
mientras que velocidad de fase, queda como;
C = (gh)1/2
o sea, C es función de la profundidad. El límite de aplicabilidad de esta aproximación, es; h/Lο <
1/20, con 5% de error.
Resumiendo, olas pueden clasificarse de acuerdo al criterio de profundidad relativa h/L como;
Clasificación h/L κh tanh(κh)

Aguas Profundas 1/2 a ∞ πa∞ ≈1
Transicionales 1/20 a 1/2 π/10 a π tanh (kh)
Someras 0 a 1/20 0 a π/10 ≈ kh
Como regla práctica, las ecuaciones completas tienen que utilizarse cuando h/L esté entre 0.5 y
0.04.
Algunas funciones útiles
Haciendo la operación C/C0 nos queda;
C/C0 = L/L0 = tanh(κh)
Si multiplicamos ambos lados de la ecuación por h/L, queda;
h/L0 = h/L tanh(κh)
Los términos h/L0 y h/L asi como otras funciones como κh = 2π h/L y tanh(κh) han sido tabuladas
por Wiegel (1954) como función de h/L0 (ver también SPM 1984, Apéndie C, Tablas C-1 y C-2).

Ejemplo 1
Una ola de período Τ = 10 seg se propaga hacia la costa sobre un fondo de pendiente uniforme
desde una profundidad de h = 200 m (656 ft) a una profundidad de h = 3 m (9.8 ft).
Encuentre las celeridades C y longitudes de la ola L correspondientes a las profundidades

anteriores.
Solución:
a) La longitud
Usando la ecuación (…)
L0 = (gΤ2/2π) = 1.56Τ2 m (5.12 Τ2 ft)
L0 = 1.56(10)2= 156 m
Para h = 200 m
h/L0 = 200/156 = 1.2821
Notese que para valores de

h/L0 > 1,
h/L0 = h/L0
Por tanto,
L = L 0 = 156 m (aguas profundas ya que h/L > 1/2)
b) La celeridad
C = L/Τ = 156/Τ =156/10=15.6 m/s
c) Para h = 3 m
h/L0 = 3/156 = 0.0192
de tablas se encuentra que
h/L = 0.05641
Por tanto,

L = 3/0.05641= 53.2 m (174 ft) (profundidad transicional ya que 1/25 < h/L < 1/2)
e) = /Τ = 53.2/10=5.32 m/s (17.4 ft/s)
Un valor aproximado para Λ puede obtenerse utilizando la ecuación
L ≈ (gΤ2/2π) (tanh(4π2h/Τ2 ))1/2
que en términos de L0 puede escribirse como;
1/2
≈ 0 ( (2π / 0))
Por tanto,
L ≈ 156 ( (2π(3)/156))1/2
L ≈ 156 ( (0.1208))1/2
L ≈ 156 (0.1202)1/2 = 54.1 (177.5 ).
que se compara may bien con el L = 53.3 que se había obtenido anteriormente. El error en este
caso es del 1.5 %. El valor de h/L pudo haberse obtenido de la Tabla C-1 (SPM 1984) o de la figura
5.
Velocidades y aceleraciones de las partículas
Las componentes de la velocidad en la vertical, ω, y en la horizontal, υ, están dadas por;
= Τ/ ( ( + )/( )) ( −ω) (22)
= Τ/ ( ( + )/( )) ( −ω) (23)
Las velocidades de las partículas de fluido bajo un tren de olas se muestra en la Figura 3-3.

Figura 3-3. Velocidades y aceleraciones locales de las partículas
Las ecuaciones expresan las componentes de la velocidad local del fluido a cualquier distancia (z +
d) sobre el fondo. Las velocidades son periódicas en x y t.
Las aceleraciones locales de las partículas de fluido se obtienen diferenciando las ecuaciones
anteriores respecto a t. Quedando,
αx= (gπH/L) [cosh[2π(z + d)/L]/cosh(2πd/L) sin θ = ∂u/∂t
αz= -(gπH/L) [sinh[2π(z + d)/L]/cosh(2πd/L) cos θ = ∂w/∂t
donde θ = (κx − ωt). Valores positivos y negativos de las aceleraciones de las partículas en la
vertical y horizontal se muestran en la figura anterior.
La figura indica que el fluido bajo la cresta se mueve en la dirección de propagación de la ola y
regresa durante el paso del valle. La Teoría lineal no predice ningún transporte de masa; por tanto,
la figura presenta solo un movimiento oscilatorio del fluido. La Figura 3-4 muestra perfiles de la
elevación superficial, η, y velocidades y aceleraciones dados por la Teoría lineal.

Figura 3.4. Perfiles de la velocidad y aceleración de las partículas (Teoría de Airy) y su relación con la
elevación superficial.
Ejemplo 2
Dada una ola de período T = 8 sec y altura H = 5.5 m (18.0 ft) en aguas donde la profundidad es d =
15 m (49 ft).
Encuentre:
Las velocidades horizontal y vertical u y w, asi como las aceleraciones and accelerations αx y αz en
z = -5 m (o z = -16.4 ft) bajo el NMM cuando θ = 2πx/L - 2πt/T = π/3 (o 600).
Solución;
0 = 1.56Τ2 =1.56(8)2=99.8
/ 0 = 15/99.8 = 0.1503

Usando la Figura II-1-5 para d/L0 = 0.1503, se encuentra que
/ = 0.1835
y
=1.742
por tanto,
= 15/0.1835=81.7μ
La evaluación de las constantes en las ecuaciones, da
( ΤΗ/2 ) 1/( ) =(5.5 × 9.8 × 8)/(2×81.7)(1/1.742)=1.515
(Η π/ ) 1/( ) =(5.5 × 9.8 × 3.1416)/(81.7)(1/1.742)=1.190
La sustitución de estos valores en las expresiones para u y w, da;
= 1.515 ( (15−5))/81.7 (600)
= 1.515 (cosh(0.7691))/81.7 (0.500)
porque
2πd/L = 0.7691 y los valores de las funciones hiperbólicas son;
cosh(0.7691)=1.3106
sinh(0.7691)=0.8472
por lo que las velocidades, son;
u = 1.515(1.3106)(0.500) = 0.99 m/s (3.26 ft/s)
w = 1.515(0.8472)(0.866) = 1.11 m/s (3.65 ft/s)

y las aceleraciones;
αx = 1.190(1.3106)(0.866) = 1.35 m/s2 (4.43 ft/s 2)
αz =1.190(0.8472)(0.500) = -0.50 m/s2 (1.65 ft/s 2)

Desplazamientos de las partículas
La integracion de las ecuaciones para u y w, da los desplazamientos horizontales y veticales de las

partículas respectivamente;
ξ = − (Η Τ2/4π ) ( ( + )/ ) θ
ζ = + (Η Τ2/4π ) ( ( + )/ ) θ
Estas expresiones pueden simplificarse mediante la utilización de la relación
(2π/Τ)2 = (2π / ) ( )
De esta manera, queda;
ξ = − (Η/2) ( ( + )/ ) θ
ζ = + (Η/2) ( ( + )/ ) θ
Elevando cada expresión al cuadrado, queda
2
θ= (ξ/ / ( + )) 2
2
θ= (ζ/ / ( + )) 2
Sumando, obtenemos;
ξ2/Α2 + ζ2/Α2 = 1
que es la ecuación de una elipse con semieje mayor A y semieje menor B, donde A y B están dados
por;
Α = Η/2 ( + )/
Β = Η/2 ( + )/
Las longitudes de A y B son una medida del desplazamiento horizontal y vertical de las partículas
de agua (ver Figura 3-5). De esta manera, la Teoría lineal predice que las partículas se moverán en
órbitas cerradas. O sea, las partículas de fluido retornan a su posición inicial después de cada ciclo
de la ola. Sin embargo, en la realidad las órbitas de las partículas no son completamente cerradas.
Esta diferencia se debe al fenómeno de transporte de masa, que discutiremos posteriormente.

( / > 1/2), sinhkd y coshkd se reducen a y, por tanto, A y B son iguales
por lo que
Α = Β = Η/2
O sea, en aguas profundas las órbitas son círculos (igual a la altura en la superficie) y disminuyen
su radio exponencialmente con la profundidad. La tasa de decremento depende del valor de
= 2π/ 0 = 2π/( Τ2/2π)
entre mayor sea el valor de κ, más rápido será el decremento con profundidad. Es decir, la
disminución es más marcada para olas de períodos cortos que para olas de períodos largos.
Para aguas someras ( / < 1/25), las ecuaciones quedan como;
Α = (Η/2) /2π y Β = (Η/2)(1 + / )
En aguas someras, las elipses son aplanadas y llegan a ser lineas rectas en toda la columna de agua.
En aguas someras casi siempre kh << 0.6, de forma tal que el diámetro orbital es mayor que la
altura. No hay dependencia en ζ, por lo que el diámetro orbital es constante desde la superficie al
fondo y la velocidad orbital υ es independiente de ζ.
En aguas de transición, las órbitas son elípticas. Entre menos profunda sea la columna, más plana
será la elipse.
Figure 3-5. Desplazamientos de las partículas de agua a partir de la posisicón media para aguas
someras, intermedias y profundas (tomado de EM 1110-2-1100)

Ejemplo 3
Dada una ola en una profundidad de d = 12 m (39.4 ft), con altura H = 3 m (9.8 ft), y período T =
10 sec. Su correspondiente altura en aguas profundas es de H0 = 3.13 m (10.27 ft). Encontrar
a) El desplazamiento horizontal y vertical máximo de la partícula de agua, a partir de su posición

media, cuando z = 0 y z = -d.
(b) El desplazamiento horizontal y vertical máximo de la partícula de agua a una profundidad de z

= -7.5 m (-24.6 ft) cuando la ola se encuentra en aguas infinitamente profundas.
(c) Muestre que para las condiciones de (b), los desplazamientos de las partículas son pequeños en
relación a la altura de la ola cuando z = -L0 /2.
Solución;
(a)
0= 1.56(10)2= 156
Para = 12
/ 0 = 12/156 = 0.0769
También,
( )=0.8306
( )=0.6389
Cuando = 0, Α y Β se reducen a;
Α = (Η/2) 1/ ( )
Β = (Η/2)
Sustituyendo valores, queda;
Α = (3/2) 1/ (0.6389) = 2.35

Β = (3/2) = 1.5
Cuando = −
Α = Η/(2 ( )) = 3/2(0.8306) = 1.81

Β=0
(b) Con H0 = 3.13 m and z = -7.5 (-24.6 ft), evaluamos el exponente de e. Notando que L = L0,
= 2π(−7.5)/156 = −0.302
de donde;
−0.302
= 0.739
Por tanto,
Α = Β = (Η0/2) = (3.13/2) 0.739 = 1.16
El desplazamiento máximo (o diámetro) de la órbita del círculo, será 2(1.16) = 2.32 (7.61 ft)
cuando z = −7.5 .
(c) A una profundidad correspondiente a la mitad de la longitud de la ola, tenemos;
ζ=− 0/2 = −156/2 = −78.0
= 2π(−78)/156 = −3.142
de donde;
−3.142
= 0.043
y
Α = Β = (Η0/2) = (3.13/2) 0.043 = 0.067
O sea, el desplazamiento máximo de la partícula es 0.067 , que es pequeño en comparación con la

altura de aguas profundas, H0 = 3.13 m (10.45 ft).
Presión Subsuperficial.
La presión bajo la superficie del mar se expresa como;
'= −ρ + ρ [( ( + ))/ ] ( −ω)
que es la suma de un componente dinámico y otro hidrostático. Aquí, ' es la presión total o
−3
absoluta, α es la presión atmosférica, y ρ es la densidad del agua de mar (ρ = 1025 ó
−3
ρ = 1000 para agua dulce, en unidades inglesas, es 1.94 slugs/ft3). Por conveniencia, es
común denotar la presión como, ,

= '− =−ρ +ρ η
que es la presión registrada por un sensor. Aqui, es el factor de respuesta a la presión, y está
dado por;
= ( + )/
Con estas relaciones, es posible obtener las elevaciones del oleaje, η, a partir de registros de
presión, como sigue;
η = [N(p + ργz)]/ρ
donde N es un factor de corrección igual a la unidad cuando se utiliza la Teoria Lineal. El valor de
Κπ para un sensor instalado en el fondo, es;
= 1/
Habíamos visto que en aguas profundas
= ω ( −ω)
= ω ( −ω)
en tanto que , queda como;

=ρ ( −ω)
O sea, está en fase con la elevación superficial a cualquier profundidad y la amplitud de la

presión, ρ , disminuye exponencialmente con la profundidad. Ni las velocidades de las
partículas ni Δp dependen de la profundidad del agua. Es decir, en aguas / > ½, Las partículas
no “sienten fondo”.
Ejemplo
Un sensor de presión instalado a 0.6m (1.97 ft) arriba del fondo marino, en aguas donde la
profundidad es = 12m (39ft), registra una fluctuación de presión p = 124 kilonewtons/m2. La
frecuencia promedio del oleaje es de φ = 0.0666 hz.
Encontrar la altura de la ola Η suponiendo que se aplica la Teoría Lineal.
Solución;
Τ = 1/ = 1/0.0666 = 15

0= 1.56Τ2 = 351 (1152 )
/ 0= 12/351 = 0.0342
Entrando en la Figura II-1-5, con d/L0, obtenemos;
/ = 0.07651
por tanto,
= 12/0.07651 = 156.8 (515 )
= 1.1178
Por tanto, el valor de queda como;
= ( + )/ = [2π(−11.4+12)/156.8]/1.1178=0.8949
Dado que η = = Η/2 cuando la presión es máxima (bajo la cresta de la ola), y N = 1.0 ya que
estamos usando la Teoría lineal, queda;
Η/2 = [N( + ρ )]/ρ = 1.0[124+(10.06)(−11.4)]/[(10.06)(0.8949)] = 1.04 (3.44 )
Por tanto,
Η = 2(1.04) = 2.08 (6.3 )
Velocidad de Grupo
Muchas veces es útil saber que tan rápido se mueve la energía de la ola. Una manera de hacerlo es
ver la velocidad de los grupos de olas, , que representan la propagación de la energía de la ola en
espacio y tiempo. Generalmente, la velocidad del grupo no es idéntica que la velocidad de las olas
individuales, , que conforman el grupo viajero.
El concepto de velocidad de grupo puede describirse considerando la interacción de dos trenes de

olas senoidales que se mueven en la misma dirección con longitudes de onda y períodos
ligeramente diferentes, como sigue;
η = η1 + η2 = Η/2 cos(2π / 1− 2π /Τ1) + Η/2 cos(2π / 2− 2π /Τ2)
donde η1 yη2 son los dos componentes. Aqui los estamos sumando porque esto es permisible
cuando la Teoría lineal se utiliza. Por simplicidad, las alturas de ambos componentes de las olas se

supone que son iguales. Dado que,  and 2, son ligeramente diferentes, para algunos valores de
a un tiempo dado, τ, los dos componentes estarán en fase y la altura de la ola será 2Η; en otros
casos, para algunos otros valores de , las dos componentes de la ola estarán fuera de fase y la
altura resultante será nula. El perfil superficial resultante está dado por la siguiente ecuación;
η = ± Η cos(π[( 2− 1)/ 1 2] − π[(Τ2 − Τ1)/Τ1Τ2] )
y se presenta en la Figura 3-6.
Figura 3-6. Características de un grupo de olas formado por la suma de dos senoides de período
diferente.
Es la velocidad de estos grupos lo que representa a la velocidad de grupo. La velocidad limitante de

los grupos de olas, cuando L1 se aproxima a L2 (y consecuentemente T1 se aproxima a T2) es la
velocidad de grupo, Cg, y puede demostrarse que esta dada por;
= 1/2(1+2 / 2 ) /Τ =
donde es el factor de someramiento;
= 1/2(1+2 / 2 )
En aguas profundas el término 2 / 2 es casi cero y = 1/2, por lo que;

= 1/2 ( 0/Τ) = 1/2
O sea, la velocidad del grupo es la mitad de la velocidad de la celeridad.
En aguas someras, ≈ , y = 1, entonces;
= /Τ = ≈( )1/2
O sea, la velocidad del grupo y la celeridad son iguales. Así, en aguas someras, como la celeridad
de la ola está determinada por la profundidad, todos las olas componentes de un tren de olas
viajarán a la misma velocidad evitando la cancelación y reforzamiento alternados de los
componentes. En aguas profundas y de transición, la celeridad depende de la longitud de la ola, por
tanto, olas ligeramente largas viajarán un poco más rápido y producirán las pequeñas diferencias de
fase que resultarán en grupos de olas. Se dice que estas olas son dispersivas o que se propagan en
un medio dispersivo (un medio donde la celeridad depende de la longitud de la ola).
La variación de los cocientes Cg/C0 and C/C0, como función de d/L0 se presenta en la Figura 3-7.
Fuera de aguas someras C > Cg. Un observador que sigue un grupo de olas a la velocidad de grupo
verá que las olas se originan al final del grupo moviéndose después hacia adelante através del grupo
viajando a la velocidad de fase, y desapareciendo posteriormente en el frente del grupo de olas.
Figura 3-7. Variación de los cocientes de la velocidad de grupo y la velocidad de fase con la velocidad
de fase en aguas profundas usando Teoría Lineal (Sarpkaya e Isaacson 1981)

El concepto de velocidad de grupo es importante porque es con esta velocidad con la que se
propaga la energía. Aunque de manera rigurosa la velocidad de grupo se obtiene de la interferencia
de dos o mas olas, el significado físico no es tan obvio como en el caso en que su derivación se
basa en el método de la energía de la ola. Eso es lo que veremos a continuación.
Energía y potencia de la ola
Si bien no hay movimiento neto del agua bajo la acción de la ola, ya que las partículas retornan a su
posición inicial, el movimiento de la ola misma constituye una transferencia de energía sobre la
superficie del mar. Es especialmente en este aspecto que las olas son importantes para las playas. El
desplazamiento de la ola sobre la superficie le da a ésta forma de energía potencial (Ε ) . Al mismo
tiempo, el movimiento orbital del agua bajo las olas constituye una energía cinética para las olas
(Ε ). La energía total Ετ por unidad de longitud de cresta de la ola para una ola definida con la
Teoría lineal, está dada por la suma de Επ y Εχ como sigue;
Ε = Ε + Ε = ∫ρ ((η + )2/2 − 2
/2) + ∫∫ρ( 2 + 2
)/2
donde los límites de integración en son entre ( , + ) y en entre (– , η) o sea, en todo un

perfil de ola. Al integrar, resulta;
Ε = Ε + Ε = ρ Η2 /16 + ρ Η2 /16 = ρ Η2 /8
que es la energía total de la ola en una longitud de onda por unidad de ancho de cresta. La energía
promedio total, denominada energía específica o densidad de energía, está dada por;
Ε’= Ε / = ρ Η2/8
Flujo de Energía
El flujo de energía de la ola es la tasa a la cual se transmite la energía en la dirección de

propagacion de la ola a travéz de un plano perpendicular a la dirección de avance de la ola y que se
extiende hasta el fondo. Suponiendo Teoría lineal, el flujo promedio de energía por unidad de
ancho de cresta de la ola transmitido a travéz de un plano perpendicular a la dirección de avance de
la ola, es;
’=1/Τ ∫∫
donde los límites de integración en son entre (–δ, η) y en τ entre (τ, τ+ρ). Al integrar, resulta;
’= Ε' = Ε'
donde ’ frecuentemente se denomina, potencia de la ola y ν ya se ha definido anteriormente.

Si existe un ángulo entre el plano a travéz del cual se transmite la energía con respecto de la
dirección de avance de la ola, entonces,
' = Ε’ θ
donde θ es el ángulo entre el plano y la dirección de avance.
Para aguas profundas y someras, ' queda como;
'0 = (Ε0 Χ0)/2
'= Ε = Ε
Dada una cierta región, atravéz de la cual las olas pasan, un balance de energía revelará que la
cantidad de energía que entra en la región es igual a la cantidad de energía que la abandona,
siempre y cuando no se añada o remueva energía. Por tanto, cuando las olas se aproximan de tal
forma que sus crestas son paralelas a los contornos de fondo, tenemos;
E’0 n0 C0 = E’nC (64)
O, dado que n0 =1/2, entonces;
(E’0C0)/2 = E’nC (65)
Cuando las crestas de la ola no son paralelas a los contornos de fondo, algunas partes viajarán a
velocidades diferentes y la ola será refractada; en este caso, la ecuación (65) ya no se aplica (see
Parts II-3 and II-4). La tasa de transmisión de energía es importante para el diseño costero y
requiere un concocimiento de Cg para determinar que tan rápido las olas se mueven hacia la costa.
La ecuación (65) establece una relación entre el cociente de la altura de la ola en una profundidad
arbitraria y la altura de aguas profundas. Este cociente se denominado: coeficiente de someramiento
(see Part II-3 for detail derivation), depende de la esbeltéz de la ola. La variación del coeficiente de
someramiento con la esbeltéz de la ola como función de la profundidad relativa δ/Λ0 se muestra en
la Figura 3-8.

Figura 3-8. Variación de ν con la esbeltéz de la ola (Sakai y Battjes 1980)
Ejemplo
Una onda oscilatoria de aguas profundas con Λ0 = 156m, Η0 = 2m y Χ0 = 15.6 m/s se mueve hacia
la costa con su cresta paralela a los contornos de fondo. Los efectos debidos a la reflexión son
ignorados. Encontrar:
a) una relación entre las alturas en aguas profundas y someras suponiendo conservación de la
energía.
b) la altura de la ola cuando la profundidad es de 3m
c) determinar la tasa a la que es transportada la energía hacia la costa, y la energía total por unidad
de ancho proporcionada a la costa en 1 hora por estas olas.
Solución:
a) Como las crestas son paralelas a los contornos, entonces Ho = H'o (no hay refracción). Entonces;
1/2 Eo Co = E C n
que, en términos de las alturas queda como;
H/H'o = [1/2n Co/C]1/2
Pero, dado que;

C/Co = L/Lo = tanh([2πd]/L)
y que
n = 1/2 (1 + [2kd]/[sinh (2kd)])
entonces;
H/H’o = [1/tanh kd] {1/[(1 + 2kh/sinh kd)]}1/2 = Ks
que es el coeficiente de "someramiento".
b) Para la ola en d = 3m,
d/Lo = 3/156 = 0.01923
H/H’o = 1.237
H = 1.237 H'o = 2.47m

c)
P’= 1/2 Eo Co = 1/8ρgH'2o Co = 39 195 N m/s
por metro de cresta de la ola. Esto representa un gasto de 14.11X107 J de energía cada hora sobre
cada metro de la playa.
Asociado con el avance de la ola hay un flujo o transmisión de moméntum. Longuet-Higgins

(1960) define el tensor de radiación como "el exceso del flujo de moméntum debido a la presencia
de las olas". Es lo que sobra de restar la presion hidrostática de la presión absoluta.
En términos sencillos, hay más flujo de moméntum en la dirección de avance de la ola porque U
está en la dirección de propagación sobre la cresta y en dirección opuesta bajo el valle. También la
presión que actúa bajo la cresta es más grande que la que actúa bajo el valle, dando por resultado un
esfuerzo neto sobre un período de la ola.
Si alineamos el eje coordenado ξ con la dirección del avance de la ola y el eje ψ paralelo a las
crestas. Entonces hay 2 componentes diferentes de cero del tensor de radiación: los flujos de
moméntum en ξ y ψ . El tensor de radiación (flujo de moméntum) a través del plano ξ = cte (
paralelo a la costa) en la dirección del avance de la ola está dado por;
Sxx = E ([2kd]/[sen h (2kd)] + 1/2) = E (2n – 1/2)
no obstante que la velocidad orbital ϖ = 0 paralelo a la cresta de la ola, hay también un flujo de
moméntum en la dirección ψ debido a que la presión se aparta de la hidrostática cuando hay olas

presentes. Este flujo de moméntum en ψ, a trav”es del plano ψ = cte, se demostró que está dado
por;
Syy = E ([kd]/[sen hkd]) = E (n – 1/2)
En aguas profundas n = 1/2, por lo que;
Sxx = E/2 ; Syy = 0
en aguas someras n = 1, quedando;
Sxx = 3/2 E; Syy = E/2
El tensor de radiación ha probado ser una herramienta bastante útil en el estudio de gran variedad
de fenómenos oceanográficos. En el contexto de los procesos litorales se ha utilizado para predecir
cambios en el nivel del mar, set up y set down, en la región cercana a la costa y para analizar la
generación de corrientes litorales.
Resumen de la Teoría Lineal
Las ecuaciones que describen el perfil de la superficie, velocidades y aceleraciones de las partículas
para la Teoría Lineal, se presentan en la Figura 3-9.

3.3 Oleaje Irregular
Anteriormente, tratamos las olas como si fueran monocromáticas. Sin embargo, imagenes de radar
y mediciones (Figura 3-9) indican que la superficie del mar se compone de olas cuyas alturas y
períodos varían y que se mueven en diferentes direcciones. Ya reconocida la variabilidad de la

superficie marina, se hace necesario tratar sus características en términos estadísticos. Lo cual
complica el análisis pero describe de manera más real la superficie del mar. El término olas
irregulares será utilizado aquí para denotar estados de mar naturales en los cuales se espera que las
caraterísticas de las olas tengan una variabilidad estadística en contraste con las olas
monocromáticas donde se supone que las propiedades son constantes.
Cuando sopla el viento y en respuesta las olas crecen, la superficie del mar tiende a ser confusa;
existe un rango amplio de olas con alturas y períodos diferentes. Dichas olas se denominan olas de
viento, o algunas veces, simplemente mar. Las olas de período largo que han viajado lejos de su
área de generación son denominadas swell
Figure 3-9. Radar Image of the sea surface in the entrance to San Francisco Bay
La superficie del océano muchas veces es una combinación de muchos componentes de oleaje.
Estos componentes individuales fueron generados por el viento en regiones del océano diferentes
porpagándose hasta el punto de observación. Si un medidor de oleaje estuviera registrando olas en
un punto fijo. El perfil registrado sería algo irregular y aleatorio. Consecuentemente, las
definiciones de altura, período y duración tienen que ser estadísticas e indicarán simplemente la
severidad de las condiciones de olas. Algunos perfiles de oleaje para diferentes condiciones se
presentan en la Figura 3-10. En la Figura 3-11 muestra un registro de oleaje típico. En las Figuras
Figura 3-12 y Figura 3-13 se presentan algunas características importantes de las olas y la
terminología que se estará utilizando.
Se ha encontrado que la altura de la ola significante es un índice bastante útil para caracterizar las
alturas de las olas de la superficie del mar. También definiremos el período pico y la dirección
promedio de la ola, los cuales caracterizan el período y dirección dominantes de las olas.

Figura 3-10. Representaciones de algunas olas oceánicas
Figura 3-11. Perfil de un estado de mar irregular

Figure 3-12. Definición de los parámetros de la ola para un estado de mar aleatorio
Figure3-13. Representación de un proceso aleatorio de oleaje (Ochi 1973)
Existen dos maneras de tratar las olas irregulares: los métodos espectrales (dominio de la
frecuencia) y el análisis ola-por-ola (o del tren de oleaje, en el dominio del tiempo)
3.3.1 Método de Cruce por Cero

El procedimiento ingenieril para el análisis en el dominio del tiempo, se denomina método de cruce
por cero. Aqui, un ola se define cuando la elevación de la superficie cruza hacia arriba la linea de
cero (o la linea de nivel medio del mar) y continua hasta el siguiente cruce hacia arriba. La altura de
la ola de cruce por cero es la diferencia en elevacion entre la cresta más alta y el valle más bajo
entre dos cruces sucesivos. A cada ola le corresponde un período. Las propiedades del campo de
oleaje se determinarán encontrando cantidades estadísticas a partir de los registros individuales.
Los registros tienen que ser suficientemente largos para que la estadística sea confiable. Con la
información obtenida, varios parámetros estadísticos de oleaje pueden calcularse.

Existen muchos parámetros candidatos que pueden utilizarse para definir la estadística de un estado
de mar irregular. Dos de los más importantes son la altura y períodos característicos. La altura
característica puede definirse de muchas maneras; estas incluyen la altura promedio, la altura
cuadrática media, la altura significante, etc. De estas, la más utilizada es la altura significante o
 El período característico podría ser el período promedio, o el período promedio de los cruces de
cero.
Parámetros de estados de mar de corto plazo.
Es bien sabido que cualquier señal periódica η(t) con una media de cero puede separarse en sus
componentes de frecuencia mediante el Análisis de Fourier. Los registors periódicos de oleaje
generalmente pueden ser tratados como procesos aleatorios regidos por leyes de la teoría de la
probabilidad. Si el registro de oleaje es una señal aleatoria, entonces utilizamos el término olas
aleatorias. Para muchos propósitos los registros de oleaje pueden considerarse como aleatorios
(Rice 1944-1945, Kinsman 1965, Phillips 1977, Price and Bishop 1974).
Las propiedades estadísticas de una señal como la del oleaje pueden obtenerse de un conjunto de
muchas observaciones simultáneas llamado ensamble o conjunto de señales {η1(t), η2(t), η3(t),
η4(t),…}, pero no de un solo registro. Para determinar las propiedades del oleaje a aprtir del
proceso η(t), tienen que hacerse ciertas suposiciones a saber:
1) el proceso debe ser estacionario. Si, en adición, η(t), se mide en diferentes localidades y sus
propiedades son invariantes, o no dependen del lugar de medición, entonces se dice que el proceso
es homogéneo.
2) se supone que el proceso η1(t), es ergódico, o sea, que cualquier registro tomado del proceso
η(t), es típico de las otras posibles realizaciones, y por tanto, el promedio de un solo registro del
ensamble es lo mismo que el promedio a través del ensamble.
Las suposiciones de ergodicidad y estacionariedad son la columna de la estadística de oleaje a

partir de mediciones. Para aplicar estos conceptos, considere un ensamble de registros que
representan un cierto estado de mar η(t), sobre un tiempo finito T. El valor medio o esperado del
estado de mar, denotado por μη ó E[η], se define como;
μη = Ε[η] = τ-1∩η(t) dt
donde los límites de integración van de –τ/2 a τ/2 y Ε denota el valor esperado de η(t).
Similarmente, la media cuadrática de η(t) corresponde al segundo momento de η(t), denotado por
Ε[η2]. La desviación estándard, ση, o raíz cuadrática media es la raiz de ese valor (la varianza). La
varianza σ2η puede expresarse como;
σ2η = Ε[η2] - μ2η

Ahora vamos a enfocarnos en definir las propiedades probabilísticas de las olas. La representación
probabilística de un cierto estado de mar es útil por dos razones: primero, permite al ingeniero
encontrar los parámetros de oleaje dentro de un cierto límite que llevará a un nivel aceptable de
riesgo. Segundo, un diseño basado en criterios probabilísticos puede resultar en ahorros

sustanciales de costos mediante la consideración de las incertidumbres en las estimaciones de
oleaje.
Para que los ingenieros usen efectivamente las olas irregulares en el diseño, las propiedades de los
registros de oleaje tienen que seguir algunas leyes de probabilidad, de forma tal que las estadísticas
puedan obtenerse analíticamente. Rice (1944-1945) desarrolló la teoría estadística de señales
aleatorias en el análisis del ruido eléctrico. Longuet-Higgins (1952) aplicó esa teoría a las olas del
mar para describir su estadística. Longuet-Higgins encontró que los parámetros de una señal de
oleaje aletorio siguen leyes de probabilidad conocidas.
Las dos funciones de probabilidad más utilizadas en el estudio de las olas, son la distribución
Gausiana (Figura 3-14) y la de Rayleigh (Figura 3-15). La distribución Gausiana se utiliza para
describir la propiedades de corto plazo de la elevación superficial η. Su densidad de probabilidad
está dada por;
π(ξ) =1/σξ(2π)1/2 ε −[(ξ−μξ)2/2σ2ξ]

donde μ es la media de ξy σξla desviación estándard. La integral de π(ξ), denotada por Π(ξ) es la
distribución de porbabilidad. No existe un forma cerrada para esta expresión, por lo que esta
distribución se presenta tabulada en función de μ yσ como;
Π(ξ)=Φ[(ξ−μξ)/σξ]
Figura 3-14. La densidad de probabilidad Gaussiana y la distribución de probabilidad acumulativa

Figura 3-15. La distribución de densidad de probabilidad de Rayleigh y la distribución de
probabilidad acumulativa
Para µx = 0 y σξ=1, π(ξ) y Φ(ξ) se reducen a;
π(ξ) =1/(2π)1/2 ε −[(ξ)2/2]

Φ(ξ) =∩ p(y) dy
donde la última integral es la función de error.
La probabilidad de excedencia Q(x) puede expresarse en términos de la probabilidad de np-

excedencia P(x) como;
Q[x(t)>x1]= 1 – P[x(t)<x1] = 1 Φ[(ξ−μξ)/σξ]
Que es la probabilidad de que x excederá x1 sobre el período de tiempo t, y se muestra como el área
sombreada al fondo de la Figura II-1-28. La probabilidad de excedancia es un parámtero de diseño
importante en proyectos basados en riesgo de diseño.
En la práctica, se trabaja generalmente con alturas en lugar de elevaciones. Longuet-Higgins (1952)

encontró a partir de teoría estadística, que las amplitudes y la alturas de las olas siguen la
distribución de Rayleigh presentada en la Figura 3-15. La función de densidad de probabilidad de
la distribución de Rayleigh asi como su probabiliad acumulativa, están dadas por;
π(ξ) =πξ/(2μ2ξ) ε −π/4[(ξ/μξ)2] para ξ≥0

Π(ξ) =1−ε −π/4[(ξ/μξ)2] para ξ≥0

donde μξ es la media.
Figura 3-16. Histogramas de las alturas normalizadas (a) y los períodos (b) de las olas junto con las
distribuciones teóricas (Chakrabarti 1987)
Distribución de Alturas
Las alturas de las olas indivduales pueden considerarse como una variable estocástica representada
por una función de distribución de probabilidad. Esta función puede obtenerse del histograma de

alturas normalizado con las alturas medias en varios registros medidos en un punto (Figura 3-16).
Thompson (1977) indicó que tan bien los registros costeros de olas siguen esta distribución de
Rayleigh. Si la energía de las olas se concentra en un rango angosto de períodos, el máximo del
perfil de la ola coincidirá con las crestas y el mínimo con los valles. Esta condición se denomina
condición de banda angosta. Bajo esta condición, las alturas están representadas por la siguiente
distribución de Rayleigh (Longuet-Higgins 1952, 1975b, 1983)
p(H) = 2H/H2rms e[-H2/H2rms]
Π(H) = 1- e[-H2/H2rms]
La altura significante H1/3 es el centroide del área para Η ≥ Η* bajo la función de densidad, donde
H > H* corresponde a las olas en el rango del tercio de las olas más grandes como se muestra en la
Figura II-1-29, o sea,
Π(Η*) = 1 - 1/3 = 1 - e[-H*2/H2rms]
de la cual, encontramos H* = 1.05Hrms. De la integración de la ecuación anterior, podemos obtener

varias estimaciones utilizando ciertas propiedades de la Funcion Error (Abramowitz and Stegun
1965). Se encuentra que;
H1/3 ≈ 4.00 (m0)1/2 = 1.416 Hrms
H1/10 = 1.27 H1/3 = 1.80 Hrms = 5.091 (m0)1/2
H1/100 = 1.67 H1/3 = 2.36 Hrms = 6.672 (m0)1/2
Hmax = 1.86 H1/3 (para 1000 ciclos de ola en el registro)
La altura máxima más probable, Hmax, en un registro que contiene N olas se relaciona con la altura
Hrms (Longuet-Higgins 1952) mediante;
Hmax = {(logN)1/2 + 0.2886/(logN)1/2 - 0.247/(logN)3/2}Hrms
Este valor de Hmax puede ser proyectado a períodos de tiempo largos ajustando el valor de N
basado en el período de cruce por cero (Tucker 1963).
Las relaciones teóricas derivadas de la distribución de Rayleigh concuerdan generalmente con los
valores determinados directamente de los registros. La distribución de probabilidad de Rayleigh se
compara con un histogramas de alturas medidas en la Figura II-1-30 (Chakrabarti 1987).
Distribución de los períodos de la ola.
Esta distribución es muy similar a la distribución normal con un período medio dado por;
T0,1 = m0/m1

La densidad de probabilidad del período de la olas está dado por (Bretschneider 1969);
p(T) = 2.7 T 3/Tavg e[-0.675τ4]
τ = T/Tavg
Longuet-Higgins (1962) presenta una distribución de densidad de probabilidad diferente;

2
p(τ) = 1/2(1 + τ )3/2
τ = (T - T0,1)/νT0,1
ν = (m0m2 - m21)/ m21
donde ν es el parámetro de ancho espectral y m0, m1, and m2 son los momentos espectrales. Esta
función de densidad de probabilidad es simétrica alrededor de τ = 0 donde es un máximo, y es
similar a la distribución normal con media igual a T0,1. Varios períodos característicos se
relacionan. De manera general, esta relación está dada por;
Tmax ≈ T1/3 ≈ ϒ Tavg
Donde el coeficiente ϒ varía entre 1.1 y 1.3.
3.3.2 Análisis Espectral.

En el período de 1950-1960, el trabajo de Rice(1944-1945) en el procesado de las señales fué
extendido al estudio de las olas (Kinsman 1965; Phillips 1977). En principio, se reconoció que la
series de tiempo de las elevaciones era similar que la de un registro de ruido. Suponiendo que esa
es una muestra discreta de un proceso contínuo, entonces, los principios del análisis de Fourier
podrían extenderse para describir el registro. En esta sección discutiremos las aproximaciones para
usar las representaciones espectrales en ingeniería, y describiremos como esa información puede
ser utilizada. La teoría estadística y las suposiciones solo serán mencionadas y los detalles de las
derivaciones solo serán referenciados.
Si consideramos la serie de tiempo de la elevación de la superficie en un solo punto tal como la de

las Figuras II-1-25, II-1-31, y II-1-32, el análisis espectral considera al registro como la variación
de la superficie alrededor de la media, y reconoce que esta variación consiste de varias
periodicidades. A diferencia del análisis en el dominio del tiempo, que busca definir olas
individuales, el análisis espectral busca definir describir las distribución de la varianza con respecto
a la frecuencia de la señal. Por convención, la distribución de la varianza con la frecuencia se
escribe como E(f) o S(f) bajo la suposición de que la función es contínua en el espacio de las
frecuencias. E(f) a menudo se denomina como el espectro de frecuencia o el espectro de energía
(suponiendo Teoría Lineal) porque la energía del campo de oleaje puede estimarse multiplicando
E(f) por ρg.

La superfice puede visualizarse no como olas individuales sino como una superficie tridimensional,
que experimenta un desplazamiento alrededor de una media y cuya varianza es periódica en tiempo
y espacio. La representación espectral más sencilla es considerar E(f,θ), que representa como se
distribuye la varianza con la frecuencia y dirección (Figura II-1-33). E(f,θ), es llamado el espectro
bi-dimensional o espectro energía direccional porque al multiplicar por ρg obtenemos la energía de
las olas. La ventaja de esta representación es que le dice al ingeniero acerca de la dirección en que
la energía se está moviendo. Un ejemplo de este tipo de espectros se muestra en la Figura II-1-34.
Figura 3-17. Surface elevation time series of a regular wave and its spectrum
(Briggs et al. 1993)

Figure 3-18. A schematic for a two-dimensional wave spectrum Ε(φ,θ)
Suposiciones
Primero, los registros de oleaje son de longitud finita (generalmente 17-68 min de largo) y son
tomados de registrar la elevación superfial en intervalos de muestro discretos (típicamente 0.5-2.0
s). También es deseable que no haya tendencia en los datos.
Figure II-1-32. Surface elevation time series of an irregular wave and its spectrum
(Briggs et al. 1993)
Figure II-1-34. A directional spectrum and its frequency and direction spectrum
(Briggs et al 1993)
Descripción del Análisis Espectral
El análisis espectral es esencialmente una transformación de la serie de tiempo de las elevaciones

del dominio del tiempo al dominio de la frecuencia. La transformación se realiza mediante una
herramienta matemática conocida como la técnica de la Transformada Rápida de Fourier ó FFT
(Cooley and Tukey 1965). Aqui, primero trataremos el análisis de un registro en un punto para
obtener su espectro y después trataremos el caso de como obtener el espectro direccional.
La densidad espectral de energía de la ola, E(f), o simplemente el espectro de energía puede

obtenerse directamente de una serie de tiempo contínua de la superficie del mar, η(t) con la ayuda
del Análisis de Fourier. Utilizando el Análisis de Fourier, el perfil original puede descomponerse en
una suma infinita de senoides de amplitud An, frecuencia ωn , y fase relativa εn, o sea,

η(t) = ∝ An cos (ωn t - εn)
= ∝ an cos nωt + bnsin nωt
Los coeficientes an y bn pueden determinarse explícitamente a partir de las propiedades de

ortogonalidad de las funciones circulares. Nótese también que a0 es la media del registro. Como las
observaciones reales son de longitud finita, se usa la tranformada discreta de Fourier y el numero
de términos, n, en la sumatoria es finito.
La covarianza de η(t) se relaciona con el espectro de energía. Esta se define en términos de los
cuadrados de las amplitudes como;
η2(t) = ∝ A2n Δφ
A2n = ½ (a2n + b2n)1/2
εn = tan-1(an/ bn)
Una estimación del espectro de energía de η(t) puede obtenerse por
E(f) = Tr -1{∝η(nΔt)ε2πφ(νΔτ) Δτ}2
donde Tr es la longitud del registro y Δt es el intervalo de muestreo.
En la práctica, la longitud total se divide en M segmentos con igual numero de puntos N. Dejando
que N sea un potencia de 2 para eficiencia computacional, el resultado es entonces promediado
sobre las M secciones. En el análisis de la FFT, las variables M, N, y t tienen que seleccionarse
independientemente, pero Tr y t se mantienen fijos de tal forma que el numero total de puntos se
obtenga de esos valores. Tradicionalmente, los valores más comunes de N van entre 512 a 2,048,
mientras que para M generalmente es 8 o mayor. Como Tr depende de N o M, y t como Tr = M N t,
entonces mayores valores de N y M darán, en general, mejor resolución y alta confianza en la
estimación del espectro (Cooley and Tukey 1965; Chakrabarti 1987).
Consideremos el perfil de la superficie con una ola de una sola amplitud y frecuencia dada por;
η(t) = a sinωt
La varianza de esta ola sobre un período de 2π, es
σ2 = <η2> = 1/2π ∩ a2 sin22πφτ d(2πφτ)

= a2/2 = 2 ∩E1(f) df = ∩E2(f) df
O sea, la cantidad a2/2 representa la contribución a la varianza asociada con el componente ω =

2πφ (Figura II-1-35). Para consistencia entre estas cantidades se requiere que el espectro no incluya

el término ρg. La diferencia entre el espectro de un lado, E1, y el de dos lados (E2) ilustrada en la
Figura II-1-36 es importante.
Figure II-1-35. Sketches of wave spectral energy and energy density (Chakrabarti 1987)
Ahora, intuitivamente podemos extender este caso para calcular la varianza de una señal aleatoria
con media cero compuesta de las contribuciones de muchas frecuencias. Usando las relaciones
anteriores, tenemos;
σ2 = ∑ an2/2 = ∩E(f) df = m0
donde m0 es el momento espectral de orden cero. Físicamente, m0 representa el área bajo la curva
de E(f). Los momentos espectrales de orden superior pueden obtenerse de;
mi = ∩ fiE(f) df ; i = 0,1,2,...
Para alturas de olas distribuidas siguiendo una distribución de Rayleigh, Hs puede aproximarse por
(Longuet-Higgins 1952)
Hs = 3.8 (m0)1/2 ≈ 4 (m0)1/2
Por lo tanto, el momento espectral de orden cero, m0, que es el área bajo el espectro de energía,
define la altura significante para un E(f) dado (Figura II-1-37).

Figure II-1-37. Energy density and frequency relationship (Chakrabarti 1987)
El espectro de oleaje y sus parámetros
Dos parámetros que se usan frecuentemente en la distribución de probabilidad de la ola son el

ancho espectral, ν, y el ancho de banda espectral, ε, los cuales se utilizan en determinar la
estrechéz del espectro. Los valores de estos parámetros están entre 0 y 1, y pueden aproximarse en
términos de los momentos espectrales como;
ν = (m0m2/m2 – 1)1/2
ε = (1 - m22/m0 m4) 1/2
Para un espectro de banda angosta, los dos ν y ε tienen que estar cercanos a cero (Figure II-1-38).
Por ejemplo, para los dos espectros más comunes, el de Pierson-Moskowitz (PM) (Pierson and
Moskowitz 1964) y el JONSWAP (Hasselmann et al. 1973), which are discussed in the next section,
ν = 0.425 y 0.389, respectivament, with ε = 1 para ambos. Las olas del oceáno tienen un espectro
de banda amplia.

Figure II-1-38. Comparison of the PM and JONSWAP spectra (Chakrabarti 1987)
Las relaciones aproximadas para los parámetros de oleaje más comunes son;
HS = 4.0; H1/10 = 5.1(m0)1/2
Tz = (m0/m2)1/2; Tc = (m2/m4)1/2
<η> = (m0)1/2
Thompson yVincent (1985) investigaron como H1/3 y Hm0 varían en aguas may someras cerca de la
rompiente. Ellos encontraron que el cociente H1/3/Hm0 varió sistemáticamente a través de la zona de
surf, aproximándose a un máximo cerca de la rompiente. Sus resultados se presentan en forma de
nomograma en la Figura II-1-40.
La cinemática y la dinámica de las olas irregulares se estiman analíticamente de tres maneras a) en

el dominio de tiempo, mediante teoría lineal); b) en el dominio de la frecuencia, mediante el
espectro y teoría lineal (Dean and Dalrymple 1991); c) mediante técnicas de simulación de olas sintéticas
aleatorias (Borgman 1990). De los tres métodos, este último produce los resultados más reales, pero es el
más complejo.

Figure II-1-40. Variation of Hs/Hmo as a function of relative depth ¯d and significant
steepness (Thompson and Vincent 1985)
3.4 Pronóstico de Oleaje
3.5 Modificaciones del Oleaje
3.5.1 Introducción
La ingeniería costera considera problemas cerca de la línea de costa, generalmente en aguas con
profundidades menores que 20 m. Los diseños de proyecto generalmente requieren del
conocimiento del campo de oleaje sobre un área de 1-10 km2 en la cual la profundidad del agua
puede variar significantemente. Adicionalmente, el estudio de los cambios de la linea de costa y de
protección de playas requieren frecuentemente del análisis de los procesos costeros sobre cendas
litorales enteras, que pueden abarcar de 10-100 km de longitud. Los datos de oleaje generalmente
no están disponibles en el lugar o profundidades requeridas. Muchas veces, el ingeniero costero
encontrará que los los datos se han colectado o pronosticado de lugares lejanos mar adentro o de
lugares cercanos a profundidades similares. Este capítulo delinea procedimientos para transformar
las olas desde esos lugares hasta el sitio de interés del ingeniero.
El entendimiento de los procesos que afectan a las olas en la costa es esencial para la ingeniería
costera. Las olas que se propagan a través de aguas someras están influenciadas fuertemente por la
batimetría y por las corrientes. Cambios batimétricos pueden provocar cambios significantes en la
altura y dirección de viaje de la ola. La magnitud de estos cambios es sensitiva particularmente al

período, dirección de la ola, y a como la energía de la ola se dispersa de acuerdo a su dirección y
frecuencia (Figura II-3-2).
La altura de la ola es, generalmente, el factor más siginificante que influencía un proyecto. El
diseño utilizando una altura que es excesivamente conservativa puede aumentar el costo de un
proyecto y puede hacerlo ineconómico. De manera contraria, subestimar la altura de la ola puede
resultar en una falla catastrófica del proyecto o costos de mantenimiento significantes. Las maneras
de transformar las olas son numerosas y difieren en complejidad y exactitud. Consecuentemente,
los estudios de transformación reuqieren de análisis cuidadoso.
La transformación de las olas sobre batimetrías irregulares es compleja. Suposiciones de

simplificación admiten aproximaciones válidas y útiles para estimar olas en la región cercana a la
costa. En la sección siguiente se da un repaso de las bases teóricas que rigen los análisis de
transformación del oleaje, seguida del desarrollo de un método sencillo para estimaciones de
refracción y someramiento. Después se discute la transformación de olas irregulares, seguida de
una discusión de los modelos de transformación usados actualmente por el CERC.
El propósito de este capítulo es el de ofrecer métodos para estimar olas en un cierto lugar dada
información de otro sitio.
Figure II-3-2. Amplification of wave height behind a shoal for waves with different
spreads of energy in frequency and direction
3.5.2 Principios que rigen la modificación del Oleaje
En esta sección se presentan los principios científicos que rigen la transformación de las olas desde
aguas profundas hasta someras. Desafortunadamente, el problemas es tan complejo que los cálculos
detallados requienen el uso de modelos numéricos complicaods cuyos antecedentes e

implementación caen fuera de los objetivos de ste curso. Este capítulo ofrece los principios d la
transformación del oleaje, una visión simplificada, asi como una introducción a tres de los modelos
más ampliamente utulizados en el CERC.
Los procesos que afectan a las olas conforme se propagan desde aguas profundas a someras,
incluyen:
(a) Refracción.
(b) Someramiento.
(c) Difracción.
(d) Disipación por fricción.
(e) Disipación por percolación.
(f) Rompimiento.
(g) Crecimiento adicional debido al viento.
(h) Interacción ola-corrientes.
(i) Interacciones ola-ola
Los primeros tres son efectos de propagación porque resultan de la convergencia o divergencia de
las olas provocada por la forma de la topografía de fondo, que influencía la dirección de avance de
la ola y provoca que la enería se concentre o se desparrame. La difracción ocurre debido a
estructuras que obstaculizan u obstruyen la propagación de las olas. Los siguientes tres efectos son
mecanismos sumidero porque remueven energía del campo de oleaje a través de la disipación. El
viento es un mecanismo fuente porque representa una adiciíon de energía para la ola. La presencia
de campos de corrientes de gran escala puede afectar la propagación y disipación de las olas. Las
interacciones ola-ola resultan del acople alineal de los componentes de la ola y resultan en una
transferencia de energía de unas olas a otras.
Los procedimentos de esta sección abarcan desde aguas profundas hasta antes de la zona de surf.
Los procesos en esta última son tratados en la sección: “Surf Zone Hydrodynamics.”
La ecuación de la transformación de la ola
El problema general de la transformación se presenta en términos del concepto de espectro

direccional, Ε( , , , ,θ) donde y θ representan un componente de frecuencia y dirección
particular, , representan un lugar geográfico en el espacio, y representa el tiempo. Se considera
que las olas se estan propagando sobre una región en que la profundidad cambia y que no hay
corrientes. La ecuación utilizada, esta dada por;
∂E/∂ t + ∇•[Cg E)] = Sw+ Sn+ SD+ SF+ SP

A B C D E F G
Aunque multidimensional, esta ecuación es fundamentalmente sencilla.
A es la tasas de cambio temporal del espectro

B representa la propagación de la energía de la ola
C representa la contribución del viento

D representa la re-distribución de energía entre diferentes componentes de ola debido a
alinealidades
E es la disipación debido a la rompiente
F son las pérdidas por fricción
G son pérdidas por percolación
Una discusión más detallada de la mecánica espectral del oleaje puede encontrarse en Leblond y
Mysak (1978), Hasselmann (1962, 1963a, 1963b), Hasselmann et al. (1973), Barnett (1968),
Phillips (1977), Resio (1981), WAMDI (1988), etc.
Los movimientos superficiales de las olas producen un campo de velocidad que se extiende a cierta
profundidad en la columna d eagua. Esta profundidad para una ola de aguas profundas, es L/2. Si la
profundidad es menor que L/2, el movimiento se extiende hasta el fondo. En los casos en que el
movimiento de la ola interactúa con el fondo, ocurren varios cambios físicos: la celeridad C y la
velocidad de grupo Cg cambian, lo mismo que L. Si las olas se están propagando sobre un fondo
variable (y suficientemente somero como para que las olas interactúen con el fondo), entonces los
cambios en velocidad cambian la dirección de viaje de la ola asi como su amplitud (refracción y
someramiento). Si los patrones de propagación llevan a una concentración fuerte de olas, energía
puede radiarse lejos de la zona de convergencia mediante difracción de la olas (Penny and Price
1944; Berkhoff 1972). La interacción de las olas con el fondo produce una capa límite, que dará
como resultado una pérdida de energía debido a fricción con el fondo (Término F) debido al
material de fondo y formas del fondo (Bagnold 1946). Si el fondo es razonablemente poroso, la
presión asociada al campo de presión puede inducir un flujo hacia dentro y fuera del piso
(Bretschneider and Reid 1953), que da como resultado pérdida de energía por percolación (Término
G). Si el piso es lodoso o viscoelástico pueden ocurrir otros procesos (Forristall and Reece 1985).
Estos términos casi no se aplican porque muchas veces no se tiene información adecuada de la
composición del material de fondo.
Tipos de Transformaciones.
Tres casos clásicos de transformación del oleaje describen la mayoría de las situaciones
encontradas en ingeniería costera:
(a) Una gran tormenta genera olas de aguas profundas que se propagan a través de aguas someras
mientras que las olas continúan creciendo debido al viento.
(b) Una gran tormenta genera vientos en un área lejana al sitio de interés y las olas se propagan a
aguas someras como swell.
(c) El viento sopla sobre un área somera y genera olas que crecen tan grandes que interactúan con
el fondo, o sea no hay propagación de olas de aguas profundas al lugar.
Todos los casos son importantes, pero el primero y el tercero son relativamente complejos y
requieren de modelos numéricos para un tratamiento razonable. El segundo caso, es clásico y ha
servido de base para muchos estudios de ingeniería costera. Muchas veces el swell es aproximado
por olas monocromáticas, y métodos sencillos de refracción y someramiento se utilizan para
estimaciones del oleaje en la zona cercana a la playa. Dado que los procesos de refracción y
someramiento son importantes en ingeniería costera, en la siguiente sección nos dedicamos a

obtener algunas aproximaciones sencilllas que ilustran la necesidad de aproximaciones más
complejas.
3.5.3 Refracción y someramiento
Para entender estos fenómenos, consideremos el caso de una ola monocromática estacionaria que se
propaga en una región en la cual la línea de costa es recta y los contornos de fondo son
equiespaciados y aquella (Figure II-3-3). Tampoco hay corrientes presentes. Si una cresta de ola
tiene inicialmente un ángulo de aproximación a la costa diferente de 0o, entonces parte de la ola
(punto A) estará en aguas más somera que las otra parte (punto B). Dado que la profundidad en A,
hA , es menosr que en la profundidad en B, hB, entonces la velocidad de la olas en A será menor que
en B porque
CA = g/ω tanh k hA < g/ω tanh k hB = CB
La velocidad diferencial a lo largo de la cresta de la ola provoca que la cresta se haga más paralela
a la costa. El problema de propagación es, entonces, el de graficar la dirección de aproximación de
la ola y calcular su altura conforme la ola se propaga de aguas profundas a someras. En el caso de
olas monocromáticas, el período se mantiene constante. En el caso de olas irregulares, el proceso de
transformación puede afectar a las olas en cada frecuencia de manera diferente; como
consecuencia, el período pico del campo de olas puede recorrerse.
El problema de propagación de las olas pueden visualizarse mediante la construcción de rayos de

oleaje. Un rayo es el camino trazado por una ortogonal a la cresta de un tren de olas, conforme
estas se propagan hacia la costa. Por tanto, un grupo de rayos mapea la trayectoria de viaje de la
cresta de las olas. Para batimetrías sencillas, estos grupos de rayos pueden construirse a mano,
aunque este es un procedimiento tedioso. También existen progamas gráficos de computadora que
automatizan este proceso (Harrison and Wilson 1964, Dobson 1967, Noda et al. 1974), pero en
gran medida esas aproximaciones se han visto suplantadas por los métodos numéricos. El análisis
de refracción y someramiento típicamente tratan de especificar la altura y dirección de la ola a lo
largo de un rayo. La Figura II-3-4 presenta gráficas idealizadas de rayos de oleaje para varias
batimetrías típicas.

Figure II-3-3. Straight shore with all depth contours evenly spaced and parallel to
the shoreline

Figure II-3-4. Idealized plots of wave rays

Batimetría con contornos rectos y paralelos
Primero, desarrollamos la ecuación que especifica como el ángulo de la ola cambia a lo largo de un
rayo, seguida por una ecuación para la altura. La derivación es para contornos rectos y paralelos sin
corrientes presentes. El componente del sistema coordenado en es ortogonal a la línea de costa;
mientras que el se toma paralelo a la costa. Esto implica que cualquier derivada en es nula
porque dh/dy es cero.
Para olas monocromáticas, la función de fase
Ω(x,y,t) = (k cos θ + k sin θ - ωt) (3)
puede utilizarse para definir el vector numero de onda Κ como
Κ = ∇Ω (4)
Como Κ es un vector, uno puede tomar el rotor de Κ
∇ × Κ=0 (5)
que es cero, porque Κ por definición es el gradiente de un escalar y el rotor de un gradiente es cero.
Haciendo la operación indicada en (5), queda
∂( θ)/∂ −∂( θ)/∂ = 0 (6)
Dado que las derivadas en son cero, y utilizando la relación de dispersión, que liga κ con (y
notando que k =2π/CT y que el período es constante), tenemos
∂( θ/ )/∂ =0 (7)
o
θ/ = . (8)
Sea C0 la celeridad de la ola en aguas profundas. En aguas profundas, θ0/ 0 es conocido si el

ángulo de aproximación es conocido, por tanto; a lo largo de un rayo, la ecuación (8) da
θ/ = θ0/ 0 (9)
que es equivalente a la ley de Snell en óptica

Figure II-3-5. Wave-height variation along a wave ray
La variación en altura de la ola a lo largo de un rayo puede estimarse si consideramos dos rayos
muy cercanos (Figura II-3-5).En aguas profundas, el flujo de energía (ECn), que también es ECg, a
una distancia 0 a través de la cresta de la ola puede estimarse mediante (ECn)0b0. A una distancia
corta a lo largo del rayo, el flujo de energía es (ECn)1b1. Dado que los rayos son ortogonales a la
cresta de la ola, entonces no debe de haber transferencia de energía a través de los rayos, y el
principio de conservación da;
(ECn)0b0= (ECn)1b1
Pero, para una ola monocromática,
E = (1/8) ρgH2
Por lo que H1 y H0 se relacionan por
H1 =H0(Cg0/Cg1)1/2 (b0/b1)1/2
Que generalmente se escribe como;
Η1 = Η0 Κ Κ
Donde Ks es el coeficiente de someramiento y Kr es el coeficiente de refracción.

Los coeficientes combinados de refracción y someramiento para olas monocromáticas sobre
contornos rectos y paralelos, se presentan en la Figura II-3-6. Esto está automatizado en el código
ACES (Leenknecht, Szuwalski, and Sherlock 1992) y en el programa NMLONG (Kraus 1991). La
figura tmabién nos da el ángulo local de la ola θ y KrKs en términos del ángulo en aguas profundas
y d/gT2. Aunque en la mayoría de los casos la batimetría es algo más complicada que esto, estos
procedimientos nos dan una estimación rápida de los ángulos de aproximación de las olas.
Ejemplo
Dada una ola con 1 m de altura y período de 15s en aguas de 500 m de profundidad, con una playa
plana de pendiente suave. Encontrar;
La altura de la ola H y el ángulo θ a profundidades de 200, 100, 90, 80, 70, 60, 50, 40, 30, 20, 10,
16, 14, 12, 10, 8, 6, y 4 m para ángulos θ0 de aguas profundas, 15o, and 45o.
Solución
La solución puede obtenerse directamente utilizando el programa ACES, o gráficamente utilizando
la figura II-3-6. La tabla II-3-1 muestra los resultados obtenidos utilizando el programa ACES.
Para una ola con una profundidad de 10 m y un ángulo inicial de 45o, la altura y el ángulo se
calculan como sigue:
Dado que 0 para una ola de 15s es
0= 1.56(15)2=351
y como 500 m es mayor que L0/2, que es la condición inicial de aguas profundas dada. La nueva
longitud de la ola a 10m tiene que estimarse de
= ( Τ2/2π) (2π / ) = Τ/ω (κ ) (II-1-10)
y resulta ser igual a 144 m (see Problem II-1-1).
Ks puede estimarse de
Ks = (Cg0/Cg1)1/2
En aguas profundas para una ola de 15s, Cg0 es;
1/2 0= 1/2 (1.56Τ) =23.4/2 = 11.7 /
En tanto que Cg está dada por

= = 1/2©1+2κ / 2κ ∏ Τ/2π (κ )
Insertando los valores de d = 10 m, L = 144 m, T = 15 s, y g = 9.8 m/s2 da = 9.05 m/s. Por lo que
Ks =(11.70/9.05)1/2 = 1.14
La solución para Kr involucra

Kr = ©(1-sin2θ0)/(1-sin2θ1)∏1/4
En aguas profundas, θ0 = 450, de la ecuación (9),
θ=( 1 θ0)/ 0
Pero, C0 = 1.56T = 23.4 m/s. En aguas de 10 m, C1 = Λ1/Τ = 144 m/15s = 9.60 m/s.
θ = (9.6 450)/23.4 = [(9.6)(7.07)]/23.4 = 0.29
Κ = ©(1−(0.707)2)/(1−(0.29)2)∏1/4=(0.5/0.91)1/4 = 0.86
Por tanto, Hi = H0 Ks Kr = 1(1.14) (0.86) = 0.98 m.
El ángulo de aproximación es arc ( θ) = 16.80. O sea, la ola de l-m, 15-s ha cambiado un 2% en

altura por 28.20 en el ángulo de aproximación.
Como puede verse en la Tabla II-3-1 las grandes diferencias provocadas por la refracción y el
someramiento se verán a profundidades bajas.
Tabla II-3-1. Resultados del Ejemplo

Prof
θ Η θ Η θ Η
500 0 1.00 15.0 1.00 45.0 1.00
400 0 1.00 15.0 1.00 45.0 1.00
300 0 1.00 15.0 1.00 45.0 1.00
200 0 1.00 15.0 1.00 45.0 1.00
100 0 0.94 14.3 0.94 42.4 0.92
90 0 0.93 14.0 0.93 41.2 0.91
80 0 0.93 13.7 0.92 30.4 0.89
70 0 0.92 13.2 0.91 38.9 0.88
60 0 0.91 12.7 0.91 37.0 0.86
50 0 0.91 12.0 0.91 34.5 0.85
40 0 0.92 11.1 0.92 31.8 0.84
30 0 0.95 9.9 0.94 28.1 0.85
20 0 1.00 8.4 0.99 23.4 0.88
18 0 1.02 8.0 1.01 22.3 0.89
16 0 1.04 7.8 1.03 21.1 0.91
14 0 1.07 7.1 1.05 19.8 0.92
12 0 1.10 6.6 1.08 18.4 0.95

10 0 1.14 6.1 1.12 16.8 0.98
8 0 1.19 5.5 1.17 15.1 1.02
6 0 1.27 4.8 1.25 13.15 1.08
4 0 1.39 3.9 1.37 10.8 1.18
Batimetría Irregular
La topografía irregular, puede provocar que las olas se refracten de manera compleja, de forma tal
que producen variaciones en la altura de la ola y en su energía. La convergencia y divergencia
producida por la refracción es importante en decidir la construcción de estructuras costeras. Un
tratamiento adecuado, requiere que se considere todo el espectro; o sea, es necesario obtener varios
diagramas de refracción que cubran todo el rango de períodos y direcciones dentro del espectro.

Figure II-3-6. Solution nomogram
Figure 3-19. Refraction coefficients from random sea waves on a coast with straight and
parallel bottom contours (Goda 1985)

Figure 3-20. Variation of predominant wave direction for random sea waves on a coast with
straight and parallel bottom contours (Goda 1985)
3.5.4 Otros procesos que contribuyen
Difracción
Es la transferencia lateral de la energía de la ola debido a variaciones en altura a lo largo de una

cresta. La difracción es una transformación específica del lugar. Las olas pueden difractarse debido
a estructuras como promontorios, escolleras, rompeolas, o por la topografía de fondo como bajos y
arrecifes. Si alguna existe alguna de estas estructuras cerca del lugar de diseño, entonces los efectos
de difracción tienen que considerarse. En el SPM (1984) se presentan algunos resultados gráficos
de refracción para olas monocromáticas a profundidad constante alrededor de espigones. Gráficos
similares para mares irregulares se presentan más adelante.
Existen modelos que combinan los efectos de refracción y difracción, estos se basan en el método
de la ecuación parabólica (Berkoff 1972).
Disipación de Energía
Es la pérdida de energía de la ola conforme se propagan hacia la costa. Este es el resultado de la

viscosidad, turbulencia, fricción con el fondo, y el movimiento de los sedimentos de fondo
inducidos por las olas. La importancia de este efecto depende del lugar. Modelos numéricos de
propagación, como aquello utilizados en el programa WIS, han incorporado este efecto.
Rompimiento de la Ola

Este es uno de los procesos de transformación más importantes en la determinación de las
corrientes litorales y el transporte de sedimento. Desafortunadamente, este proceso no está bien
entendido. Modelos tales como los desarrollados a partir del trabajo de Longuet-Higgins y Cokelet
(1976) son prometedores en la cuantificación de la cinemática de las olas cerca de la rompiente. Sin
embargo, esos modelos no han evolucionado al nivel de utilizarse en el diseño rutinario. Se han
desarrollado otros modelos que se basan en el flujo de moméntum o de la energía (Peregrine y
Svendsen 1978; Thornton y Guza 1983; Dally, Dean, y Dalrymple 1984). El modelo de Dally está
incorporado en el RCPWAVE.
Uno de las prácticas bastante utilizadas para estimar condiciones de olas rompientes lo constituye el
uso de curvas empíricas (ver Figuras 3-21 y 3-22). Las olas rompientes a menudo son clasificadas
como spilling, plunging, or surging (Figura 3-23). La figura 3-21 muestra las condiciones para las
cuales ocurren este tipo de rompientes.
Figure 3-21. Breaker height dependency on wave steepness and bottom slope

Figure 3-22. Breaker depth dependency on wave steepness and bottom slope
Figure 3-23. Breaking wave profiles

3.5.5 Transformación de olas irregulares
La sección precedente enfatizó la refracción y someramiento de olas monocromáticas. Cuando este
proceso se aplica a una ola inicial de altura y período significante, se denomina análisis de la ola
significante. Este análisis puede ser adecuado cuando las condiciones de oleaje tienen un espectro
de ancho espectral amplio, o picos espectrales múltiples. Casos donde el ánalisis de la ola
significante es adecuado envuelven swell de espectro de ancho de banda angosta. This section
outlines differences that may be expected between the application of significant wave analyses and
application of an irregular wave approach.
La dinámica espectral surge porque olas de diferentes longitudes y esbeltéz se propagan a través y
junto con otras olas. Estas olas pueden intercambiar energía (transferencias alineales) y su
superposición puede conducir a la disipación debido a que rompen. El análisis de miles de registros
de oleaje (Bouws et al. 1987; Bouws, Gunther, and Vincent 1985; Miller and Vincent 1990) indica
que un espectro de aguas profundas adquiere una forma característica que es diferente a la de un
espectro de aguas someras. Como resultado, el nivel de energía de someramiento de una ola
irregular tiende a ser menor que el predicho para una ola monocromática, especialmente cerca de la
zona de surf. Smith y Vincent (1992) también indican que el someramiento y rompiente de olas
irregulares con dos picos espectrales puede diferir sustancialmente del caso en que las olas son
monocromáticas. Más aun, el espectro obtenido después de la refracción y someramiento puede
tener un período pico espectral diferente que el inicial. Aunque aun no se cuenta con una
explicación saisfactoria de este fenómeno, su impacto es importante, con diferencias de hasta 30-
40% en comparación con la aproximación de análisis significante.
Como ya se indicó anteriormente, conforme las olas se propagan, continuan creciendo debido a la
acción del viento, o pueden perder energía. Estos efectos no pueden ser incorporados
realisticamente en cálculos manuales. Los cálculos que envuelven rayos son tediosos y sujetos a
muchas inexactitudes. Se han hecho muchos avances en el cálculo de la transformación del oleaje;
estos ya se indicaron en las secciones precedentes. Aqui describimos tres de los códigos
computacionales disponibles y en uso por el CERC. Cada programa es complicado y requiere algo
de esfuerzo para su uso apropiado. Los tres modelos son modelos de estado estacionario. Existen
modelos de ecuaciones de aguas someras que son dependientes del tiempo (Jensen et al. 1987;
Demirbilek and Webster 1992a, 1992b), pero no se discuten aqui porque su uso requiere juegos
extensos de datos meteorológicos y no pueden ser aplicados tan fácilmente. A continuación, se
presentan las características básicas de los tres modelos mencionados:
RCPWAVE
RCPWAVE es un modelo de olas estacionario y lineal basado en la ecuación de pendiente suave e

incluye olas rompientes. Es aplicable a áreas costeras abiertas sin estructuras. Básicamente es un
modelo para olas monocromáticas.
REFDIF1
REFDIF es un modelo de estado estacionario que incluye olas rompeintes, atenuación y algunos
efectos alineales. Aunque es un modelo primariamente para olas monocromáticas, existe una

versión espectral disponible. El modelo puede simular aspectos de propagación asociados con
corrientes sencillas y puede incluir estruturas.
STWAVE
STWAVE es un modelo lineal estacionario que calcula la evolución del espectro direccional en el
espacio (Equation II-3-1). El modelo incluye rompiente, fricción con el fondo, percolación, y
forzamiento por viento, y resuelve para las transferencias alineales de energía dentro del espectro
de las olas. Tiene dos modos para tratar la difracción de la energía y el dominio computacional
puede incluir estructuras simples. El modelo puede manejar aspectos de propagación asociados con
corrinentes simples.

Mecánica de Olas

Cargado por

Copyright:

Formatos disponibles

Mecánica de Olas

Cargado por

Información del documento

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

Mecánica de Olas

Cargado por

Copyright:

Formatos disponibles

TRANSPORTE LITORAL COSTERO

CAPITULO 3. MECANICA DE OLAS

La descripción detallada de las ecuaciones básicas de la mecánica de olas puede consultarse en

carlos r torres iio/uabc 4/18/2008 1

Teorías más Importantes

i) Teoría Lineal de Airy: olas de amplitud pequeña en aguas profundas

carlos r torres iio/uabc 4/18/2008 2

3.2 Oleaje Regular

3.2.2 Definición de parámetros de la ola

Otros parámetros incluyen la frecuencia angular (ω = 2π/Τ), el número de onda κ = 2π/L, la

carlos r torres iio/uabc 4/18/2008 3

3.2.3 Teoría Lineal o de Airy

u = ∂ϕ/∂x; w = ∂ϕ/∂z (2.1)

Si la densidad del agua de mar, ρ, es constante, entonces;

∂u/∂x + ∂w/∂z = 0 (2.2)

o, en términos de ϕ(x, z, t), queda como;

carlos r torres iio/uabc 4/18/2008 4

p/ρ = ∂ϕ/∂t – (u2 + w2)/2 − gz = 0 (2.4)

η( x, t) = Η/2 cos(κx − ωt) (2.5)

Un resultado fundamental obtenido en la solución, lo constituye la relación de dispersión;

Sustituyendo los equivalentes de κ y ω, queda;

L = (gΤ2/2π) tanh (2πh/L ) = gΤ/ω tanh(κh) (II-1-10)

L ≈ (gΤ2/2π)(tanh (4π2 h /gΤ2))1/2

y su velocidad de fase, será ;

carlos r torres iio/uabc 4/18/2008 5

mientras que velocidad de fase, queda como;

Clasificación h/L κh tanh(κh)

Algunas funciones útiles

Haciendo la operación C/C0 nos queda;

C/C0 = L/L0 = tanh(κh)

Si multiplicamos ambos lados de la ecuación por h/L, queda;

h/L0 = h/L tanh(κh)

carlos r torres iio/uabc 4/18/2008 6

Encuentre las celeridades C y longitudes de la ola L correspondientes a las profundidades

Usando la ecuación (…)

L0 = (gΤ2/2π) = 1.56Τ2 m (5.12 Τ2 ft)

h/L0 = 200/156 = 1.2821

Notese que para valores de

L = L 0 = 156 m (aguas profundas ya que h/L > 1/2)

C = L/Τ = 156/Τ =156/10=15.6 m/s

h/L0 = 3/156 = 0.0192

de tablas se encuentra que

carlos r torres iio/uabc 4/18/2008 7

e) = /Τ = 53.2/10=5.32 m/s (17.4 ft/s)

Un valor aproximado para Λ puede obtenerse utilizando la ecuación

L ≈ (gΤ2/2π) (tanh(4π2h/Τ2 ))1/2

que en términos de L0 puede escribirse como;

L ≈ 156 (0.1202)1/2 = 54.1 (177.5 ).

Velocidades y aceleraciones de las partículas

Las componentes de la velocidad en la vertical, ω, y en la horizontal, υ, están dadas por;

= Τ/ ( ( + )/( )) ( −ω) (22)

= Τ/ ( ( + )/( )) ( −ω) (23)

carlos r torres iio/uabc 4/18/2008 8

αx= (gπH/L) [cosh[2π(z + d)/L]/cosh(2πd/L) sin θ = ∂u/∂t

αz= -(gπH/L) [sinh[2π(z + d)/L]/cosh(2πd/L) cos θ = ∂w/∂t

carlos r torres iio/uabc 4/18/2008 9

carlos r torres iio/uabc 4/18/2008 10

La evaluación de las constantes en las ecuaciones, da

( ΤΗ/2 ) 1/( ) =(5.5 × 9.8 × 8)/(2×81.7)(1/1.742)=1.515

(Η π/ ) 1/( ) =(5.5 × 9.8 × 3.1416)/(81.7)(1/1.742)=1.190