Model order reduction for compressible flows solved using the discontinuous Galerkin methods

Abstract

Projection-based reduced order models (ROM) based on the weak form and the strong form of the discontinuous Galerkin (DG) method are proposed and compared for shock-dominated problems. The incorporation of dissipation components of DG in a consistent manner, including the upwinding flux and the localized artificial viscosity model, is employed to enhance stability of the ROM. To ensure efficiency, the discrete empirical interpolation method (DEIM) is adopted to enable hyper-reduction, for which the upwinding flux is decomposed into the central part and the dissipation part. The maximum local wave speed in the upwinding dissipation part is compressed and approximated using the DEIM approach, and the same strategy is applied to the artificial viscosity. Energy stability is proved with the strong-form-based ROM prior to hyper-reduction for the one-dimensional scalar cases. Eigenvalue spectrum is analyzed to verify and compare the stability properties of the two proposed ROMs. Several benchmark cases are conducted to test the performance of the proposed models. Results show that stable computations with reasonable acceleration for shock-dominated cases can be achieved with the ROM built on the strong form.