The surrogate matrix methodology: Accelerating isogeometric analysis of waves

Daniel Drzisga, Brendan Keith, Barbara Wohlmuth

Research output: Contribution to journalArticlepeer-review

6 Scopus citations

Abstract

The surrogate matrix methodology delivers low-cost approximations of matrices (i.e., surrogate matrices) which are normally computed in Galerkin methods via element-scale quadrature formulas. In this paper, the methodology is applied to a number of model problems in wave mechanics treated in the Galerkin isogeometric setting. Herein, the resulting surrogate methods are shown to significantly reduce the assembly time in high frequency wave propagation problems. In particular, the assembly time is reduced with negligible loss in solution accuracy. This paper also extends the scope of previous articles in its series by considering multi-patch discretizations of time-harmonic, transient, and nonlinear PDEs as particular use cases of the methodology. Our a priori error analysis for the Helmholtz equation demonstrates that the additional consistency error introduced by the presence of surrogate matrices is independent of the wave number. In addition, our floating point analysis establishes that the computational complexity of the methodology compares favorably to other contemporary fast assembly techniques for isogeometric methods. Our numerical experiments demonstrate clear performance gains for time-harmonic problems, both with and without the presence of perfectly matched layers. Notable speed-ups are also presented for a transient problem with a compressible neo-Hookean material.

Original languageEnglish
Article number113322
JournalComputer Methods in Applied Mechanics and Engineering
Volume372
DOIs
StatePublished - 1 Dec 2020

Keywords

  • Helmholtz equation
  • Hyperelasticity
  • Isogeometric analysis
  • Linear elasticity
  • Matrix assembly
  • Surrogate numerical methods

Fingerprint

Dive into the research topics of 'The surrogate matrix methodology: Accelerating isogeometric analysis of waves'. Together they form a unique fingerprint.

Cite this