Structure preserving model order reduction of shallow water equations

Karasözen, Bülent
In this paper, we present two different approaches for constructing reduced-order models (ROMs) for the two-dimensional shallow water equation (SWE). The first one is based on the noncanonical Hamiltonian/Poisson form of the SWE. After integration in time by the fully implicit average vector field method, ROMs are constructed with proper orthogonal decomposition(POD)/discrete empirical interpolation method that preserves the Hamiltonian structure. In the second approach, the SWE as a partial differential equation with quadratic nonlinearity is integrated in time by the linearly implicit Kahan's method, and ROMs are constructed with the tensorial POD that preserves the linear-quadratic structure of the SWE. We show that in both approaches, the invariants of the SWE such as the energy, enstrophy, mass and circulation are preserved over a long period of time, leading to stable solutions. We conclude by demonstrating the accuracy and the computational efficiency of the reduced solutions by a numerical test problem.


Energy preserving model order reduction of the nonlinear Schrodinger equation
Karasözen, Bülent (2018-12-01)
An energy preserving reduced order model is developed for two dimensional nonlinear Schrodinger equation (NLSE) with plane wave solutions and with an external potential. The NLSE is discretized in space by the symmetric interior penalty discontinuous Galerkin (SIPG) method. The resulting system of Hamiltonian ordinary differential equations are integrated in time by the energy preserving average vector field (AVF) method. The mass and energy preserving reduced order model (ROM) is constructed by proper orth...
Reduced order modelling of nonlinear cross-diffusion systems
Karasözen, Bülent; Uzunca, Murat; Yıldız, Süleyman (2021-07-15)
In this work, we present reduced-order models (ROMs) for a nonlinear cross-diffusion problem from population dynamics, the Shigesada-Kawasaki-Teramoto (SKT) equation with Lotka-Volterra kinetics. The formation of the patterns of the SKT equation consists of a fast transient phase and a long stationary phase. Reduced order solutions are computed by separating the time into two time-intervals. In numerical experiments, we show for one- and two-dimensional SKT equations with pattern formation, the reduced-orde...
Model order reduction for nonlinear Schrodinger equation
Karasözen, Bülent; Uzunca, Murat (2015-05-01)
We apply the proper orthogonal decomposition (POD) to the nonlinear Schrodinger (NLS) equation to derive a reduced order model. The NLS equation is discretized in space by finite differences and is solved in time by structure preserving symplectic mid-point rule. A priori error estimates are derived for the POD reduced dynamical system. Numerical results for one and two dimensional NLS equations, coupled NLS equation with soliton solutions show that the low-dimensional approximations obtained by POD reprodu...
Structure-preserving reduced-order modeling of Korteweg–de Vries equation
Uzunca, Murat; Karasözen, Bülent; Yıldız, Süleyman (2021-10-01)
Computationally efficient, structure-preserving reduced-order methods are developed for the Korteweg–de Vries (KdV) equations in Hamiltonian form. The semi-discretization in space by finite differences is based on the Hamiltonian structure. The resulting skew-gradient system of ordinary differential equations (ODEs) is integrated with the linearly implicit Kahan's method, which preserves the Hamiltonian approximately. We have shown, using proper orthogonal decomposition (POD), the Hamiltonian structure of t...
Reduced order optimal control of the convective FitzHugh-Nagumo equations
Karasözen, Bülent; KÜÇÜKSEYHAN, TUĞBA (2020-02-15)
In this paper, we compare three model order reduction methods: the proper orthogonal decomposition (POD), discrete empirical interpolation method (DEIM) and dynamic mode decomposition (DMD) for the optimal control of the convective FitzHugh-Nagumo (FHN) equations. The convective FHN equations consist of the semi-linear activator and the linear inhibitor equations, modeling blood coagulation in moving excitable media. The semilinear activator equation leads to a non-convex optimal control problem (OCP). The ...
Citation Formats
B. Karasözen and M. UZUNCA, “Structure preserving model order reduction of shallow water equations,” MATHEMATICAL METHODS IN THE APPLIED SCIENCES, pp. 0–0, 2020, Accessed: 00, 2020. [Online]. Available: