A Geometric Approach to Dynamical Model-Order Reduction

Abstract

Any model order reduced dynamical system that evolves a modal decomposition to approximate the discretized solution of a stochastic PDE can be related to a vector field tangent to the manifold of fixed rank matrices. The Dynamically Orthogonal (DO) approximation is the canonical reduced order model for which the corresponding vector field is the orthogonal projection of the original system dynamics onto the tangent spaces of this manifold. The embedded geometry of the fixed rank matrix manifold is thoroughly analyzed. The curvature of the manifold is characterized and related to the smallest singular value through the study of the Weingarten map. Differentiability results for the orthogonal projection onto embedded manifolds are reviewed and used to derive an explicit dynamical system for tracking the truncated Singular Value Decomposition (SVD) of a time-dependent matrix. It is demonstrated that the error made by the DO approximation remains controlled under the minimal condition that the original solution stays close to the low rank manifold, which translates into an explicit dependence of this error on the gap between singular values. The DO approximation is also justified as the dynamical system that applies instantaneously the SVD truncation to optimally constrain the rank of the reduced solution. Riemannian matrix optimization is investigated in this extrinsic framework to provide algorithms that adaptively update the best low rank approximation of a smoothly varying matrix. The related gradient flow provides a dynamical system that converges to the truncated SVD of an input matrix for almost every initial data.

Cited by in corpus (10)

• Variational approach to closure of nonlinear dynamical systems: Autonomous case
• Neural Closure Models for Dynamical Systems
• Model order reduction for stochastic dynamical systems with continuous symmetries
• On Asymptotic Linear Convergence of Projected Gradient Descent for Constrained Least Squares
• Implicit low-rank Riemannian schemes for the time integration of stiff partial differential equations
• Bayesian Learning of Coupled Biogeochemical-Physical Models
• Sensitivity of low-rank matrix recovery
• Stability properties of a projector-splitting scheme for dynamical low rank approximation of random parabolic equations
• Low-rank optimization methods based on projected projected-gradient descent that accumulate at Bouligand stationary points
• Model order reduction via Lie groups