STRUCTURE-PRESERVING REDUCED BASIS METHODS FOR POISSON SYSTEMS

We develop structure-preserving reduced basis methods for a large class of nondissipative problems by resorting to their formulation as Hamiltonian dynamical systems. With this perspective, the phase space is naturally endowed with a Poisson manifold structure which encodes the physical properties, symmetries, and conservation laws of the dynamics. The goal is to design reduced basis methods for the general state-dependent degenerate Poisson structure based on a two-step approach. First, via a local approximation of the Poisson tensor, we split the Hamiltonian dynamics into an "almost symplectic" part and the trivial evolution of the Casimir invariants. Second, canonically symplectic reduced basis techniques are applied to the nontrivial component of the dynamics, preserving the local Poisson tensor kernel exactly. The global Poisson structure and the conservation properties of the phase flow are retained by the reduced model in the constant-valued case and up to errors in the Poisson tensor approximation in the state-dependent case. A priori error estimates for the solution of the reduced system are established. A set of numerical simulations is presented to corroborate the theoretical findings.