Exact operator inference with minimal data

TL;DR

This paper presents a new snapshot data generation method for operator inference that guarantees exact reconstruction of reduced-order models with minimal data, improving stability and preserving properties like symmetry.

Contribution

The authors introduce a guaranteed, minimal snapshot data generation technique for operator inference that ensures exact ROM reconstruction without heuristics or re-projection.

Findings

Exact ROM reconstruction achieved with minimal snapshots.

Method preserves properties like symmetry in operators.

Numerical tests confirm stability and accuracy.

Abstract

This work introduces a novel method to generate snapshot data for operator inference that guarantees the exact reconstruction of intrusive projection-based reduced-order models (ROMs). To ensure exact reconstruction, the operator inference least squares matrix must have full rank, without regularization. Existing works have achieved this full rank using heuristic strategies to generate snapshot data and a-posteriori checks on full rank, but without a guarantee of success. Our novel snapshot data generation method provides this guarantee thanks to two key ingredients: first we identify ROM states that induce full rank, then we generate snapshots corresponding to exactly these states by simulating multiple trajectories for only a single time step. This way, the number of required snapshots is minimal and orders of magnitude lower than typically reported with existing methods. The method avoids non-Markovian terms and does not require re-projection. Since the number of snapshots is minimal, the least squares problem simplifies to a linear system that is numerically more stable. In addition, because the inferred operators are exact, properties of the intrusive ROM operators such as symmetry or skew-symmetry are preserved. Numerical results for differential equations involving 2nd, 3rd and 8th order polynomials demonstrate that the novel snapshot data generation method leads to exact reconstruction of the intrusive reduced order models.