A 3D-ACA accelerated time domain boundary element method for elastodynamics using FMM and $\mathcal{H}$-matrix techniques

Abstract

The time-domain Boundary Element Method (BEM) for linear elastodynamics with vanishing initial conditions is considered. Spatial discretization uses standard low-order boundary elements, while temporal discretization employs the generalized Convolution Quadrature (gCQ) method. The gCQ framework requires evaluating BEM matrices in the Laplace domain at several complex frequencies along a chosen contour, producing a three-dimensional tensor with one spatial matrix slice per frequency. To reduce storage and computational cost, a low rank approximation of the tensor is computed using 3D-Adaptive Cross Approximation (3D-ACA), extending the classical ACA to handle both the additional frequency dimension and the tensorial structure of elastodynamics. Within each frequency slice, the BEM matrices are further compressed using either the classical ACA algorithm using the $\mathcal{H}$-matrix approach or a Chebyshev interpolation based Fast Multipole Method (FMM). A comparative study of all proposed methods is carried out using two academic examples, and the structural vibration of an induction machine is analyzed.