Premium
A subspace iteration eigensolver based on Cauchy integrals for vibroacoustic problems in unbounded domains
Author(s) -
Baydoun Suhaib Koji,
Voigt Matthias,
Goderbauer Benedikt,
Jelich Christopher,
Marburg Steffen
Publication year - 2021
Publication title -
international journal for numerical methods in engineering
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.421
H-Index - 168
eISSN - 1097-0207
pISSN - 0029-5981
DOI - 10.1002/nme.6701
Subject(s) - subspace topology , finite element method , boundary element method , projection (relational algebra) , mathematics , krylov subspace , generalized minimal residual method , eigenvalues and eigenvectors , projector , methods of contour integration , iterative method , computer science , mathematical optimization , algorithm , mathematical analysis , physics , artificial intelligence , quantum mechanics , thermodynamics
Despite the potential and the increasing popularity of the boundary element method (BEM), modal analyses based on BEM are not yet put into engineering practice, mainly due to the lack of efficient solvers for the underlying nonlinear eigenvalue problem (EVP). In this article, we review a subspace iteration method based on FEAST for the solution of vibroacoustic EVPs involving the finite element method (FEM) and BEM. The subspace is obtained by applying a spectral projector and is computed by contour integration, whereas the contour is also used to subsequently solve the projected EVP by rational approximation. The computation of the projection matrices is addressed by two approaches. In the case of heavy fluid loading, we solve the underlying coupled linear systems by an iterative block Krylov method. In the case of light fluid loading, we exploit the fact that the coupled system admits accurate model order reduction solely based on the structural subsystem. Applications to a spherical shell and to a musical bell indicate that only a few contour points are required for an accurate solution without inducing spurious eigenvalues. The results are compared with those of a contour integral method and illustrate the efficiency of the proposed eigensolver.