<strong>Predicción</strong> <strong><strong>de</strong>l</strong> <strong>tiempo</strong> <strong>de</strong> <strong>reverberación</strong> <strong>en</strong> <strong>salas</strong> mediante procesos <strong>de</strong> difusión <strong>de</strong> la <strong>en</strong>ergía acústica∂w(r, t)– D 1∇ 2 w(r, t) = P(r, t) in V 1,∂t∂w(r, t)– D 2∇ 2 w(r, t) = 0 in V 2,∂t(10)las cuales se usan junto con la ecuación (9) <strong>de</strong> las condiciones<strong>de</strong> contorno. Es importante resaltar que este métodopue<strong>de</strong> ampliarse para su uso <strong>en</strong> un número arbitrario<strong>de</strong> espacios acoplados o con un número arbitrario <strong>de</strong>fu<strong>en</strong>tes sonoras sin añadir más suposiciones. Comoejemplo <strong>de</strong> una respuesta <strong>de</strong> doble p<strong>en</strong>di<strong>en</strong>te, se mo<strong><strong>de</strong>l</strong>arondos <strong>salas</strong> acopladas a través <strong>de</strong> una puerta <strong>de</strong> dim<strong>en</strong>siones0,2 m <strong>de</strong> ancho y 3,9 m <strong>de</strong> alto. La primerasala está recubierta <strong>de</strong> un material absorb<strong>en</strong>te con coefici<strong>en</strong>tealto igual a 0,6, mi<strong>en</strong>tras que la segunda sala sepue<strong>de</strong> consi<strong>de</strong>rar reflectante pues su coefici<strong>en</strong>te <strong>de</strong> absorciónes 0,1. La fu<strong>en</strong>te sonora está situada <strong>en</strong> el primerrecinto <strong>de</strong> mayores dim<strong>en</strong>siones. Al com<strong>en</strong>zar la emisión<strong>de</strong> la fu<strong>en</strong>te, parte <strong>de</strong> la <strong>en</strong>ergía acústica atraviesa laapertura alcanzando el recinto contiguo. La Fig. 7 muestrala respuesta impusional <strong>en</strong>ergética <strong>de</strong> este sistemaacoplado formado por dos recintos. Al cesar <strong>de</strong> emitir lafu<strong>en</strong>te sonora la <strong>en</strong>ergía acústica <strong>en</strong> la sala más absorb<strong>en</strong>tese at<strong>en</strong>úa a mayor velocidad que <strong>en</strong> la sala másreflectante, creándose la primera p<strong>en</strong>di<strong>en</strong>te <strong>de</strong> la caída<strong>en</strong>ergética, tal y como se observa <strong>en</strong> la Fig. 7 (b). Seguidam<strong>en</strong>te,se alcanza un instante don<strong>de</strong> la <strong>en</strong>ergía reman<strong>en</strong>te<strong>en</strong> la sala reflectante es mayor que <strong>en</strong> la sala absorb<strong>en</strong>tey parte <strong>de</strong> esta <strong>en</strong>ergía regresa a la salaemisora apareci<strong>en</strong>do una segunda p<strong>en</strong>di<strong>en</strong>te <strong>de</strong> caída,ver Fig. 7 (d). El instante don<strong>de</strong> se produce este cambio<strong>de</strong> p<strong>en</strong>di<strong>en</strong>te es cuando el vector <strong>de</strong> flujo <strong>de</strong> <strong>en</strong>ergíaacústica, o int<strong>en</strong>sidad, se anula <strong>en</strong> las inmediaciones <strong><strong>de</strong>l</strong>a apertura y se le conoce como punto <strong>de</strong> inflexión.4 ConclusionesLa estimación <strong><strong>de</strong>l</strong> <strong>tiempo</strong> <strong>de</strong> <strong>reverberación</strong> es una herrami<strong>en</strong>tanecesaria <strong>en</strong> el campo <strong>de</strong> la acústica recintos. A día<strong>de</strong> hoy, las formulas clásicas <strong><strong>de</strong>l</strong> <strong>tiempo</strong> <strong>de</strong> <strong>reverberación</strong> <strong><strong>de</strong>l</strong>a acústica estadística se sigu<strong>en</strong> usando ampliam<strong>en</strong>te, olvidándosea m<strong>en</strong>udo el hecho <strong>de</strong> que esta teoría está basada<strong>en</strong> suposiciones que pue<strong>de</strong>n limitar su aplicación.En este artículo se realiza una revisión <strong><strong>de</strong>l</strong> mo<strong><strong>de</strong>l</strong>o <strong><strong>de</strong>l</strong>a ecuación <strong>de</strong> difusión acústica como método alternativopara pre<strong>de</strong>cir el <strong>tiempo</strong> <strong>de</strong> <strong>reverberación</strong> <strong>en</strong> <strong>salas</strong>complejas y con bajo coste computacional. Se ha mostradoque el mo<strong><strong>de</strong>l</strong>o <strong>de</strong> difusión pue<strong>de</strong> consi<strong>de</strong>rarse unaaproximación geométrica-estadística <strong><strong>de</strong>l</strong> campo sonoro<strong>en</strong> <strong>salas</strong> y que ti<strong>en</strong>es sus fundam<strong>en</strong>tos <strong>en</strong> una ecuaciónconocida como la ecuación <strong>de</strong> transfer<strong>en</strong>cia radiativa <strong>de</strong><strong>en</strong>ergía. Las difer<strong>en</strong>tes ext<strong>en</strong>siones que se han <strong>de</strong>sarrollado<strong><strong>de</strong>l</strong> mo<strong><strong>de</strong>l</strong>o <strong>de</strong> la ecuación <strong>de</strong> difusión le permit<strong>en</strong>simular <strong>en</strong>tornos complejos y realistas, don<strong>de</strong> la absorción<strong><strong>de</strong>l</strong> aire es importante y los acoples <strong>en</strong>tre varios volúm<strong>en</strong>espue<strong>de</strong>n afectar a la respuesta acústica <strong>de</strong> lasala. A<strong>de</strong>más, como se ha mostrado <strong>en</strong> el ejemplo <strong>de</strong>una sala alargada, los resultados obt<strong>en</strong>idos por este métodoson <strong>de</strong>p<strong>en</strong>di<strong>en</strong>tes <strong>de</strong> la posición fu<strong>en</strong>te-receptor y<strong>de</strong> la distribución <strong>de</strong> materiales <strong>en</strong> las pare<strong>de</strong>s.Sin embargo, quedan todavía algunos <strong>de</strong>safíos queresolver <strong>en</strong> un futuro <strong>en</strong> forma <strong>de</strong> limitaciones <strong><strong>de</strong>l</strong> mo<strong><strong>de</strong>l</strong>o<strong>de</strong> difusión acústica. Por ejemplo, <strong>en</strong> <strong>salas</strong> con pare<strong>de</strong>scuyas reflexiones sean principalm<strong>en</strong>te especulares, lasestimaciones <strong>de</strong> <strong>tiempo</strong> <strong>de</strong> <strong>reverberación</strong> realizadas porel mo<strong><strong>de</strong>l</strong>o <strong>de</strong> difusión no son consist<strong>en</strong>tes físicam<strong>en</strong>te.A<strong>de</strong>más la ecuación no pue<strong>de</strong> mo<strong><strong>de</strong>l</strong>ar el f<strong>en</strong>óm<strong>en</strong>o <strong>de</strong> ladifracción que sufr<strong>en</strong> las ondas <strong>en</strong> algunos obs táculos.5 Agra<strong>de</strong>cimi<strong>en</strong>tosEste trabajo ha sido subv<strong>en</strong>cionado por el Ministerio<strong>de</strong> Economía <strong>de</strong> España bajo el proyecto con refer<strong>en</strong>ciaTEC2009-14414-C03-01 y el proyecto con refer<strong>en</strong>ciaTEC2012-37945-C02-01, y por la Universidad CatólicaSan Antonio bajo el proyecto con refer<strong>en</strong>cia PMAFI-PI-13/11.6 Refer<strong>en</strong>cias[1] L. L. Beranek, “Acoustics and the Concert Hall,”. J.Acoust. Soc. Am. 57:1258-1262, 1975.[2] H. Ste<strong>en</strong>ek<strong>en</strong>, T. Houtgast, “A Physical Method forMeasuring Speech Transmission Quality,” J. Acoust.Soc. Am., 67:318-326, 1980.[3] T. Hidaka, N. Nishihara, “Objective evaluation ofchamber-music halls in Europe and Japan,” J.Acoust. Soc. Am., 116:357-372, 2004.[4] W. C. Sabine, “Collected Papers on Acoustics; Paper1 - “Reverberation,” Cambridge: Harvard Press, 1927(Reimpreso <strong>en</strong> 1964 por Dover).[5] C. F. Eyring, ‘‘Reverberation time in ‘‘<strong>de</strong>ad’’ rooms,’’ J.Acoust. Soc. Am. 1:217- 241, 1930.[6] G. Millington, ‘A modified formula for reverberation,’’J. Acoust. Soc. Am. 4:69–82, 1932.[7] H. <strong>Arau</strong>-Pucha<strong>de</strong>s, ‘‘An improved reverberation formula,’’Acustica 65:163–180, 1988.[8] M. Hodgson, “Wh<strong>en</strong> is diffuse-field theory accurate?”In Proc. Wallace Clem<strong>en</strong>t Sabine C<strong>en</strong>t<strong>en</strong>nial Symposium,157–160, Cambridge, MA, 1994.[9] P. M. Morse, K. U. Ingard, Theoretical Acoustics. Mc-Graw-Hill, New York, 1968.[10 W. B. Joyce, “Sabine’s reverberation time and ergodicauditoriums,” J. Acoust. Soc. Am., 58:643-655, 1975.[18]revista <strong>de</strong> acústica | Vol. 44 | N. os 1 y 2
<strong>Predicción</strong> <strong><strong>de</strong>l</strong> <strong>tiempo</strong> <strong>de</strong> <strong>reverberación</strong> <strong>en</strong> <strong>salas</strong> mediante procesos <strong>de</strong> difusión <strong>de</strong> la <strong>en</strong>ergía acústica[11] A. Krokstad, S. Strom, S. Sorsdal, “Calculating theacoustical room response by the use of a ray tracingtechnique,” J. Sound Vibration, 8(1): 118-125, 1968.[12] J. B. All<strong>en</strong>, D. A. Berkley, “Image method for effici<strong>en</strong>tlysimulating small-room acoustics,” J. Acoust.Soc. Am., 65(4):943-950, 1979.[13] J. H. Rin<strong><strong>de</strong>l</strong>, “The Use of Computer Mo<strong><strong>de</strong>l</strong>ing inRoom Acoustics,” Journal of Vibro<strong>en</strong>gineering,3:219-224, 2000.[14] M. Vorlän<strong>de</strong>r, “International round robin on roomacoustical computer simulation,” In Proceedings of15th International Congress on Acoustics, 577-580,Trondheim, Norway, 1995.[15] H. Lehnert, “Systematic errors of the ray-tracing algorithm,”Applied Acoustics, 38:207-221, 1993.[16] J. M. Navarro, F. Jacobs<strong>en</strong>, J. Escolano, J. J. Lopez,“A theoretical approach to room acoustic mo<strong><strong>de</strong>l</strong>based on a radiative transfer mo<strong><strong>de</strong>l</strong>,” Acta AcusticaUnited with Acustica, 96:1078– 1089, 2010.[17] J. Picaut, L. Simon, J. D. Polack, “A mathematicalmo<strong><strong>de</strong>l</strong> of diffuse sound field based on a diffusionequation,” Acustica, 83:614-621, 1997.[18] V. Valeau, J. Picaut, M. Hodgson, “On the use of adiffusion equation for room-acoustic prediction,” J.Acoust. Soc. Am., 119:1504-1513, 2006.[19] L. Savioja. Mo<strong><strong>de</strong>l</strong>ing Techniques for Virtual Acoustics.PhD thesis, Helsinki University of Technology,Telecommunications Software and Multimedia Laboratory,Espoo, Finland, 1999.[20] H. E. Bass, H. J. Bauer, L. B. Evans, “Atmosphericabsorption of sound: Analytical expressions,” J.Acoust. Soc. Am., 52(3):821-825, 1972.[21] Y. Jing, N. Xiang, “On boundary conditions for thediffusion equation in room acoustic predictions:Theory, simulations, and experim<strong>en</strong>ts,” J. Acoust.Soc. Am., 123:145-153, 2008.[22] M. R. Schroe<strong>de</strong>r, “New method of measuring reverberationtime,” J. Acoust. Soc. Am., 37:409-412, 1965.[23] J. M. Navarro, J. Escolano, J. J. Lopez, “Implem<strong>en</strong>tationand evaluation of a diffusion equation mo<strong><strong>de</strong>l</strong>based on finite differ<strong>en</strong>ce schemes for sound fieldprediction in rooms,” Applied Acoustics, 73:659-665, 2012.[24] ISO/DIS 3382, Acoustics – Measurem<strong>en</strong>t of roomacoustic parameters – Part 1: Performance spaces,2009.[25] A. Billon, V. Valeau, A. Sakout, and J. Picaut, “On theuse of a diffusion mo<strong><strong>de</strong>l</strong> for acoustically coupledrooms,” J. Acoust. Soc. Am., 120:2043-2054,2006.Envíese a: Revista Española <strong>de</strong> <strong>Acústica</strong> - SEAe-mail: secretaria@sea-acustica.eshttp://www.sea-acustica.esApellidos:Estoy interesado <strong>en</strong>:• ASOCIARME A LA SOCIEDAD ESPAÑOLA DE ACÚSTICA• SUSCRIBIRME A LA REVISTA DE ACÚSTICADirección para correspon<strong>de</strong>ncia:Nombre:C.P.: Ciudad: Provincia:Tel.: Fax: e-mail:C<strong>en</strong>tro <strong>de</strong> trabajo:Puesto <strong>de</strong> trabajo:Dirección:C.P.: Ciudad: Provincia:Tel.: Fax: e-mail:revista <strong>de</strong> acústica | Vol. 44 | N. os 1 y 2 [ 19]