StabOp: A Data-Driven Stabilization Operator for Reduced Order Modeling

TL;DR

This paper introduces StabOp, a data-driven stabilization operator for reduced order models that improves accuracy and robustness over traditional filter-based methods, especially in convection-dominated flows.

Contribution

The paper proposes a novel data-driven stabilization operator (StabOp) that replaces traditional spatial filters in ROMs, optimized via PDE-constrained learning for enhanced stability and accuracy.

Findings

StabOp-L-ROM significantly outperforms classical L-ROM in accuracy.

The new method effectively stabilizes ROMs across various flow regimes.

StabOp provides a different smoothing mechanism than traditional filters.

Abstract

Spatial filters have played a central role in large eddy simulation and, more recently, in reduced order model (ROM) stabilization for convection-dominated flows. Nevertheless, important open questions remain: in under-resolved regimes, which filter is most suitable for a given stabilization or closure model? Moreover, once a filter is selected, how should its parameters, such as the filter radius, be determined? Addressing these questions is essential for the reliable design and performance of filter-based stabilization strategies. To answer these questions, we propose a novel strategy that differs fundamentally from current filter-based approaches: we replace traditional spatial filters with a data-driven stabilization operator (StabOp) that yields accurate results for a given resolution, quantity of interest, and stabilization strategy. Although the new StabOp can be used for both classical discretizations and ROMs, and for different types of filter-based stabilization or closure, for clarity we focus on ROMs and the Leray stabilization. To build the new StabOp, we postulate its model form as a linear, quadratic, or nonlinear mapping, and then solve a PDE-constrained optimization problem to minimize a given loss function. Using the resulting StabOp in the Leray ROM (L-ROM) yields a new stabilized ROM, StabOp-L-ROM. To assess the StabOp-L-ROM, we compare it with the L-ROM and the standard ROM in numerical simulations of four flows: 2D flow past a cylinder at Re=500, lid-driven cavity at Re=10000, 3D flow past a hemisphere at Re=2200, and minimal channel flow at Re=5000. Our numerical results show that the StabOp-L-ROM can be orders of magnitude more accurate than the classical L-ROM tuned with an optimal filter radius in the predictive regime. Furthermore, while the new StabOp smooths the input flow fields, its smoothing mechanism differs from that of classical spatial filters.