One of the main challenges in time domain wave-based acoustics is the accurate simulation of both boundary conditions and barriers capable of reflecting and transmitting energy. Such scattering structures are generally of irregular geometry and characterised in terms of frequency-dependent reflectances and transittances. Conditions for numerical stability can be difficult to obtain in either case. Immersed boundary methods, which are heavily used in computational fluid dynamics applications, replace boundaries by discrete driving terms, avoiding volumetric meshing and staircasing approaches altogether. The main contribution of this article is a unified numerical treatment of both impedance boundary conditions and barriers capable of transmitting energy and suitable for use in the setting of wave-based acoustics. It is framed in terms of the immersed boundary method within a finite difference time domain scheme, using a dual set of matched discrete driving terms in both the conservation of mass and momentum equations that can be tuned against a desired reflectance or transmittance. Numerical results in three dimensions are presented, illustrating non-porous barriers and impedance boundary conditions, and highlight important features such as spurious leakage through an immersed boundary. A brief demonstration of conditions for numerical stability of the immersed boundary method in this context is provided in an appendix.
Read full abstract