Time domain room acoustic simulations using the spectral element method.

This paper presents a wave-based numerical scheme based on a spectral element method, coupled with an implicit-explicit Runge-Kutta time stepping method, for simulating room acoustics in the time domain. The scheme has certain features which make it highly attractive for room acoustic simulations, namely (a) its low dispersion and dissipation properties due to a high-order spatio-temporal discretization; (b) a high degree of geometric flexibility, where adaptive, unstructured meshes with curvilinear mesh elements are supported; and (c) its suitability for parallel implementation on modern many-core computer hardware. A method for modelling locally reacting, frequency dependent impedance boundary conditions within the scheme is developed, in which the boundary impedance is mapped to a multipole rational function and formulated in differential form. Various numerical experiments are presented, which reveal the accuracy and cost-efficiency of the proposed numerical scheme.

[1]  M. Vorländer.  Computer simulations in room acoustics: concepts and uncertainties. , 2013, The Journal of the Acoustical Society of America.

[2]  Chi-Wang Shu,et al.  High order finite difference and finite volume WENO schemes and discontinuous Galerkin methods for CFD , 2001 .

[3]  Toru Otsuru,et al.  A finite-element method using dispersion reduced spline elements for room acoustics simulation , 2014 .

[4]  Nikunj Raghuvanshi,et al.  Numerical wave simulation for interactive audio-visual applications , 2016 .

[5]  Heewon Lee,et al.  An efficient algorithm for the image model technique , 1988 .

[6]  Vesa Välimäki,et al.  Fifty Years of Artificial Reverberation , 2012, IEEE Transactions on Audio, Speech, and Language Processing.

[7]  Stéphanie Salmon,et al.  Arbitrary High-Order Finite Element Schemes and High-Order Mass Lumping , 2007, Int. J. Appl. Math. Comput. Sci..

[8]  Shinichi Sakamoto,et al.  Phase-error analysis of high-order finite difference time domain scheme and its influence on calculation results of impulse response in closed sound field , 2007 .

[9]  S. P. Oliveira.  On Multiple Modes of Propagation of High-Order Finite Element Methods for the Acoustic Wave Equation , 2017 .

[10]  Lauri Savioja,et al.  REAL-TIME 3D FINITE-DIFFERENCE TIME-DOMAIN SIMULATION OF LOW- AND MID-FREQUENCY ROOM ACOUSTICS , 2010 .

[11]  Brian Hamilton,et al.  Audibility of dispersion error in room acoustic finite-difference time-domain simulation as a function of simulation distance. , 2016, The Journal of the Acoustical Society of America.

[12]  Gregor Gassner,et al.  A Comparison of the Dispersion and Dissipation Errors of Gauss and Gauss-Lobatto Discontinuous Galerkin Spectral Element Methods , 2011, SIAM J. Sci. Comput..

[13]  Thomas Krijnen,et al.  openPSTD: The open source pseudospectral time-domain method for acoustic propagation , 2016, Comput. Phys. Commun..

[14]  Dinesh Manocha,et al.  WAVE: Interactive Wave-based Sound Propagation for Virtual Environments , 2015, IEEE Transactions on Visualization and Computer Graphics.

[15]  A. Kulowski.  Algorithmic representation of the ray tracing technique , 1985 .

[16]  J. Hesthaven,et al.  Multiscale modelling of sound propagation through the lung parenchyma , 2014 .

[17]  Spencer J. Sherwin,et al.  Dispersion Analysis of the Continuous and Discontinuous Galerkin Formulations , 2000 .

[18]  Marie-Annick Galland,et al.  Broadband liner impedance eduction for multimodal acoustic propagation in the presence of a mean flow , 2016 .

[19]  Isaac Fried,et al.  Finite element mass matrix lumping by numerical integration with no convergence rate loss , 1975 .

[20]  David A. Kopriva,et al.  Implementing Spectral Methods for Partial Differential Equations , 2009 .

[21]  Peter Møller Juhl,et al.  Fundamentals of General Linear Acoustics , 2013 .

[22]  A. Semlyen,et al.  Rational approximation of frequency domain responses by vector fitting , 1999 .

[23]  Michael Vorlnder,et al.  Auralization: Fundamentals of Acoustics, Modelling, Simulation, Algorithms and Acoustic Virtual Reality , 2020 .

[24]  Stefan Bilbao,et al.  Finite Volume Time Domain Room Acoustics Simulation under General Impedance Boundary Conditions , 2016, IEEE/ACM Transactions on Audio, Speech, and Language Processing.

[25]  Moshe Dubiner.  Spectral methods on triangles and other domains , 1991 .

[26]  Mark Ainsworth,et al.  Dispersive and Dissipative Behavior of the Spectral Element Method , 2009, SIAM J. Numer. Anal..

[27]  Ivo Babuška,et al.  The h, p and h-p version of the finite element method: basis theory and applications , 1992 .

[28]  Toru Otsuru,et al.  An explicit time-domain finite element method for room acoustics simulations: Comparison of the performance with implicit methods , 2016 .

[29]  Jing Xia,et al.  Effects of reverberation and noise on speech intelligibility in normal-hearing and aided hearing-impaired listeners. , 2018, The Journal of the Acoustical Society of America.

[30]  Michael Vorländer,et al.  Combined wave and ray based room acoustic simulations of audio systems in car passenger compartments, Part I: Boundary and source data , 2014 .

[31]  Y. Lam.  Issues for computer modelling of room acoustics in non-concert hall settings , 2005 .

[32]  Maarten Hornikx,et al.  Hybrid Fourier pseudospectral/discontinuous Galerkin time-domain method for wave propagation , 2017, J. Comput. Phys..

[33]  Claes Eskilsson,et al.  A stabilised nodal spectral element method for fully nonlinear water waves , 2015, J. Comput. Phys..

[34]  Carlos A. Felippa,et al.  Mass Matrix Templates: General Description and 1D Examples , 2015 .

[35]  M. Y. Hussaini,et al.  An Analysis of the Discontinuous Galerkin Method for Wave Propagation Problems , 1999 .

[36]  Damian T. Murphy,et al.  Room Impulse Response Synthesis and Validation Using a Hybrid Acoustic Model , 2013, IEEE Transactions on Audio, Speech, and Language Processing.

[37]  Cheol-Ho Jeong,et al.  Estimation of surface impedance at oblique incidence based on sparse array processing. , 2017, The Journal of the Acoustical Society of America.

[38]  Per-Olof Persson,et al.  A Simple Mesh Generator in MATLAB , 2004, SIAM Rev..

[39]  Ronan Guivarch,et al.  Parallelization of Spectral Element Methods , 2002, VECPAR.

[40]  Jian Yu,et al.  Effects of artificial viscosity and upwinding on spectral properties of the discontinuous Galerkin method , 2018, Computers & Fluids.

[41]  S. P. Oliveira,et al.  DFT MODAL ANALYSIS OF SPECTRAL ELEMENT METHODS FOR ACOUSTIC WAVE PROPAGATION , 2008 .

[42]  D. Botteldooren.  Finite‐difference time‐domain simulation of low‐frequency room acoustic problems , 1995 .

[43]  Géza Seriani,et al.  A Parallel Spectral Element Method for Acoustic Wave Modeling , 1997 .

[44]  Y. Miki.  Acoustical Properties of porous materials : Modifications of Delany-Bazley models , 1990 .

[45]  Stefan Bilbao,et al.  FDTD Methods for 3-D Room Acoustics Simulation With High-Order Accuracy in Space and Time , 2017, IEEE/ACM Transactions on Audio, Speech, and Language Processing.

[46]  Sven‐Ingvar Thomasson.  Reflection of waves from a point source by an impedance boundary , 1976 .

[47]  Joe LoVetri,et al.  Modeling of the seat dip effect using the finite‐difference time‐domain method , 1996 .

[48]  Damian Murphy,et al.  Explicit Higher-Order FDTD Schemes for 3D Room Acoustic Simulation , 2014, IEEE/ACM Transactions on Audio, Speech, and Language Processing.

[49]  Dinesh Manocha,et al.  Wave-based sound propagation in large open scenes using an equivalent source formulation , 2013, TOGS.

[50]  Chi-Wang Shu.  Essentially non-oscillatory and weighted essentially non-oscillatory schemes for hyperbolic conservation laws , 1998 .

[51]  Trevor J Cox,et al.  A transient boundary element method model of Schroeder diffuser scattering using well mouth impedance. , 2008, The Journal of the Acoustical Society of America.

[52]  Kirill V. Horoshenkov,et al.  Acoustic Absorbers and Diffusers: Theory, Design and Application, third ed., Trevor J. Cox, Peter D’Antonio. CRC Press, Boca Raton (2017). ISBN: 9781498740999 , 2017 .

[53]  S. P. Oliveira,et al.  Optimal blended spectral-element operators for acoustic wave modeling , 2007 .

[54]  Mcj Maarten Hornikx.  Ten questions concerning computational urban acoustics , 2016 .

[55]  Lauri Savioja,et al.  Overview of geometrical room acoustic modeling techniques. , 2015, The Journal of the Acoustical Society of America.

[56]  H. Kreiss,et al.  Comparison of accurate methods for the integration of hyperbolic equations , 1972 .

[57]  M. Carpenter,et al.  Additive Runge-Kutta Schemes for Convection-Diffusion-Reaction Equations , 2003 .

[58]  M. Schroeder,et al.  On Frequency Response Curves in Rooms. Comparison of Experimental, Theoretical, and Monte Carlo Results for the Average Frequency Spacing between Maxima , 1962 .

[59]  A. Krokstad,et al.  Calculating the acoustical room response by the use of a ray tracing technique , 1968 .

[60]  Samuli Laine,et al.  Accelerated beam tracing algorithm , 2009 .

[61]  Stefan Bilbao,et al.  Modeling of Complex Geometries and Boundary Conditions in Finite Difference/Finite Volume Time Domain Room Acoustics Simulation , 2013, IEEE Transactions on Audio, Speech, and Language Processing.

[62]  M. Vercammen,et al.  Sound concentration caused by curved surfaces , 2013 .

[63]  A. Patera.  A spectral element method for fluid dynamics: Laminar flow in a channel expansion , 1984 .

[64]  Yazhou Wang,et al.  Space-time spectral element method solution for the acoustic wave equation and its dispersion analysis , 2017 .