Universidade Técnica De Lisboa Instituto Superior
UNIVERSIDADE TÉCNICA DE LISBOA
INSTITUTO SUPERIOR TÉCNICO
ACOUSTIC MODELLING FOR VIRTUAL SPACES
Diogo Gonçalo Franco Falcão Osório de Alarcão
(Licenciado)
Dissertação para obtenção do Grau de Doutor em
Engenharia Electrotécnica e de Computadores
Orientador: Doutor José Luís Bento Coelho
Júri
Presidente: Reitor da Universidade Técnica de Lisboa
Vogais: Doutor Luís Manuel Braga da Costa Campos Doutor Jorge Viçoso Patrício Doutor José Luís Bento Coelho Doutor João Manuel Domingues Perdigão Doutor Afonso Manuel dos Santos Barbosa
Novembro 2005 UNIVERSIDADE TÉCNICA DE LISBOA
INSTITUTO SUPERIOR TÉCNICO
ACOUSTIC MODELLING FOR VIRTUAL SPACES
Diogo Gonçalo Franco Falcão Osório de Alarcão
(Licenciado)
Dissertação para obtenção do Grau de Doutor em
Engenharia Electrotécnica e de Computadores
Orientador: Doutor José Luís Bento Coelho
Júri
Presidente: Reitor da Universidade Técnica de Lisboa
Vogais: Doutor Luís Manuel Braga da Costa Campos Doutor Jorge Viçoso Patrício Doutor José Luís Bento Coelho Doutor João Manuel Domingues Perdigão Doutor Afonso Manuel dos Santos Barbosa
Novembro 2005
ii
To my son Rodrigo and to my daughter Inês
iii Abstract
This thesis presents a new model for the sound energy propagation inside enclosed spaces. The new model is developed starting from the physical first principles and from rigorous mathematical definitions. The theoretical foundation allows the rigorous definition of the energy-based quantities used for describing sound energy propagation inside enclosures and the definition of a general energy balance equation for sound particles, applicable to any room. In addition, new theoretical tools for the study of sound energy fields in rooms are introduced, which are based on the language of functional analysis. The governing equations are thus expressed in terms of linear operators with very convenient properties, which are mathematically detailed.
In addition, a new combined method that solves the derived equations is presented. This new combined method uses an extended mirror image source method solving for the propagation of the specularly reflected energy components inside enclosures, and a time-dependent hierarchical radiosity method solving for the propagation of the diffusely reflected energy components. New algorithmic refinements are introduced in the computer implementation of this combined method. Finally, it is shown that the new proposed method is fast, flexible and accurate enough to be applied as an alternative method for room acoustics simulation.
Key words: room acoustics, sound particles, linear integral operators, diffuse and specular reflections, extended image source method, time-dependent hierarchical radiosity.
iv Resumo
A presente Tese descreve um novo modelo da propagação da energia sonora em recintos fechados. O novo modelo é desenvolvido a partir de primeiros princípios físicos e de definições matemáticas rigorosas. As fundações teóricas obtidas permitem definir de forma rigorosa as quantidades energéticas usadas na descrição da propagação de energia sonora dentro de salas bem como a definição de uma equação geral de balanço de energia para partículas sonoras, aplicável a uma sala qualquer. São introduzidas novas ferramentas teóricas para o estudo de campos de energia sonora em salas, baseadas na linguagem da análise funcional. As equação que governam a propagação de energia sonora são pois expressas em termos de operadores lineares com propriedades convenientes.
Adicionalmente, é apresentado um novo método combinado para a resolução das equações derivadas. Este método combinado utiliza o método das imagens estendido para a solução da propagação das componentes reflectidas especularmente e o método de radiosidade hiérarquica, dependente do tempo, para obter a solução da propagação das componentes reflectidas difusamente. São ainda apresentados vários aperfeiçoamentos utilizados nos algoritmos implementados do método combinado.
Finalmente, mostra-se que o novo método combinado é de computação rápida, fléxivel e suficientemente preciso estabelecendo-se num método alternativo para a simulação de acústica de salas.
Palavras chave: acústica de salas, partículas sonoras, operadores integrais lineares, reflexões difusas e especulares, método das fontes imagem estendido, radiosidade hierárquica dependente do tempo.
v ACKNOWLEDGMENTS
I would like to sincerely thank my advisor Professor J. L. Bento Coelho for giving me the opportunity of conducting this work under his supervison. I thank him for the continuous encouragement and unhesitating support during this endeavour and for the the help in the preparation of this thesis.
I must also thank my family, specially my wife Marta, for her love, understanding and patience at every step along the way. A very special thanks to my parents, especially to my mother, who was the first person that pressed me to start such a difficult task, giving me confidence and encouragement. Finally, to my young son and daughter, I thank them for their great patience when, many times, their father was wholly absorbed in his research and could not give them all the attention they deserve.
This work was partially financially supported by FCT – Portuguese Foundation for Science and Technology under the III QCA of the EU.
v i TABLE OF CONTENTS
Chapter 1 – Introduction...... 1 Chapter 2 – Basic Theory Concepts...... 7 2.1 Introduction...... 7 2.2 Physical Descriptors...... 8 2.2.1 Characterisation of Acoustic Disturbances...... 9 2.2.2 The Wave Equation...... 11 2.2.3 Intensity and Energy Density ...... 16 2.2.4 Harmonic and Non-harmonic Sound Waves...... 20 2.3 Sound Sources and Receivers...... 22 2.4 Sound Radiation...... 24 2.5 Sound Receivers...... 25 2.6 Propagation of Sound Waves...... 26 2.6.1 The Sound Field in Front of a Wall...... 27 2.6.2 Diffraction of Sound Waves ...... 39 2.6.3 The Attenuation of Sound in Air...... 41 2.6.4 Sound Absorbers...... 43 2.6.4.1 Membrane Absorbers...... 43 2.6.4.2 Porous Absorbers ...... 45 2.6.4.3 Resonant Absorbers ...... 46 2.6.4.4 People and Furniture...... 47 2.7 Sound Perception ...... 47 2.7.1 Fundamental Properties of the Human Ear ...... 48 2.7.1.1 Intensity Perception of the Auditory System. Thresholds...... 49 2.7.1.2 Equal Loudness Contours...... 50 2.7.1.3 Frequency Perception of the Auditory System...... 51 2.7.1.4 Critical Bandwidths and Masking...... 52 2.7.1.5 Time Perception...... 53 2.7.1.6 Spatial Sound Perception...... 53 Chapter 3 – Sound Fields in Enclosures...... 57 3.1 Introduction...... 57 3.2 The Wave Equation and Boundary Conditions...... 58 3.3 Natural Modes for the Rectangular Room ...... 59 3.4 Steady-state and Transient Sound Fields inside Enclosures...... 65 3.5 Reverberation Time...... 69 3.5.1 Time Distribution of Reflections...... 71 3.5.2 Spatial-Directional Characteristics of Reverberation. Diffuse Sound Fields. .. 72 3.5.3 Quantitative Measures for Sound Fields inside Enclosures...... 73 3.5.4 Subjective Evaluation of Room Acoustics...... 78 Chapter 4 – Acoustic Modelling of the Sound Field in Enclosures...... 81 4.1 Wave-based Methods...... 82 4.1.1 Analytic Methods ...... 84 4.1.1.1 Modal Decomposition ...... 84 4.1.1.2 Integral-based Methods. The Helmoltz-Huygens-Kirchhoff Theorem....85
v ii 4.1.2 Finite Element Formulations...... 86 4.1.2.1 The Finite Element Method (FEM)...... 86 4.1.2.2 The Boundary Element Method (BEM) ...... 87 4.1.2.3 Finite Difference Methods (FDM)...... 88 4.2 Geometrical Acoustics Based Methods...... 88 4.2.1 The Diffuse Field Theory...... 90 4.2.2 The Mirror Image Source Method (MISM)...... 95 4.2.3 The Ray Tracing Method...... 99 4.2.4 Beam, Cone and Pyramid Tracing ...... 102 4.2.5 Finite Sound Ray Integration Method (FSRIM) ...... 105 4.2.6 The Radiant Exchange Method. The Radiosity Method...... 105 4.2.7 Hybrid Methods ...... 109 Chapter 5 – A New Model for Sound Energy Propagation Inside Enclosures...... 111 5.1 Introduction...... 111 5.2 Sound Particles and Sound Energy ...... 114 5.2.1 Measure Spaces, Measures and Domains...... 114 5.2.2 Sound Particles and the Phase Space ...... 118 5.2.3 Particle Measures and the Phase Space Density...... 120 5.2.4 The Trajectory Space and Particle Events...... 123 5.2.5 The Trajectory Space Flux...... 126 5.3 Transport Equations...... 128 5.3.1 A Transport Equation for Particles...... 129 5.3.2 Boundary Conditions ...... 135 5.3.3 The Integral Version of the Transport Equation for Particles...... 137 5.3.4 Definition of Energy Quantities for describing Sound Fields...... 143 5.3.5 A General Transport Equation for Sound Energy...... 150 5.3.6 Steady-state Regime ...... 155 5.3.7 Time-dependent Regime...... 155 5.4 Functional Analysis Formalism...... 156 5.4.1 Linear Operator Formalism of the Governing Equations...... 158 5.4.1.1 Steady-state Regime ...... 161 5.4.1.2 Time-dependent Regime...... 162 5.4.2 Properties of the Operators ...... 162 5.4.2.1 Norms of the Operators...... 163 5.5 Approximate Solution Methods for the Boundary Integral Equation ...... 171 5.5.1 Truncated Neumann Series and Iterated Kernels...... 171 5.5.2 Successive Approximations...... 173 5.5.3 The Nyström Method ...... 173 5.5.4 Finite Basis Methods ...... 174 5.5.4.1 The Point Collocation Method...... 175 5.5.4.2 The Least Squares Method...... 175 5.5.4.3 The Galerkin Method...... 176 5.6 The Diffuse Reflection Case...... 177 5.6.1 Kuttruff’s Integral Equation ...... 180 Chapter 6 – A New Combined Method for the Prediction of Room Impulse Responses ...... 183 6.1 Equation of Motion for the Acoustic Energy inside an Enclosure ...... 184 6.2 The Transport and Reflection Operator...... 185
viii 6.3 Solution of the Specular Operator Equation...... 192 6.4 Solution of the Diffuse Operator Equation ...... 193 6.5 The Combined Method...... 197 Chapter 7 – Implementation of the New Combined Method ...... 199 7.1 Introduction...... 199 7.2 Geometric Representation of an Enclosure...... 200 7.3 The Combined Prediction Method...... 202 7.3.1 The Extended Mirror Image Source Method...... 203 7.3.1.1 Geometric Calculation of the Mirror Image Sources...... 203 7.3.1.2 Accelerating Techniques...... 207 7.3.1.3 Extension to High Orders of Reflection...... 211 7.3.1.4 Specular Room Impulse Responses ...... 216 7.3.2 The Time-dependent Hierarchical Radiosity Method...... 217 7.3.2.1 Hierarchy Construction and Representation ...... 218 7.3.2.2 Solution of the Time-dependent Hierarchical Radiosity Method...... 224 7.3.2.3 Obtaining Diffuse Room Impulse Responses...... 229 7.3.3 Combination of the Results from Both Methods...... 231 Chapter 8 – Validation and Application of the Combined Method...... 235 8.1 Measurement Setup...... 235 8.2 Case Studies: Acoustic Measurements vs. Simulation Results ...... 238 8.2.1 Room 1: VA2 of IST...... 238 8.2.1.1 Acoustic Measurements in Room VA2 ...... 239 8.2.1.2 Simulation Results for Room VA2...... 243 8.2.2 Room 2: classroom V007 of IST ...... 250 8.2.2.1 Acoustic Measurements in Room V007...... 252 8.2.2.2 Simulation Results for Classroom V007...... 255 8.2.3 Room 3: Meeting Room 01.1 of IST...... 262 8.2.3.1 Acoustic Measurements in Room 01.1 ...... 263 8.2.3.2 Simulation Results for Room 01.1...... 266 8.2.4 Room 4: Room C9 of IST...... 273 8.2.4.1 Acoustic Measurements in Room C9...... 274 8.2.4.2 Simulation Results for Classroom C9 ...... 277 8.2.5 Room 5: Congress Centre Auditorium of IST ...... 284 8.2.5.1 Acoustic Measurements in Congress Centre Auditorium...... 285 8.2.5.2 Simulation Results for Congress Centre Auditorium...... 290 Chapter 9 – Conclusions...... 299 9.1 Summary and General Conclusions...... 299 9.2 Final Conclusions ...... 301 9.3 Future Work ...... 302 References ...... 304 Appendix ...... 312
ix LIST OF FIGURES
Number Page 1. Relative attenuation of spherical wave above impedance plane (grass) compared to free field (dependence on frequency) 34 2. Relative attenuation of spherical wave above impedance plane (grass) compared to free field (dependence on source-receiver distance) 34 3. Modulus of scattering factor for a 10 cm x 10 cm patch with specific impedance 10 over a plane boundary with specific impedance 1, plotted in function of incidence angle and frequency 36 4. Squared modulus of scattering factor for a 10 cm2 rough patch 38 5. Three-dimensional representation of the directional scattering characteristics of a wall 39 6. Notation used for the diffraction of a plane sound wave by the edge of a rigid half plane 40 7. Mean squared pressure of diffracted plane wave from a rigid half-plane. 40 8. Air absorption coefficient plotted in function of sound frequency 43 9. Example of absorption coefficient for a panel with 80 Kg/m2 density, plotted as function of frequency 45 10. Example of absorption coefficient for a porous absorber 46 11. Example of absorption coefficient for a resonant absorber 47 12. Sketch of the human ear and cross section of the cochlea duct 48 13. Equal Loudness Contours 50 14. Example of a reverberation decay curve of a room 70 15. Echogram with temporal structure of the room energy impulse response 72 16. Example of geometrical construction of mirror image sources of an arbitrary room 96 17. Finite element formulation with constant basis functions 194 18. Validity and visibility test for potential mirror image sources 204 19. Example of the validity and visibility test for potential mirror image sources 205 20. Back-face Culling accelerating technique 208 21. Example of two impossible polygon combinations 208
x 22. Example of a view frustum for discarding subsequent higher order images 209 23. Example of clustering of input polygons into one single “parent” polygon 210 24. Example of least squares fit of the number of visible images 213 25. Example of least squares fit of the distance of visible images to a receiving point 214 26. Example of specular echograms calculated with the EMISM with statistical extension to higher orders 217 27. Polygons are substructured and interact at an appropriate level. The corresponding quadtrees are shown 219 28. Hierarchical subdivision and links at various levels 221 29. Example of ray-casting for determining the occluded form factor between two polygons 222 30. Example of the subdivision of a convex polygon with six sides in six child polygons 223 31. Input polygons for a room 229 32. Mesh of the room, after the hierarchical refinement 229 33. Specular energy impulse response obtained with the implemented EMISM 232 34. Diffuse energy impulse response obtained by the implemented time-dependent hierarchical method 233 35. Total energy impulse response obtained by the combined method 233 36. Acoustic measurements setup 236 37. Anechoic chamber measurements 237 38. Auditorium VA2, looking at the backside of the room 238 39. Auditorium VA2, looking at the front side of the room 239 40. Wireframe drawing of the model of Auditorium VA2 240 41. Coloured legend for the materials used in the model of room VA2 244 42. 1000 Hz Specular energy impulse response for combination S1-A, room VA2 (5 ms integrated; linear scale) 245 43. 1000 Hz Diffuse energy impulse response for combination S1-A, room VA2 (5 ms integrated; linear scale) 245 44. 1000 Hz Total energy impulse response for combination S1-A, room VA2 (5 ms integrated; linear scale) 246
xi 45. 1000 Hz Total energy impulse response for combination S1-A, room VA2 (5 ms integrated; logarithmic scale 246 46. 1000 Hz Schröder backwards total energy impulse response for combination S1-A, room VA2 (logarithmic scale). 247 47. View of classroom V007 251 48. Another view of classroom V007 251 49. Wireframe drawing of the model of classroom V007 252 50. Coloured legend for the materials used in the model of room V007 256 51. 1000 Hz Specular energy impulse response for combination S1-C, room V007 (5 ms integrated; linear scale) 257 52. 1000 Hz Diffuse energy impulse response for combination S1-C, room V007 (5 ms integrated; linear scale) 257 53. 1000 Hz Total energy impulse response for combination S1-C, room V007 (5 ms integrated; linear scale) 258 54. 1000 Hz Total energy impulse response for combination S1-C, room V007 (5 ms integrated; logarithmic scale) 258 55. 1000 Hz Schröder backwards total energy impulse response for combination S1-C, room V007 (logarithmic scale) 259 56. View of meeting room 01.1 towards the front 262 57. View of meeting room 01.1 towards the back 263 58. Wireframe drawing of the model of meeting room 01.1 264 59. Coloured legend for the materials used in the model of meeting room 01.1 267 60. 1000 Hz Specular energy impulse response for combination S1-C, room 01.1 (1/225 ms integrated; linear scale) 268 61. 1000 Hz Diffuse energy impulse response for combination S1-C, room 01.1 (1/225 ms integrated; linear scale) 268 62. 1000 Hz Total energy impulse response for combination S1-C, room 01.1 (1/225 ms integrated; linear scale) 269 63. 1000 Hz Total energy impulse response for combination S1-C, room 01.1 (1/225 ms integrated; logarithmic scale) 269 64. 1000 Hz Schröder backwards total energy impulse response for combination S1-C, room 01.1 (logarithmic scale) 270
xii 65. Classroom C9, looking at the frontside of the room 273 66. Classroom C9, looking at the backside of the room 274 67. Wireframe drawing of the model of classroom C9 275 68. Coloured legend for the materials used in the model of classroom C9 277 69. 1000 Hz Specular energy impulse response for combination S1-A, classroom C9 (1/150 ms integrated; linear scale) 279 70. 1000 Hz Diffuse energy impulse response for combination S1-A, classroom C9 (1/150 ms integrated; linear scale) 279 71. 1000 Hz Total energy impulse response for combination S1-A, classroom C9 (1/150 ms integrated; linear scale) 280 72. 1000 Hz Total energy impulse response for combination S1-A, classroom C9 (1/150 ms integrated; logarithmic scale) 280 73. 1000 Hz Schröder backwards total energy impulse response for combination S1-A, classroom C9 (logarithmic scale) 281 74. Congress Centre Auditorium, looking at the frontside of the room 284 75. Congress Centre Auditorium, looking at the backside of the room 285 76. Wireframe drawing of the model of the Congress Centre Auditorium 286 77. Coloured legend for the materials used in the model of the Congress Centre Auditorium 291 78. 1000 Hz Specular energy impulse response for combination S1-E, Congress Centre Auditorium (5 ms integrated; linear scale) 292 79. 1000 Hz Diffuse energy impulse response for combination S1-E, Congress Centre Auditorium (5 ms integrated; linear scale) 292 80. 1000 Hz Total energy impulse response for combination S1-E, Congress Centre Auditorium (5 ms integrated; linear scale) 293 81. 1000 Hz Total energy impulse response for combination S1-E, Congress Centre Auditorium (5 ms integrated; logarithmic scale) 293 82. 1000 Hz Schröder backwards total energy impulse response for combination S1-E, Congress Centre Auditorium (logarithmic scale) 294
xiii LIST OF TABLES
Number Page 1. Audible ranges of the human ear 49 2. Standard centre, lower and upper frequencies for octave bands 52 3. Intensity levels measured in anechoic chamber for the Meyer Sound UPM-1 loudspeaker 237 4. Acoustic power of the sound source UPM-1 238
5. Measured T30 values - room VA2 241 6. Measured EDT values - room VA2 241
7. Measured Definition values D50 - room VA2 242
8. Measured Clarity values C80 - room VA2 242
9. Measured Lp values - room VA2 243
10. Reverberation times T30 predicted by the combined method - room VA2 247 11. Early decay times EDT predicted by the combined method - room VA2 248
12. Definition values D50 predicted by the combined method - room VA2 248
13. Clarity values C80 predicted by the combined method - room VA2 248
14. Steady-state sound pressure level Lp predicted by the combined method - room VA2 248
15. Difference between predicted and measured values of T30 - room VA2 249 16. Difference between predicted and measured values of EDT - room VA2 249
17. Difference between predicted and measured values of D50 - room VA2 249
18. Difference between predicted and measured values of C80 - room VA2 250
19. Difference between predicted and measured values of Lp - room VA2 250
20. Measured T30 values - room V007 253 21. Measured EDT values - room VA2 253
22. Measured Definition-D50 values - room V007 254
23. Measured Clarity-C80 values – room V007 254
24. Measured Lp values - room V007 255
25. Reverberation times T30 predicted by the combined method - room V007 259 26. Early decay times EDT predicted by the combined method - room V007 259
27. Definition values D50 predicted by the combined method - room V007 260
xi v 28. Clarity values C80 predicted by the combined method - room V007 260
29. Steady-state sound pressure level Lp predicted by the combined method - room V007 260
30. Difference between predicted and measured values of T30 - room V007 261 31. Difference between predicted and measured values of EDT - room V007 261
32. Difference between predicted and measured values of D50 - room V007 261
33. Difference between predicted and measured values of C80 - room V007 261
34. Difference between predicted and measured values of Lp - room V007 262
35. Measured T30 values - room 01.1 264 36. Measured EDT values - room 01.1 265
37. Measured Definition-D50 values - room 01.1 265
38. Measured Clarity-C80 values – room 01.1 265
39. Measured Lp values - room 01.1 266
40. Reverberation times T30 predicted by the combined method - room 01.1 270 41. Early decay times EDT predicted by the combined method – room 01.1 271
42. Definition values D50 predicted by the combined method – room 01.1 271
43. Clarity values C80 predicted by the combined method – room 01.1 271
44. Steady-state sound pressure level Lp predicted by the combined method – room 01.1 271
45. Difference between predicted and measured values of T30 – room 01.1 272 46. Difference between predicted and measured values of EDT – room 01.1 272
47. Difference between predicted and measured values of D50 – room 01.1 272
48. Difference between predicted and measured values of C80 – room 01.1 272
49. Difference between predicted and measured values of Lp – room 01.1 273
50. Measured T30 values - room C9 275 51. Measured EDT values – room C9 276
52. Measured Definition-D50 values – room C9 276
53. Measured Clarity-C80 values – room C9 276
54. Measured Lp values – room C9 277
55. Reverberation times T30 predicted by the combined method – room C9 281 56. Early decay times EDT predicted by the combined method – room C9 282
57. Definition values D50 predicted by the combined method – room C9 282
58. Clarity values C80 predicted by the combined method – room C9 282
59. Steady-state sound pressure level Lp predicted by the combined method – room C9 282
xv 60. Difference between predicted and measured values of T30 – room C9 283 61. Difference between predicted and measured values of EDT – room C9 283
62. Difference between predicted and measured values of D50 – room C9 283
63. Difference between predicted and measured values of C80 – room C9 283
64. Difference between predicted and measured values of Lp – room C9 284
65. Measured T30 values – Congress Centre Auditorium 287 66. Measured EDT values – Congress Centre Auditorium 288
67. Measured Definition-D50 values – Congress Centre Auditorium 288
68. Measured Clarity-C80 values – Congress Centre Auditorium 289
69. Measured Lp values – Congress Centre Auditorium 289
70. Reverberation times T30 predicted by the combined method – Congress Centre Auditorium 294 71. Early decay times EDT predicted by the combined method – Congress Centre Auditorium 295
72. Definition values D50 predicted by the combined method – Congress Centre Auditorium 295
73. Clarity values C80 predicted by the combined method – Congress Centre Auditorium 295
74. Steady-state sound pressure level Lp predicted by the combined method – Congress Centre Auditorium 296
75. Difference between predicted and measured values of T30 – Congress Centre Auditorium 296 76. Difference between predicted and measured values of EDT – Congress Centre Auditorium 296
77. Difference between predicted and measured values of D50 – Congress Centre Auditorium 297
78. Difference between predicted and measured values of C80 – Congress Centre Auditorium 297
79. Difference between predicted and measured values of Lp – Congress Centre Auditorium 297 80. Octave bands Sabine absorption coefficients 312
xvi LIST OF SYMBOLS
Symbol Page
p0 Ambient pressure 7 p Acoustic pressure 8
pinst. Instantaneous pressure 8
Lp Sound pressure level 8
pref Standard reference pressure 8 i Vector 8 v Particle velocity 9 ρ Fluid density 9
ρ0 Ambient fluid density 9 V Volume 9
SS Surface (area) 9 ∇i Divergence 10 c Speed of sound 11 ∇2 Laplacian 12 Φ Velocity potential 12 k Wave number 13 f Sound frequency 13 ω Angular frequency 13 λ Wavelength 13 T Period 13 ˆi Complex quantity 14 r Radial coordinate 15 I Acoustic intensity 16 w Acoustic energy density 16 n Unit normal vector 17
xvii i Expected value; time average 18
i Absolute value 18
∇ Gradient in direction of increasing r 18 r i* Conjugate 18 Re[i] Real part 19 e Unit vector in direction of increasing r 19 r Im[i] Imaginary part 21 S()ω Fourier spectrum 21 Ψ()ω Fourier phase 21 Π()ω Power; Power spectrum 21 θ,ϕ Angles 24 Γ(,,θ ϕω ) Directivity function 24 J (,,θ ϕω ) Radiation pattern function 25
Zˆ(,ω Ω ) Wall impedance at frequency ω for sound waves incident from direction Ω 27 R(,ω Ω ) Wall resistance at frequency ω for sound waves incident from direction Ω 27 X (,ω Ω ) Wall reactance at frequency ω for sound waves incident from direction Ω 27 ˆ ˆ ζω(,Ω ) Specific wall impedance Z(,ω Ω )ρ0c 28
Rˆ(,ω Ω ) Complex plane wave pressure reflection factor 28
θir,θ Incident direction, reflected direction (polar representation) 28 α(,ω Ω ) Absorption coefficient, at frequency ω , for waves incident from direction Ω 29 erfc Complementary error function 33 σ Effective flow resistivity 33 ˆ Ψ(,θiiϕθϕ ,,) Scattering factor 35
βˆ Specific admittance 35
22 Wξ ξ Correlation length of random fluctuation 36
Iir, I Absolute value of the incident intensity, of the reflected intensity 37
xviii i Norm 37
h Relative humidity 42 I Absolute value of the intensity 42
mT(,ω Celsius ,) h Air absorption coefficient (abbreviated: m()ω ) 42
att(,ω TCelsius ,) h Air attenuation (in dB Km ) 42
M p Panel mass density 44 b Plate thickness 44 Ξ Flow resistivity 45 d Width of enclosed space 46 ωˆ =+ωδi Complex angular frequency 58 δ Damping constant 59 κγˆ =+ki Eigenvalue (characteristic value) 60
Ls Length (coordinate s ) 60
Nfmodes () Number of eigenmodes up to frequency f 64 dN modes ()f Eigenmodes density at frequency f 65 df
δ NM Kronecker symbol 65 Q Volume velocity strength 66 δ ()rr− Dirac delta function 66 0 G(,rr ,ω ) Green’s function 66 0
fSch Schröder frequency 67
δ Mean of damping factors 67 gt(,rr ,) Impulse response (pressure) 67 0 Ψ(,rr ,)t Signal response (pressure) 68 0
T60 Reverberation time 69 D()e Directional energy density 72
T15 Reverberation time (measured from -5dB to -20dB) 73
xix T20 Reverberation time (measured from -5dB to -25dB) 73
T30 Reverberation time (measured from -5dB to -35dB) 73 Et() Decay curve 74 ht() Energy impulse response 74 EDT Early decay time 74
D50 Definition 74
C50 Speech clarity index 75
C80 Music clarity index 75 G Relative sound level 75
tI Initial-time-delay gap 75
tS Centre-of-gravity time 76 BR Bass ratio 76 LEF Lateral energy fraction 76 IACC Interaural cross correlation coefficient 76 STI Speech transmission index 77 ST Stage parameters 77 A()ω Absorption at frequency ω 91 α()ω Mean absorption coefficient at frequency ω 93
αi ()ω Area-averaged random-incidence energy absorption coefficient at frequency ω 93 M Number of walls 97 K Reflection order 97 ∆()ω Diffusivity coefficient 98 rr; Source’s position vector; Receiver’s position vector 98 SR NR Number of rays 100 B “Irradiation strength” (after [Kuttruff, 1979]) 106 (,Γ P ,)ρ Measure space 114 i Measure 114 M Two-dimensional manifold; Boundary of enclosure 115 σ Surface area measure (Lebesgue) 115
xx V Volume measure (Lebesgue) 115 Ω = (,θ ϕ ) Unit length vector in R3 , representing directions 116 S2 Unit sphere 116 χ Subset of the unit sphere 116