Intrusive and data-driven reduced order modelling of the rotating thermal shallow water equation

• Published in: Applied mathematics and computation 2022-05, Vol.421, p.126924, Article 126924
• Main Authors: Karasözen, Bülent, Yıldız, Süleyman, Uzunca, Murat
• Format: Article
• Language: English
• Subjects: Finite differences; Fluids; Hamiltonian systems; Least-squares; Model order reduction
• ISSN: 0096-3003; 1873-5649
• DOI: 10.1016/j.amc.2022.126924

•Intrusive and data-driven reduced-order solutions are equal up to numerical errors.•Incorporating the knowledge of parameter dependency, reduced-order solutions are predicted for new parameter values without interpolation.•The overall system behaviour is captured accurately in the training and prediction phases.•Preservation of polynomial invariants ensure the long-term stability of the reduced solution.•Operator inference with re-projection is tailored to learn physics-informed low-order models. In this paper, we investigate projection-based intrusive and data-driven model order reduction in numerical simulation of rotating thermal shallow water equation (RTSWE) in parametric and non-parametric form. Discretization of the RTSWE in space with centered finite differences leads to Hamiltonian system of ordinary differential equations with linear and quadratic terms. The full-order model (FOM) is obtained by applying linearly implicit Kahan’s method in time. Applying proper orthogonal decomposition with Galerkin projection (POD-G), we construct the intrusive reduced-order model (ROM). We apply operator inference (OpInf) with re-projection as data-driven ROM. In the parametric case, we make use of the parameter dependency at the level of the PDE without interpolating between the reduced operators. The least-squares problem of the OpInf is regularized with the minimum norm solution. Both ROMs behave similarly and are able to accurately predict the in the test and training data and capture system behaviour in the prediction phase with several orders of magnitude in computational speed-up over the FOM. The preservation of system physics such as the conserved quantities of the RTSWE by both ROMs enable that the models fit better to data and stable solutions are obtained in long-term predictions which are robust to parameter changes.