V. Talasila

Hamiltonian discretization of boundary control systems

A fundamental problem in the simulation and control of complex physical systems containing distributed-parameter components concerns finite-dimensional approximation. Numerical methods for partial differential equations (PDEs) usually assume the boundary conditions to be given, while more often than not the interaction of the distributed-parameter components with the other components takes place precisely via the boundary. On the other hand, finite-dimensional approximation methods for infinite-dimensional input-output systems (e.g., in semi-group format) are not easily relatable to numerical techniques for solving PDEs, and are mainly confined to linear PDEs. In this paper we take a new view on this problem by proposing a method for spatial discretization of boundary control systems based on a particular type of mixed finite elements, resulting in a finite-dimensional input-output system. The approach is based on formulating the distributed-parameter component as an infinite-dimensional port-Hamiltonian system, and exploiting the geometric structure of this representation for the choice of appropriate mixed finite elements. The spatially discretized system is again a port-Hamiltonian system, which can be treated as an approximating lumped-parameter physical system of the same type. In the current paper this program is carried out for the case of an ideal transmission line described by the telegraphers equations, and for the two-dimensional wave equation.

Approximation of the Telegrapher's equations

The problem of approximating a distributed parameter system with free boundary conditions is solved for the 1-dimensional Telegraphers equations. The Telegraphers equations are described using an infinite-dimensional port-Hamiltonian model, and we derive a finite dimensional port-Hamiltonian model using a mixed finite-element procedure. We show that energy conservation, passivity and some dynamic invariants are preserved in the discretization.

DISCRETE PORT HAMILTONIAN SYSTEMS

Either from a control theoretic viewpoint or from an analysis viewpoint it is necessary to convert smooth systems to discrete systems, which can then be implemented on computers for numerical simulations. Discrete models can be obtained either by discretizing a smooth model, or by directly modeling at the discrete level itself. One of the goals of this paper is to model port-Hamiltonian systems at the discrete level. We also show that the dynamics of the discrete models we obtain exactly correspond to the dynamics obtained via a usual discretization procedure. In this sense we offer an alternative to the usual procedure of modeling (at the smooth level) and discretization.