Stabilized Time-Discontinuous Galerkin Methods With Applications to Structural Acoustics

The time-discontinuous Galerkin (TDG) method possesses high-order accuracy and desirable C-and L-stability for second-order hyperbolic systems including structural acoustics. C- and L-stability provide asymptotic annihilation of high frequency response due to spurious resolution of small scales. These non-physical responses are due to limitations in spatial discretization level for large-complex systems. In order to retain the high-order accuracy of the parent TDG method for high temporal approximation orders within an efficient multi-pass iterative solution algorithm which maintains stability, generalized gradients of residuals of the equations of motion expressed in state-space form are added to the TDG variational formulation. The resultant algorithm is shown to belong to a family of Pade approximations for the exponential solution to the spatially discrete hyperbolic equation system. The final form of the algorithm uses only a few iteration passes to reach the order of accuracy of the parent solution. Analysis of the multi-pass algorithm shows that the first iteration pass belongs to the family of (p+1)-stage stiff accurate Singly-Diagonal-Implicit-Runge-Kutta (SDIRK) method. The methods developed can be viewed as a generalization to the SDIRK method, retaining the desirable features of efficiency and stability, now extended to high-order accuracy. An example of a transient solution to the scalar wave equation demonstrates the efficiency and accuracy of the multi-pass algorithms over standard second-order accurate single-step/single-solve (SS/SS) methods.