<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ústicapor profesionales <strong>de</strong> la acústica para pre<strong>de</strong>cir el camposonoro <strong>en</strong> <strong>salas</strong> <strong>de</strong> todo tipo. Sin embargo, la teoría <strong><strong>de</strong>l</strong>campo difuso ti<strong>en</strong>e limitaciones importantes <strong>de</strong> aplicaciónque suel<strong>en</strong> olvidarse. Fundam<strong>en</strong>talm<strong>en</strong>te, sus prediccionespier<strong>de</strong>n precisión cuando aum<strong>en</strong>ta la irregularida<strong>de</strong>n la geometría, por ejemplo <strong>salas</strong> alargadas, ycuando aum<strong>en</strong>ta la no uniformidad <strong>de</strong> distribución <strong>de</strong> laabsorción [8]. A<strong>de</strong>más, las ecuaciones <strong>de</strong> la acústica estadísticadan como resultado un valor único para toda lasala, el cual no <strong>de</strong>p<strong>en</strong><strong>de</strong> <strong>de</strong> la posición <strong>de</strong> la fu<strong>en</strong>te, <strong><strong>de</strong>l</strong>a posición <strong><strong>de</strong>l</strong> receptor, ni la distribución <strong>de</strong> los materialesabsorb<strong>en</strong>tes perdiéndose la distribución espacial <strong><strong>de</strong>l</strong><strong>tiempo</strong> <strong>de</strong> <strong>reverberación</strong>.Aunque estos mo<strong><strong>de</strong>l</strong>os estadísticos pue<strong>de</strong>n versecomo casos límite <strong>de</strong> la ecuación <strong>de</strong> ondas [9], tales mo<strong><strong>de</strong>l</strong>osse <strong>en</strong>ti<strong>en</strong><strong>de</strong> más rigurosam<strong>en</strong>te como casos límite<strong>de</strong> la acústica geométrica [10] <strong>en</strong> la cual el sonido semo<strong><strong>de</strong>l</strong>a como un conjunto <strong>de</strong> partículas sonoras sin interacción.Los métodos geométricos más importantes quese han aplicado a acústica son los métodos <strong><strong>de</strong>l</strong> trazado<strong>de</strong> rayos [11] y <strong>de</strong> la fu<strong>en</strong>te imag<strong>en</strong> [12]. Estos métodoshan crecido <strong>en</strong> popularidad <strong>de</strong>bido a su capacidad <strong>de</strong>consi<strong>de</strong>rar geometrías complejas y difer<strong>en</strong>tes condicionesacústicas interiores, lo que se ha plasmado <strong>en</strong> numerososprogramas comerciales [13]. G<strong>en</strong>eralm<strong>en</strong>te,ofrec<strong>en</strong> un grado <strong>de</strong> precisión aceptable <strong>en</strong> la predicción<strong>de</strong> difer<strong>en</strong>tes parámetros <strong>de</strong> acústica <strong>de</strong> <strong>salas</strong> [14]. Sinembargo, algunos caminos <strong>de</strong> reflexión importantes pue<strong>de</strong>nllegar a per<strong>de</strong>rse al simular <strong>salas</strong> <strong>de</strong> gran<strong>de</strong>s dim<strong>en</strong>sionesy con múltiples espacios conectados [15]. A<strong>de</strong>más,el aum<strong>en</strong>to <strong>en</strong> el número <strong>de</strong> rayos lanzados paraalcanzar resultados <strong>de</strong> confianza pue<strong>de</strong> llevar <strong>tiempo</strong>s <strong>de</strong>computación prohibitivos.Reci<strong>en</strong>tem<strong>en</strong>te, se ha propuesto un mo<strong><strong>de</strong>l</strong>o teóricobasado <strong>en</strong> la ecuación <strong>de</strong> transfer<strong>en</strong>cia radiativa <strong>de</strong> la<strong>en</strong>ergía para g<strong>en</strong>eralizar los algoritmos <strong>de</strong> predictivos basados<strong>en</strong> la propagación <strong>de</strong> partículas sonoras [16]. Estateoría g<strong>en</strong>eral permite s<strong>en</strong>tar las bases <strong>de</strong> los mo<strong><strong>de</strong>l</strong>osgeométricos y <strong>en</strong>globar <strong>en</strong> ellos el mo<strong><strong>de</strong>l</strong>o <strong>de</strong> la ecuación<strong>de</strong> difusión. El uso <strong><strong>de</strong>l</strong> proceso <strong>de</strong> difusión <strong>de</strong> la <strong>en</strong>ergíaacústica [17] se pue<strong>de</strong> utilizar para pre<strong>de</strong>cir, con precisióny con un coste computacional bajo, el campo sonoro<strong>en</strong> <strong>salas</strong> <strong>de</strong> forma arbitraría y distribución <strong>de</strong> la absorciónno uniforme, principalm<strong>en</strong>te <strong>en</strong> <strong>salas</strong> pocoabsorb<strong>en</strong>tes [18].En este trabajo se repasan los difer<strong>en</strong>tes <strong>de</strong>sarrollosrealizados <strong>en</strong> el mo<strong><strong>de</strong>l</strong>o <strong>de</strong> la ecuación <strong>de</strong> difusión parasu aplicación <strong>en</strong> <strong>salas</strong> <strong>de</strong> formas complejas mostrandosus v<strong>en</strong>tajas fr<strong>en</strong>te a los mo<strong><strong>de</strong>l</strong>os estadísticos. Este artículoestá dividido <strong>en</strong> las sigui<strong>en</strong>tes secciones. En primerlugar se pres<strong>en</strong>tan los principios <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 repasando el mo<strong><strong>de</strong>l</strong>o <strong>de</strong> transfer<strong>en</strong>ciaradiativa <strong>de</strong> la <strong>en</strong>ergía. Seguidam<strong>en</strong>te se explican las difer<strong>en</strong>tesexpansiones que se han ido <strong>de</strong>sarrollando sobrela ecuación <strong>de</strong> difusión que permit<strong>en</strong> mo<strong><strong>de</strong>l</strong>ar los f<strong>en</strong>óm<strong>en</strong>os<strong>de</strong> propagación pres<strong>en</strong>tes <strong>en</strong> <strong>salas</strong> complejas.Finalm<strong>en</strong>te, se pres<strong>en</strong>tan las conclusiones <strong><strong>de</strong>l</strong> trabajo.2 Principios <strong><strong>de</strong>l</strong> mo<strong><strong>de</strong>l</strong>o2.1 El mo<strong><strong>de</strong>l</strong>o <strong>de</strong> transfer<strong>en</strong>cia radiativa <strong>de</strong> <strong>en</strong>ergíaacústica2.1.1 IntroducciónEl sonido, como la luz, es un f<strong>en</strong>óm<strong>en</strong>o <strong>de</strong> ondas.Aunque exist<strong>en</strong> varias difer<strong>en</strong>cias <strong>en</strong>tre luz y sonido,usando las similitu<strong>de</strong>s <strong>en</strong>tre ambos f<strong>en</strong>óm<strong>en</strong>os <strong>de</strong> propagaciónfue posible <strong>de</strong>sarrollar técnicas <strong>de</strong> mo<strong><strong>de</strong>l</strong>ado<strong>de</strong> acústica <strong>de</strong> <strong>salas</strong> conocidas como acústica geométrica.Exist<strong>en</strong> un gran número <strong>de</strong> métodos ampliam<strong>en</strong>tereconocidos basados <strong>en</strong> estas técnicas aunque cadauno <strong>de</strong> ellos suele <strong>de</strong>rivarse por separado [19]. Por lotanto, existía la necesidad <strong>de</strong> establecer unos fundam<strong>en</strong>tosg<strong>en</strong>erales para el mo<strong><strong>de</strong>l</strong>ado geométrico <strong>de</strong> la acústica<strong>de</strong> <strong>salas</strong>, más aun cuando reci<strong>en</strong>tem<strong>en</strong>te la comunidadci<strong>en</strong>tífica ha prestado especial at<strong>en</strong>ción a unmo<strong><strong>de</strong>l</strong>o <strong>en</strong>ergético que hace uso <strong>de</strong> la ecuación <strong>de</strong> difusión[17]. Es posible usar la teoría clásica <strong>de</strong> la transfer<strong>en</strong>ciaradiativa <strong>de</strong> la óptica para proponer un mo<strong><strong>de</strong>l</strong>og<strong>en</strong>eral <strong>de</strong> propagación <strong>de</strong> <strong>en</strong>ergía acústica que permitaunificar los fundam<strong>en</strong>tos <strong>de</strong> una amplia variedad <strong>de</strong> métodosgeométricos acústicos incluy<strong>en</strong>do el mo<strong><strong>de</strong>l</strong>o <strong>de</strong> laecuación <strong>de</strong> difusión [16].2.1.2 Ecuación <strong>de</strong> transfer<strong>en</strong>cia radiativa acústicaEn analogía con la óptica, <strong>en</strong> este mo<strong><strong>de</strong>l</strong>o se suponeque las ondas se comportan como partículas conocidascomo fonones. La interacción <strong>de</strong> las partículas sonorascon el medio se pue<strong>de</strong> <strong>de</strong>scribir mediante una ecuaciónintegro-difer<strong>en</strong>cial. La teoría <strong>de</strong> la transfer<strong>en</strong>cia trabajacon la <strong>en</strong>ergía que se propaga a través o es emitida poruna <strong>de</strong>terminada superficie dA y la proporción <strong>de</strong> partículasque la atraviesan <strong>de</strong>fini<strong>en</strong>do la radiancia sonoraL(r,ŝ ,t) don<strong>de</strong> r repres<strong>en</strong>ta un vector <strong>de</strong> posición,ŝes elvector <strong>de</strong> dirección y t correspon<strong>de</strong> con el <strong>tiempo</strong> y cuyasunida<strong>de</strong>s son Wm –2 sr –1 . La Figura 1 muestra un esquema<strong><strong>de</strong>l</strong> flujo <strong>de</strong> <strong>en</strong>ergía a través <strong>de</strong> un elem<strong>en</strong>to difer<strong>en</strong>cial<strong>de</strong> área dA <strong>de</strong>ntro <strong>de</strong> un elem<strong>en</strong>to difer<strong>en</strong>cial <strong>de</strong>ángulo sólido dΩ.La ecuación <strong>de</strong> transfer<strong>en</strong>cia radiativa acústica se <strong>de</strong>rivausando el principio <strong>de</strong> conservación <strong>de</strong> <strong>en</strong>ergía. Definiremosun rayo sonoro como un grupo <strong>de</strong> partículassonoras que se muev<strong>en</strong> <strong>en</strong> una direcciónŝ<strong>de</strong>ntro <strong>de</strong> unángulo solido dΩ y se suele repres<strong>en</strong>tar como un haz tal[12]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µ s–(1 – a). (3)lPara t<strong>en</strong>er <strong>en</strong> cu<strong>en</strong>ta la absorción <strong><strong>de</strong>l</strong> medio y <strong>de</strong> losobjetos, la fracción <strong>de</strong> radiación inci<strong>de</strong>nte absorbida <strong>en</strong>un camino libre medio se <strong>de</strong>signa como µ a,µ a= m a l , (4)don<strong>de</strong> m es la at<strong>en</strong>uación <strong>de</strong>bida a la absorción <strong><strong>de</strong>l</strong> aire[20] y todas las variables ti<strong>en</strong>e una dim<strong>en</strong>sión <strong>de</strong> m –1 .La ecuación (1) es la ecuación <strong>de</strong> transfer<strong>en</strong>cia radiativaacústica. A difer<strong>en</strong>cia <strong>de</strong> lo que suele ocurrir <strong>en</strong> elcampo <strong>de</strong> la óptica, esta última ecuación <strong>de</strong>p<strong>en</strong><strong>de</strong> <strong><strong>de</strong>l</strong><strong>tiempo</strong> <strong>de</strong>bido a que las partículas sonoras se propagana una velocidad finita.2.1.3 Condiciones <strong>de</strong> contornoFigura 1. Esquema <strong><strong>de</strong>l</strong> flujo <strong>de</strong> <strong>en</strong>ergía a través <strong>de</strong> un elem<strong>en</strong>todifer<strong>en</strong>cial <strong>de</strong> área dA <strong>de</strong>ntro <strong>de</strong> un elem<strong>en</strong>to difer<strong>en</strong>cial <strong>de</strong> ángulosólido dΩ.y como muestra <strong>en</strong> la Fig. 1. Según la teoría <strong>de</strong> la transfer<strong>en</strong>cia,un rayo pier<strong>de</strong> <strong>en</strong>ergía a través <strong>de</strong> la diverg<strong>en</strong>ciay la at<strong>en</strong>uación, incluy<strong>en</strong>do la absorción y la dispersiónfuera <strong><strong>de</strong>l</strong> haz, mi<strong>en</strong>tras que gana <strong>en</strong>ergía <strong>de</strong> la dispersiónque <strong>en</strong>tra <strong>en</strong> el haz <strong>de</strong>s<strong>de</strong> el exterior y <strong>de</strong> las fu<strong>en</strong>tessonoras <strong>en</strong> el medio. Por lo tanto, la ecuación <strong>de</strong> transfer<strong>en</strong>ciaes una expresión matemática <strong>de</strong> los cambiosque sufre la radiación sonora <strong>en</strong> función <strong><strong>de</strong>l</strong> <strong>tiempo</strong> <strong>de</strong>bidoa un balance <strong>de</strong> contribuciones <strong>de</strong> <strong>en</strong>ergía positiva ynegativa:1 ∂L(r, ŝ , t)= –ŝ · L(r, ŝ , t) – m tL(r, ŝ , t)c ∂t+ m s∫ Ω´P(ŝ ´, ŝ )L(r,ŝ ´, t)∂Ω´ + L(r, ŝ ´, t), <strong>en</strong> V (1)don<strong>de</strong> m tes el coefici<strong>en</strong>te <strong>de</strong> at<strong>en</strong>uación como la suma<strong><strong>de</strong>l</strong> coefici<strong>en</strong>te <strong>de</strong> absorción m ay el coefici<strong>en</strong>te <strong>de</strong> dispersiónm s. Estos coefici<strong>en</strong>tes <strong>de</strong>fin<strong>en</strong> las propieda<strong>de</strong>s <strong>de</strong>at<strong>en</strong>uación <strong><strong>de</strong>l</strong> medio mediante absorción y dispersión,si<strong>en</strong>do proporcionales al camino libre medio recorrido.Para una sala con un volum<strong>en</strong> V y área superficial S t, elcamino libre medio <strong>de</strong> la sala l pue<strong>de</strong> calcularse con unaecuación analítica s<strong>en</strong>cilla [10]:l = 4VS t(2)Por lo tanto, la fracción <strong>de</strong> radiación inci<strong>de</strong>nte dispersadapor objetos <strong>en</strong> un camino libre medio se <strong>de</strong>signacomo µ s,La ecuación (1) solo es valida para puntos <strong>en</strong> el interior<strong><strong>de</strong>l</strong> recinto, no sirve para los límites <strong><strong>de</strong>l</strong> mismo. Laradiación sonora que sale <strong>de</strong> las superficies <strong><strong>de</strong>l</strong> recinto<strong>de</strong>be <strong>de</strong>terminarse mediante la resolución <strong>de</strong> unas condiciones<strong>de</strong> contorno. Unas condiciones implícitas o reflectantesintroduc<strong>en</strong> una interpretación física <strong>de</strong> dispersióny absorción <strong>en</strong> los límites <strong>de</strong> la sala, o sea, laspartículas creadas <strong>en</strong> la superficie son el resultado <strong>de</strong> lasreflexiones <strong>de</strong> las partículas que impactan <strong>en</strong> ella.∫ Ω +L(r b, ŝ , t)(ŝ · nˆ)∂Ω == ∫ Ω –R F(r b; ŝ ´, ŝ )L(r b, ŝ ´, t)(ŝ ´ · –nˆ)∂Ω´ <strong>en</strong> ∂V (5)En la última ecuación, R Frepres<strong>en</strong>ta el factor <strong>de</strong> reflexión<strong>de</strong> la superficie con unida<strong>de</strong>s sr –1 <strong>de</strong>finida como laprobabilidad <strong>de</strong> que una partícula <strong>en</strong> el punto r bmoviéndose<strong>en</strong> la dirección ŝ ´ sea reflejada hacia la nueva direcciónŝ , tal y como se muestran <strong>en</strong> el esquema <strong>de</strong> laFig. 2. Finalm<strong>en</strong>te, la ecuación (1), para el interior <strong><strong>de</strong>l</strong> recinto,junto con la ecuación (5), para los límites <strong><strong>de</strong>l</strong> mismo,compon<strong>en</strong> el sistema <strong>de</strong> ecuaciones para el transporte<strong>de</strong> partículas sonoras <strong>en</strong> el mo<strong><strong>de</strong>l</strong>o <strong>de</strong> transfer<strong>en</strong>ciaradiativa acústica.2.2 El mo<strong><strong>de</strong>l</strong>o <strong>de</strong> la ecuación <strong>de</strong> difusiónacústica2.2.1 Ecuación <strong>de</strong> difusiónLa ecuación <strong>de</strong> transfer<strong>en</strong>cia radiativa es difícil <strong>de</strong> resolverdado que ti<strong>en</strong>e seis variables in<strong>de</strong>p<strong>en</strong>di<strong>en</strong>tes. Elnúmero <strong>de</strong> variables pue<strong>de</strong> reducirse realizando unassuposiciones a<strong>de</strong>cuadas sobre el comportami<strong>en</strong>to <strong>de</strong> laspartículas sonoras <strong>en</strong> un medio dispersivo con el objeto<strong>de</strong> obt<strong>en</strong>er un proceso <strong>de</strong> difusión <strong>de</strong> la <strong>en</strong>ergía. En g<strong>en</strong>eral,la ecuación <strong>de</strong> difusión <strong>de</strong>scribe situaciones físicasrevista <strong>de</strong> acústica | Vol. 44 | N. os 1 y 2 [ 13]