Resonador de Helmontz
Resonador de Helmontz
Resonador de Helmontz
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. INTRODUCCIN
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 armnicos con diferentes amplitudes. Las teoras 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 fsica, la psicologa y la msica. Rudolph Koenig, el fabricante del resonador que se muestra en la figura 1, hizo suya la profesin 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 presin 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 presin. Debido a la inercia del aire en el cuello una pequea porcin de aire sale al exterior de modo que la presin de la cavidad disminuye, alcanzando un valor ligeramente menor que la externa y esto hace que el aire de nuevo fluya hacia el interior. Este proceso tiende a repetirse de manera continua, sin embargo, los cambios de presin van decreciendo. Este comportamiento es similar a las oscilaciones libres amortiguadas de una masa atada al extremo de un resorte. El aire atrapado en la cmara acta como el
Lat. Am. J. Phys. Educ. Vol. 3, No. 1, Jan. 2009 127
Los cambios en las dimensiones de la cmara modifican las propiedades elsticas del resorte: una cmara grande tendr una constante menor, y viceversa. El aire en el cuello es la masa del oscilador. El aire en el cuello tiene una velocidad mayor que el aire en la cmara. Si el cuello
http://www.journal.lapen.org.mx
es ms largo su masa ser ms grande, y viceversa. Si el cuello tiene una seccin muy pequea comparada con la seccin de la cmara la velocidad en el cuello ser muy grande mientras que si su seccin aumenta la velocidad se reducir. En la seccin 2, vamos a mostrar cmo Helmholtz se vali de un modelo basado en los razonamientos de los prrafos anteriores para calcular la frecuencia de resonancia de una cavidad acstica que tiene la forma de una botella. Se puede tambin encontrar la frecuencia de resonancia de la botella resolviendo la ecuacin de ondas en su interior, con las condiciones de borde adecuadas [4]. Este procedimiento se presenta en la seccin 3. Con este mtodo, sin embargo, se enmascara el elegante razonamiento fsico realizado por Helmholtz [1] [5]. El objetivo de este trabajo es mostrar como, a partir de la solucin analtica de la ecuacin de ondas, se puede recuperar el modelo utilizado por Helmholtz y determinar las condiciones de su validez. En la seccin 4 se muestra que cuando la longitud del sistema es pequea comparada con la longitud de onda y la seccin del cuello pequea comparada con la seccin de la cmara, la energa potencial del sistema se concentra en la cmara y la energa cintica en el cuello. La energa total del sistema se mantiene constante, alternndose de potencial a cintica de la misma forma que en un resorte. Utilizando las expresiones analticas de estas energas se calcula la constante elstica efectiva de la cmara. La solucin de la ecuacin de ondas muestra que la botella, adems del modo Helmholtz, tiene infinitas frecuencias de resonancias. El modo Helmholtz tiene lugar solamente cuando la impedancia acstica de las dos partes del sistema (cuello y cmara) son muy diferentes. Si el radio del cuello aumenta, la frecuencia del modo Helmholtz aumenta y tiende a la frecuencia ms baja del tubo cerrado abierto de seccin constante que se obtiene cuando la seccin del cuello se hace igual a la seccin de la cmara.
de llevar a esta masa de aire a su posicin de equilibrio. Esta fuerza de restitucin Fx se debe a la diferencia de presin p , entre la presin interna y la presin externa, y viene dada por:
Fx = pS1 ,
(1)
k m
x
La variacin de presin p se relaciona con la variacin de volumen V a travs del mdulo de compresibilidad B del gas [6]: V p = B . (2) V2 La variacin de volumen viene dada por V = S1 x , entonces: S2x Fx = B 1 = keq x , (3) V2 donde S2 keq = B 1 , (4) V2 es la constante elstica equivalente del aire contenido en el volumen principal. La masa del aire en el cuello de la botella es m = 0 S1 L1 . Por lo tanto, la frecuencia angular viene dada por k S1 H = eq = c , (5) m V2 L1 donde c = B / 0 es la velocidad del sonido en el aire. Para una botella abierta a la atmsfera, c 340 m / s , y de dimensiones, V2 = 103 m3 , S1 = 104 m 2 , L1 = 5 102 m , se tiene H 477.15 s 1 o sea f 75.98 Hz . Si se sopla suave y sostenidamente la boca de la botella se obtiene un sonido de esta frecuencia aproximadamente.
http://www.journal.lapen.org.mx
vx ( x, t ) = vx ( x)eit .
(11)
Suponemos que la seccin transversal del tubo es circular y que sus dimensiones y/o el intervalo de frecuencias exploradas son tales que solamente se excita el modo longitudinal de ms baja frecuencia. Esto se cumple si la frecuencia de la onda f es menor que la frecuencia de corte f c del modo transversal superior ms cercano que viene dado por f c 0.92c / a donde c es la velocidad del sonido y a el dimetro del tubo [7]. Simbolizamos con p , y vx las pequeas variaciones de la presin, densidad y la componente x de la velocidad de un fluido, que se encuentra en el interior de un tubo de seccin constante, respecto a los valores constantes de la presin p0 , de la densidad 0 y de la velocidad v0 x = 0 . Nota: p es ahora la variacin de presin que en la seccin anterior designamos con p . Si retenemos solamente los trminos de primer orden y si se supone que los movimientos del fluido son adiabticos, la ecuacin de conservacin de la masa y la ecuacin de movimiento de Euler toman la forma [8]:
v p = B x , t x vx 1 p , = 0 x t (a )
Si reemplazamos (11) y (10) en la ecuacin (6b) encontramos 1 ( Ae i x Bei x ) , vx = (12) Z0 donde Z 0 = 0 c es la impedancia caracterstica del fluido.
(6)
(b)
donde B es el mdulo de compresibilidad adiabtico del fluido. De estas ecuaciones se deduce que p (y vx ) satisface la ecuacin de ondas:
2 p 1 2 p =0. x 2 c 2 t 2
(7)
FIGURA 3. Botella de seccin circular. La seccin S2 del volumen principal cambia abruptamente a S1 en el cuello.
En la regin 1, la presin (si no se tiene en cuenta el efecto de borde) debe ser nula en x = L1 . La solucin armnica de la ecuacin (7) que satisface esta condicin de borde viene dada por
p1 ( x) = A1 sin ( L1 x) .
(8)
(13)
(9)
v1 ( x) = i
A1 cos ( L1 x) . Z0
(14)
En la regin 2, la presin debe ser mxima en x = L2 . La solucin armnica de la onda de presin que satisface esta condicin viene dada por
p2 ( x) = A2 cos ( L2 + x) .
(10)
(15)
La velocidad del fluido asociada a esta onda de presin viene dada por
Lat. Am. J. Phys. Educ. Vol. 3, No. 1, Jan. 2009 129
v2 ( x) = i
A2 sin ( L2 + x) . Z0
(16)
Para empalmar las soluciones en x = 0 debemos establecer las condiciones de contorno en esa seccin. Como en una onda sonora p p0 y 0 , el aire se comprime muy poco, entonces podemos suponer que el volumen se conserva. Si despreciamos el flujo de aire en la direccin radial, que tiene lugar en las proximidades de x = 0 , la condicin de contorno viene dada por [7]
S1v1 (0) = S 2 v2 (0) .
de parmetros concentrados [7] como supuso Helmholtz para construir su modelo. Si S1 / S 2 = 1 , el sistema resulta un tubo uniforme de longitud L1 + L2 . En esta situacin, la ecuacin (22) se reduce a
tan L1 = cot L2 ,
Como, 1 tan tan = (tan + tan ) cot( + ) , esta ecuacin implica que cot(kL1 + kL2 ) = 0 , lo que nos da
(17)
L1 + L2 = / 2, 3 / 2, 5 / 2, ,
que corresponden a los modos de un tubo de longitud L1 + L2 cerrado en un extremo y abierto en el otro. Si se mantiene S 2 constante mientras S1 tiende a cero, se obtiene una cavidad de longitud L2 cerrada en sus dos extremos. La ecuacin (22) se transforma en tan L2 = 0 cuyas soluciones son L2 = , 2 , . Estos autovalores corresponden a los modos de un cilindro de longitud L2 cerrado en ambos extremos. En estos dos casos, S1 / S 2 1 y S1 / S2 0 , la longitud de onda del modo ms bajo es comparable a la dimensin del sistema (la cuarta parte y la mitad, respectivamente). Estos modos, si bien son los que tienen la frecuencia ms baja, no son modos Helmholtz, ya que sus longitudes de onda son comparables con las dimensiones del sistema.
(18)
en efecto, si p1 (0) p2 (0) , la masa infinitesimal que se encuentra en un pequeo entorno alrededor de x = 0 tendra una aceleracin infinita. Las condiciones (17) y (18) nos proporcionan las ecuaciones
A1 A cos L1 2 sin L2 = 0 , Z a1 Za2 A1 sin L1 A2 cos L2 = 0 ,
(19)
(20)
donde Z a1 = Z 0 / S1 y Z a 2 = Z 0 / S2 son las impedancias acsticas de los medios 1 y 2, respectivamente. Los autovalores del sistema se encuentran anulando el discriminante del sistema, lo que nos da
Z a1 sin kL1 sin kL2 Z a 2 cos kL1 cos kL2 = 0 .
(21)
En el caso general, para calcular los autovalores debemos encontrar, grfica o numricamente, las races de la ecuacin (22). Sin embargo, en lugar de encontrar los ceros de la ecuacin
tan L1 tan L2 S1 =0, S2
(22)
Si L1
1 y L2
2 L1 L2 =
S1 , S2
U= (23)
1 , tan L1 tan L2 S1 / S 2 +
(25)
donde es un nmero positivo muy pequeo. En la figura se muestra la grfica de log10 (U ) en funcin de f = / 2 , para una botella que tiene las
(24)
dimensiones,
3
= 106 . La frecuencia ms baja, f H = 75.98 Hz es la frecuencia del modo Helmholtz. La frecuencia que se
http://www.journal.lapen.org.mx
obtiene a partir de la frmula de Helmholtz (ecuacin (5)) nos da f H = 76.59 Hz . La diferencia, que es del orden del 1% se debe a la aproximacin de tan kL por kL.
6
0,5
75.98 Hz 1148 Hz
4
0 -0.15
-0.1
-0.05
0.05
Log10 (U)
-2 0
500
1000
1500
x (m)
Frecuencia (Hz)
FIGURA 4. Log10(U) en funcin de las frecuencias. Los picos se encuentran en la frecuencia de los modos normales. En la figura se observa el pico correspondiente al modo Helmholtz (75.98 Hz) y el de un modo de frecuencia ms alta (1148 Hz).
FIGURA 5. Ondas de presin (arriba) y velocidad (abajo), que se calculan con las ecuaciones (28) y (29) para k = 2 f H / c , donde f H = 75.98 Hz .
El modo Helmholtz tiene una longitud de onda. = c / f 4.5 m que es mucho mayor que la longitud total del sistema, L = L1 + L2 = 0.20 m . Sin embargo, la longitud de onda del segundo modo (0.30 m) es comparable con las dimensiones del sistema.
En la Fig. 5 se grafican las ondas de presin y de velocidad (en unidades arbitrarias) en funcin de x, para una botella que tiene las dimensiones dadas en el prrafo anterior, cuando est oscilando en su modo Helmholtz. La onda de presin tiene una amplitud casi constante en el interior de 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 ms grande que la velocidad del aire en la cavidad principal .La razn entre estas velocidades es S 2 / S1 (del orden de 68 para este ejemplo).
1
(26)
presin (u.a.)
sin L1 . cos L2
(27)
-1 -0.15
-0.1
-0.05
0.05
Para encontrar las ondas estacionarias de presin y velocidad se multiplican las ecuaciones (13) a (16) por eit . Considerando (26) y (27) y tomando la parte real se obtiene p2 ( x, t ) = pm cos ( x + L2 ) cos(t + ) , (28a)
p v2 ( x, t ) = m sin ( x + L2 ) sin(t + ) , Z0 p1 ( x, t ) = pm v1 ( x, t ) =
(28b)
x (m)
FIGURA 6. Ondas de presin (arriba) y velocidad (abajo), en el interior de la botella, para el segundo modo cuya autofrecuencia es f = 1148 Hz.
Este comportamiento de las ondas de presin y de velocidad es coherente con el modelo que elabor Helmholtz para encontrar la frecuencia propia de este modo pues de acuerdo a estas grficas en la cavidad
http://www.journal.lapen.org.mx
principal la energa va a resultar puramente potencial elstica y en el cuello puramente cintica. En la Fig. 6 se muestran las amplitudes de las ondas de presin y de velocidad, en el interior del tubo, para el segundo modo, cuya autofrecuencia es 1148 Hz. La curvatura que tiene este modo es una evidencia que su longitud de onda es comparable con las dimensiones del sistema. Por el contrario, las grficas sin curvatura de la Fig. 5 confirman que, para el modo Helmholtz, la longitud de onda es muy grande comparada con las dimensiones del sistema.
de energa potencial u p 2 ; esto es, uc 2 h2 2 u p 2 . En el cuello, sin embargo, la energa potencial elstica es un infinitsimo de segundo orden con respecto a la densidad de energa cintica; esto es, u p1 h12 uc1 . Por lo tanto, si despreciamos los trminos de segundo orden en h1 y h2 , la densidad de energa es puramente potencial en el volumen principal y puramente cintica en el cuello. En esta aproximacin la energa potencial total del sistema E p se obtiene integrando la densidad de energa potencial u p 2 sobre el volumen principal, lo que da:
E p = E pm cos 2 (H t + ) ,
(34a)
(34b)
De la misma manera la energa cintica del sistema Ec se obtiene integrando la densidad de energa cintica uc1 sobre el volumen del cuello, as se tiene: Ec = Ecm sin 2 (H t + ) , (35a)
( 30a) ( 30b)
donde Ecm es la energa cintica mxima que viene dada por p 2 (S L ) 1 Ecm = 0 m 12 1 2 . (35b) 2 ( L1 ) Z 0 Si en (35b) reemplazamos por H / c , ecuacin (24), se puede demostrar que
( 31a)
Ecm = E pm
v1
(31b)
La densidad de energa total u , en cada una de las regiones, es la suma de la densidad de energa potencial elstica [7] up = p2 , 2Z 0 c (32)
Esto es, la energa total del resonador se mantiene constante. Entonces, hay una transformacin permanente de energa cintica en energa potencial y viceversa, del mismo modo que en un oscilador masa resorte en ausencia de rozamiento.
(33)
De las ecuaciones (30a) y (30b) se deduce que la densidad de energa cintica en el volumen principal uc 2 es un infinitsimo de segundo orden con respecto a la densidad
Lat. Am. J. Phys. Educ. Vol. 3, No. 1, Jan. 2009 132
Para calcular la constante elstica equivalente de la cavidad principal tenemos que relacionar su energa potencial con la amplitud de desplazamiento de la masa, que designaremos como 1m . De la ecuacin (31b) obtenemos
http://www.journal.lapen.org.mx
v1m =
pm , Z 0 h1
(36a)
(36b)
(37)
1 2 keq 1 m, 2
(38)
de una botella con los de un sistema masa resorte, no se limita slo a su modo ms bajo, o modo Helmholtz, sino tambin a sus modos de frecuencias superiores. La utilizacin de resonadores Helmholtz es cada vez ms frecuente en los laboratorios del nivel medio [11] [12] y universitario [13] [14], tanto como un prototipo de un oscilador o como un dispositivo para filtrar sonido en una dada gama de frecuencias, pues son muy fciles de construir y las mediciones se realizan mediante micrfonos, parlantes y detectores de seales que se encuentran normalmente en los laboratorios de estos establecimientos. Son tambin muy empleados [15] [16] como amortiguadores de ruido en ductos, vehculos y mquinas o en nuevos dispositivos experimentales [17] para simular el comportamiento de metamateriales. Es por lo tanto un tema relevante para incluirlo en los cursos de vibraciones y ondas o de ecuaciones diferenciales en las carreras de ingeniera y de ciencias fsicas.
REFERENCIAS
[1] von Helmholtz, H., On the Sensations of Tone as a Physiological basis of the Theory of Music, (Dover Publications, New York, 1954). [2] Boring, E. G., Sensation and Perception in the History of Experimental Psychology, (AppletonCentury, New York, 1942). [3] Boyle, R., New Experiments Physico-Mechanicall, Touching the Spring of the Air and its Effects, (Oxford, 1660). [4] Crawford, F. S., Lowest modes of a bottle, Am. J. Phys. 56, 702 (1988). [5] Rayleigh, J. W. S., The Theory of Sound, (Dover Publications, New York, 1954). [6] French, A. P., Vibration and Waves, (W.W. Norton, New York, 1971). [7] Kinsler, L. E. , Frey, A. R., Coppens, A. J. and Sanders, J. V., Fundamentals of Acoustics, (Wiley, New York, 1982). [8] Granger, R. A., Fluids Mechanics, (Dover Publications, New York, 1995). [9] Pask, C., Mathematics and the science of analogies, Am. J. Phys. 71, 526-534 (2003). [9] Khonder, A. N., Rezwan Khan, N., Abwar, A. F. M., Transmission line analogy of resonance tunneling phenomena: The generalized impedance concept, J. Appl. Phys. 63, 5191 - 5193 (1988). [11] Silverman, M. P. and Worthy, E. R. Musical mastery of a Coke bottle: Physical modeling by analogy, Phys. Teach. 36, 7074 (1998). [12] Gluck, P., Ben-Sultan, S., Dinur, T., Resonance in Flasks and Pipes, Phys. Teach. 44, 10-15 (2006). [13] Guiguet, A., Welti, R., Supresin de modos de vibracin acsticos con un resonador Helmholtz, Rev. Bras. Ens. Fis. 25, 287-293 (2003). [14] Moloney, M. J., Quality factors and conductances in Helmholtz resonators, Am. J. Phys. 72, 10351039 (2004). [15] Pierce, A. D., Acoustics: An Introduction To Its Physical Principles And Applications, (McGraw-Hill Book Co., New York, 1981).
133 http://www.journal.lapen.org.mx
VI. CONCLUSIONES
Hemos mostrado que, partiendo de la solucin analtica de la ecuacin de ondas acstica, se puede recuperar la fsica del modelo utilizado por Helmholtz para encontrar la frecuencia del modo que lleva su nombre. Calcular la frecuencia del modo Helmholtz de esta manera es un poco ms laborioso. Sin embargo, la deduccin que se obtiene utilizando esta matemtica ms dura puede generalizarse a otros sistemas fsicos. En efecto, las analogas matemticas basadas en las ecuaciones diferenciales, en muchos casos, son ms fructferas que las analogas basadas en los mecanismos [9]. Existe una gran variedad de fenmenos fsicos que se describen por medio de una ecuacin de ondas lineal: ondas electromagnticas, acsticas, ondas mecnicas en medios elsticos, etc. En todos estos casos se puede introducir el concepto de impedancia [7], an para las ondas cunticas [10]. En cualquiera de estos fenmenos puede aparecer, por lo tanto, un modo tipo Helmholtz. La condicin que se debe cumplir es que se tenga un salto grande en la impedancia de los dos medios y que sus dimensiones sean pequeas comparada con la longitud de onda asociada. Aunque habitualmente no se lo subraya, el modo de oscilacin del sistema clsico masa resorte es en verdad un modo Helmholtz. Este sistema tiene infinitos modos de frecuencias superiores, que en la prctica se los confunde con los modos normales de un resorte fijo en sus dos extremos, pues el nodo de estos modos est muy prximo a la masa del oscilador. Sin embargo, si la masa es comparable con la masa del resorte, el nodo del segundo modo est separado de la masa y se pueden observar oscilaciones de amplitud relativamente grandes en esta segunda resonancia. Por lo tanto, la analoga de los modos
Lat. Am. J. Phys. Educ. Vol. 3, No. 1, Jan. 2009
[16] Cameron, T. y Russell, D., Laboratory Instruction in Acoustics and Vibration, (Proceedings of the American Society for Engineering Education, ASSE, Washington D. C., 1996).
[17] Fang, N., Xi, D., Xu, J., Ambati, M., Srituravanich, W., Sun, C., Zhang, X., Ultrasonic metamaterials with negative modulus, Nature Materials 5, 452 456, (2006).
134
http://www.journal.lapen.org.mx