Dirac structures and boundary control systems associated with skew-symmetric differential operators

Associated with a skew-symmetric linear operator on the spatial domain $[a,b]$ we define a Dirac structure which includes the port variables on the boundary of this spatial domain. This Dirac structure is a subspace of a Hilbert space. Naturally, associated to this Dirac structure is infinite dimensional system. We parameterize the boundary port variables for which the $C_{0}$-semigroup associated to this system is contractive or unitary. Furthermore, this parameterization is used to split the boundary port variables into inputs and outputs. Similarly, we define a linear port controlled Hamiltonian system associated with the previously defined Dirac structure and a symmetric positive operator defining the energy of the system. We illustrate this theory on the example of the Timoshenko Beam. \u

Port Hamiltonian formulation of infinite dimensional systems I. Modeling

In this paper, some new results concerning the modeling of distributed parameter systems in port Hamiltonian form are presented. The classical finite dimensional port Hamiltonian formulation of a dynamical system is generalized in order to cope with the distributed parameter and multivariable case. The resulting class of infinite dimensional systems is quite general, thus allowing the description of several physical phenomena, such as heat conduction, piezoelectricity and elasticity. Furthermore, classical PDEs can be rewritten within this framework. The key point is the generalization of the notion of finite dimensional Dirac structure in order to deal with an infinite dimensional space of power variables

Boundary control for a class of dissipative differential operators including diffusion systems

In this paper we study a class of partial differential equations (PDE's), which includes Sturm-Liouville systems and diffusion equations. From this class of PDE's we define systems with control and observation through the boundary of the spatial domain. That is, we describe how to select boundary conditions, such that the resulting system has inputs and outputs acting through the boundary. Furthermore, these boundary conditions are chosen in a way that the resulting system has a nonincreasing energy.\u

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

Differential operator Dirac structures

As shown before, skew-adjoint linear differential operators, mapping efforts into flows, give rise to Dirac structures on a bounded spatial domain by a proper definition of boundary variables. In the present paper this is extended to pairs of linear differential operators defining a formally skew-adjoint relation between flows and efforts. Furthermore it is shown how the underlying repeated integration by parts operation is streamlined by the use of two-variable polynomial calculus. Dirac structures defined by formally skew adjoint operators together with differential operator effort constraints are treated within the same framework. Finally it is sketched how the approach can be also used for Lagrangian subspaces on bounded domains

Linear Boundary Port-Hamiltonian Systems with Implicitly Defined Energy

In this paper we extend the previously introduced class of boundary port-Hamiltonian systems to boundary control systems where the variational derivative of the Hamiltonian functional is replaced by a pair of reciprocal differential operators. In physical systems modelling, these differential operators naturally represent the constitutive relations associated with the implicitly defined energy of the system and obey Maxwell's reciprocity conditions. On top of the boundary variables associated with the Stokes-Dirac structure, this leads to additional boundary port variables and to the new notion of a Stokes-Lagrange subspace. This extended class of boundary port-Hamiltonian systems is illustrated by a number of examples in the modelling of elastic rods with local and non-local elasticity relations. Finally it shown how a Hamiltonian functional on an extended state space can be associated with the Stokes-Lagrange subspace, and how this leads to an energy balance equation involving the boundary variables of the Stokes-Dirac structure as well as of the Stokes-Lagrange subspace.Comment: 23 page

Coupling between hyperbolic and diffusive systems: A port-Hamiltonian formulation

The aim of this paper is to study a conservative wave equation coupled to a diffusion equation : this coupled system naturally arises in musical acoustics when viscous and thermal effects at the wall of the duct of a wind instrument are taken into account. The resulting equation, known as Webster-Lokshin model, has variable coefficients in space, and a fractional derivative in time. The port-Hamiltonian formalism proves adequate to reformulate this coupled system, and could enable another well-posedness analysis, using classical results from port-Hamiltonian systems theory. First, an equivalent formulation of fractional derivatives is obtained thanks to so-called diffusive representations: this is the reason why we first concentrate on rewriting these diffusive representations into the port-Hamiltonian formalism; two cases must be studied separately, the fractional integral operator as a low-pass filter, and the fractional derivative operator as a high-pass filter. Second, a standard finite-dimensional mechanical oscillator coupled to both types of dampings, either low-pass or high-pass, is studied as a coupled pHs. The more general PDE system of a wave equation coupled with the diffusion equation is then found to have the same structure as before, but in an appropriate infinite-dimensional setting, which is fully detailled

Pseudo-spectral methods for the spatial symplectic reduction of open systems of conservation laws

International audienceA reduction method is presented for systems of conservation laws with boundary energy flow. It is stated as a generalized pseudo-spectral method which performs exact differentiation by using simultaneously several approximation spaces generated by polynomials bases and suitable choices of port-variables. The symplecticity of this spatial reduction method is proved when used for the reduction of both closed and open systems of conservation laws, for any choice of collocation points (i.e. for any polynomial bases). The symplecticity of some more usual collocation schemes is discussed and finally their accuracy on approximation of the spectrum, on the example of the ideal transmission line, is discussed in comparison with the suggested reduction scheme

Reduction of Stokes-Dirac structures and gauge symmetry in port-Hamiltonian systems

Stokes-Dirac structures are infinite-dimensional Dirac structures defined in terms of differential forms on a smooth manifold with boundary. These Dirac structures lay down a geometric framework for the formulation of Hamiltonian systems with a nonzero boundary energy flow. Simplicial triangulation of the underlaying manifold leads to the so-called simplicial Dirac structures, discrete analogues of Stokes-Dirac structures, and thus provides a natural framework for deriving finite-dimensional port-Hamiltonian systems that emulate their infinite-dimensional counterparts. The port-Hamiltonian systems defined with respect to Stokes-Dirac and simplicial Dirac structures exhibit gauge and a discrete gauge symmetry, respectively. In this paper, employing Poisson reduction we offer a unified technique for the symmetry reduction of a generalized canonical infinite-dimensional Dirac structure to the Poisson structure associated with Stokes-Dirac structures and of a fine-dimensional Dirac structure to simplicial Dirac structures. We demonstrate this Poisson scheme on a physical example of the vibrating string