Weak Form of Stokes-Dirac Structures and Geometric Discretization of Port-Hamiltonian Systems

We present the mixed Galerkin discretization of distributed parameter port-Hamiltonian systems. On the prototypical example of hyperbolic systems of two conservation laws in arbitrary spatial dimension, we derive the main contributions: (i) A weak formulation of the underlying geometric (Stokes-Dirac) structure with a segmented boundary according to the causality of the boundary ports. (ii) The geometric approximation of the Stokes-Dirac structure by a finite-dimensional Dirac structure is realized using a mixed Galerkin approach and power-preserving linear maps, which define minimal discrete power variables. (iii) With a consistent approximation of the Hamiltonian, we obtain finite-dimensional port-Hamiltonian state space models. By the degrees of freedom in the power-preserving maps, the resulting family of structure-preserving schemes allows for trade-offs between centered approximations and upwinding. We illustrate the method on the example of Whitney finite elements on a 2D simplicial triangulation and compare the eigenvalue approximation in 1D with a related approach.Comment: Copyright 2018. This manuscript version is made available under the CC-BY-NC-ND 4.0 license http://creativecommons.org/licenses/by-nc-nd/4.0

Script geometry

In this paper we describe the foundation of a new kind of discrete geometry and calculus called Script Geometry. It allows to work with more general meshes than classic simplicial complexes. We provide the basic denitions as well as several examples, like the Klein bottle and the projective plane. Furthermore, we also introduce the corresponding Dirac and Laplace operators which should lay the groundwork for the development of the corresponding discrete function theory

Port-Hamiltonian discretization for open channel flows

A finite-dimensional Port-Hamiltonian formulation for the dynamics of smooth open channel flows is presented. A numerical scheme based on this formulation is developed for both the linear and nonlinear shallow water equations. The scheme is verified against exact solutions and has the advantage of conservation of mass and energy to the discrete level

Lattice calculations on the spectrum of Dirac and Dirac-K\"ahler operators

We present a matrix technique to obtain the spectrum and the analytical index of some elliptic operators defined on compact Riemannian manifolds. The method uses matrix representations of the derivative which yield exact values for the derivative of a trigonometric polynomial. These matrices can be used to find the exact spectrum of an elliptic operator in particular cases and in general, to give insight into the properties of the solution of the spectral problem. As examples, the analytical index and the eigenvalues of the Dirac operator on the torus and on the sphere are obtained and as an application of this technique, the spectrum of the Dirac-Kahler operator on the sphere is explored.Comment: 11 page

Port-Hamiltonian systems on discrete manifolds

This paper offers a geometric framework for modeling port-Hamiltonian systems on discrete manifolds. The simplicial Dirac structure, capturing the topological laws of the system, is defined in terms of primal and dual cochains related by the coboundary operators. This finite-dimensional Dirac structure, as discrete analogue of the canonical Stokes-Dirac structure, allows for the formulation of finite-dimensional port-Hamiltonian systems that emulate the behaviour of the open distributed-parameter systems with Hamiltonian dynamics.Comment: MATHMOD 2012 - 7th Vienna International Conference on Mathematical Modellin