The boundary element method (BEM) is an often used tool for numerically solving acoustic radiation and reflection problems. Most of the time, a formulation in the frequency domain can be used, however, for short impulses  or when the acoustic simulation is coupled with a non-linear behaviour caused by structure deformation, a formulation in the time domain is necessary.


The boundary integral equations and the fundamental solution necessary for the BEM in the time domain are derived by inverse Fourier transformation of the corresponding formulations in the frequency domain. These equations are then discretized using the Galerkin method in the spatial dimensions and the collocation method in the time dimension. The MOT (Marching-On-in-Time) method is used to solve the resulting system of equations. The well known stability problem of the MOT-method is handled by using the Burton-Miller approach in combination with the Galerkin method in the spatial discretization and high order temporal interpolations. It is well known that these measures enhance the stability of MOT.

Additionally it is planned to enhance the efficiency of the method by using a modified plane wave time decomposition (PWTD) algorithm.