FDTD/K-DWM simulation of 3D room acoustics on general purpose graphics hardware using compute unified device architecture (CUDA)

The growing demand for reliable prediction of sound fields in rooms have resulted in adaptation of various approaches for physical modeling, including the Finite Difference Time Domain (FDTD) and the Digital Waveguide Mesh (DWM). Whilst considered versatile and attractive methods, they suffer from dispersion errors that increase with frequency and vary with direction of propagation, thus imposing a high frequency calculation limit. Attempts have been made to reduce such errors by considering different mesh topologies, by spatial interpolation, or by simply oversampling the grid. As the latter approach is computationally expensive, its application to three-dimensional problems has often been avoided. In this paper, we propose an implementation of the FDTD on general purpose graphics hardware, allowing for high sampling rates whilst maintaining reasonable calculation times. Dispersion errors are consequently reduced and the high frequency limit is increased. A range of graphics processors are evaluated and compared with traditional CPUs in terms of accuracy, calculation time and memory requirements

Low-frequency sound source localization as a function of closed acoustic spaces

Further development of an emerging generalized theory of low-frequency sound localization in closed listening spaces is presented that aims to resolve the ambiguities inherent in previous research. The approach takes a robust set of equations based on source/listener location, reverberation time and room dimensions and tests them against a set of evaluation procedures to explore image location against theoretical expectations. Phantom imaging is germane to the methodology and its match within the theoretical framework is investigated. Binaural recordings are used to inspect a range of closed environments for localization clues each with a range of source-listener placements. A complementary series of small-scale listening tests are included for perceptual validation

Acoustic modeling using the digital waveguide mesh

The digital waveguide mesh has been an active area of music acoustics research for over ten years. Although founded in 1-D digital waveguide modeling, the principles on which it is based are not new to researchers grounded in numerical simulation, FDTD methods, electromagnetic simulation, etc. This article has attempted to provide a considerable review of how the DWM has been applied to acoustic modeling and sound synthesis problems, including new 2-D object synthesis and an overview of recent research activities in articulatory vocal tract modeling, RIR synthesis, and reverberation simulation. The extensive, although not by any means exhaustive, list of references indicates that though the DWM may have parallels in other disciplines, it still offers something new in the field of acoustic simulation and sound synth

Contributions to discrete-time methods for room acoustic simulation

The sound field distribution in a room is the consequence of the acoustic properties of radiating sources and the position, geometry and absorbing characteristics of the surrounding boundaries in an enclosure (boundary conditions). Despite there existing a consolidated acoustic wave theory, it is very difficult, nearly impossible, to find an analytical expression of the sound variables distribution in a real room, as a function of time and position. This scenario represents as an inhomogeneous boundary value problem, where the complexity of source properties and boundary conditions make that problem extremely hard to solve. Room acoustic simulation, as treated in this thesis, comprises the algebraical approach to solve the wave equation, and the way to define the boundary conditions and source modeling of the scenario under analysis. Numerical methods provide accurate algorithms for this purpose and among the different possibilities, the use of discrete-time methods arises as a suitable solution for solving those partial differential equations, particularized by some specific constrains. Together with the constant growth of computer power, those methods are increasing their suitability for room acoustic simulation. However, there exists an important lack of accuracy in the definition of some of these conditions so far: current frequency-dependent boundary conditions do not comply with any physical model, and directive sources in discrete-time methods have been hardly treated. This thesis discusses about the current state-of-the-art of the boundary conditions and source modeling in discrete-time methods for room acoustic simulation, and it contributes some algorithms to enhance boundary condition formulation, in a locally reacting impedance sense, and source modelling in terms of directive sources under a defined radiation pattern. These algorithms have been particularized to some discrete-time methods such as the Finite Difference Time Domain and the Digital Waveguide Mesh.Escolano Carrasco, J. (2008). Contributions to discrete-time methods for room acoustic simulation [Tesis doctoral no publicada]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/8309Palanci

Source excitation strategies for obtaining impulse responses in finite difference time domain room acoustics simulation

This paper considers source excitation strategies in finite difference time domain room acoustics simulations for auralization purposes. We demonstrate that FDTD simulations can be conducted to obtain impulse responses based on unit impulse excitation, this being the shortest, simplest and most efficiently implemented signal that might be applied. Single, rather than double, precision accuracy simulations might be implemented where memory use is critical but the consequence is a remarkably increased noise floor. Hard source excitation introduces a discontinuity in the simulated acoustic field resulting in a shift of resonant modes from expected values. Additive sources do not introduce such discontinuities, but instead result in a broadband offset across the frequency spectrum. Transparent sources address both of these issues and with unit impulse excitation the calculation of the compensation filters required to implement transparency is also simplified. However, both transparent and additive source excitation demonstrate solution growth problems for a bounded space. Any of these approaches might be used if the consequences are understood and compensated for, however, for room acoustics simulation the hard source is the least favourable due to the fundamental changes it imparts on the underlying geometry. These methods are further tested through the implementation of a directional sound source based on multiple omnidirectional point sources