A Fast Multi-Frequency Iterative Acoustic Boundary Element Method Solver Based on a Preconditioned Accelerated Recycling Strategy

Abstract

Common multi-frequency solution strategies for acoustic systems are either associated with high computational costs or resort to model reduction strategies, which potentially induce considerable approximation errors. In that context, this work proposes a fast multi-frequency solver for acoustic Boundary Element Method (BEM) analyses that combines conventional preconditioners with an accelerated recycling strategy, achieving in that way a fast and accurate iterative solution. Specifically, by approximating both the system and the related preconditioner by individual affine expressions, an upfront Galerkin projection on a global deflation basis is facilitated, therefore enabling the deployment of the accelerated recycling scheme. The latter acts by accelerating the construction of the related deflation projectors, therefore leading to a significantly faster recycling procedure. Finally, the efficiency of the global deflation basis is guaranteed through an Automatic Krylov subspaces Recycling for Deflation (AKR-D) [D. Panagiotopoulos, W. Desmet and E. Deckers, An accelerated subspaces recycling strategy for the deflation of parametric linear systems based on model order reduction, Comput. Methods Appl. Mech. Eng. 403 (2023) 115765] algorithm and thus, the cumulative number of iterations required by an iterative solver within a frequency sweep is drastically decreased. The proposed multi-frequency solver is tested on two industrially relevant examples of an interior and an exterior scattering problem, benchmarked against the traditional solution strategies.