Fundamentos de Oscilaciones, ondas y óptica

Cinemática de las ondas

Felipe Valencia Hernandez fvalenciah@[Link]

Departamento de física, Universidad Nacional de Colombia

[Link]

Septiembre de 2014

F. Valencia 1/35 1 / 35
Repaso

• La clase anterior estudiamos el movimiento de una cuerda

con extremos fijos.
• Ya fuera tomando el límite continuo de una cadena de
osciladores, o escribiendo las ecuaciones de Newton para
cada segmento de la cuerda, encontramos que el
movimiento obedece la ecuación de D’ Alambert:
∂2 ∂2
2
y(x, t) = v 2 2 y(x, t)
∂t ∂x

F. Valencia 2/35 2 / 35
Forma de las soluciones

Las soluciones de esta ecuación de D’Alambert, para la cuerda

con extremos fijos, pueden escribirse como superposición de
modos normales:
X
y(x, t) = Cn sin(kn x) cos(ωn t)
n
con
nπ
kn =
L

r
T
ωn = vkn v=
λ

F. Valencia 3/35 3 / 35
Modos normales

play

F. Valencia 4/35 4 / 35
Superposición de modos

10 r
X nπ T nπ
y(x, t) = sin( x) cos( )
L λ L
n=0

play

F. Valencia 5/35 5 / 35
Ondas viajeras

En todos los fenómenos ondulatorios en una dimensión, las

soluciones se describen como combinaciones de pulsos que
viajan hacia la derecha o hacia la izquierda:

y(x, t) = f (x − vt) + g(x + vt)

play

F. Valencia 6/35 6 / 35
Ondas estacionarias
En el caso particular de los modos normales de la cuerda con
extremos fijos, tendríamos:
yn (x, t) ∝ sin (kn (x − vt)) + sin (kn (x + vt))
De hecho, este tipo de superposición de ondas viajeras
siempre da origen a perturbaciones con puntos estacionarios
en todo tiempo (nodos):
sin(k(x − vt) + φ) + sin(k(x + vt) + φ) = 2 sin(kx + φ) cos(ωt)

nπ − φ
= 0 si x = ∀t
k
Este tipo de ondas se conoce como ondas estacionarias.

play
F. Valencia 7/35 7 / 35
Ondas en general

• Este tipo de perturbaciones que viajan en el espacio con

una velocidad que depende del medio y del tipo de
perturbación, es lo que consideramos una onda.
• En el caso de sistemas unidimensionales, la
descomposición del movimiento como suma de ondas
viajeras con velocidad fija, implica la ecuación de onda de
D’ Alambert.
∂2y ∂2y
y(x, y) = g(x + vt) + h(x − vt) → 2
= v2 2
∂t ∂x

F. Valencia 8/35 8 / 35
Cinemática de las ondas (1d)

Entonces, estamos estudiando las propiedades de las ondas,

asociadas a la ecuación de D’Alambert:
∂2 2 ∂
2
y(x, t) = v y(x, t)
∂t2 ∂x2
Ya no sólo para la cuerda con extremos fijos, sino en cualquier
escenario en que la encontremos.

F. Valencia 9/35 9 / 35
Solución general, ondas viajeras

Siempre la solución será suma de ondas viajeras, pues este

resultado no depende de ninguna de nuestras suposiciones de
la clase anterior:En efecto, sin en lugar de las variables x, t
tomamos como variables independientes u1 = x − vt, y
u2 = x + vt, tenemos:
∂y ∂y ∂y ∂2y ∂2y ∂2y ∂2y
= + → = + + 2
∂x ∂u1 ∂u2 ∂x2 ∂u21 ∂u22 ∂u1 ∂u2

∂y ∂y ∂y ∂2y 2
2∂ y
2
2∂ y
2
2 ∂ y
= −v +v → = v + v − 2v
∂t ∂u1 ∂u2 ∂x2 ∂u21 ∂u22 ∂u1 ∂u2

∂2 2 ∂
2 ∂2y
y = v → = 0 → y = g(u1 ) + h(u2 )
∂t2 ∂x2 ∂u1 ∂u2

F. Valencia 10/35 10 / 35
Solución general, modos normales
También podemos escribir siempre las soluciones en términos
de los modos normales, que corresponden, nuevamente, a
formas de moverse en las que cada parte del sistema se
mueve con la misma frecuencia y en fase con las demás:

yn (x, t) = Fn (x) cos(ωn t + φn )

Que corresponde a una separación de variables:

∂2y ∂2y
y(x, t) = Fn (x)G(t) → = Fn00 (x)G(t) = Fn (x)G̈(t)
∂x2 ∂t2
que implica, en la ecuación de movimiento:

Fn00 (x) G̈(t)

Fn (x)G̈(t) = v 2 Fn00 (x)G(t) → v 2 = = −ω 2
Fn (x) G(t)

F. Valencia 11/35 11 / 35
Modos normales, fronteras

Tenemos, entonces, que los modos normales satisfacen:

d2 d2 ωn2
−ωn2 F (x) = v 2 Fn (x) → Fn (x) = − F (x)
dx2 dx2 v2
Que es la ecuación de un oscilador armónico en x:
ωn
Fn (x) = A cos(kn x + ψn ) kn =
v
Los valores aceptables para kn y ψn dependerán de las
condiciones en la frontera.

F. Valencia 12/35 12 / 35
Algunos casos comunes

• Sistema sin fronteras (propagación libre en el vacío:kn

puede tomar cualquier valor real, k, y hay modos normales
sin(kx) y cos(kx).
• Extremos fijos , y(0) = y(L) = 0,como en nuestra clase
nπ
anterior: kn = L , los modos solo pueden ser del tipo
sin(kn x)
• Condiciones periódicas y(0, t) = y(L, t) 6= 0,kn = 2nπ
L ,y
hay modos normales sin(kx) y cos(kx).
∂y ∂y
• Problemita Frontera libre: ∂x (0, t) = ∂x (L, t) = 0

F. Valencia 13/35 13 / 35
Ortogonalidad de los modos

Para el caso de osciladores discretos vimos que los vectores

que describían los modos normales eran ortogonales:
X
yin = ain cos(ωn t + φn ) → ain aim = δmn
i

llevando al caso continuo, tenemos, entonces:

yn (x, t) = Fn (x) cos(ωn t + φn )

X Z
→ aiν → Fn (x) ain aim → Fn (x)Fm (x)dx = δmn
i

F. Valencia 14/35 14 / 35
Ortogonalidad de los modos...

Por ejemplo, para los modos de la cuerda con extremos fijos

nπ
Fn (x) ∝ sin( x)
L

nπ mπ
sin( x) sin( x) = sin(a + b) sin(a − b)
L L
kn +km kn −km
donde hemos definido a = 2 x b= 2 x entonces:

sin(kn x) sin(km x) = (sin a cos b + cos a sin b)(sin a cos b − cos a sin b)

= sin2 a cos2 b−cos2 a sin2 b = (1−cos2 a) cos2 b−cos2 a sin2 b = cos2 b−cos2 a

1 + cos(2b) 1 + cos(2a) cos(2b) − cos(2a)

=( − )=
2 2 2

F. Valencia 15/35 15 / 35
Ortogonalidad de los modos...
Entonces:
cos( n−m n+m
Z L L
L πx) − cos( L πx)
Z
sin(kn x) sin(km x)dx =
0 0 2

entonces, si n 6= m
Z L π π
L sin((n − m)x)  L sin((n + m)x) 
sin(kn x) sin(km x)dx = − =0
0 2(n − m)π 0 2(n + m)π 0

y sí n = m
Z L
L
sin(kn x) sin(km x) =
0 2
así que, con la normalizacion apropiada, escribiríamos :
r
2
fn (x) = sin(kn (x))
L
(nuevamente hemos usado las minúsculas para la forma
normalizada.)
F. Valencia 16/35 16 / 35
Descomposición en modos normales
Entonces, ya que sabemos que cualquier solución con las
condiciones de frontera propuestas tiene la forma:
X
y(x, t) = Cn fn (x) cos(ωn t + φn )
n

∂y X
(x, t) = − ωn Cn fn (x) sin(ωn t + φn )
∂t n

podemos encontrar las amplitudes de cada modo Cn y sus

fases de forma análoga al caso discreto, proyectando las
posiciones y velocidades un tiempo cualquiera sobre modos
normales:
Z
Cn cos(ωn to + φn ) = y(x, to )fn (x)dx
Z
∂y
−ωn Cn sin(ωn to + φn ) = (x, to )fn (x)dx
∂t
F. Valencia 17/35 17 / 35
 Saturday Morning
Breakfast Cereal
audacity

F. Valencia 18/35 18 / 35
Ondas armónicas

Una onda viajera cuya dependencia con el tiempo y el espacio

es armónica, se llama una onda armónica. Por simplicidad,
escribimos esas ondas como:

yn (x) = Akn e−i(kn x−ωn )t ωn = vkn

Pero recordando que:

• Las cantidades físicas son siempre números
reales:tomamos combinaciones de las exponenciales
complejas que nos den números reales.

F. Valencia 19/35 19 / 35
Solución numérica de la ecuación

La ecuación de D’Alambert
∂ 2 y(x, t) 2
2 ∂ y(x, t)
= v
∂t2 ∂x2
se presta para una solución numérica explícita relativamente
[Link] aproximamos las derivadas con diferencias
centradas:
∂ y (xi , tj )
 
1 y(xi , tj+1 ) − y(xi , tj ) y(xi , tj ) − y(xi , tj−1 )
= −
∂t2 ∆t ∆t ∆t

y(xi , tj+1 ) − 2y(xi , tj ) + y(xi , tj−1 )

=
∆t2

F. Valencia 20/35 20 / 35
Solución numérica de la ecuación

y, análogamente, para las derivadas espaciales:

∂ y (xi , tj ) y(xi+1 , tj ) − 2y(xi , tj ) + y(xi−1 , tj )
2
=
∂x ∆x2
Entonces, la ecuación diría:

y(xi , tj+1 ) − 2y(xi , tj ) + y(xi , tj−1 ) y(xi+1 , tj ) − 2y(xi , tj ) + y(xi−1 , tj )

2
= v2
∆t ∆x2
y podemos despejar:
y(xi , tj+1 ) =
∆t2 ∆t2
−y(xi , tj−1 )+2y(xi , tj )∗(1−v 2 2
)+v 2 (y(xi+1 , tj )+y(xi−1 , tj ))
∆x ∆x2

F. Valencia 21/35 21 / 35
Pulsos y trenes de onda

play

play

F. Valencia 22/35 22 / 35
Superposición

Consideremos, para empezar, dos ondas armónicas viajando

en el espacio libre, con diferentes números de onda, por
ejemplo:

y1 (x, t) = A sin(k1 (x − vt)) y2 (x, t) = A sin(k2 (x − vt)

por supuesto, la combinación de las dos:

y(x, t) = A (sin(k1 (x − vt)) + sin(k2 (x − vt))

es una posible onda, que también viaja con velocidad v hacia la

derecha. cómo se ve este tipo de superposición?.

F. Valencia 23/35 23 / 35
Superposición: modulación

Para ondas de igual magnitud, la superposición se verá como

una onda armónica viajera modulada:
cos(x*(k1−k2)/2)
2

1.5

0.5

-0.5

-1

-1.5

-2
-20 -15 -10 -5 0 5 10 15 20

play

F. Valencia 24/35 24 / 35
Superposición; modulación

En efecto, tenemos que, usando u = x − vt, ∆k = (k1 − k2)/2,

k̄ = (k1 + k2)/2; k1 = k̄ + ∆k, k2 = k̄ − ∆k

sin(k1 u) = sin(k̄u) cos(∆ku) + cos(k̄u) sin(∆ku)

sin(k2 u) = sin(k̄u) cos(∆ku) − cos(k̄u) sin(∆ku)

entonces:

sin(k1 u) + sin(k2 u) = 2 cos(∆ku) sin(k̄u)

F. Valencia 25/35 25 / 35
Superposición de pulsos simétricos

play

Pregunta: cómo pueden reaparecer los pulsos después de

que se han anulado sus amplitudes?

F. Valencia 26/35 26 / 35
Velocidad de cada punto

Si consideramos la velocidad de cada punto para ondas

viajeras (y(x, t) = y(x ± vt)):

∂y(x, t) ∂y(x, t)
ẏ(x) = = ±v
∂t ∂x
Cómo son las velocidades en el ejemplo anterior en el
momento en que se anulan las amplitudes?

F. Valencia 27/35 27 / 35
Superposición; interferencia

Considere ahora la superposición de dos ondas armónicas con

las mismas frecuencias y amplitudes pero alguna diferencia de
fase:

y(x, t) = A(sin(kn x − ωn t) + sin(kn x − ωn t + ))

Si la diferencia de fase  = 0 las amplitudes se sumarán, se

dice que hay interferencia completamente [Link] el
contrario, si  = π las ondas se anularán, se dice entonces que
hay interfefencia completamente destructiva.

F. Valencia 28/35 28 / 35
Interferencia

En los casos intermedios, la resultante de ésta superposición

será una onda armónica con una amplitud diferente. En efecto,
puede mostrar como problema en el cuaderno que:
 
y(x, t) = 2A cos( ) sin(kn x − ωn t + )
2 2

F. Valencia 29/35 29 / 35
Densidades de energía

Para un segmento pequeño de cuerda (o su equivalente para

otras ondas) la energía cinética estaría dada por:
1 ∂y ∆K 1 ∂y
∆K = ∆x λ( )2 → = λ( )2
2 ∂t ∆x 2 ∂t
y la energía potencial sería el trabajo efectuado contra la
tensión para producir la deformación:
p 1 ∂y
∆U = T (ds − dx) = T ( dx2 + dy 2 − dx) ∼ T ( )2 dx
2 ∂x

F. Valencia 30/35 30 / 35
Densidades de energía..

Así que tenemos:

∆U 1 ∂y
= T ( )2
∆x 2 ∂x
Como teníamos que para las ondas viajeras:

∂y(x, t) ∂y(x, t)
= ±v
∂t ∂x
p
con v = T /λ, tendríamos

∆U ∆K
=
∆x ∆x

F. Valencia 31/35 31 / 35
Energía de los modos normales

Nuevamente, recordemos que el movimiento en general se

puede descomponer en modos normales:
X
y(x, t) = Ak ei(kx−ωk t) + c.c.
k

La densidad de energía por cada modo sería de la forma:

∆Ek 1
= λ|Ak |2 ωk2
∆x 2

F. Valencia 32/35 32 / 35
Transporte de energía

Entonces, si consideramos una onda armónica, la cantidad de

energía que se transfiere de un punto a los vecinos por unidad
de tiempo es:
∆Ek ∆Ek ∆x 1
P = = = λ|Ak |2 ωk2 v
∆t ∆x ∆t 2

F. Valencia 33/35 33 / 35
Ondas transversales y longitudinales

En las ondas de deformación en la cuerda que discutimos

anteriormente, el movimiento de las partes del sistema es
perpendicular a la propagación de la onda, estas ondas se
llaman transversales. Por supuesto, en algunos medios es
posible tener también ondas longitudinales, en las que el
desplazamiento sea paralelo a la propagación.

F. Valencia 34/35 34 / 35
Final, final, final!

Es todo por hoy.

Gracias por su atención!

F. Valencia 35/35 35 / 35