@katharinafuchs

PSD-Greedy Basis Generation for Structure-Preserving Model Order Reduction of Hamiltonian Systems

, , and . Proceedings of the Conference Algoritmy 2020, page 151--160. Vydavateľstvo SPEKTRUM, (August 2020)

Abstract

Hamiltonian systems are central in the formulation of non-dissipative physical systems. They are characterized by a phase-space, a symplectic form and a Hamiltonian function. In numerical simulations of Hamiltonian systems, algorithms show improved accuracy when the symplectic structure is preserved 10. For structure-preserving model order reduction (MOR) of Hamiltonian systems, symplectic MOR 17, 14, 11, 3, 16 can be used. It is based on a reduced-order basis that is symplectic, which requires symplectic basis generation techniques. In our work, we discuss greedy algorithms for symplectic basis generation. We complement the procedure presented in 14 with ideas of the POD-greedy 9, which results in a new greedy symplectic basis generation technique, the PSD-greedy. Inspired by POD-greedy, we use compression techniques in the greedy iterations to enrich the basis iteratively. We prove that this algorithm computes a symplectic basis when symplectic techniques are used for compression. In the numerical experiments, we compare the discussed methods for a linear elasticity problem. The results show that improvements of up to one order of magnitude in the relative reduction error are achievable with the new basis generation technique compared to the existing greedy approach from 14.

Links and resources

Tags

community