The surrogate matrix methodology: Accelerating isogeometric analysis of waves

by   Daniel Drzisga, et al.

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 isogeometic 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.


page 9

page 20

page 23


Fast assembly of Galerkin matrices for 3D solid laminated composites using finite element and isogeometric discretizations

This work presents a novel methodology for speeding up the assembly of s...

The surrogate matrix methodology: A reference implementation for low-cost assembly in isogeometric analysis

A reference implementation of a new method in isogeometric analysis (IGA...

Boundary Element Methods for the Wave Equation based on Hierarchical Matrices and Adaptive Cross Approximation

Time-domain Boundary Element Methods (BEM) have been successfully used i...

A numerical study of the pollution error and DPG adaptivity for long waveguide simulations

High-frequency wave propagation has many important applications in acous...

A hierarchical preconditioner for wave problems in quasilinear complexity

The paper introduces a novel, hierarchical preconditioner based on neste...

Fast formation and assembly of isogeometric Galerkin matrices for trimmed patches

This work explores the application of the fast assembly and formation st...