Application of a Krylov subspace method for an efficient solution of acoustic transfer functions

Abstract

Solving acoustic radiation problems, arising from systems including fluid–structure interaction, is of interest in many engineering applications. Computing frequency response functions over a large frequency range is a concern in such applications. A method which solves the Helmholtz equation for multiple frequencies in one step is the matrix-Padé-via-Lanczos connection for unsymmetric systems, as presented by Wagner et al. [1]. The present work is based on Ref. [1] and presents a method for efficiently computing frequency responses over a frequency range for coupled structural-acoustic problems, where the structure and the acoustic near field are discretized with finite elements and an analytical Dirichlet-to-Neumann map approximates the far field. The method is based on a Krylov-subspace projection technique which derives a matrix-valued Padé approximation for a restricted area in the near field and the pressure field on a spherical boundary. On the spherical boundary, where the finite domain is truncated, the non-local modified Dirichlet-to-Neumann operator is applied as a low-rank update matrix. The present contribution extends this method and incorporates new techniques for a more stable model reduction through the Lanczos algorithm and a novel weighted adaptive windowing technique. Further, structural damping is incorporated, for computing the acoustic radiation of a harmonically excited plate. These computed results are compared with acoustic measurements in an anechoic chamber and verified with computational results obtained with a commercial code that uses the perfectly matched layer method.