Conservative model reduction for finite-volume models

Abstract

This work proposes a method for model reduction of finite-volume models that guarantees the resulting reduced-order model is conservative, thereby preserving the structure intrinsic to finite-volume discretizations. The proposed reduced-order models associate with optimization problems characterized by a minimum-residual objective function and nonlinear equality constraints that explicitly enforce conservation over subdomains. Conservative Galerkin projection arises from formulating this optimization problem at the time-continuous level, while conservative least-squares Petrov-Galerkin (LSPG) projection associates with a time-discrete formulation. We equip these approaches with hyper-reduction techniques in the case of nonlinear flux and source terms, and also provide approaches for handling infeasibility. In addition, we perform analyses that include deriving conditions under which conservative Galerkin and conservative LSPG are equivalent, as well as deriving a posteriori error bounds. Numerical experiments performed on a parameterized quasi-1D Euler equation demonstrate the ability of the proposed method to ensure not only global conservation, but also significantly lower state-space errors than nonconservative reduced-order models such as standard Galerkin and LSPG projection. (d) a posteriori bounds (Section 5.3) for the error in the quantities conserved over subdomains (Theorem 5.3), in the null space (Lemma 5.1) and row space (Lemma 5.2) of the constraints, in the full state (Theorem 5.2), and in the conserved quantities (Lemma 5.3 and Theorem 5.3).3. Numerical experiments on a parameterized quasi-1D Euler equation associated with modeling inviscid compressible flow in a converging-diverging nozzle (Section 6). These experiments demonstrate the merits of the proposed method and illustrate the importance of ensuring reduced-order models are globally conservative.We remark that this work was first presented publically at the "Recent Developments in Numerical Methods for Model Reduction" workshop at the Institut Henri Poincaré on November 10, 2016.Other works have also explored formulating reduced-order models that associate with constrained optimization problems. Zimmermann et al.[59] equip equality 'aerodynamic constraints' to ROMs applied to steady-state external flows, where the constraints associate with matching experimental data or target performance metrics in a design setting. Recently, Reddy et al. [45] propose equipping the time-discrete Galerkin ROM with inequality constraints that enforce solution positivity or a bound on the gas-void fraction. Relatedly, Fick et al. [26] proposed a modified Galerkin optimization problem applicable to the incompressible Navier-Stokes equations, where the inequality constraints associate with bounds on the generalized coordinates; these bounds correspond to the extreme values of the generalized coordinates arising during the training simulations.The remainder of this paper is organized as follows. Section 2 describes finite-volume discretizations of conservation ...

Physically constrained data‐driven correction for reduced‐order modeling of fluid flows

“…More recently, physical constraints have started to be used in standard (ie, without ROM) LES closure modeling (see, eg, the works of Duraisamy et al and Wang et al). Finally, physical constraints have also been used in standard ROM (ie, without closure modeling) . The CDDC‐ROM proposed in this paper uses physical constraints to improve the physical accuracy of the ROM closure model (ie, the Correction term in the DDC‐ROM).…”

“…Finally, physical constraints have also been used in standard ROM (ie, without closure modeling). [45][46][47][48][49][50][51][52] The CDDC-ROM proposed in this paper uses physical constraints to improve the physical accuracy of the ROM closure model (ie, the Correction term in the DDC-ROM).…”

Deep Conservation: A latent dynamics model for exact satisfaction of physical conservation laws [Report]

“…Further, many projection-based dynamics learning methods perform projection using a minimum-residual formulation [7] that does not preclude the violation of important physical properties such as conservation. To mitigate this issue, recent work has proposed a projection technique that explicitly enforces conservation over subdomains by adopting a constrained least-squares formulation to define the projection [8].…”

A locally conservative reduced flux reconstruction for elliptic problems

“…In the context of model order reduction of elliptic problems (1) by the reduced Basis (RB) method, it is sometimes desirable to obtain a locally conservative flux (4), for instance for a posteriori error estimation [1] or in the context of flow problems [2]. However, this is not straightforward for RB solutions (5), even if the underlying scheme (2) yields locally conservative approximations.…”

“…However, this is not straightforward for RB solutions (5), even if the underlying scheme (2) yields locally conservative approximations. While [3] ensure global mass conservation and [4] ensure mass conservation with respect to few subdomains, we are not aware of any work ensuring local mass conservation with respect to the underlying grid in the RB context. The present work closes this gap.…”