Quantized local reduced-order modeling in time (ql-ROM)

TL;DR

This paper introduces a quantized local reduced-order modeling (ql-ROM) approach that partitions complex chaotic systems into clusters, enabling more stable and accurate local models for nonlinear PDEs like Kuramoto-Sivashinsky and Navier-Stokes.

Contribution

The paper proposes a novel ql-ROM framework that improves stability and accuracy of reduced-order models by local clustering and modeling, addressing limitations of global ROMs.

Findings

Enhanced numerical stability over global ROMs

Improved short-term prediction accuracy

Accurate long-term statistical predictions

Abstract

Spatiotemporally chaotic systems, such as the solutions of some nonlinear partial differential equations, are dynamical systems that evolve toward a lower dimensional manifold. This manifold has an intricate geometry with heterogeneous density, which makes the design of a single (global) nonlinear reduced-order model (ROM) challenging. In this paper, we turn this around. Instead of modeling the manifold with one single model, we partition the manifold into clusters within which the dynamics are locally modeled. This results in a quantized local reduced-order model (ql-ROM), which consists of (i) quantizing the manifold via unsupervised clustering; (ii) constructing intrusive ROMs for each cluster; and (iii) seamlessly patch the local models with a change of basis and assignment functions. We test the method on two nonlinear partial differential equations, i.e., the Kuramoto-Sivashinsky and 2D Navier-Stokes equations (Kolmogorov flow), across bursting, chaotic, quasiperiodic, and turbulent regimes. The local models are built via Galerkin projection onto the local principal directions, which are centered on the cluster centroids. The dynamics are modeled by switching a local ROM based on the cluster proximity. The proposed ql-ROM framework has three advantages over global ROMs (g-ROMs): (i) numerical stability, (ii) improved short-term prediction accuracy in time, and (iii) accurate prediction of long-term statistics, such as energy spectra and probability distributions. The computational overhead is minimal with respect to g-ROMs. The proposed framework retains the interpretability and simplicity of intrusive projection-based ROMs, whilst overcoming their limitations in modeling complex, high-dimensional, nonlinear dynamics.