Stable model order reduction for time-domain exterior vibro-acoustic finite element simulations

Abstract This paper presents a novel method that enables model order reduction of a fully-coupled, exterior vibro-acoustic finite element model for time domain simulations. The method preserves the stability of the full model and reduces the amount of degrees of freedom significantly, with only a moderate amount of calculation complexity. Infinite elements are used on the finite element boundary to satisfy the Sommerfeld radiation condition. Two different strategies to calculate the reduced order model are compared. The first strategy works with a split reduced basis and can be applied on any fully stable model. The second strategy starts from a modified Everstine formulation and directly builds a reduced basis from the full model, leading to more compact reduced order models. Furthermore, a method is derived to perform explicit time integration on the reduced system, while avoiding the inversion of the mass matrix, which might not be possible due to the presence of the infinite elements. Also this method is shown to preserve the stability of the model and a computationally efficient way for implementation of the method is discussed. The effectiveness of the novel methodology is demonstrated with two numerical models.

[1]  Ren-Cang Li,et al.  Structure-Preserving Model Reduction Using a Krylov Subspace Projection Formulation , 2005 .

[2]  Jean-Pierre Berenger,et al.  A perfectly matched layer for the absorption of electromagnetic waves , 1994 .

[3]  Barbara Kaltenbacher,et al.  A modified and stable version of a perfectly matched layer technique for the 3-d second order wave equation in time domain with an application to aeroacoustics , 2013, J. Comput. Phys..

[4]  R. J. Astley,et al.  Three-dimensional wave-envelope elements of variable order for acoustic radiation and scattering , 1998 .

[5]  B. Moore.  Principal component analysis in linear systems: Controllability, observability, and model reduction , 1981 .

[6]  Zhaojun Bai,et al.  SOAR: A Second-order Arnoldi Method for the Solution of the Quadratic Eigenvalue Problem , 2005, SIAM J. Matrix Anal. Appl..

[7]  J. Keller,et al.  Exact non-reflecting boundary conditions , 1989 .

[8]  C. Wallace.  Radiation Resistance of a Rectangular Panel , 1972 .

[9]  B. Lohmann,et al.  Structure Preserving Order Reduction of Large Scale Second Order Systems , 2004 .

[10]  G. Duan.  Analysis and Design of Descriptor Linear Systems , 2010 .

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

[12]  A. Majda,et al.  Absorbing boundary conditions for the numerical simulation of waves , 1977 .

[13]  David S. Burnett,et al.  A three‐dimensional acoustic infinite element based on a prolate spheroidal multipole expansion , 1994 .

[14]  J. A. Hamilton,et al.  The stability of infinite element schemes for transient wave problems , 2006 .

[15]  Zhaojun Bai,et al.  Dimension Reduction of Large-Scale Second-Order Dynamical Systems via a Second-Order Arnoldi Method , 2005, SIAM J. Sci. Comput..

[16]  G. C. Everstine.  A symmetric potential formulation for fluid-structure interaction , 1981 .

[17]  Serkan Gugercin,et al.  Model Reduction of Descriptor Systems by Interpolatory Projection Methods , 2013, SIAM J. Sci. Comput..

[18]  Roland W. Freund,et al.  Structure-Preserving Model Order Reduction of RCL Circuit Equations , 2008 .

[19]  Awad H. Al-Mohy,et al.  A New Scaling and Squaring Algorithm for the Matrix Exponential , 2009, SIAM J. Matrix Anal. Appl..

[20]  Daniel Dreyer,et al.  Efficient Infinite Elements based on Jacobi Polynomials , 2008 .

[21]  Wim Desmet,et al.  Stability‐preserving model order reduction for time‐domain simulation of vibro‐acoustic FE models , 2017 .

[22]  R. J. Astley,et al.  Finite element formulations for acoustical radiation , 1983 .

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