Implicit finite difference schemes for the 3-D wave equation using a 27-point stencil on the cubic grid are presented, for use in room acoustics modelling and artificial reverberation. The system of equations that arises from the implicit formulation is solved using the Jacobi iterative method. Numerical dispersion is analysed and computational efficiency is compared to second-order accurate 27-point explicit schemes. Timing results from GPU implementations demonstrate that the proposed algorithms scale over their explicit counterparts as expected: by a factor of M + 2, where M is a fixed number of Jacobi iterations (eight can be sufficient in single precision). Thus, the accuracy of the approximation can be improved over explicit counterparts with only a linear increase in computational costs, rather than the quartic (in operations) and cubic (in memory) increases incurred when oversampling the grid. These implicit schemes are advantageous in situations where less than 1% dispersion error is desired.

Download

Collisions play a major role in various models of musical instruments; one particularly interesting case is that of the guitar fretboard, the subject of this paper. Here, the string is modelled including effects of tension modulation, and the distributed collision both with the fretboard and individual frets, and including both effects of free string vibration, and under finger-stopped conditions, requiring an additional collision model. In order to handle multiple distributed nonlinearities simultaneously, a finite difference time domain method is developed, with a penalty potential allowing for a convenient model of collision within a Hamiltonian framework, allowing for the construction of stable energy-conserving methods. Implementation details are discussed, and simulation results are presented illustrating a variety of features of such a model.

Download

In this paper, a physics-based model for a snare drum will be discussed, along with its finite difference simulation. The interactions between a mallet and the membrane and between the snares and the membrane will be described as perfectly elastic collisions. A novel numerical scheme for the implementation of collisions will be presented, which allows a complete energy analysis for the whole system. Viscothermal losses will be added to the equation for the 3D wave propagation. Results from simulations and sound examples will be presented.

Download