A practical examination of the errors arising in the direct collocation boundary element method for acoustic scattering

Abstract

a b s t r a c tFor many engineers and acousticians, the boundary element method (BEM) provides an invaluable tool in the analysis of complex problems. It is particularly well suited for the examination of acoustical problems within large domains. Unsurprisingly, the widespread application of the BEM continues to produce an academic interest in the methodology. New algorithms and techniques are still being proposed, to extend the functionality of the BEM, and to compute the required numerical tasks with greater accuracy and efficiency. However, for a given global error constraint, the actual computational accuracy that is required from the various numerical procedures is not often discussed. Within this context, this paper presents an investigation into the discretisation and computational errors that arise in the BEM for acoustic scattering. First, accurate routines to compute regular, weakly singular, and nearly weakly singular integral kernels are examined. These are then used to illustrate the effect of the requisite boundary discretisation on the global error. The effects of geometric and impedance singularities are also considered. Subsequently, the actual integration accuracy required to maintain a given global error constraint is established. Several regular and irregular scattering examples are investigated, and empirical parameter guidelines are provided.

k-Wave: MATLAB toolbox for the simulation and reconstruction of photoacoustic wave fields

“…There are two significant advantages to using Fourier series: first, the amplitudes of the Fourier components can be calculated efficiently using the FFT, and second, the basis functions are sinusoidal, so only two nodes per wavelength are required, rather than the six to ten required in other methods. 20,21 While the PS method improves efficiency in the spatial domain, conventional FD schemes are still necessary to calculate the gradients ͑rates of change͒ in the time domain. The FD approximation introduces instability into the numerical simulation that can only be controlled by limiting the size of the time-step.…”

Novel basis functions for the partition of unity boundary element method for Helmholtz problems

“…The integrals within (11) can be evaluated using quadrature, taking care to use an appropriate scheme (e.g. [23,24]) when considering the singular integrals that arise.…”

Survey on predicting and controlling of structure-borne noise from rail transit bridges

“…20 60 Statistical Energy Analysis [32] Janssens & Thompson [8] TWINS NORBERT Bewes Thompson [33] [34] C. [36] 1 2 3 4 [37] ANSYS [24] ANSYS SYSNOISE D. [25,38,39] [40] [41] Helmholtz [39,[42][43][44][45] [ 39,[46][47] Helmholtz Helmholtz Galerkin [39] [42] [26,45] [ 44,48] [ [42][43][44][45]48] [49] [50] Ansys Sysnoise…”