In previous work, a set of time-domain equations for sound propagation in a rigid porous medium, including viscous and thermal dissipation effects, was derived [Ostashev et al., J. Acoust. Soc. Am. 115, 2624 (2004)]. From those equations, time-domain counterparts to frequency-domain impedance boundary conditions (BCs) were also derived. In this paper, several approaches to computational implementation of time-domain acoustic interactions with rigid porous surfaces are discussed that are based on these equations. Most of the approaches involve convolutions between relaxation or transformed impedance functions and the acoustic wavefield variables. For these approaches, when the relaxation times for viscous and thermal diffusion in the pores (which are inversely proportional to the static flow resistivity) are large, the convolution integrals must be evaluated over many periods of the acoustic disturbance and therefore become very computationally demanding. Some example finite-difference time-domain (FDTD) calculations show the significantly increased attenuation within the porous medium resulting from the convolution terms. Computational instabilities are observed when the wave is strongly attenuated over the FDTD spatial grid interval. Alternative approaches, based on Pade and series approximations to the time-domain impedance BC, are also considered that provide more compact numerical solutions.
Read full abstract