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

The modeling of diffuse boundaries in the 2-D digital waveguide mesh

The digital waveguide mesh can be used to simulate the propagation of sound waves in an acoustic system. The accurate simulation of the acoustic characteristics of boundaries within such a system is an important part of the model. One significant property of an acoustic boundary is its diffusivity. Previous approaches to simulating diffuse boundaries in a digital waveguide mesh are effective but exhibit limitations and have not been analyzed in detail. An improved technique is presented here that simulates diffusion at boundaries and offers a high degree of control and consistency. This technique works by rotating wavefronts as they pass through a special diffusing layer adjacent to the boundary. The waves are rotated randomly according to a chosen probability function and the model is lossless. This diffusion model is analyzed in detail, and its diffusivity is quantified in the form of frequency dependent diffusion coefficients. The approach used to measuring boundary diffusion is described here in detail for the 2-D digital waveguide mesh and can readily be extended for the 3-D case

Efficient Synthesis of Room Acoustics via Scattering Delay Networks

An acoustic reverberator consisting of a network of delay lines connected via scattering junctions is proposed. All parameters of the reverberator are derived from physical properties of the enclosure it simulates. It allows for simulation of unequal and frequency-dependent wall absorption, as well as directional sources and microphones. The reverberator renders the first-order reflections exactly, while making progressively coarser approximations of higher-order reflections. The rate of energy decay is close to that obtained with the image method (IM) and consistent with the predictions of Sabine and Eyring equations. The time evolution of the normalized echo density, which was previously shown to be correlated with the perceived texture of reverberation, is also close to that of IM. However, its computational complexity is one to two orders of magnitude lower, comparable to the computational complexity of a feedback delay network (FDN), and its memory requirements are negligible

The KW-boundary hybrid digital waveguide mesh for room acoustics applications

The digital waveguide mesh is a discrete-time simulation used to model acoustic wave propagation through a bounded medium. It can be applied to the simulation of the acoustics of rooms through the generation of impulse responses suitable for auralization purposes. However, large-scale three-dimensional mesh structures are required for high quality results. These structures must therefore be efficient and also capable of flexible boundary implementation in terms of both geometrical layout and the possibility for improved mesh termination algorithms. The general one-dimensional N-port boundary termination is investigated, where N depends on the geometry of the modeled domain and the mesh topology used. The equivalence between physical variable Kirchoff-model, and scattering-based wave-model boundary formulations is proved. This leads to the KW-hybrid one-dimensional N-port boundary-node termination, which is shown to be equivalent to the Kirchoff- and wave-model cases. The KW-hybrid boundary-node is implemented as part of a new hybrid two-dimensional triangular digital waveguide mesh. This is shown to offer the possibility for large-scale, computationally efficient mesh structures for more complex shapes. It proves more accurate than a similar rectilinear mesh in terms of geometrical fit, and offers significant savings in processing time and memory use over a standard wave-based model. The new hybrid mesh also has the potential for improved real-world room boundary simulations through the inclusion of additional mixed modeling algorithms

Modelling of Path Arrival Rate for In-Room Radio Channels with Directive Antennas

We analyze the path arrival rate for an inroom radio channel with directive antennas. The impulse response of this channel exhibits a transition from early separate components followed by a diffuse reverberation tail. Under the assumption that the transmitter's (or receiver's) position and orientation are picked uniformly at random we derive an exact expression of the mean arrival rate for a rectangular room predicted by the mirror source theory. The rate is quadratic in delay, inversely proportional to the room volume, and proportional to the product of beam coverage fractions of the transmitter and receiver antennas. Making use of the exact formula, we characterize the onset of the diffuse tail by defining a "mixing time" as the point in time where the arrival rate exceeds one component per transmit pulse duration. We also give an approximation for the power-delay spectrum. It turns out that the power-delay spectrum is unaffected by the antenna directivity. However, Monte Carlo simulations show that antenna directivity does indeed play an important role for the distribution of instantaneous mean delay and rms delay spreadComment: Submitted to IEEE Trans. Antennas and Propagatio

Real-time dynamic articulations in the 2-D waveguide mesh vocal tract model

Time domain articulatory vocal tract modeling in one-dimensional (1-D) is well established. Previous studies into two-dimensional (2-D) simulation of wave propagation in the vocal tract have shown it to present accurate static vowel synthesis. However, little has been done to demonstrate how such a model might accommodate the dynamic tract shape changes necessary in modeling speech. Two methods of applying the area function to the 2-D digital waveguide mesh vocal tract model are presented here. First, a method based on mapping the cross-sectional area onto the number of waveguides across the mesh, termed a widthwise mapping approach is detailed. Discontinuity problems associated with the dynamic manipulation of the model are highlighted. Second, a new method is examined that uses a static-shaped rectangular mesh with the area function translated into an impedance map which is then applied to each waveguide. Two approaches for constructing such a map are demonstrated; one using a linear impedance increase to model a constriction to the tract and another using a raised cosine function. Recommendations are made towards the use of the cosine method as it allows for a wider central propagational channel. It is also shown that this impedance mapping approach allows for stable dynamic shape changes and also permits a reduction in sampling frequency leading to real-time interaction with the model

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