Non-intrusive Balancing Transformation of Highly Stiff Systems with Lightly-damped Impulse Response

Abstract

Balanced truncation (BT) is a model reduction method that utilizes a coordinate transformation to retain eigen-directions that are highly observable and reachable. To address realizability and scalability of BT applied to highly stiff and lightly-damped systems, a non-intrusive data-driven method is developed for balancing discrete-time systems via the eigensystem realization algorithm (ERA). The advantage of ERA for balancing transformation makes full-state outputs tractable. Further, ERA enables balancing despite stiffness, by eliminating computation of balancing modes and adjoint simulations. As a demonstrative example, we create balanced ROMs for a one-dimensional reactive flow with pressure forcing, where the stiffness introduced by the chemical source term is extreme (condition number $10^{13}$), preventing analytical implementation of BT. We investigate the performance of ROMs in prediction of dynamics with unseen forcing inputs and demonstrate stability and accuracy of balanced ROMs in truly predictive scenarios whereas without ERA, POD-Galerkin and Least-squares Petrov-Galerkin projections fail to represent the true dynamics. We show that after the initial transients under unit impulse forcing, the system undergoes lightly-damped oscillations, which magnifies the influence of sampling properties on predictive performance of the balanced ROMs. We propose an output domain decomposition approach and couple it with tangential interpolation to resolve sharp gradients at reduced computational costs.

References in corpus (12)

• Direct numerical simulation of turbulent channel flow up to $R{e}_{\tau }\approx 5200$
• Proper Orthogonal Decomposition Closure Models For Turbulent Flows: A Numerical Comparison
• Reduced order models for control of fluids using the Eigensystem Realization Algorithm
• On closures for reduced order models $-$ A spectrum of first-principle to machine-learned avenues
• Adaptive $h$-refinement for reduced-order models
• Data-driven reduced-order models via regularized operator inference for a single-injector combustion process
• Minimal subspace rotation on the Stiefel manifold for stabilization and enhancement of projection-based reduced order models for the compressible Navier-Stokes equations
• Data-driven Discovery of Closure Models
• A Priori Estimation Of Memory Effects In Coarse-Grained Nonlinear Systems Using The Mori-Zwanzig Formalism
• Entropy stable reduced order modeling of nonlinear conservation laws
• Projection-free approximate balanced truncation of large unstable systems
• Invariant Galerkin Ansatz Spaces and Davison-Maki Methods for the Numerical Solution of Differential Riccati Equations