Structure-Preserving Finite Volume Method for 2D Linear and Non-Linear Port-Hamiltonian Systems

Abstract

In this work we extend the results of a high order finite volume semi-discretization for port-Hamiltonian system 1D linear case (Kotyczka (2016)) to the 2D linear case, worked on the wave equation. The existing pHs discretization methods deal only with the geometric part, in this paper we perform an adapted symplectic time stepping to get the fully discrete scheme in order to preserve both the geometrical properties and the energy aspects. We also show that staggered finite volume method carry over to a non-linear problem, the 2D irrotational shallow water equations. However, due to the non linearity and the non separability of the Hamiltonian, some difficulties arise both for the high order accuracy in the spatial discretization, and also for the symplecticity of the time integration.