Capturando la física de los resonadores Helmholtz

con la ecuación de ondas acústica

Maricel Matar, Reinaldo Welti

Laboratorio de Vibraciones y Ondas, Departamento de Física, FCEIA, Universidad
Nacional de Rosario, Avenida Pellegrini 250, CP 2000, Rosario, Argentina.

E-mail: [email protected]

(Recibido el 2 de Diciembre; aceptado el 10 de Enero de 2009)

Resumen
En este trabajo se muestra cómo a partir de la solución analítica de la ecuación de ondas se puede recuperar el
ingenioso modelo utilizado por Helmholtz para determinar la frecuencia de resonancia más baja de una cavidad
acústica que tiene forma de botella.

Palabras clave: Resonador Helmholtz, ecuación de ondas, elasticidad del aire.

Abstract
In this paper we show how from the analytic solution of the wave equation we can recover the ingenious model used
by Helmholtz to determine the lowest resonance frequency of an acoustic cavity that is shaped like a bottle.

Keywords: Helmholtz resonador, waves equations, air elasticity

PACS: 43.10.Sv, 43.20. Ks, 43.20, Wd ISSN 1870-9095

I. INTRODUCCIÓN resorte. El aire, como es compresible, tiene una constante

elástica [3].
Los primeros resonadores Helmholtz [1] fueron
construidos por H. von Helmholtz en 1860, quien los
utilizó, entre otras cosas, para demostrar que cualquier
sonido musical complejo puede ser considerado como una
mezcla de la frecuencia fundamental y sus armónicos con
diferentes amplitudes. Las teorías de Helmholtz sobre el
sonido han tenido un impacto considerable en su época y
han estimulado a muchos investigadores a realizar estudios
sobre las relaciones entre la física, la psicología y la
música. Rudolph Koenig, el fabricante del resonador que
se muestra en la figura 1, hizo suya la profesión de
“materializar las ideas de Helmholtz en aparatos” [2].
Un ejemplo de la resonancia Helmholtz es el sonido
que se crea cuando uno sopla transversalmente sobre el
cuello de una botella. Cuando se fuerza al aire a
introducirse en una cavidad, la presión interior de la misma
crece. Una vez que el agente externo que fuerza al aire
hacia el interior de la cavidad desaparece, el aire
comprimido del interior fluye hacia afuera. Este flujo de
aire tiende a sobre-compensar la diferencia de presión.
Debido a la inercia del aire en el cuello una pequeña
porción de aire sale al exterior de modo que la presión de
la cavidad disminuye, alcanzando un valor ligeramente FIGURA 1. Resonador Helmholtz construido por R. Koening en
menor que la externa y esto hace que el aire de nuevo fluya el año 1865.
hacia el interior. Este proceso tiende a repetirse de manera
continua, sin embargo, los cambios de presión van Los cambios en las dimensiones de la cámara modifican
decreciendo. las propiedades elásticas del “resorte”: una cámara grande
Este comportamiento es similar a las oscilaciones tendrá una constante menor, y viceversa. El aire en el
libres amortiguadas de una masa atada al extremo de un cuello es la masa del oscilador. El aire en el cuello tiene
resorte. El aire atrapado en la cámara actúa como el una velocidad mayor que el aire en la cámara. Si el cuello
Lat. Am. J. Phys. Educ. Vol. 3, No. 1, Jan. 2009 127 http://www.journal.lapen.org.mx
Maricel Matar y Reinaldo Welti
es más largo su masa será más grande, y viceversa. Si el de llevar a esta masa de aire a su posición de equilibrio.
cuello tiene una sección muy pequeña comparada con la Esta fuerza de restitución Fx se debe a la diferencia de
sección de la cámara la velocidad en el cuello será muy presión δ p , entre la presión interna y la presión externa, y
grande mientras que si su sección aumenta la velocidad se
viene dada por:
reducirá.
En la sección 2, vamos a mostrar cómo Helmholtz se
valió de un modelo basado en los razonamientos de los Fx = δ pS1 , (1)
párrafos anteriores para calcular la frecuencia de
resonancia de una cavidad acústica que tiene la forma de donde S1 es la sección del cuello.
una botella.
Se puede también encontrar la frecuencia de
resonancia de la botella resolviendo la ecuación de ondas
en su interior, con las condiciones de borde adecuadas [4].
Este procedimiento se presenta en la sección 3. Con este
método, sin embargo, se enmascara el elegante
razonamiento físico realizado por Helmholtz [1] [5]. El k
objetivo de este trabajo es mostrar como, a partir de la m
solución analítica de la ecuación de ondas, se puede
recuperar el modelo utilizado por Helmholtz y determinar
las condiciones de su validez. En la sección 4 se muestra x
que cuando la longitud del sistema es pequeña comparada
con la longitud de onda y la sección del cuello pequeña
comparada con la sección de la cámara, la energía
potencial del sistema se concentra en la cámara y la
energía cinética en el cuello. La energía total del sistema FIGURA 2. Modelo para encontrar la frecuencia más baja de
se mantiene constante, alternándose de potencial a cinética oscilación de la botella.
de la misma forma que en un resorte. Utilizando las
expresiones analíticas de estas energías se calcula la La variación de presión δ p se relaciona con la variación
constante elástica efectiva de la cámara. de volumen δ V a través del módulo de compresibilidad
La solución de la ecuación de ondas muestra que la B del gas [6]:
botella, además del modo Helmholtz, tiene infinitas δV
frecuencias de resonancias. El modo Helmholtz tiene lugar δ p = −B . (2)
solamente cuando la impedancia acústica de las dos partes V2
del sistema (cuello y cámara) son muy diferentes. Si el La variación de volumen viene dada por δ V = S1 x ,
radio del cuello aumenta, la frecuencia del modo entonces:
Helmholtz aumenta y tiende a la frecuencia más baja del S2x
tubo cerrado – abierto de sección constante que se obtiene Fx = − B 1 = −keq x , (3)
cuando la sección del cuello se hace igual a la sección de V2
la cámara. donde
S2
keq = B 1 , (4)
V2
II. LA FÓRMULA DE HELMHOLTZ es la constante elástica equivalente del aire contenido en el
volumen principal.
Consideremos una botella como la que se muestra en la
La masa del aire en el cuello de la botella es
figura 2 que tiene un gran volumen V2 y un cuello
m = ρ 0 S1 L1 . Por lo tanto, la frecuencia angular ω viene
estrecho de longitud L1 y sección S1 . El volumen del pico dada por
V1 = S1 L1 es muy pequeño comparado con V2 . La botella k S1
está abierta a la atmósfera que tiene una densidad ρ 0 y ωH = eq = c , (5)
m V2 L1
presión p0 .
Helmholtz [1] [5] para encontrar la frecuencia de una donde c = B / ρ 0 es la velocidad del sonido en el aire.
de las oscilaciones libres que puede realizar este sistema,
supuso que el aire en el cuello se mueve como un pistón Para una botella abierta a la atmósfera, c ≈ 340 m / s , y
sólido mientras que el aire en el volumen principal de la de dimensiones, V2 = 10−3 m3 , S1 = 10−4 m 2 , L1 = 5 × 10−2 m ,
botella se comprime y expande de manera alternativa
como si fuera un resorte. se tiene ωH ≈ 477.15 s −1 o sea f ≈ 75.98 Hz . Si se sopla
Si la masa de aire m en el interior del cuello se suave y sostenidamente la boca de la botella se obtiene un
desplaza, en un cierto instante una distancia x hacia la sonido de esta frecuencia aproximadamente.
derecha como se muestra en la figura 2, la presión interna
desciende y como resultado se obtiene una fuerza que trata
Lat. Am. J. Phys. Educ. Vol. 3, No. 1, Jan. 2009 128 http://www.journal.lapen.org.mx
 Capturando la física del resonador Helmholtz con la ecuación de ondas acústica
III SOLUCIÓN DE LA ECUACIÓN DE ONDAS EN vx ( x, t ) = vx ( x)eiωt . (11)
EL INTERIOR DE UN TUBO DE SECCIÓN
CONSTANTE Si reemplazamos (11) y (10) en la ecuación (6b)
encontramos
Suponemos que la sección transversal del tubo es circular
1
y que sus dimensiones y/o el intervalo de frecuencias vx = ( Ae − i β x − Bei β x ) , (12)
exploradas son tales que solamente se excita el modo Z0
longitudinal de más baja frecuencia. Esto se cumple si la donde Z 0 = ρ0 c es la impedancia característica del fluido.
frecuencia de la onda f es menor que la frecuencia de
corte f c del modo transversal superior más cercano que
viene dado por f c ≈ 0.92c /π a donde c es la velocidad del IV SOLUCIÓN DE LA ECUACIÓN DE ONDAS
sonido y a el diámetro del tubo [7]. EN EL INTERIOR DE UNA BOTELLA
Simbolizamos con p , ρ y vx las pequeñas
Para simplificar el problema vamos a suponer que la
variaciones de la presión, densidad y la componente x de botella tiene una sección circular que abruptamente cambia
la velocidad de un fluido, que se encuentra en el interior de de radio, al pasar del cuello (región 1) a su volumen
un tubo de sección constante, respecto a los valores principal (región 2) como se muestra en la Fig. 3. La
constantes de la presión p0 , de la densidad ρ 0 y de la sección y longitud del cuello son S1 y L1 , mientras que la
velocidad v0 x = 0 . Nota: p es ahora la variación de sección y longitud del volumen principal son S 2 y L2
presión que en la sección anterior designamos con δ p . respectivamente. Como en el cuello y en el volumen
Si retenemos solamente los términos de primer orden principal la sección es constante, en cada una de estas
y si se supone que los movimientos del fluido son regiones la presión debe satisfacer la ecuación de ondas
adiabáticos, la ecuación de conservación de la masa y la (7).
ecuación de movimiento de Euler toman la forma [8]:

∂p ∂v
= −B x , (a )
∂t ∂x
(6)
∂vx 1 ∂p
=− , (b)
∂t ρ0 ∂x

donde B es el módulo de compresibilidad adiabático del

fluido.
De estas ecuaciones se deduce que p (y vx ) satisface
la ecuación de ondas:
FIGURA 3. Botella de sección circular. La sección S2 del
∂2 p 1 ∂2 p volumen principal cambia abruptamente a S1 en el cuello.
− =0. (7)
∂x 2 c 2 ∂t 2
En la región 1, la presión (si no se tiene en cuenta el efecto
Para encontrar la solución armónica de esta ecuación de borde) debe ser nula en x = L1 . La solución armónica
hacemos de la ecuación (7) que satisface esta condición de borde
viene dada por
p( x, t ) = p( x)eiωt . (8)
p1 ( x) = A1 sin β ( L1 − x) . (13)
Remplazando (8) en (7) encontramos
La onda de velocidad del fluido asociada es
2
d p
+ β2p = 0 , (9) A1
dx 2 v1 ( x) = −i cos β ( L1 − x) . (14)
Z0
donde β = ω / c .
La solución de esta ecuación diferencial es inmediata En la región 2, la presión debe ser máxima en x = − L2 . La
y viene dada por solución armónica de la onda de presión que satisface esta
condición viene dada por
p = Ae − iβ x + Bei β x . (10)
p2 ( x) = A2 cos β ( L2 + x) . (15)
La velocidad del fluido asociada a esta onda de presión
viene dada por La onda de velocidad asociada es
Lat. Am. J. Phys. Educ. Vol. 3, No. 1, Jan. 2009 129 http://www.journal.lapen.org.mx
Maricel Matar y Reinaldo Welti
A2 de parámetros concentrados [7] como supuso Helmholtz
v2 ( x) = −i sin β ( L2 + x) . (16) para construir su modelo.
Z0
Si S1 / S 2 = 1 , el sistema resulta un tubo uniforme de
Para empalmar las soluciones en x = 0 debemos establecer longitud L1 + L2 . En esta situación, la ecuación (22) se
las condiciones de contorno en esa sección. Como en una reduce a
onda sonora p p0 y ρ ρ 0 , el aire se comprime muy
poco, entonces podemos suponer que el volumen se tan β L1 = cot β L2 ,
conserva. Si despreciamos el flujo de aire en la dirección
radial, que tiene lugar en las proximidades de x = 0 , la Como, 1 − tan α tan β = (tan α + tan β ) cot(α + β ) , esta
condición de contorno viene dada por [7] ecuación implica que cot(kL1 + kL2 ) = 0 , lo que nos da

S1v1 (0) = S 2 v2 (0) . (17)

β L1 + β L2 = π / 2, 3π / 2, 5π / 2,… ,
La otra condición de contorno es la continuidad de la
que corresponden a los modos de un tubo de longitud
presión,
L1 + L2 cerrado en un extremo y abierto en el otro.
p1 (0) = p2 (0) , (18) Si se mantiene S 2 constante mientras S1 tiende a cero,
se obtiene una “cavidad” de longitud L2 cerrada en sus
en efecto, si p1 (0) ≠ p2 (0) , la masa infinitesimal que se dos extremos. La ecuación (22) se transforma en
encuentra en un pequeño entorno alrededor de tan β L2 = 0 cuyas soluciones son β L2 = π , 2π ,… . Estos
x = 0 tendría una aceleración infinita. autovalores corresponden a los modos de un cilindro de
Las condiciones (17) y (18) nos proporcionan las longitud L2 cerrado en ambos extremos.
ecuaciones En estos dos casos, S1 / S 2 → 1 y S1 / S2 → 0 , la
longitud de onda del modo más bajo es comparable a la
A1 A
cos β L1 − 2 sin β L2 = 0 , (19) dimensión del sistema (la cuarta parte y la mitad,
Z a1 Za2 respectivamente). Estos modos, si bien son los que tienen
la frecuencia más baja, no son modos Helmholtz, ya que
A1 sin β L1 − A2 cos β L2 = 0 , (20) sus longitudes de onda son comparables con las
dimensiones del sistema.
donde Z a1 = Z 0 / S1 y Z a 2 = Z 0 / S2 son las impedancias
acústicas de los medios 1 y 2, respectivamente. A. Cálculo de los autovalores
Los autovalores del sistema se encuentran anulando el
discriminante del sistema, lo que nos da En el caso general, para calcular los autovalores debemos
encontrar, gráfica o numéricamente, las raíces de la
Z a1 sin kL1 sin kL2 − Z a 2 cos kL1 cos kL2 = 0 . (21) ecuación (22). Sin embargo, en lugar de encontrar los
ceros de la ecuación
Esta ecuación puede escribirse de la forma
S1
Z S tan β L1 tan β L2 − =0,
tan β L1 tan β L2 = a 2 = 1 . (22) S2
Z a1 S 2
es más conveniente encontrar los máximos de la función
Si β L1 1 y β L2 1 , la ecuación (22) se reduce a
1
U= , (25)
S1 tan β L1 tan β L2 − S1 / S 2 + ε
β 2 L1 L2 = , (23)
S2
donde ε es un número positivo muy pequeño.
que nos lleva a la fórmula de Helmholtz: En la figura se muestra la gráfica de log10 (U ) en
función de f = ω / 2π , para una botella que tiene las
ωH
S1
β= = (24) dimensiones, L1 = 0.05m , L2 = 0.15m , S1 = 10−4 m 2 y
c V2 L1 −3
S 2 = 6.8 × 10 m . Se supuso que la velocidad del sonido c
2

Como β L1 1 y β L2 1 , la longitud de onda λ es es igual a 343.5 m/s (T = 20º C), ρ0 = 1.2 Kg / m3 y que
mucho mayor que L1 y L2 y el sistema puede considerarse ε = 10−6 . La frecuencia más baja, f H = 75.98 Hz es la
frecuencia del modo Helmholtz. La frecuencia que se

Lat. Am. J. Phys. Educ. Vol. 3, No. 1, Jan. 2009 130 http://www.journal.lapen.org.mx
 Capturando la física del resonador Helmholtz con la ecuación de ondas acústica
obtiene a partir de la fórmula de Helmholtz (ecuación (5)) 1
nos da f H = 76.59 Hz . La diferencia, que es del orden del

presión (u. a.)

1% se debe a la aproximación de tan kL por kL.
0,5

75.98 Hz
0
1148 Hz -0.15 -0.1 -0.05 0 0.05
4 x (m)

15
Log10 (U)

velocidad (u. a.)

2
10

5
0

0
-0.15 -0.1 -0.05 0 0.05
-2
0 500 1000 1500 x (m)
Frecuencia (Hz)

FIGURA 5. Ondas de presión (arriba) y velocidad (abajo), que se

FIGURA 4. Log10(U) en función de las frecuencias. Los picos se calculan con las ecuaciones (28) y (29) para k = 2π f H / c ,
encuentran en la frecuencia de los modos normales. En la figura
se observa el pico correspondiente al modo Helmholtz (75.98 Hz) donde f H = 75.98 Hz .
y el de un modo de frecuencia más alta (1148 Hz).
En la Fig. 5 se grafican las ondas de presión y de velocidad
El modo Helmholtz tiene una longitud de onda. (en unidades arbitrarias) en función de x, para una botella
λ = c / f 4.5 m que es mucho mayor que la longitud total que tiene las dimensiones dadas en el párrafo anterior,
del sistema, L = L1 + L2 = 0.20 m . Sin embargo, la cuando está oscilando en su modo Helmholtz. La onda de
presión tiene una amplitud casi constante en el interior de
longitud de onda del segundo modo (0.30 m) es
comparable con las dimensiones del sistema. la cavidad principal y decrece casi linealmente, en el
cuello, hasta tomar el valor cero en su extremo. La
velocidad del aire en el cuello es mucho más grande que la
velocidad del aire en la cavidad principal .La razón entre
B. Análisis de las ondas estacionarias de presión y estas velocidades es S 2 / S1 (del orden de 68 para este
velocidad
ejemplo).
De la ecuación (20) obtenemos
1

sin β L1
presión (u.a.)

A2 = A1 . (26)
cos β L2 0

Como en general A1 es complejo, hacemos

-1
A1 = A1 eiα . (27) -0.15 -0.1 -0.05 0 0.05
x (m)
Para encontrar las ondas estacionarias de presión y 1
velocidad (u.a.)

velocidad se multiplican las ecuaciones (13) a (16) por

eiωt . Considerando (26) y (27) y tomando la parte real se 0
obtiene
p2 ( x, t ) = pm cos β ( x + L2 ) cos(ωt + α ) , (28a) -1
p
v2 ( x, t ) = m sin β ( x + L2 ) sin(ωt + α ) , (28b) -0.15 -0.1 -0.05 0 0.05
Z0
x (m)
cos β L2
p1 ( x, t ) = pm sin β ( L1 − x) cos(ωt + α ) , (29a) FIGURA 6. Ondas de presión (arriba) y velocidad (abajo), en el
sin β L1 interior de la botella, para el segundo modo cuya autofrecuencia
pm cos β L2 es f = 1148 Hz.
v1 ( x, t ) = cos β ( L1 − x) sin(ωt + α ) , (29b)
Z 0 sin β L1
Este comportamiento de las ondas de presión y de
donde pm = A2 es la amplitud de la onda de presión en el velocidad es coherente con el modelo que elaboró
medio 2. Helmholtz para encontrar la frecuencia propia de este
modo pues de acuerdo a estas gráficas en la cavidad

Lat. Am. J. Phys. Educ. Vol. 3, No. 1, Jan. 2009 131 http://www.journal.lapen.org.mx
Maricel Matar y Reinaldo Welti
principal la energía va a resultar puramente potencial de energía potencial u p 2 ; esto es, uc 2 ≤ h2 2 u p 2 . En el
elástica y en el cuello puramente cinética.
cuello, sin embargo, la energía potencial elástica es un
En la Fig. 6 se muestran las amplitudes de las ondas de
infinitésimo de segundo orden con respecto a la densidad
presión y de velocidad, en el interior del tubo, para el
segundo modo, cuya autofrecuencia es 1148 Hz. La de energía cinética; esto es, u p1 ≤ h12 uc1 .
curvatura que tiene este modo es una evidencia que su Por lo tanto, si despreciamos los términos de segundo
longitud de onda es comparable con las dimensiones del orden en h1 y h2 , la densidad de energía es puramente
sistema. Por el contrario, las gráficas sin curvatura de la potencial en el volumen principal y puramente cinética en
Fig. 5 confirman que, para el modo Helmholtz, la longitud el cuello. En esta aproximación la energía potencial total
de onda es muy grande comparada con las dimensiones del del sistema E p se obtiene integrando la densidad de
sistema.
energía potencial u p 2 sobre el volumen principal, lo que
da:
V. CÁLCULO DE LAS ENERGÍAS CINÉTICA
Y POTENCIAL DEL MODO HELMHOLTZ E p = E pm cos 2 (ωH t + α ) , (34a)

Para el modo Helmholtz, β L1 1 y β L2 1 , por lo donde E pm es la energía potencial máxima que viene dada
tanto, podemos hacer las siguientes simplificaciones: por
cos β ( x + L2 ) 1 , para cualquier x en el intervalo pm2 V2
(− L2 , 0) , sin β ( x + L2 ) β ( x + L2 ) , para cualquier x en el E pm = . (34b)
2Z 0 c
intervalo (− L2 , 0) , cos β ( x + L1 ) 1 , para cualquier x en
el intervalo (0, L1 ) , sin β ( x + L1 ) β ( x + L1 ) , para De la misma manera la energía cinética del sistema Ec se
cualquier x en el intervalo (0, L1 ) y obtiene integrando la densidad de energía cinética uc1
cos β L2 / sin β L1 1/ β L1 y su frecuencia es la sobre el volumen del cuello, así se tiene:
correspondiente al autovalor dado por la ecuación (24).
Utilizando estas aproximaciones, la presión y la Ec = Ecm sin 2 (ωH t + α ) , (35a)
velocidad en el volumen principal y en el cuello vienen
dadas por: donde Ecm es la energía cinética máxima que viene dada
por
p2 ≈ pm cos(ωH t + α ) , ( 30a)
1 p 2 (S L )
Ecm = ρ 0 m 12 1 2 . (35b)
2 ( β L1 ) Z 0
pm
v2 ≈ h2 (1 + x /L2 ) sin(ωH t + α ) , ( 30b)
Z0
Si en (35b) reemplazamos β por ωH / c , ecuación (24), se
puede demostrar que
L1 − x
p1 ≈ pm cos(ωH t + α ) , ( 31a)
L1 Ecm = E pm

pm 1 La energía total del sistema es entonces:

v1 ≈ sin(ωH t + α ) . (31b)
Z 0 h1
E = E pm cos 2 (ωH t + α ) + Ecm sin 2 (ωH t + α ) = E pm = Ecm
La densidad de energía total u , en cada una de las
regiones, es la suma de la densidad de energía potencial Esto es, la energía total del resonador se mantiene
elástica [7] constante. Entonces, hay una transformación permanente
de energía cinética en energía potencial y viceversa, del
p2 mismo modo que en un oscilador masa – resorte en
up = , (32) ausencia de rozamiento.
2Z 0 c

y la densidad de energía cinética A. Cálculo de la constante elástica equivalente de la

cavidad principal
1
uc = ρ0 v 2 . (33)
2 Para calcular la constante elástica equivalente de la
cavidad principal tenemos que relacionar su energía
De las ecuaciones (30a) y (30b) se deduce que la densidad potencial con la amplitud de desplazamiento de la masa,
de energía cinética en el volumen principal uc 2 es un que designaremos como Φ1m .
infinitésimo de segundo orden con respecto a la densidad De la ecuación (31b) obtenemos
Lat. Am. J. Phys. Educ. Vol. 3, No. 1, Jan. 2009 132 http://www.journal.lapen.org.mx
 Capturando la física del resonador Helmholtz con la ecuación de ondas acústica
pm de una botella con los de un sistema masa – resorte, no se
v1m = , (36a) limita sólo a su modo más bajo, o modo Helmholtz, sino
Z 0 h1
también a sus modos de frecuencias superiores.
La utilización de resonadores Helmholtz es cada vez
donde v1m es la máxima amplitud de la velocidad en el más frecuente en los laboratorios del nivel medio [11] [12]
cuello y como v = ∂Φ / ∂t , se tiene y universitario [13] [14], tanto como un prototipo de un
oscilador o como un dispositivo para filtrar sonido en una
v1m = ωH Φ1m . (36b) dada gama de frecuencias, pues son muy fáciles de
construir y las mediciones se realizan mediante
A partir de (36a) y (36b) tenemos: micrófonos, parlantes y detectores de señales que se
encuentran normalmente en los laboratorios de estos
pm = Z 0 h1ωH Φ1m . (37) establecimientos. Son también muy empleados [15] [16]
como amortiguadores de ruido en ductos, vehículos y
máquinas o en nuevos dispositivos experimentales [17]
Reemplazando (37) en (34b) nos queda: para simular el comportamiento de metamateriales. Es por
lo tanto un tema relevante para incluirlo en los cursos de
1 vibraciones y ondas o de ecuaciones diferenciales en las
E pm = keq Φ12m , (38)
2 carreras de ingeniería y de ciencias físicas.

donde keq , la constante elástica equivalente de la cavidad,

es REFERENCIAS
2
S
keq = B
, 1
(39) [1] von Helmholtz, H., On the Sensations of Tone as a
V2 Physiological basis of the Theory of Music, (Dover
resultado que coincide con la expresión encontrada en la Publications, New York, 1954).
sección II. [2] Boring, E. G., Sensation and Perception in the History
of Experimental Psychology, (Appleton–Century, New
York, 1942).
VI. CONCLUSIONES [3] Boyle, R., New Experiments Physico-Mechanicall,
Touching the Spring of the Air and its Effects, (Oxford,
Hemos mostrado que, partiendo de la solución analítica de 1660).
la ecuación de ondas acústica, se puede recuperar la física [4] Crawford, F. S., Lowest modes of a bottle, Am. J. Phys.
del modelo utilizado por Helmholtz para encontrar la 56, 702 (1988).
frecuencia del modo que lleva su nombre. Calcular la [5] Rayleigh, J. W. S., The Theory of Sound, (Dover
frecuencia del modo Helmholtz de esta manera es un poco Publications, New York, 1954).
más laborioso. Sin embargo, la deducción que se obtiene [6] French, A. P., Vibration and Waves, (W.W. Norton,
utilizando esta matemática más dura puede generalizarse a New York, 1971).
otros sistemas físicos. En efecto, las analogías matemáticas [7] Kinsler, L. E. , Frey, A. R., Coppens, A. J. and
basadas en las ecuaciones diferenciales, en muchos casos, Sanders, J. V., Fundamentals of Acoustics, (Wiley, New
son más fructíferas que las analogías basadas en los York, 1982).
mecanismos [9]. Existe una gran variedad de fenómenos [8] Granger, R. A., Fluids Mechanics, (Dover
físicos que se describen por medio de una ecuación de Publications, New York, 1995).
ondas lineal: ondas electromagnéticas, acústicas, ondas [9] Pask, C., Mathematics and the science of analogies,
mecánicas en medios elásticos, etc. En todos estos casos se Am. J. Phys. 71, 526-534 (2003).
puede introducir el concepto de impedancia [7], aún para [9] Khonder, A. N., Rezwan Khan, N., Abwar, A. F. M.,
las ondas cuánticas [10]. En cualquiera de estos fenómenos Transmission line analogy of resonance tunneling
puede aparecer, por lo tanto, un modo tipo Helmholtz. La phenomena: The generalized impedance concept, J. Appl.
condición que se debe cumplir es que se tenga un salto Phys. 63, 5191 - 5193 (1988).
grande en la impedancia de los dos medios y que sus [11] Silverman, M. P. and Worthy, E. R. Musical mastery
dimensiones sean pequeñas comparada con la longitud de of a Coke bottle: Physical modeling by analogy, Phys.
onda asociada. Aunque habitualmente no se lo subraya, el Teach. 36, 70–74 (1998).
modo de oscilación del sistema clásico masa – resorte es [12] Gluck, P., Ben-Sultan, S., Dinur, T., Resonance in
en verdad un modo Helmholtz. Este sistema tiene infinitos Flasks and Pipes, Phys. Teach. 44, 10-15 (2006).
modos de frecuencias superiores, que en la práctica se los [13] Guiguet, A., Welti, R., Supresión de modos de
confunde con los modos normales de un resorte fijo en sus vibración acústicos con un resonador Helmholtz, Rev.
dos extremos, pues el nodo de estos modos está muy Bras. Ens. Fis. 25, 287-293 (2003).
próximo a la masa del oscilador. Sin embargo, si la masa [14] Moloney, M. J., Quality factors and conductances in
es comparable con la masa del resorte, el nodo del segundo Helmholtz resonators, Am. J. Phys. 72, 1035–1039 (2004).
modo está separado de la masa y se pueden observar [15] Pierce, A. D., Acoustics: An Introduction To Its
oscilaciones de amplitud relativamente grandes en esta Physical Principles And Applications, (McGraw-Hill Book
segunda resonancia. Por lo tanto, la analogía de los modos Co., New York, 1981).
Lat. Am. J. Phys. Educ. Vol. 3, No. 1, Jan. 2009 133 http://www.journal.lapen.org.mx
Maricel Matar y Reinaldo Welti
[16] Cameron, T. y Russell, D., Laboratory Instruction in [17] Fang, N., Xi, D., Xu, J., Ambati, M., Srituravanich,
Acoustics and Vibration, (Proceedings of the American W., Sun, C., Zhang, X., Ultrasonic metamaterials with
Society for Engineering Education, ASSE, Washington D. negative modulus, Nature Materials 5, 452 – 456, (2006).
C., 1996).

Lat. Am. J. Phys. Educ. Vol. 3, No. 1, Jan. 2009 134 http://www.journal.lapen.org.mx