A fast multipole boundary element method for acoustics in viscothermal fluids

Standard numerical models in acoustics rely on the isentropic Helmholtz equation. Its derivation assumes adiabatic and reversible, i.e., dissipation-free, wave propagation. Sound waves in fluids are, however, subject to viscous and thermal losses. These losses originate from viscous friction and heat conduction, leading to the formation of acoustic boundary layers. Considering these effects becomes significant in setups with acoustic cavities of similar dimension as the boundary layers. Recently, boundary element methods (BEM) accounting for the viscothermal dissipation have been proposed. These methods limit the discretization to the surface of the fluid domain and require significantly fewer degrees of freedom than comparable finite element models. However, the BEM coefficient matrices are fully populated, resulting in high computational costs and storage requirements. This study develops a new BEM formulation for acoustics in viscothermal fluids that uses the fast multipole method - a low-rank approximation technique based on a hierarchical subdivision of the computational domain - to alleviate this shortcoming. It is shown that the fast viscothermal BEM improves the algorithmic complexity over the conventional formulation, reducing both the memory requirements and solution time. The results indicate good scalability of the method making it feasible for larger applications such as acoustic metamaterials.