A-priori sparsification of Galerkin-based reduced order models

Abstract

A methodology to generate sparse Galerkin models of chaotic/unsteady fluid flows containing a minimal number of active triadic interactions is proposed. The key idea is to find an appropriate set of basis functions for the projection representing elementary flow structures that interact minimally one with the other and thus result in a triadic interaction coefficient tensor with sparse structure. Interpretable and computationally efficient Galerkin models can be thus obtained, since a reduced number of triadic interactions needs to be computed to evaluate the right hand side of the model. To find the basis functions, a subspace rotation technique is used, whereby a set of Proper Orthogonal Decomposition (POD) modes is rotated into a POD subspace of larger dimension using coordinates associated to low-energy dissipative scales to alter energy paths and the structure of the triadic interaction coefficient tensor. This rotation is obtained as the solution of a non-convex optimisation problem that maximises the energy captured by the new basis, promotes sparsity and ensures long-term temporal stability of the sparse Galerkin system. We demonstrate the approach on two-dimensional lid-driven cavity flow at $Re=2\times 10^4$ where the motion is chaotic. We show that the procedure generates Galerkin models with a reduced set of active triadic interactions, distributed in modal space according to established knowledge of scale interactions in two-dimensional flows. This property, however, is only observed if long-term temporal stability is explicitly included in the formulation, indicating that a dynamical constraint is necessary to obtain a physically consistent sparsification.

References in corpus (8)

• Discovering governing equations from data: Sparse identification of nonlinear dynamical systems
• Sparsity-promoting dynamic mode decomposition
• A Novel Model Order Reduction Approach for Navier-Stokes Equations at High Reynolds Number
• Study of dynamics in post-transient flows using Koopman mode decomposition
• Minimal subspace rotation on the Stiefel manifold for stabilization and enhancement of projection-based reduced order models for the compressible Navier-Stokes equations
• Faster identification of optimal contraction sequences for tensor networks
• Energy transfer in turbulent channel flows and implications for resolvent modelling
• Energy transfer mechanisms and resolvent analysis in the cylinder wake