Teorı́a de ondas cortas de pequeña

amplitud en superficie libre

Notas del curso Hidráulica Marı́tima

Ismael Piedra-Cueva Agosto 2019

IMFIA - Facultad de Ingenierı́a Universidad de la República
Capı́tulo 1

Introducción

En este capitulo se presentan las ecuaciones de ondas de pequeña amplitud en

flujos a superficie libre cuando se propagan en un medio en el cual la profundidad
es constante en el espacio.
Las ecuaciones de onda en flujo irrotacional, están gobernadas por la ecuación
de Laplace con las correspondientes condiciones de contorno.
Un problema principal para la obtención de una solución, es que las condiciones
de borde que se aplican en la superficie libre, se aplican en un lugar del espacio
fı́sico desconocido a priori, por lo cual el problema no tiene solución exacta desde
un punto de vista formal.
Adicionalmente, las condiciones de borde en la superficie libre son no lineales, lo
cual es una complejidad adicional.
La teorı́a lineal fue desarrollada por Airy en 1845, y Stokes en 1947 desarrollo so-
luciones de 1er orden y ordenes superiores. Esta aproximación se basa en asumir que
la amplitud de la onda es pequeña en relación a la longitud de onda  = a/L << 1,
de forma tal que la solución no se aparta demasiado de la posición de reposo, de
donde proviene el nombre de ecuaciones de onda de pequeña amplitud. Esta sim-
plificación permite aplicar las condiciones de borde en z = 0 o en la superficie del
reposo, la cual es conocida, y en consecuencia es posible una solución. Adicional-
mente se puede demostrar que cumpliendo esta hipótesis se pueden despreciar los
términos no lineales. Por otra parte, aunque menos explicitado, se debe verificar que
la amplitud sea también pequeña en relación a la profundidad a/h << 1.
Si bien estas simplificaciones limitan el rango de validez de la solución, en la
práctica se ha verificado que su solución puede ser aplicada con muy buena aproxi-
mación aun no verificando estrictamente las restricciones anteriores.
Adicionalmente, esta teorı́a es la base para comprender las soluciones de ordenes
superiores y otras teorı́as de propagación de ondas que son frecuentemente utilizadas.

2
Draft Notas de curso. Teoria de ondas cortas de pequeña amplitud en superficie libre

Muchos estudios de esta ecuación se han realizado desde entonces.

Capı́tulo 1 IPC 3
Capı́tulo 2

Teorı́a de olas en fondo plano

En este capitulo se presenta el problema 2-D y ecuaciones a resolver junto a las

condiciones de borde correspondientes.
La hipótesis de base de las teorı́a de olas en su conjunto, es que el flujo es
irrotacional, y en consecuencia se desprecian los efectos viscosos y se asume que
el flujo puede ser representado por un potencial de velocidades, esto es el flujo es
potencial.
Sin embargo, no es inmediato aceptar esta hipótesis, ya que en los flujos a super-
ficie libre en general son turbulentos, con un desarrollo de capa limite bien definido.
En estos flujos, se entiende que la capa limite alcanza la superficie libre, o sea abar-
ca toda la profundidad del flujo. Por otra parte, esta capa limite permite distinguir
distintas zonas, que no es necesario detallar ahora, pero en las proximidades de la
pared predominan los efectos viscosos. Dependiendo de la rugosidad de fondo, puede
existir o no la sub-capa viscosa.
El flujo inducido por el oleaje es un movimiento periódico de periodo T , al igual
que el inducido por las mareas. Pero el hecho peculiar de las olas, es que el el periodo
de la oscilación es sustancialmente más corto que el de las mareas (10-20 segundos
las olas, 12 horas las mareas).
El estudio de la capa limite laminar oscilante en el plan horizontal, estudiado
hace muchos años por Schlichting (1960), ver esto, permite extraer dos conclusiones
importantes. La 1ra, es que la altura de la capa limite de fondo δ depende, ademas
de la viscosidad, de la frecuencia de oscilación, y puede ser estimada como δ ≈
p p
5 2ν/σ = 5 T ν/π. Por lo tanto, cuando la frecuencia es muy baja, o sea el
periodo el grande, el espesor de la capa limite crece. Para periodo pequeños, δ es
pequeño. Por ejemplo para un flujo laminar y T=10 segundos, δ ≈ 8,9mm. Si la capa
limite fuera turbulenta, la viscosidad turbulenta serı́a mayor que la laminar, aunque
la presencia de la pared limita el desarrollo de la turbulencia y en consecuencia no

4
Draft Notas de curso. Teoria de ondas cortas de pequeña amplitud en superficie libre

es demasiado elevada, por lo cual δ permanece muy pequeña, en general menor a

10cm.
El segundo elemento, se refiere a la vorticidad generada en la capa limite.
En la capa lı́mite se generan vórtices. Durante la oscilación, en los momentos
de inversión del flujo, se produce el desprendimiento de vórtices que son ejectados
hacia la columna de agua. Sin embargo, esta ejección no es continua, solo ocurre
en los momentos próximos a la inversión del flujo, por lo cual una gran parte de
los vórtices se destruirán entre una ejección y la siguiente. Por lo tanto, solo una
parte muy pequeña de los vórtices generados se transportan y permanecen en el
flujo externo a la capa limite, y se puede concluir que la capa lı́mite casi no afecta
al resto de la columna de agua.
En consecuencia, fuera de la capa lı́mite, el flujo se puede asumir como irrotacio-
nal: las fuerzas viscosas son despreciables y las fuerzas externas son esencialmente
conservativas, ya que la fuerza de gravitación es la dominante (inicialmente se des-
precia la capilaridad y Coriolis). Esto significa que si consideramos olas más largas
que unos pocos centı́metros (de forma que la capilaridad no sea relevante) y más
cortas que algunos kilómetros (de forma que Coriolis no sea relevante), podemos
asumir que las fuerzas externas son conservativas.
La conclusión es que el flujo ambiental (las ondas) con buena precisión pueden
describirse como un flujo potencial.

2.1. Ecuaciones de la onda

Se considera que las olas se propagan en el plano x-z en el cual la profundidad h

es constante en tiempo y espacio.
En la Figura (2.1) se presenta un esquema de la onda y principales variables.
La profundidad del cuerpo de agua es h. El origen de coordenadas se ubica en
la superficie libre en reposo (z = 0), y el fondo se ubica en z = −h.
La onda en superficie se describe por una función η(x, t).
Las velocidades horizontal y vertical en el cuerpo de agua son (u, w).
Se buscaran soluciones periódicas de perı́odo T , cuya frecuencia es σ = 2π/T .
La longitud de onda es L y el numero de onda k = 2π/L.
La semi-amplitud de la onda es a, y la altura de la ola H = 2a.
La celeridad de la onda o velocidad de la fase, es C = L/T = σ/k.

Capı́tulo 2 IPC 5
Notas de curso. Teoria de ondas cortas de pequeña amplitud en superficie libre Draft

Figura 2.1: Esquema onda superficie libre.

2.1.1. Ecuaciones de movimiento

Las ecuaciones de partida son la ecuación de conservación de la masa y de can-
tidad de movimiento para un fluido no viscoso, esto es las ecuaciones de Euler.
La ecuación de conservación de la masa se expresa como:
∂u ∂w
+ =0 (2.1)
∂x ∂z
La ecuación de conservación de la cantidad de movimiento se expresa como:
∂u ∇p
+ (u.∇)u = − +g (2.2)
∂t ρ

donde u es el vector velocidad y g vector gravedad.

La fuerza de la gravedad es conservativa, por lo cual puede expresarse como
−∇(gz).
Por otra parte, la parte convectiva de la velocidad puede expresarse de la siguiente
forma:

1
(u.∇)u = ∇( |u|2 ) + w × u
2
donde w = ∇ x u es el rotor de la velocidad, y |u|2 = (u2 + w2 ).
Sustituyendo:

∂u (u2 + w2 ) ∇p
+ ∇( )+w×u=− − ∇(gz)
∂t 2 ρ
∂u (u2 + w2 ) ∇p
+ ∇( )+ + ∇(gz) = −w × u = u × w (2.3)
∂t 2 ρ
Si el flujo es irrotacional, w = ∇ x u = 0, por lo cual el término de la derecha
se anula, y se obtiene la ecuación de flujo irrotacional:

∂u (u2 + w2 ) ∇p
+ ∇( )+ + ∇(gz) = 0 (2.4)
∂t 2 ρ

6 IPC Capı́tulo 2
Draft Notas de curso. Teoria de ondas cortas de pequeña amplitud en superficie libre

Matemáticamente, la condición necesaria para que el rotor sea nulo es que la

velocidad derive de un potencial, o sea que es posible encontrar un potencial de
velocidades que verifique u = −∇φ, siendo φ el potencial de velocidades, función
escalar φ(x, z, t).
En consecuencia es sencillo verificar que el producto vectorial es nulo:
w = ∇ x u = ∇ x ∇φ = 0.
 
i j k
 
 
∂x ∂y ∂z  = 0
 
 
 φx φy φz 

i(φyz − φzy ) + j(φzx − φxz ) + k(φxy − φyx ) = 0

Lo expresado anteriormente, significa entonces que existe un potencial de velo-

cidades u = −∇φ(x, z, t).
La ecuación de cantidad de movimiento puede expresarse como:

∂φ (u2 + w2 ) p
∇ (− + + + gz) = 0
∂t 2 ρ
Una primera integración permite obtener la ecuación de Bernoulli:

∂φ (u2 + w2 ) p
− + + + gz = C(t) (2.5)
∂t 2 ρ

El valor de la constante C(t) depende de la ubicación del sistema coordenado.

Para un sistema coordenado como el que nosotros utilizamos, en el cual el origen
del sistema se ubica en la superficie libre, esto es en el reposo z = 0, para que se
verifique que en la superficie libre la presión p(z = 0) = 0 , la constante debe valer
C = 0.
Por el contrario, si el origen del sistema coordenado fuera el fondo, en el cual la
superficie libre en reposo se ubica en z = h, entonces para verificar la presión nula
p(z = h) = 0 , C debe valer C = gh.
La ecuación de Bernoulli queda expresada entonces como:
∂φ p 1 ∂φ ∂φ
− + gz + + [( )2 + ( )2 ] = 0 (2.6)
∂t ρ 2 ∂x ∂z
Sustituyendo ahora el potencial de velocidades se en la ecuación de conservación
de la masa obtiene la ecuación de Laplace:

∂ 2φ ∂ 2φ
∇2 φ(x, z, t) = + =0 (2.7)
∂x2 ∂z 2
la cual permite calcular el potencial de velocidades en todo el dominio de solución.

Capı́tulo 2 IPC 7
Notas de curso. Teoria de ondas cortas de pequeña amplitud en superficie libre Draft

2.1.2. Energı́a y Flujo de energı́a

La energı́a del movimiento es la suma de la energı́a cinética mas la energı́a
potencial.
La energı́a total en un volumen V esta dada por:
Z
1
E= ( ρu2 + ρgz)dV (2.8)
V 2
Donde V es un volumen arbitrario, que en este caso tomamos como fijo en el
tiempo.
El flujo de energı́a esta dado por definición por a variación de la energı́a en el
tiempo.
Z
dE d 1
F = = ( ρu2 + ρgz)dV
dt dt V 2
y como el volumen V es fijo:
Z
∂ 1 2
F = ( ρu + ρgz)dV
V ∂t 2
La energı́a potencial tiene derivada parcial nula, por lo tanto
Z
∂ 1 2
F = ( ρu )dV (2.9)
V ∂t 2

Ahora podemos utilizar la ecuación de movimiento (2.3). Multiplicando por la

velocidad u.
∂u u2 ∇p
u. + u.∇( ) + u. + u.∇(gz) = 0
∂t 2 ρ
donde se aplica la propiedad de que (∇ × u).u = 0 ya que (u es ortogonal a
w), y su producto escalar es nulo, y lo que se verifica automáticamente si el flujo es
irrotacional w = ∇ × u = 0 . Por lo tanto
∂ 1 2 ∂u u2
ρ ( u ) = ρu = −ρu.∇( ) − u.∇p − u.∇(ρgz)
∂t 2 ∂t 2
Teniendo en cuenta que

∇(u2 .u) = u.∇(u2 ) + u2 ∇.u = u.∇(u2 )

ya que ∇.u = 0 por conservación de la masa.

∂u (u2 .u)
ρu = −ρ∇( ) − ∇pu − ∇(ρgzu)
∂t 2
El flujo de energı́a se expresa entonces

(u2 .u)
Z
dE
F = =− (ρ∇( ) + ∇pu + ∇(ρgzu)) dV
dt V 2

8 IPC Capı́tulo 2
Draft Notas de curso. Teoria de ondas cortas de pequeña amplitud en superficie libre

Aplicando teorema de Gauss:

ρu2
Z
dE
F = =− ( + p + ρgz) u.n dS (2.10)
dt S 2

donde n es la normal exterior a la superficie S.

El flujo de energı́a por unidad de superficie se define entonces como:

ρu2
P =( + p + ρgz) u (2.11)
2
lo cual representa la suma del flujo de energı́a cinética mas la potencia de las
presiones mas el flujo de energı́a potencial.

2.2. Condiciones de borde

La ecuación de Laplace describe el movimiento de todos los flujos potenciales, por
lo cual el movimiento en cada caso particular esta determinado por las condiciones
de borde.
Hay dos tipos de condiciones de borde: cinemática y dinámica. En el fondo, se
aplica una condición cinemática. En la superficie libre se debe aplicar una condición
cinemática y otra condición dinámica.
Condición de borde de fondo z = −h. En el fondo, la condición cinematica es que
la velocidad del flujo normal al fondo sea nula. Si el fondo es inclinado la condición
serı́a que u.n = 0, donde n es el vector normal al fondo.
Para el caso de fondo horizontal, entonces w(z = −h) = 0, o ∂φ/∂z|−h = 0 .
En la superficie libre z = η(x, t) se debe verificar que la presión es nula, p(x, η, t) =
0, por tanto

∂φ 1 ∂φ ∂φ
− + gη + [( )2 + ( )2 ] = 0 (2.12)
∂t 2 ∂x ∂z
Adicionalmente, en la superficie libre se debe verificar la condición de borde
cinemática.
Se asume que la superficie libre F (x, z, t) = 0 puede expresar en forma semi-
implicita como una función F (x, z, t) = z − η(x, t) = 0.
Un observador que se desplaza adherido a la superficie libre no observa flujo a
través de la misma, por lo cual la derivada total o Lagragiana debe ser nula:

d
F (x, z, t) = 0
dt
d
(z − η(x, t)) = 0
dt

Capı́tulo 2 IPC 9
Notas de curso. Teoria de ondas cortas de pequeña amplitud en superficie libre Draft

dz d
= η(x, t)
dt dt
Aplicando la definición de velocidad,

dz d
= w(η) = η(x, t)
dt dt

∂η(x, t) ∂η(x, t)
w(η) = + u(η)
∂t ∂x
que es la condición cinemática de la superficie libre.

2.3. Resumen ecuaciones de la onda

En resumen, el problema a resolver es:

∂ 2φ ∂ 2φ
∇2 φ = + =0 − h <= z <= η(x, t) (2.13)
∂x2 ∂z 2

∂φ/∂z = 0 z = −h (2.14)

∂φ 1 ∂φ ∂φ
− + gη + [( )2 + ( )2 ] = 0 z = η(x, t) (2.15)
∂t 2 ∂x ∂z
∂η(x, t) ∂η(x, t)
w(η) = + u(η) z = η(x, t) (2.16)
∂t ∂x
Este sistema de ecuaciones no tiene solución por dos razones.
1.- Las 2 condiciones de borde que se aplican en la superficie libre, se deben
aplicar en una superficie η(x, t) que es desconocida su forma y posición.
2.- Ambas condiciones de borde son no lineales.
Por estas razones, es necesario efectuar aproximaciones de algún tipo.

10 IPC Capı́tulo 2
Capı́tulo 3

Teorı́a lineal del oleaje

En este capı́tulo se desarrollan las ecuaciones correspondientes a la teorı́a lineal

de ondas, y se obtiene su solución, lo cual permitirá estudiar posteriormente la
cinemática y distribución de presiones, entre otras cosas.
La derivación de la ecuación lineal de la ola comienza con la ecuación de Laplace.
El punto de partida es conocer una solución exacta del potencial de velocida-
des, que en este caso es el reposo, y esta descrito por la solución φo (x, z, t) = 0 y
ηo (x, z, t) = 0.
La solución general es no lineal, por lo cual para resolverlas, utilizamos una
expansión en serie en un parámetro pequeño  llamado de Poincare, parámetro no
definido a priori, pero  << 1.

3.1. Teorı́a lineal del oleaje 1er orden

Buscamos soluciones próximas a la solución exacta que se conoce, o sea soluciones
de la forma

φ = φo + φ1 (x, z, t) + 2 φ2 (x, z, t) + 3 φ3 (x, z, t) + ... (3.1)

η = ηo + η1 (x, z, t) + 2 η2 (x, z, t) + 3 η3 (x, z, t) + ... (3.2)

O sea un desarrollo en serie de potencia entera en .

Ahora sustituimos este desarrollo en serie en la ecuación de Laplace. El Lapla-
ciano o las derivadas segundas son operadores lineales, por lo cual se puede expresar
como:

∇2 φ = ∇2 φo + ∇2 φ1 (x, z, t) + 2 ∇2 φ2 (x, z, t) + 3 ∇2 φ3 (x, z, t) + ... = 0 (3.3)

donde sabemos que ∇2 φo (x, z, t) = 0 por ser el reposo.

11
Notas de curso. Teoria de ondas cortas de pequeña amplitud en superficie libre Draft

Al orden 1 se tiene ∇2 φ1 (x, z, t) = 0, y al orden 2 ∇2 φ2 (x, z, t) = 0, etc.

Ahora sustituimos en las condiciones de borde.
En el fondo z = −h

∂φ ∂φo ∂φ1 ∂φ2 ∂φ3

= + + 2 + 3 + ... = 0 (3.4)
∂z ∂z ∂z ∂z ∂z
Al 1er orden
∂φ1
=0 z = −h (3.5)
∂z
Al 2do orden
∂φ2
= 0 z = −h (3.6)
∂z
Ahora aplicamos las condiciones de borde cinemática en la superficie libre.

∂φ ∂η(x, t) ∂φ ∂η(x, t)
− = − z = η(x, t) (3.7)
∂z ∂t ∂x ∂x

∂φo ∂φ1 2 ∂φ2 ∂ηo ∂η1 2 ∂η2 ∂φo ∂φ1 2 ∂φ2 ∂ηo ∂η1 2 ∂η2
− − − = + + −[ + + ][ + + ]
∂z ∂z ∂z ∂t ∂t ∂t ∂x ∂x ∂x ∂x ∂x ∂x
Al 1er orden
∂φ1 ∂η1
− = z = η(x, t) (3.8)
∂z ∂t
Al 2do orden

∂φ2 ∂η2 ∂φ1 ∂η1

−= − z = η(x, t) (3.9)
∂z ∂t ∂x ∂x
Ahora aplicamos la condición dinámica:

∂φo ∂φ1 ∂φ2

− − − 2 + g[ηo + η1 + 2 η2 + 3 η3 ]
∂t ∂t ∂t
1 ∂φo ∂φ1 ∂φ2 2 ∂φo ∂φ1 ∂φ2
+ [[ + + 2 ] +[ + + 2 ]] = 0
2 ∂x ∂x ∂x ∂z ∂z ∂z
Al 1er orden

∂φ1
− + gη1 = 0 z = η(x, t) (3.10)
∂t
Al 2do orden

∂φ2 1 ∂φ1 2 ∂φ1 2

− + gη2 + [( ) +( ) ] = 0 z = η(x, t) (3.11)
∂t 2 ∂x ∂z
Las condiciones de borde en la superficie libre deben aplicarse en z = η que es
desconocida.
Para resolver este aspecto, la aproximación a seguir es efectuar un desarrollo en
serie de Taylor tomando como referencia la coordenada z = 0.

12 IPC Capı́tulo 3
Draft Notas de curso. Teoria de ondas cortas de pequeña amplitud en superficie libre

∂φ1 ∂φ1 ∂ 2 φ1
|z=η = |z=0 + η 2 |z=0 + ....
∂z ∂z ∂z

∂η1 ∂η1 ∂ 2 η1
|z=η = |z=0 + η |z=0 + ....
∂t ∂t ∂t∂z

∂φ1 ∂φ1 ∂ 2 φ1
|z=η = |z=0 + η |z=0 + ....
∂t ∂t ∂t∂z
Reteniendo solo los términos de 1er orden y sustituyendo en la condiciones de
borde se obtiene finalmente al 1er orden:

∂φ1 ∂η1
− = z=0 (3.12)
∂z ∂t
∂φ1
− + gη1 = 0 z=0 (3.13)
∂t
Combinando ambas condiciones de borde en superficie libre se tiene:

∂ 2 φ1 ∂φ1
2
+g z=0 (3.14)
∂t ∂z
Finalmente, se debe verificar que la onda sea periódica en el espacio y en el
tiempo, para lo cual se debe verificar que:
φ(x, z, t) = φ(x + L, z, t)
φ(x, z, t) = φ(x, z, t + T )

3.2. Resumen ecuaciones lineales 1er orden

En resumen las ecuaciones lineales al 1er orden son las siguientes. La Figura (3.2)
muestra un esquema.

∂ 2 φ1 ∂ 2 φ1
+ =0 − h <= z <= 0 (3.15)
∂x2 ∂z 2
∂φ1
=0 z = −h (3.16)
∂z

∂ 2 φ1 ∂φ1
+ g z=0 (3.17)
∂t2 ∂z

φ(x, z, t) = φ(x + L, z, t) (3.18)

φ(x, z, t) = φ(x, z, t + T ) (3.19)

Capı́tulo 3 IPC 13
Notas de curso. Teoria de ondas cortas de pequeña amplitud en superficie libre Draft

Figura 3.1: Esquema onda lineal superficie libre.

3.3. Solución ecuaciones lineales 1er orden

La ecuación diferencial del Laplaciano se resuelve en forma clásica utilizando el
método de separación de variables, esto es se admite que la solución del potencial
esta compuesto por dos funciones independientes entre si:

φ(x, z, t) = N (t)Z(z)X(x) (3.20)

Sustituyendo en el laplaciano se tiene:

X 00 (x) Z 00 (z)
=− = C1 (3.21)
X(x) Z(z)
La condición de borde de la superficie libre (3.17)

N 00 (t) Z 0 (z)
= −g = C2 z=0 (3.22)
N (t) Z(z)

En este caso, la constante C2 no puede ser positiva, porque las soluciones de

N (t) son exponenciales en el tiempo, y no son soluciones periódicas. Por lo cual solo
consideramos soluciones C2 < 0, C2 = −σ 2 , donde σ es real, y desconocida a priori.
Las soluciones de N(t) serán de la forma N (t) = (sin σt, cos σt) .
En cuanto a C1 hay dos posibilidades: que sea positiva o negativa.

Caso C1 < 0

Caso 1: C1 < 0, C1 = −k 2 , k real, y también desconocido a priori.

La ecuación (3.21) para la componente Z(z)

14 IPC Capı́tulo 3
Draft Notas de curso. Teoria de ondas cortas de pequeña amplitud en superficie libre

Z 00 (z) = k 2 Z(z)

cuya solución es de la forma

Z(z) = A cosh k(h + z) + B sinh k(h + z)

Aplicando la condición de borde de fondo Z 0 (z = −h) = 0

Z 0 (z = −h) = Ak sinh k(h − h) + Bk cosh k(h − h) = 0

por lo cual se verifica solo si B = 0.

Z(z) = A cosh k(h + z)

Para la componente en x se tiene:

X 00 (x) = −k 2 X(x) (3.23)

Las soluciones de X(x) serán de la forma X(x) = (sin kx, cos kx) .
En resumen, para el caso C1 < 0, C1 = −k 2 , la solución de la ecuación de
Laplace es una combinación lineal de:

φ(x, z, t) = (sin kx, cos kx) cosh k(h + z)(sin σt, cos σt)

o en forma más explicita el potencial puede tener cualquiera de las siguientes

expresiones:
φ(x, z, t) = A cosh k(h + z) cos kx cos σt (3.24)

φ(x, z, t) = A cosh k(h + z) sin kx sin σt (3.25)

φ(x, z, t) = A cosh k(h + z) sin kx cos σt (3.26)

φ(x, z, t) = A cosh k(h + z) cos kx sin σt (3.27)

Por otra parte, se debe verificar la condición de borde en la superficie libre

Z 0 (z)
−g = −σ 2 z=0
Z(z)

g Z 0 (z) = σ 2 Z(z) z=0

gkA sinh kh = σ 2 A sinh kh

σ 2 = gk tanh kh (3.28)

denominada ecuación de dispersión de la ola.

Capı́tulo 3 IPC 15
Notas de curso. Teoria de ondas cortas de pequeña amplitud en superficie libre Draft

3.3.1. Ecuación de dispersión

La ecuación (3.28) denominada ecuación de dispersión, es una de las ecuaciones
mas trascendentes de la teorı́a del oleaje, lineal o no lineal, ya que relaciona la
profundidad h, el perı́odo T o frecuencia de la ola σ, y la longitud de onda L o
número de onda k.
Si bien la ecuación permite calcular en forma directa la frecuencia σ a partir de
los datos de profundidad h y número de onda k, no es lo usual porque en general el
problema es el inverso, conocido el perı́odo T de la onda, se quiere calcular k.
La solución de k debe resolverse en forma iterativa. La ecuación de dispersión
tiene solución única. En efecto, puede expresarse como

σ2h
= tanh kh
g kh

que son dos curvas que se interceptan en una sola solución para kh.
Pero además, permite calcular la celeridad o velocidad de fase de la onda, como
C = σ/k.
p p
C= g/k tanh kh = gL/2π tanh kh (3.29)

Se observa que la celeridad de las ondas, dependen de su perı́odo, por lo cual

ondas de mayor periodo tiene mayor celeridad que las de menor periodo.
En consecuencia, un tren de ondas que originalmente tiene componentes de dife-
rente periodo, al propagarse se dispersan, ya que las ondas de mayor periodo viajan
mas rápido que la de menor perı́odo. De aquı́ su nombre ecuación de dispersión.
Este resultado es sustancialmente diferente al que se obtiene en la teorı́a de
ondas en aguas poco profundas o de la celeridad de las perturbaciones en las cuales
√
C = gh.
En esta expresión, la celeridad de la onda solo depende de la profundidad, inde-
pendientemente del periodo. Las ondas de marea, de diferente perı́odo, todas viajas
√
a la misma velocidad, dada por gh. Son ondas no-dispersivas.
La celeridad de onda dada por (3.29) es exacta, y es la mejor solución conocida
al presente para calcular la velocidad de propagación de una onda.
La bondad de esta solución, y a pesar de las diferentes aproximaciones realizadas
en la solución, surge del hecho de que resuelve el problema de la forma mas completa
posible en el plano vertical x-z.
Las diferentes teorı́as de olas de la familia del tipo Serre, Boussinesq, KdV, sufren
el problema que no reproducen correctamente la celeridad de propagación de la onda
dada por (3.29). Numerosas alternativas de solución se han buscado para mitigar
este aspecto.

16 IPC Capı́tulo 3
Draft Notas de curso. Teoria de ondas cortas de pequeña amplitud en superficie libre

En la teorı́a de oleaje es de interés analizar dos casos limites de la teorı́a. Uno

se refiere al caso de aguas profundas, en la cual la profundidad h >> L, y en
consecuencia kh → ∞. El otro extremo es el caso de aguas poco profundas, en el
cual h << L, y en consecuencia kh → 0.

Caso aguas profundas kh → ∞

En este caso la tanh kh → 1 y la celeridad en aguas profundas, denominada

p
Co = g/k. O sea la celeridad de la onda solo depende del número de onda k.
Por otra parte, como C = L/T , podemos calcular la longitud de onda corres-
pondiente a aguas profundas como (L/T )2 = g/k = gL/2π, de donde

g T2
Lo = ≈ 1,56 T 2 (3.30)
2π
La longitud de onda en aguas profundas Lo varia con el cuadrado del perı́odo.
Por ejemplo una ola de T=10 s tendrá una longitud de onda de 156m. Su número
de onda ko = 2π/Lo = 0,04rad/m.
Se puede inferir entonces que la longitud de onda da la dimensión espacial hori-
zontal y en vertical del espacio afectado por la onda. Para profundidades mayores a
156m, estas olas no están afectadas por la profundidad.
En general esta condición se cumple cuando h/L > 1/2.
gT
p
La celeridad en aguas profundas esta dada entonces Co = g/ko = Lo /T = 2π
.
Ondas de mayor perı́odo tiene mayor longitud de onda y mayor celeridad, por lo
tanto viajan más rápido.
Por otra parte, la longitud de onda se puede expresar en función de Lo

L2 gL
C2 = 2
= tanh kh
T 2π
Dividiendo por esta misma ecuación en aguas profundas:

gT 2
L= tanh kh = Lo tanh kh
2π
La longitud de onda tiene su máximo valor en aguas profundas y decrece con la
tanh kh.

Caso aguas poco profundas kh → 0

En este caso, la tanh kh = sinh kh/ cosh kh → kh y la celeridad en aguas poco

√
profundas C = gh. O sea la celeridad de la onda solo depende de la profundidad
h, siendo independiente del perı́odo. Estas ondas son no-dispersivas, todas viajan a
la misma velocidad.

Capı́tulo 3 IPC 17
Notas de curso. Teoria de ondas cortas de pequeña amplitud en superficie libre Draft

√
Adicionalmente se se observa la longitud de onda L = T gh.
En consecuencia, L → 0 cuando h → 0.
Este resultado, esperable, introduce una limitación muy fuerte al rango de validez
de la teorı́a lineal de pequeña amplitud, basada en la hipótesis de que  es un
parámetro pequeño  = a/L << 1.
Cuando L → 0,  → ∞ y el desarrollo en serie sobre el cual se fundamenta la
teorı́a deja de ser válido.
Por esta razón, se han desarrollado otras teorı́as basadas en otras aproximaciones
que incorporan un segundo parámetro a/h, que describen mejor el comportamiento
de las ondas en aguas poco profundas. Las ecuaciones de Serre, Boussinesq y KdV
son ejemplos de estas aproximaciones.

3.4. Ondas progresivas

En la sección anterior, se calculó el potencial de velocidad que verifica la ecuación
de Laplace y las condiciones de borde.
La ecuación de Laplace es una ecuación lineal, por lo cual es posible superponer
las diferentes soluciones para encontrar distintos casos de interes.
El caso de la onda progresiva se puede obtener combinando dos soluciones del
potencial de velocidades:

φ(x, z, t) = A cosh k(h + z) sin kx cos σt − A cosh k(h + z) cos kx sin σt

φ(x, z, t) = A cosh k(h + z) sin(kx − σt)

Esta ecuación debe verificar la condición

∂φ
− + gη = 0 z=0
∂t
1
η = (−σ)A cosh kh cos(kx − σt) = a cos(kx − σt)
g
Por lo cual
−ga
A=
σ cosh kh
El potencial de velocidades de una onda progresiva es

ga cosh k(h + z)
φ(x, z, t) = − sin(kx − σt) (3.31)
σ cosh kh
Finalmente, se debe verificar que la onda sea periódica en el espacio y en el
tiempo, para lo cual se debe verificar que:

18 IPC Capı́tulo 3
Draft Notas de curso. Teoria de ondas cortas de pequeña amplitud en superficie libre

φ(x, z, t) = φ(x + L, z, t)
Para que esto se cumpla, se debe verificar que
sin kx = sin k(x + L) sin kx = sin(kx) cos(kL) + sin(kL) cos(kx)
lo cual se cumple solo si sin(kL) = 0 y cos(kL) = 1.
Por lo tanto, kL=0,2π. Por lo cual k = 2π/L. Es el numero de onda.
De la misma forma, para que φ(x, z, t) = φ(x, z, t + T ),
se debe verificar que
sin σt = sin σ(t + T )
sin σt = sin(σt) cos(σT ) + sin(σT ) cos(σt)
lo cual se cumple solo si sen(σT ) = 0 y cos(σT ) = 1.
Por lo tanto, σT = 0, 2π.
Por lo cual σ = 2π/T , es la frecuencia angular.
La ecuación de la superficie libre esta dada por

η(x, t) = a cos(kx − σt) (3.32)

Al argumento de la función cos se lo denomina fase de la onda S(x, t) = (kx−σt).

Para un observador que viaja con la onda, la derivada total de S(x, t) es nula,
por lo que cual dS(x, t)/dt = kdx/dt − σdt/dt = 0.
dx/dt = σ/k = C. La celeridad de la onda es positiva, la onda progresiva viaja
hacia el x>0.
Una función de fase S(x, t) = (kx + σt), da una celeridad C = −σ/k que repre-
senta una onda progresiva que viaja en sentido de los x negativos.
Por otra parte, también se observa que ∇S(x, t) = k , o sea que el vector
número de onda es el gradiente de la fase; y su derivada temporal es la frecuen-
cia ∂S(x, t)/∂t = −σ.
Derivando respecto del tiempo
∂(∇S)/∂t = ∂k/∂t.
∇∂S/∂t = ∂k/∂t
−∇σ = ∂k/∂t

∂k/∂t + ∇σ = 0 (3.33)

Esta ecuación expresa que una variación espacial de la frecuencia de la onda debe
tener acompañada la variación del número de onda o de la longitud de onda en un
punto fijo en el tiempo.
Como la onda tiene un periodo fijo (y una frecuencia), en un punto x fijo su
longitud de onda dada por la ecuación de dispersión no varı́a en el tiempo. Por lo
tanto, cuando la onda se propaga en el espacio, su frecuencia (y en consecuencia su

Capı́tulo 3 IPC 19
Notas de curso. Teoria de ondas cortas de pequeña amplitud en superficie libre Draft

periodo) tampoco puede variar durante la propagación, ∇σ = 0. En otras palabras,

una onda con un periodo dado en aguas profundas, lo mantiene constante durante
su propagación hacia la costa.

3.5. Velocidad de las partı́culas de agua

El campo de velocidades y presiones se puede obtener directamente del potencias
de velocidades calculado anteriormente (3.31).
Las componentes de la velocidad de las partı́culas fluidas, llamada velocidad
orbital, se calculan en forma directa.

∂φ gka cosh k(h + z)

u(x, z, t) = − = cos(kx − σt)
∂x σ cosh kh
cosh k(h + z)
u(x, z, t) = σa cos(kx − σt) (3.34)
sinh kh
σ2
donde se utilizo la ecuación de dispersión gk = tanh kh

∂φ gka sinh k(h + z)

w(x, z, t) = − = sin(kx − σt)
∂z σ cosh kh
sinh k(h + z)
w(x, z, t) = σa sin(kx − σt) (3.35)
sinh kh
La superficie libre esta dada por

η(x, z, t) = a cos(kx − σt) (3.36)

Por lo cual se observa que la velocidad horizontal esta en fase con la superficie
libre. Las máximas velocidades horizontales se producen en el momento de la cresta
de la ola, y las mı́nimas en el valle, donde es negativa. En el momento que pasa la
cresta de la ola, la velocidad u se orienta hacia el sentido de propagación de la onda,
mientras que cuando pasa el valle de la onda la velocidad apunta hacia donde viene
la ola.
Por otra parte, se observa que la máxima velocidad ocurre en la superficie libre,
esto es z = 0, y que decrece con el cosh k(h + z), siendo z negativo hacia abajo,
alcanzando el menor valor en z = −h.
Por otra parte, la componte vertical w esta π/2 fuera de fase. Cuando u toma
un máximo/mı́nimo w es nula, y cuando u es nula, w es máxima/mı́nima.
Las Figuras 3.5 y 3.5 muestran un esquema de la distribución de velocidades, y
la Figura 3.5 muestra un esquema del movimiento orbital de las partı́culas durante
el pasaje de la onda.

20 IPC Capı́tulo 3
Draft Notas de curso. Teoria de ondas cortas de pequeña amplitud en superficie libre

Figura 3.2: Esquema distribución de velocidades.

Figura 3.3: Esquema distribución de velocidades. Tomado de Andersen and Frigaard

2011

La distribución de presiones esta dada por la ecuación de Bernoulli (2.6 ) linea-

lizada

Capı́tulo 3 IPC 21
Notas de curso. Teoria de ondas cortas de pequeña amplitud en superficie libre Draft

Figura 3.4: Movimiento orbital de las partı́culas de una ola que viaja hacia la derecha.
Tomado de Holthuijsen 2007.

p ∂φ
= − gz (3.37)
ρ ∂t
Reeplazando el potencial de velocidades

ga cosh k(h + z)
p(x, z, t) = σ cos(kx − σt) − ρgz
σ cosh kh

cosh k(h + z)
p(x, z, t) = ρga cos(kx − σt) − ρgz
cosh kh
La presión en tiene una componente hidrostática −ρgz, recordando que z es
negativo hacia abajo, y una componente dinámica pd . La componente dinámica esta
en fase con η ambas oscilan como cos(kx − σt). Además la presión dinámica decrece
con la profundidad.
Se debe notar, que el dominio de validez del potencial de velocidades es −h <
z < 0, por lo cual la zona comprendida entre la cresta y z=0 y entre z=0 y el valle,
están fuera del dominio de solución.
En z=0, la presión dinámica pd vale pd (x, 0, t) = ρga cos(kx − σt), lo cual
representa la presión hidrostática del fluido que se encuentra por encima de z=0.
Sin embargo, como varı́a esta presión en z=0 y la presión nula p = 0 en z = η no esta
definida. Una aproximación hidrostática es en algunos casos utilizada. La Figura 3.5
muestra un esquema de la distribución de presiones.

22 IPC Capı́tulo 3
Draft Notas de curso. Teoria de ondas cortas de pequeña amplitud en superficie libre

Figura 3.5: Distribución de presiones.

3.6. Trayectorias de las partı́culas de agua

Las componentes de velocidad permiten calcular las trayectorias de la partı́culas,
o sus órbitas.
Por definición, la trayectoria es la integral del campo de velocidades siguiendo
el movimiento de la partı́cula fluida. Por definición de velocidad u = dx/dt, la
trayectoria se calcula como

Z Z Z Z
dx(t) dz(t)
x(t) = dt = u(x(t), z(t), t)dt z(t) = dt = w(x(t), z(t), t)dt
t dt t t dt t

La integración no puede efectuarse porque las expresiones del campo de velocida-

des no permiten despejar las variables x(t), z(t), por lo cual se debe efectuar alguna
aproximación.
La velocidad que va variando en función del tiempo a lo largo de la trayectoria, la
aproximamos como la velocidad de un punto fijo (x,y), no dependiente del tiempo,
más una variación. En definitiva hacemos un desarrollo en serie de Taylor en 2
variables:

∂u ∂u
u(x(t), z(t), t) = u(x, z, t) + ∆x |x,z + ∆z |x,z + ....
∂x ∂z

Capı́tulo 3 IPC 23
Notas de curso. Teoria de ondas cortas de pequeña amplitud en superficie libre Draft

∂w ∂w
w(x(t), z(t), t) = w(x, z, t) + ∆x|x,z + ∆z |x,z + ....
∂x ∂z
Denominando (ξ, ζ) las coordenada (x,z) de la trayectoria se tiene entonces
Z
∂u ∂u
ξ(t) = [u(x, z, t) + ∆x |x,z + ∆z |x,z ] dt
t ∂x ∂z
Z
∂w ∂w
ζ(t) = [w(x, z, t) + ∆x |x,z + ∆z |x,z ] dt
t ∂x ∂z
Al 1er orden de aproximación, se puede despreciar el 2do término de la integral,
y entonces aproximar la ecuación de la trayectoria por
Z Z
cosh k(h + z)
ξ(t) = u(x, z, t) dt = [σa cos(kx − σt)] dt
t t sinh kh
Z Z
sinh k(h + z)
ζ(t) = w(x, z, t) dt = [σa sin(kx − σt)] dt
t t sinh kh
donde (x,z) son un punto fijo del espacio, por lo cual el integrando solo depende
del tiempo. Efectuando la integral:

cosh k(h + z)
ξ(t) = −a sin(kx − σt) = −A(x, z) sin(kx − σt)
sinh kh
sin k(h + z)
ζ(t) = a cos(kx − σt) = B(x, z) cos(kx − σt)
sinh kh
donde
cosh k(h + z) sinh k(h + z)
A(x, z) = a B(x, z) = a (3.38)
sinh kh sinh kh
Elevando al cuadrado y sumando ambas componentes:

ξ(t)2 /A(x, z)2 ) = sin2 (kx − σt),

ζ(t)2 /B 2 (x, z)) = cos2 (kx − σt)

ξ(t)2 ζ(t)2
+ =0 (3.39)
A(x, z)2 B 2 (x, z)
Es la ecuación de la elipse, por lo cual las trayectoria al 1er orden de aproximación
son elipses, con orbitas cerradas, en la cual A el el radio mayor y B el radio menor.
La Figura 3.6 presenta un esquema de las órbitas o trayectorias de las particulas
de agua al 1er orden de aproximación. Se observa que las trayectorias son elipses
cerradas, luego de un periodo t=T las articulas vuelven a su punto de salida t=0,
lo cual confirma que en la terı́a lineal de 1er orden, en 1 periodo de onda no hay
transporte de masa.
En la superficie libre z = 0 el radio menor B = a, lo cual es compatible con el
desplazamiento de la superficie libre, y en el fondo z = −h el radio menor B = 0, lo

24 IPC Capı́tulo 3
Draft Notas de curso. Teoria de ondas cortas de pequeña amplitud en superficie libre

cual es compatible con la condición de borde de fondo que la velocidad vertical es

nula en z = −h.
La amplitud del eje mayor en el fondo es A = a/ sinh(kh). Si se compara esta
expresión con el valor del modulo de la velocidad orbital horizontal ûb en el fondo
z = −h, ûb = σ a/ sinh(kh) = σ A. O sea la velocidad de fondo se vincula con la
amplitud de la órbita multiplicado or la frecuencia.

Figura 3.6: Esquema trayectorias al 1er orden. Tomado de Andersen and Frigaard
2011

Caso aguas profundas kh → ∞

Para el caso de aguas profundas, la función hiperbólica se puede expresar como

cosh k(h + z) ekh ekz + e−kh e−kz

A(x, z) = a =a
sinh kh ekh − e−kh

como e−kh → 0
ekh ekz + e−kh e−kz
A(x, z) = a = a (ekz + e−2kh e−kz ) ≈ a ekz
ekh
siendo z negativo hacia abajo, De la misma forma

B(x, z) ≈ a ekz

Se observa que las amplitudes de la elipse son iguales, o sea las órbitas son
circulares con igual radio.

Capı́tulo 3 IPC 25
Notas de curso. Teoria de ondas cortas de pequeña amplitud en superficie libre Draft

A la profundidad z = −L/2,
A = B = a e(2π/L)(−L/2) = a e−π = 0,072 a = 0,036 H,
con lo cual a mayores profundidades que L/2 la penetración del oleaje es des-
preciable.
Por esta razón, la condición de aguas profundas se toma h/L > 1/2.

Caso aguas poco profundas kh → 0

Para el caso de aguas poco profundas, la función hiperbólica se puede expresar

como
cosh k(h + z) a
A(x, z) = a =
sinh kh kh
La amplitud horizontal de la órbita es constante en la vertical.

sinh k(h + z) k(h + z)

B(x, z) = a =a = a(1 + z/h)
sinh kh kh
La amplitud vertical de la órbita varia desde B = 0 en el fondo z = −h y crece
linealmente hasta el valor B = a en la superficie libre z = 0.
La Figura 3.6 presenta un esquema de las órbitas en los casos considerados.

Figura 3.7: Esquema trayectorias al 1er orden.

3.7. Energı́a de la ola

En la sección anterior se estudio la cinemática de la ola, esto es su campo de
velocidades y de presiones. En esta sección se estudia la energı́a de la ola y la
velocidad con que se transmite la energı́a.
Cuando hablamos de energı́a ondulatoria, normalmente pensamos en el contenido
de energı́a mecánica, es decir, energı́a cinética y potencial. La energı́a cinética se
origina en el movimiento de las partı́culas de agua y la energı́a potencial se origina
por el desplazamiento de la superficie del agua desde una superficie plana horizontal.

26 IPC Capı́tulo 3
Draft Notas de curso. Teoria de ondas cortas de pequeña amplitud en superficie libre

No obstante hay otras componentes de la energı́a, en particular la energı́a perdida

por fricción con el fondo o la energı́a disipada por la rotura del oleaje, temas que se
trataran en otros capı́tulos.
Se puede mostrar que la energı́a en la onda se propaga en la dirección de propa-
gación de la onda, o de otra manera, la onda se propaga en la dirección en la que se
propaga la energı́a.

3.7.1. Energı́a promedio de la ola

La energı́a de la ola esta compuesta por su energı́a cinética y potencial.
La energı́a cinética se expresa como c (x, z, t) = 1/2ρ(u2 + w2 ), lo cual expresa
la energı́a cinética por unidad de volumen en función del punto y tiempo (x,z,t). En
teorı́a de oleaje, es muy utilizada la energı́a de la ola integrada en toda la columna
de agua y promediada en el perı́odo de la ola T o en su longitud de onda L. De esta
forma, se obtiene una energı́a media, que localmente, a escala de la longitud de onda
o del perı́odo, es constante.
Por simplicidad llamamos θ = kx − σt.

Z η
1 σa 2
Ec = ρ( ) [cosh2 k(h + z) cos2 θ + sinh2 k(h + z) sin2 θ] dz
−h 2 sinh kh

Z η
1 σa 2
Ec = ρ ( ) [cosh2 k(h + z) cos2 θ + sinh2 k(h + z) sin2 θ] dz
2 sinh kh −h

En primera aproximación la integral puede calcularse entre z=-h y z=0, ya que

los términos del integrando son cuadráticos:
Z 0
1 σa 2
Ec = ρ ( ) [cos2 θ + sinh2 k(h + z)] dz
2 sinh kh −h

donde se ha utilizado la relación cosh2 k(h + z) − sinh2 k(h + z) = 1.

Resolviendo la integral:

1 σa 2 sinh k(h + z) cosh k(h + z) k(h + z) 0

Ec = ρ ( ) [cos2 θ z + ( − )]−h
2 sinh kh 2k 2k

1 σa 2 sinh kh cosh kh kh
Ec = ρ ( ) [cos2 θ h + − ]
2 sinh kh 2k 2k
Utilizando la ecuación de dispersión σ 2 = gk tanh kh
1 gk 1 sinh kh cosh kh
Ec = ρa2 [(cos2 θ − ) h + ]
2 sinh kh cos kh 2 2k

Capı́tulo 3 IPC 27
Notas de curso. Teoria de ondas cortas de pequeña amplitud en superficie libre Draft

1 k 1 1
Ec = ρga2 [ (cos2 θ − ) h + ]
2 sinh kh cos kh 2 2
Ahora se promedia en una longitud de onda. Como cos2 θ = 1/2(1 + cos 2θ), el
promedio es 1/2, y la energı́a cinética promedio en una longitud de onda es

1
Ec = ρga2 (3.40)
4

Para el cálculo de la energı́a potencial, se supone que el fluido es incompresible y

se deprecia nuevamente la tensión superficial, toda la energı́a potencial se origina en
las fuerzas gravitacionales. Solo se considera la energı́a causada por el desplazamiento
de la superficie del agua desde el nivel medio del agua. Con estos supuestos podemos
escribir el valor instantáneo de la energı́a potencial p (θ) por unidad de área en el
plano horizontal como:

Z η Z 0 Z η
Ep = ρgzdz − ρgzdz = ρgzdz
−h −h 0

1 1 1 1
Ep = ρgη 2 = ρga2 cos2 θ = ρga2 (1 + cos 2θ)
2 2 2 2
Promediando en una longitud de onda se tiene la energı́a potencial promedio en
la columna de agua
1
Ep = ρga2 (3.41)
4
En consecuencia, la energı́a total, cinética y potencial promedio en una longitud
de onda esta dada por

1
E = Ec + Ep = ρga2 (3.42)
2
Se hace notar que la energı́a potencial y cinética de se distribuyen en partes
iguales. Esta energı́a es una energı́a por unidad de area N.m/m2 .

3.7.2. Flujo de Energı́a de la ola

En la subsección anterior se calculo la energı́a o densidad de energı́a promedio en
una longitud de onda. Ahora debemos considerar la energı́a promedio (durante un
perı́odo) que se transporta a través de una sección vertical fija integrada a lo largo
de la profundidad. Esta sección vertical es paralela a los frentes de onda y tiene un
ancho de 1 m, se denomina flujo de energı́a transportado medio o simplemente el
flujo de energı́a FE .

28 IPC Capı́tulo 3
Draft Notas de curso. Teoria de ondas cortas de pequeña amplitud en superficie libre

Flujo de Energı́a: Potencia fuerzas externas

El flujo de energı́a se puede calcular utilizando la expresión (2.11).

En forma simple, el flujo de energı́a se puede calcular como el trabajo producido
por fuerzas externas: en un elemento vertical dz actúa la fuerza de presión horizontal
pdz. Durante el intervalo de tiempo dt, el elemento se mueve la distancia udt hacia
la derecha. El trabajo producido por unidad de ancho es fuerza x distancia: ∆E =
pu∆t. Entonces el flujo de energı́a se puede calcular como FE = ∆E/∆t = p u.
La presión esta dada por la suma de la parte dinámica más la componente
hidrostática p = pd − ρgz.
Si se desprecia el flujo de energı́a cinética (velocidades al cubo), que es de un
orden mayor al de las presiones, el flujo puede calcularse como:
Z 0
FE = pd (x, z, t) u(x, z, t) dz
−h

Sustituyendo la expresión de la componente dinámica de la presión y de la velocidad:

Z 0
cosh k(h + z) cosh k(h + z)
FE = ρga cos θ σa cos θ dz
−h cosh kh sinh kh
0
ρga2 σ
Z
FE = cos2 θ cosh2 k(h + z) dz
cosh kh sinh kh −h

ρga2 σ sinh kh cosh kh kh

FE = cos2 θ( + )
cosh kh sinh kh 2k 2k
ρga2 σ kh
FE = cos2 θ(1 + )
2k cosh kh sinh kh
ρga2 σ 1 2kh
FE = (1 + cos 2θ)(1 + )
2k 2 sinh 2kh
Promediando nuevamente en 1 longitud de onda , utilizando la definición de energı́a
E y C = σ/k
C 2kh
FE = E (1 + ) = E Cg (3.43)
2 sinh 2kh

C 2kh
Cg = (1 + ) (3.44)
2 sinh 2kh
Cg es la velocidad con la cual se transfiere la densidad de energı́a E y se la
denomina celeridad de grupo .
En aguas profundas, Cg = C/2, o sea la energı́a se transfiere a la mitad de la
celeridad de la onda. Mientras que en aguas poco profundas Cg = C.
En un tren de ondas establecido, la forma de la onda viaja a celeridad C pero la
energı́a se trasfiere a velocidad Cg .

Capı́tulo 3 IPC 29
Notas de curso. Teoria de ondas cortas de pequeña amplitud en superficie libre Draft

Pero si el tren de ondas avanza en un medio no perturbado, en reposo, el frente de

onda debe viajar a la velocidad Cg que es la velocidad con que se transfiere energı́a
desde la zona perturbada a la zona en reposo.
Por esta razón, el frente de onda de un tren de ondas, en general avanza mas
lentamente que las ondas en la zona perturbada. Un observador que viaja con el
frente de ondas, ve llegar las ondas, que se van desintegrando cuando llegan a su
frente, justamente transfiriendo ası́ su energı́a a la zona no perturbada.
En aguas poco profundas, Cg = C la velocidad del grupo es igual a la velocidad
de fase, por lo que las ondas individuales viajan a la misma velocidad que la energı́a.
Esto significa que, en aguas poco profundas no desaparecen olas en el frente de onda.

3.7.3. Grupo de ondas

Una forma sencilla de analizar el agrupamiento de olas, se puede realizar si
sumamos dos ondas armónicas (η1 y η1 , ver Figura 3.7.3, con frecuencias ligeramente
diferentes, viajando en la misma dirección.
Estas dos ondas se sumarán entre sı́, y en un momento cuando están en fase, es
decir, cuando las crestas de las dos ondas componentes coinciden, la amplitud será
el doble de la onda incidente, pero se cancelan entre sı́ en otro momento (cuando
están desfasadas 180 grados, es decir, cuando la cresta de una onda coincide con la
depresión de la otra). Esto se repetirá una y otra vez, en otras palabras, se genera
una serie de grupos de ondas.

η = η1 + η2 = a cos(k1 x − σ1 t) + a cos(k2 x − σ2 t)
(k1 x − σ1 t) − (k2 x − σ2 t (k1 x − σ1 t) + (k2 x − σ2 t
η = 2a cos( ) cos( )
2 2
(k1 − k2 )x − (σ1 − σ2 )t (k1 + k2 )x − (σ1 + σ2 )t
η = 2a cos( ) cos( )
2 2
Se observan que hay dos ondas. Por un lado
(k1 + k2 )x − (σ1 + σ2 )t
cos( )
2
representa una onda cuya frecuencia y número de onda es el promedio de la frecuen-
cia y numero de onda de cada componente. Esta es la onda rápida.
Por otro lado, hay otra onda
(k1 − k2 )x − (σ1 − σ2 )t
2a cos( )
2
cuya frecuencia y número de onda es la diferencia entre las frecuencias y número de
onda de cada componente. Como ∆k = k1 − k2 y ∆σ = σ1 − σ2 . La amplitud de la

30 IPC Capı́tulo 3
Draft Notas de curso. Teoria de ondas cortas de pequeña amplitud en superficie libre

onda esta modulada por esta función.

∆k x − ∆σ t
2a cos( )
2

cuya celeridad es
∆σ ∂σ
= = Cg
∆k ∂k
en la cual σ esta dada por la ecuación de dispersión. En consecuencia, la velocidad
con la cual se propaga la modulación es la celeridad de grupo.

Figura 3.8: Modulación de 2 ondas. Grupo de ondas que viaja hacia la derecha.
Tomado de Holthuijsen 2007.

3.8. Ondas estacionarias. Reflexión.

En la sección anterior, se calculó el potencial de velocidad de una onda progresiva
que verifica la ecuación de Laplace y las condiciones de borde.
La ecuación de Laplace es una ecuación lineal, por lo cual es posible superponer
las diferentes soluciones para encontrar distintos casos de interés.
El caso de la onda progresiva el potencial de velocidades es:

ga cosh k(h + z)
φ(x, z, t) = − sin(kx − σt)
σ cosh kh

η = a cos(kx − σt)

Capı́tulo 3 IPC 31
Notas de curso. Teoria de ondas cortas de pequeña amplitud en superficie libre Draft

Ahora podemos superponer una onda progresiva que viaja en un sentido de

amplitud a con otra otra onda progresiva que viaja en sentido contrario con amplitud
Ra, siendo R un coeficiente, denominado coeficiente de Reflexión. Si R=1 la reflexión
es perfecta. En general 0 < R < 1.

η = a cos(kx − σt) + Ra cos(kx + σt + )

η = aei(kx−σt) + Rae−i(kx+σt+)

η = a(1 + Re−i(2kx++) )ei(kx−σt)

η/a = |A| ei(kx−σt++θ)

donde tan θ = Imag(A)/Real(A) es la diferencia de fase.

Asumiendo que la diferencia de fase en la reflexión es nula,  = 0, se tiene:

A = 1 + Re−i(2kx+) = (1 + R cos 2kx) − iR sin 2kx

|A|2 = AA∗ = (1 + R cos 2kx)2 + (R sin 2kx)2 = (1 + R2 + 2R cos 2kx)

El valor máximo de A se produce cuando cos 2kx = 1, o sea 2kx = 0, 2π, por lo
cual se produce en kx = 0, π. |A| = (1 + R). Si el coeficiente de reflexión vale R=1,
|A| = 2.
El valor mı́nimo de A se produce cuando cos 2kx = −1, o sea 2kx = π, 3π, por
lo cual se produce en kx = π/2, 3π/2. |A| = (1 − R). Si el coeficiente de reflexión
vale R=1, |A| = 0.
La ecuación de la onda es entonces

η = a(1 + R2 + 2R cos 2kx) cos i(kx − σt +  + θ)

La Figura 3.8 muestra un esquema del perfil de la onda y del campo de veloci-
dades.
Considerando R = 1, en los puntos kx = π/2, 3π/2 la amplitud es cero |A| = 0,
y en consecuencia η = 0 para todo tiempo t. Estos puntos se denominan nodos de
la onda, permanecen fijos en el espacio, no se desplazan.
En los puntos kx = 0, π, |A| = 2, por lo cual la amplitud de la onda es el doble
de la amplitud de la onda incidente. Estos puntos se denominan antinodos.
Las velocidades se pueden calcular de la misma forma, a través de la suma de 2
potenciales de velocidad

ga cosh k(h + z) gRa cosh k(h + z)

φ(x, z, t) = − sin(kx − σt) − sin(−kx − σt)
σ cosh kh σ cosh kh

32 IPC Capı́tulo 3
Draft Notas de curso. Teoria de ondas cortas de pequeña amplitud en superficie libre

Figura 3.9: Onda estacionaria. Reflexión. Tomado de Holthuijsen 2007.

Las velocidades están dadas por u − ∇phi:

gka cosh k(h + z)
u(x, z, t) = [cos(kx − σt) − R cos(kx + σt)]
σ cosh kh
gka cosh k(h + z)
u(x, z, t) = [cos kx cos σt+sin kx sin σt−R(cos kx cos σt−sin kx sin σt)]
σ cosh kh
Para R = 1, y utilizando la ecuación de dipersión:
cosh k(h + z)
u(x, z, t) = 2aσ sin kx sin σt
sinh kh
sinh k(h + z)
w(x, z, t) = −2aσ cos kx sin σt
sinh kh
Se observa que en kx = 0, 2π, la velocidad u = 0, o sea en los antinodos la
velocidad horizontal es siempre nula, mientras que la velocidad vertical w es máxi-
ma/mı́nima.
Se observa que en los antinodos, la componente horizontal de la velocidad es
siempre nula, mientras que la vertical toma los valores máximos/mı́nimos. Por esta
razón, fı́sicamente es como si en los antinodos se ubicara una pared impermeable.
Normalmente en situación reales la reflexión no es perfecta, y el coeficiente de re-
flexión R < 1, con lo cual la onda no es perfectamente estacionaria, y se la denomina
parcialmente estacionaria.
La Figura 3.8 muestra un esquema del perfil de la onda parcialmente estacionaria
y del campo de velocidades.

Capı́tulo 3 IPC 33
Notas de curso. Teoria de ondas cortas de pequeña amplitud en superficie libre Draft

Figura 3.10: Onda parcialmente estacionaria. Reflexión. Tomado de Holthuijsen

2007.

Se debe notar, que la costa, aunque no sea una estructura vertical, es un elemento
reflejante, cuyo coeficiente de reflexión depende de las caracterı́sticas de la ola y de
la pendiente de fondo. Asimismo, durane la reflexión puede existir un desfasaje entre
la onda incidente y la reflejada, por lo cual  6= 0.
El coeficiente de reflexión puede estimarse en función del número de Iribarren ζ:

tan α 2
R = 0,1ζ 2 = 0,1( p ) (3.45)
H/Lo

34 IPC Capı́tulo 3
Bibliografı́a
Andersen, T. L., and Frigaard, P. (2011). Lecture Notes for the Course in Water
Wave Mechanics. Aalborg: Department of Civil Engineering, Aalborg University.
DCE Lecture notes, No. 24.
Jean Pierre Germain. Notas de curso. Teorı́a de la ola. Universidad Joseph Fou-
rier, Grenoble I. Francia.
Dean, R.G. and Dalrymple, R.A., 1991. Water wave mechanics for engineers and
scientists. Second printing with correction, World Scientific Publishing Co. Pte. Ltd.,
Singapore.
Holthuijsen, L. H., 2007. Waves in Oceanic and Coastal Waters. Cambridge Uni-
versity Press, UK.

35