Time domain boundary element method for room acoustics

This thesis is about improving the suitability of the time domain Boundary Element Method (BEM) for predicting the scattering from surface treatments used to improve the acoustics of rooms. The discretised integral equations are typically solved by marching on in time from initial silence; however, this being iterative has potential for divergence. Such instability and high computational cost have prohibited the time domain BEM from widespread use. The underlying integral equation is known to not possess unique solutions at certain frequencies, physically interpreted as cavity resonances, and these manifest as resonant poles, all excited and potentially divergent due to numerical error. This has been addressed by others using the combined field integral equation; an approach built upon in this thesis. Accuracy and stability may also be compromised by poor discretisation and integration accuracy. The latter is investigated on real-world surfaces, demonstrating that the popular Gaussian integration schemes are not suitable in some circumstances. Instead a contour integration scheme capable of resolving the integrands‟ singular nature is developed. Schroeder diffusers are Room Acoustic treatments which comprise wells separated by thin fins. The algorithm is extended to model such surfaces, applying the combined field integral equation to the body and an open surface model to the fins. It is shown that this improves stability over an all open surface model. A new model for compliant surfaces is developed, comparable to the surface impedance model used in the frequency domain. This is implemented for surfaces with welled and absorbing sections, permitting modelling of a Schroeder diffuser as a box with surface impedances that simulate the delayed reflections caused by the wells. A Binary Amplitude Diffuser - a partially absorbing diffuser - is also modelled. These new models achieve good accuracy but not universal stability and avenues of future research are proposed to address the latter issue.

[1]  O. D. Kellogg Foundations of potential theory , 1934 .

[2]  Eric Michielssen,et al.  Fast transient analysis of acoustic wave scattering from rigid bodies using a two-level plane wave time domain algorithm , 1999 .

[3]  Gérard C. Herman,et al.  Scattering of transient acoustic waves by an inhomogeneous obstacle , 1981 .

[4]  M. Salazar-Palma,et al.  Solving time domain electric field Integral equation without the time variable , 2006, IEEE Transactions on Antennas and Propagation.

[5]  Ergin,et al.  Fast analysis of transient acoustic wave scattering from rigid bodies using the multilevel plane wave time domain algorithm , 2000, The Journal of the Acoustical Society of America.

[6]  S. Gupta,et al.  Delta Function , 1964 .

[7]  M. B. Friedman,et al.  Diffraction of Pulses by Cylindrical Obstacles of Arbitrary Cross Section , 1962 .

[8]  Jian-Ming Jin,et al.  A Fast Fourier Transform Accelerated Marching-on-in-Time Algorithm for Electromagnetic Analysis , 2001 .

[9]  T. W. Wu A direct boundary element method for acoustic radiation and scattering from mixed regular and thin bodies , 1995 .

[10]  Paul H. L. Groenenboom Wave propagation phenomena , 1983 .

[11]  H. Saunders,et al.  Acoustics: An Introduction to Its Physical Principles and Applications , 1984 .

[12]  P. Smith,et al.  Instabilities in Time Marching Methods for Scattering: Cause and Rectification , 1990 .

[13]  Y. F. Wang,et al.  Transient scattering by arbitrary axisymmetric surfaces , 1978 .

[14]  Sadasiva M. Rao,et al.  Transient scattering from dielectric cylinders - E-field, H-field, and combined field solutions , 1992 .

[15]  R. P. Shaw Transient scattering by a circular cylinder , 1975 .

[16]  Joseph B. Keller,et al.  Diffraction and reflection of pulses by wedges and corners , 1951 .

[17]  E. Michielssen,et al.  Analysis of transient electromagnetic scattering from closed surfaces using a combined field integral equation , 2000 .

[18]  T. Cox,et al.  Prediction and evaluation of the scattering from quadratic residue diffusers , 1994 .

[19]  S. Amini,et al.  Multi-level fast multipole solution of the scattering problem , 2003 .

[20]  Peter D'Antonio,et al.  Two Dimensional Binary Amplitude Diffusers , 1999 .

[21]  T. Ha-Duong,et al.  On Retarded Potential Boundary Integral Equations and their Discretisation , 2003 .

[22]  M. Salazar-Palma,et al.  Solution of time domain electric field Integral equation using the Laguerre polynomials , 2004, IEEE Transactions on Antennas and Propagation.

[23]  Arif Ahmet Ergin Plane -Wave Time -Domain Algorithms for Efficient Analysis of Three-Dimensional Transient Wave Phenomena , 2000 .

[24]  E. Michielssen,et al.  Analysis of transient wave scattering from rigid bodies using a Burton–Miller approach , 1999 .

[25]  D. Henwood,et al.  Stability analysis of a collocation method for solving the retarded potential integral equation , 2005 .

[26]  R. Martinez,et al.  The thin‐shape breakdown (TSB) of the Helmholtz integral equation , 1990 .

[27]  M. A. Stalzer,et al.  A prescription for the multilevel Helmholtz FMM , 1998 .

[28]  M. Schroeder Diffuse sound reflection by maximum−length sequences , 1975 .

[29]  S. P. Walker,et al.  Hybridization of curvilinear time-domain integral equation and time-domain optical methods for electromagnetic scattering analysis , 1998 .

[30]  Hongbin Ju,et al.  Impedance and Its Time-Domain Extensions , 2000 .

[31]  Kenneth M. Mitzner Numerical Solution for Transient Scattering from a Hard Surface of Arbitrary Shape , 1966 .

[32]  Chi Hou Chan,et al.  Improved temporal basis function for time domain electric field integral equation method , 1999 .

[33]  Richard Barakat Transient Diffraction of Scalar Waves by a Fixed Sphere , 1960 .

[34]  T. Ha-Duong,et al.  A Galerkin BEM for transient acoustic scattering by an absorbing obstacle , 2003 .

[35]  S. Amini,et al.  Multiwavelet Galerkin boundary element solution of Laplace's equation , 2006 .

[36]  Christopher K. W. Tam,et al.  Time-Domain Impedance Boundary Conditions for Computational Aeroacoustics , 1996 .

[37]  Yasuhito Kawai,et al.  A numerical method for the calculation of transient acoustic scattering from thin rigid plates , 1990 .

[38]  M. Schroeder Binaural dissimilarity and optimum ceilings for concert halls: More lateral sound diffusion , 1979 .

[39]  Richard Paul Shaw Diffraction of Acoustic Pulses by Obstacles of Arbitrary Shape with a Robin Boundary Condition , 1967 .

[40]  B. P. Rynne,et al.  Stability of Time Marching Algorithms for the Electric Field Integral Equation , 1990 .

[41]  Dj Chappell,et al.  A stable boundary element method for modeling transient acoustic radiation , 2006 .

[42]  G. F. Miller,et al.  The application of integral equation methods to the numerical solution of some exterior boundary-value problems , 1971, Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences.

[43]  T. Terai A note on the paper “on calculation of sound fields around three dimensional objects by integral equation methods” , 1980 .

[44]  S. Wandzurat,et al.  Symmetric quadrature rules on a triangle , 2003 .

[45]  Mingyu Lu,et al.  Fast Evaluation of Two-Dimensional Transient Wave Fields , 2000 .

[46]  A. Forestier,et al.  A Galerkin scheme for the time domain integral equation of acoustic scattering from a hard surface , 1989 .

[47]  Yuan Xu,et al.  A new temporal basis function for the time-domain integral equation method , 2001, IEEE Microwave and Wireless Components Letters.

[48]  Peter M. van den Berg,et al.  A least‐square iterative technique for solving time‐domain scattering problems , 1982 .

[49]  C. T Dyka,et al.  Transient fluid-structure interaction in naval applications using the retarded potential method , 1998 .

[50]  Hongbin Ju,et al.  Broadband Time-Domain Impedance Models , 2001 .

[51]  Agostino Monorchio,et al.  A space-time discretization criterion for a stable time-marching solution of the electric field integral equation , 1997 .

[52]  G. Maier,et al.  Symmetric Galerkin Boundary Element Methods , 1998 .

[53]  Yehuda Leviatan,et al.  On the use of spatio-temporal multiresolution analysis in method of moments solutions of transient electromagnetic scattering , 2001 .

[54]  Richard Paul Shaw,et al.  Transient acoustic scattering by a free (pressure release) sphere , 1972 .

[55]  Raymond D. Mindlin,et al.  Response of an Elastic Cylindrical Shell to a Transverse, Step Shock Wave , 1989 .

[56]  Hongbin Ju,et al.  Time-domain Impedance Boundary Conditions for Computational Acoustics and Aeroacoustics , 2004 .

[57]  Simon P. Walker,et al.  ANALYSIS OF THREE-DIMENSIONAL TRANSIENT ACOUSTIC WAVE PROPAGATION USING THE BOUNDARY INTEGRAL EQUATION METHOD , 1996 .

[58]  Bryan P. Rynne,et al.  INSTABILITIES IN TIME MARCHING METHODS FOR SCATTERING PROBLEMS , 1986 .

[59]  M. Abramowitz,et al.  Handbook of Mathematical Functions With Formulas, Graphs and Mathematical Tables (National Bureau of Standards Applied Mathematics Series No. 55) , 1965 .

[60]  Richard Paul Shaw Comments on “Numerical Solution for Transient Scattering from a Hard Surface of Arbitrary Shape—Retarded Potential Technique“ [K. M. Mitzner, J. Acoust. Soc. Am. 42, 391–397 (1967)] , 1968 .

[61]  Trevor J. Cox Predicting the scattering from reflectors and diffusers using two‐dimensional boundary element methods , 1994 .

[62]  H. A. Schenck Improved Integral Formulation for Acoustic Radiation Problems , 1968 .