Energy-conserving hyper-reduction and temporal localization for reduced order models of the incompressible Navier-Stokes equations

Journal Article (2024)
Author(s)

R. B. Klein (Centrum Wiskunde & Informatica (CWI), TU Delft - Fluid Mechanics)

Benjamin Sanderse (Centrum Wiskunde & Informatica (CWI))

Research Group
Fluid Mechanics
Copyright
© 2024 R.B. Klein, B. Sanderse
DOI related publication
https://doi.org/10.1016/j.jcp.2023.112697
More Info
expand_more
Publication Year
2024
Language
English
Copyright
© 2024 R.B. Klein, B. Sanderse
Research Group
Fluid Mechanics
Volume number
499
Reuse Rights

Other than for strictly personal use, it is not permitted to download, forward or distribute the text or part of it, without the consent of the author(s) and/or copyright holder(s), unless the work is under an open content license such as Creative Commons.

Abstract

A novel hyper-reduction method is proposed that conserves kinetic energy and momentum for reduced order models of the incompressible Navier-Stokes equations. The main advantage of conservation of kinetic energy is that it endows the hyper-reduced order model (hROM) with a nonlinear stability property. The new method poses the discrete empirical interpolation method (DEIM) as a minimization problem and subsequently imposes constraints to conserve kinetic energy. Two methods are proposed to improve the robustness of the new method against error accumulation: oversampling and Mahalanobis regularization. Mahalanobis regularization has the benefit of not requiring additional measurement points. Furthermore, a novel method is proposed to perform energy- and momentum-conserving temporal localization with the principle interval decomposition: new interface conditions are derived such that energy and momentum are conserved for a full time-integration instead of only during separate intervals. The performance of the new energy- and momentum-conserving hyper-reduction methods and the energy- and momentum-conserving temporal localization method is analysed using three convection-dominated test cases; a shear-layer roll-up, two-dimensional homogeneous isotropic turbulence and a time-periodic inviscid flow consisting of a vortex in a uniform background flow. Our main finding is that energy conservation in combination with oversampling or regularization leads to a robust method with excellent long time stability properties. When any of these two ingredients is missing, accuracy and/or stability is significantly impaired.