Well‐scaled, a‐posteriori error estimation for model order reduction of large second‐order mechanical systems

• Published in: Zeitschrift für angewandte Mathematik und Mechanik 2020-08, Vol.100 (8), p.n/a
• Main Authors: Grunert, Dennis, Fehr, Jörg, Haasdonk, Bernard
• Format: Article
• Language: English
• Subjects: Algorithms; a‐posteriori error estimation; Computer simulation; Errors; Galerkin method; mechanical system; Mechanical systems; model order reduction; Model reduction; Polynomials; Power series; Series expansion; spectral theorem; Upper bounds
• ISSN: 0044-2267; 1521-4001
• DOI: 10.1002/zamm.201900186

Model Order Reduction is used to vastly speed up simulations but it also introduces an error to the simulation results, which needs to be controlled. The a‐posteriori error estimator of Ruiner et al. for second‐order systems, which is based on the residual, has the advantage of having provable upper bounds and being usable independently of the reduction method. Nevertheless a bottleneck is found in the offline phase, making it unusable for larger models. We use the spectral theorem, power series expansions, monotonicity properties, and self‐tailored algorithms to largely speed up the offline phase by one polynomial order both in terms of computation time as well as storage complexity. All properties are proven rigorously. This eliminates the aforementioned bottleneck. Hence, the error estimator of Ruiner et al. can finally be used for large, linear, second‐order mechanical systems reduced by any model reduction method based on Petrov–Galerkin reduction. The examples show speedups of up to 28.000 and the ability to compute much larger systems with a fixed amount of memory. Model Order Reduction is used to vastly speed up simulations but it also introduces an error to the simulation results, which needs to be controlled. The a‐posteriori error estimator of Ruiner et al. for second‐order systems, which is based on the residual, has the advantage of having provable upper bounds and being usable independently of the reduction method. Nevertheless a bottleneck is found in the offline phase, making it unusable for larger models….