EFFICIENT SOUND POWER COMPUTATION OF OPEN STRUCTURES WITH INFINITE/FINITE ELEMENTS AND BY MEANS OF THE PADÉ-VIA-LANCZOS ALGORITHM

The solution of time-harmonic acoustic problems suffers from a frequency dependency which usually requires to solve the systems for many discrete frequencies independently. Among others, the Pade-via-Lanczos approximation provides an efficient solution. Herein, this method is applied to the external acoustic problem by using a finite and infinite element formulation. The state-space system matrix is first reduced to a condensed transfer matrix which is approximated in a power series, i.e. Pade approximation. The coefficients of this power series are reconstructed by using a Krylov basis, i.e. Lanczos method. Based on this method, the authors provide a formulation for the radiated sound power. The necessary integration over the boundary for the sound power is done in the first vectors of the Lanczos method. Thus the scalar sound power can be approximated in a certain frequency band by means of a transfer matrix, which is several orders of magnitude smaller than the total numbers of degrees of freedom of the overall problem. The error of this additional approximation can be estimated. The method is tested in an example of an open cavity including two loadcases, one representing a cavity with an open window, the other one being a radiating obstacle.

[1]  C. Brebbia,et al.  Boundary Elements IX , 1987 .

[2]  Roland W. Freund,et al.  Reduced-Order Modeling of Large Linear Subcircuits via a Block Lanczos Algorithm , 1995, 32nd Design Automation Conference.

[3]  Michel Rochette,et al.  Calculation of vibro-acoustic frequency response functions using a single frequency boundary element solution and a Pade expansion , 1999 .

[4]  Zhaojun Bai,et al.  Error Estimation of the Padé Approximation of Transfer Functions via the Lanczos Process , 1998 .

[5]  Steffen Marburg Normal modes in external acoustics. Part III: Sound power evaluation based on superposition of frequency-independent modes , 2006 .

[6]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[7]  Omar Laghrouche,et al.  Short wave modelling using special finite elements , 2000 .

[8]  M. Malhotra EFFICIENT COMPUTATION OF MULTI-FREQUENCY FAR-FIELD SOLUTIONS OF THE HELMHOLTZ EQUATION USING PADÉ APPROXIMATION , 2000 .

[9]  K. Gerdes The conjugated vs. the unconjugated infinite element method for the Helmholtz equation in exterior domains , 1998 .

[10]  Thomas J. R. Hughes,et al.  Recent developments in finite element methods for structural acoustics , 1996 .

[11]  F. Ihlenburg Finite Element Analysis of Acoustic Scattering , 1998 .

[12]  L. Thompson A review of finite-element methods for time-harmonic acoustics , 2006 .

[13]  R. Freund Model reduction methods based on Krylov subspaces , 2003, Acta Numerica.

[14]  Henrik Wallén,et al.  Broadband multilevel fast multipole algorithm for acoustic scattering problems , 2006 .

[15]  R. Astley Infinite elements for wave problems: a review of current formulations and an assessment of accuracy , 2000 .

[16]  R. J. Astley,et al.  Three-dimensional wave-envelope elements of variable order for acoustic radiation and scattering. Part I. Formulation in the frequency domain , 1998 .

[17]  Sabine Langer,et al.  SCALED BOUNDARY FINITE ELEMENT METHOD FOR ACOUSTICS , 2006 .

[18]  James P. Tuck-Lee,et al.  Adaptive frequency windowing for multifrequency solutions in structural acoustics based on the matrix Padé‐via‐Lanczos algorithm , 2008 .

[19]  Gene H. Golub,et al.  Matrix computations (3rd ed.) , 1996 .

[20]  Stefan Schneider,et al.  Application of Fast Methods for Acoustic Scattering and Radiation Problems , 2003 .

[21]  P. Pinsky,et al.  Application of Padé via Lanczos approximations for efficient multifrequency solution of Helmholtz problems. , 2003, The Journal of the Acoustical Society of America.

[22]  O. von Estorff,et al.  Improved conditioning of infinite elements for exterior acoustics , 2003 .

[23]  G. Alexits Approximation theory , 1983 .

[24]  T. Wu,et al.  A DIRECT MIXED-BODY BOUNDARY ELEMENT METHOD FOR MUFFLERS WITH INTERNAL THIN COMPONENTS COVERED BY LINING AND A PERFORATED PANEL , 2007 .

[25]  Tuomo Rossi,et al.  Fast direct solution of the Helmholtz equation with a perfectly matched layer or an absorbing boundary condition , 2003 .

[26]  Roland W. Freund,et al.  A Lanczos-type method for multiple starting vectors , 2000, Math. Comput..

[27]  A. Bermúdez,et al.  VALIDATION OF ACOUSTIC MODELS FOR TIME-HARMONIC DISSIPATIVE SCATTERING PROBLEMS , 2007 .

[28]  Isaac Harari,et al.  A survey of finite element methods for time-harmonic acoustics , 2006 .

[29]  S. Marburg,et al.  Performance of iterative solvers for acoustic problems. Part I. Solvers and effect of diagonal preconditioning , 2003 .

[30]  R. J. Astley,et al.  Conditioning of infinite element schemes for wave problems , 2001 .

[31]  R. J. Astley,et al.  Mapped Wave Envelope Elements for Acoustical Radiation and Scattering , 1994 .

[32]  Tetsuya Sakuma,et al.  A TECHNIQUE FOR PLANE-SYMMETRIC SOUND FIELD ANALYSIS IN THE FAST MULTIPOLE BOUNDARY ELEMENT METHOD , 2005 .

[33]  L. Gaul,et al.  APPLICATION OF THE FAST MULTIPOLE BEM FOR STRUCTURAL–ACOUSTIC SIMULATIONS , 2005 .

[34]  Karl Meerbergen,et al.  Theory and software for Krylov methods for the computation of the frequency response of large acoustic finite element models , 2002 .

[35]  R. W. Freund,et al.  Passive reduced-order modeling via Krylov-subspace methods , 2000, CACSD. Conference Proceedings. IEEE International Symposium on Computer-Aided Control System Design (Cat. No.00TH8537).

[36]  Steffen Marburg,et al.  Normal modes in external acoustics. Part II: Eigenvalues and eigenvectors in 2D , 2006 .

[37]  S. Marburg Developments in structural-acoustic optimization for passive noise control , 2002 .

[38]  R. Freund Krylov-subspace methods for reduced-order modeling in circuit simulation , 2000 .