Premium
An adaptive sparse grid rational Arnoldi method for uncertainty quantification of dynamical systems in the frequency domain
Author(s) -
Römer Ulrich,
Bollhöfer Matthias,
Sreekumar Harikrishnan,
Blech Christopher,
Christine Langer Sabine
Publication year - 2021
Publication title -
international journal for numerical methods in engineering
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.421
H-Index - 168
eISSN - 1097-0207
pISSN - 0029-5981
DOI - 10.1002/nme.6761
Subject(s) - sparse grid , discretization , monte carlo method , frequency domain , computer science , algorithm , interpolation (computer graphics) , grid , sampling (signal processing) , uncertainty quantification , mathematics , mathematical optimization , statistics , filter (signal processing) , artificial intelligence , motion (physics) , mathematical analysis , geometry , computer vision
In this paper, we address discrete linear systems in the frequency domain, where both frequency and random parameters are considered. Sampling such a system many times is computationally challenging if the size is too large, which is often the case for a finite element discretization of partial differential equations. We propose an adaptive method, steered by dual error indicators, which combines rational Arnoldi model order reduction and sparse grid interpolation with hierarchical Leja nodes. At each Leja node, a reduced order model (ROM) is constructed such that the ROM‐error in the frequency domain is balanced with the sparse grid error in the random parameter domain. Both the ROM basis and the sparse grid set are enlarged in a hierarchical manner to achieve a prescribed accuracy in statistical moments. Moreover, parameter sensitivities over the frequency range can be easily extracted from the combined reduced order‐surrogate model. In the numerical tests considered in the paper, the method employs sampling sets, which are reduced by at least an order of magnitude, compared to Monte Carlo simulation. Additionally, for an example from vibroacoustics, building a ROM reduces the system size by roughly a factor of 100, while still providing an acceptable accuracy.