Discrete gradient methods for irreversible port-Hamiltonian systems

Alexandre Anahory Simoes,David Martín de Diego,Bernhard Maschke

2023-03-15

Abstract:In this paper we introduce discrete gradient methods to discretize irreversible port-Hamiltonian systems showing that the main qualitative properties of the continuous system are preserved using this kind discretizations methods.

Numerical Analysis,Mathematical Physics,Dynamical Systems

What problem does this paper attempt to address?

Energy-Preserving and Passivity-Consistent Numerical Discretization of Port-Hamiltonian Systems

Elena Celledoni,Eirik Hoel Høiseth

DOI: https://doi.org/10.48550/arXiv.1706.08621

2017-06-26

Numerical Analysis

Abstract:In this paper we design discrete port-Hamiltonian systems systematically in two different ways, by applying discrete gradient methods and splitting methods respectively. The discrete port-Hamiltonian systems we get satisfy a discrete notion of passivity, which lets us, by choosing the input appropriately, make them globally asymptotically stable with respect to an equilibrium point. We test methods designed using the discrete gradient approach in numerical experiments, and the results are encouraging when compared to relevant existing integrators of identical order.
Structure-Preserving Time Discretization of Port-Hamiltonian Systems via Discrete Gradient Pairs

Philipp Schulze

2023-11-01

Abstract:We discuss structure-preserving time discretization for nonlinear port-Hamiltonian systems with state-dependent mass matrix. Such systems occur, for instance, in the context of structure-preserving nonlinear model order reduction for port-Hamiltonian systems and, in this context, structure-preserving time discretization is crucial for preserving some of the properties of the time-continuous reduced-order model. For this purpose, we introduce a new class of time discretization schemes which is based on so-called discrete gradient pairs and leads to an exact power balance on the time-discrete level. Moreover, for the special case of a pointwise symmetric and positive definite mass matrix, we present an explicit construction of a discrete gradient pair. Finally, we illustrate the theoretical findings by means of a numerical example, where the time-continuous system is a nonlinear reduced-order model for an advection-diffusion problem.

Numerical Analysis
Discrete nonlinear elastodynamics in a port-Hamiltonian framework

Philipp L. Kinon,Tobias Thoma,Peter Betsch,Paul Kotyczka

DOI: https://doi.org/10.1002/pamm.202300144

2023-06-30

Abstract:We provide a fully nonlinear port-Hamiltonian formulation for discrete elastodynamical systems as well as a structure-preserving time discretization. The governing equations are obtained in a variational manner and represent index-1 differential algebraic equations. Performing an index reduction one obtains the port-Hamiltonian state space model, which features the nonlinear strains as an independent state next to position and velocity. Moreover, hyperelastic material behavior is captured in terms of a nonlinear stored energy function. The model exhibits passivity and losslessness and has an underlying symmetry yielding the conservation of angular momentum. We perform temporal discretization using the midpoint discrete gradient, such that the beneficial properties are inherited by the developed time stepping scheme in a discrete sense. The numerical results obtained in a representative example are demonstrated to validate the findings.

Dynamical Systems,Computational Engineering, Finance, and Science,Systems and Control,Numerical Analysis
Port maps of Irreversible Port Hamiltonian Systems

Bernhard Maschke,Jonas Kirchhoff

DOI: https://doi.org/10.48550/arXiv.2302.09026

2023-02-18

Abstract:Irreversible Port Hamiltonian Systems are departure of Port Hamiltonian Systems as they are generated not only by a Hamiltonian function but also by an entropy function and defined with respect to a quasi-Poisson bracket which embeds the definition of the irreversible phenomena taking place in the system. However the port map, consisting in the input map and the output map were poorly justified and lacked any physical consistency. In this paper, we suggest a novel definition of the port maps which allows to recover not only the energy balance equation (when the Hamiltonian equals the total energy of the system) but also a entropy balance equation including the irreversible entropy creation at the interface (the port) of the system in addition to the entropy creation term due to internal irreversible phenomena.

Dynamical Systems
A geometric framework for discrete time port-Hamiltonian systems

Karim Cherifi,Hannes Gernandt,Dorothea Hinsen,Volker Mehrmann

2023-08-18

Abstract:Port-Hamiltonian systems provide an energy-based formulation with a model class that is closed under structure preserving interconnection. For continuous-time systems these interconnections are constructed by geometric objects called Dirac structures. In this paper, we derive this geometric formulation and the interconnection properties for scattering passive discrete-time port-Hamiltonian systems.

Optimization and Control,Dynamical Systems
Conformal variational discretisation of infinite dimensional Hamiltonian systems with gradient flow dissipation

Damiano Lombardi,Cecilia Pagliantini

2024-12-09

Abstract:Nonconservative evolution problems describe irreversible processes and dissipative effects in a broad variety of phenomena. Such problems are often characterised by a conservative part, which can be modelled as a Hamiltonian term, and a nonconservative part, in the form of gradient flow dissipation. Traditional numerical approximations of this class of problem typically fail to retain the separation into conservative and nonconservative parts hence leading to unphysical solutions. In this work we propose a mixed variational method that gives a semi-discrete problem with the same geometric structure as the infinite-dimensional problem. As a consequence the conservation laws and the dissipative terms are retained. A priori convergence estimates on the solution are established. Numerical tests of the Korteweg-de Vries equation and of the two-dimensional Navier-Stokes equations on the torus and on the sphere are presented to corroborate the theoretical findings.

Numerical Analysis
On discrete-time dissipative port-Hamiltonian (descriptor) systems

Karim Cherifi,Hannes Gernandt,Dorothea Hinsen,Volker Mehrmann

2023-10-17

Abstract:Port-Hamiltonian (pH) systems have been studied extensively for linear continuous-time dynamical systems. This manuscript presents a discrete-time pH descriptor formulation for linear, completely causal, scattering passive dynamical systems based on the system coefficients. The relation of this formulation to positive and bounded real systems and the characterization via positive semidefinite solutions of Kalman-Yakubovich-Popov inequalities is also studied.

Optimization and Control,Dynamical Systems
Discrete-Time Models for Implicit Port-Hamiltonian Systems

Fernando Castaños,Hannah Michalska,Dmitry Gromov,Vincent Hayward

DOI: https://doi.org/10.48550/arXiv.1501.05097

2015-01-21

Abstract:Implicit representations of finite-dimensional port-Hamiltonian systems are studied from the perspective of their use in numerical simulation and control design. Implicit representations arise when a system is modeled in Cartesian coordinates and when the system constraints are applied in the form of additional algebraic equations (the system model is in a DAE form). Such representations lend themselves better to sample-data approximations. An implicit representation of a port-Hamiltonian system is given and it is shown how to construct a sampled-data model that preserves the port-Hamiltonian structure under sample and hold.

Systems and Control,Dynamical Systems,Numerical Analysis,Optimization and Control
Trends in ferromagnetism, hole localization, and acceptor level depth for Mn substitution in

Priya Mahadevan,Alex Zunger

DOI: https://doi.org/10.1063/1.1799245

2004-09-13

Abstract:We examine the intrinsic mechanism of ferromagnetism in dilute magnetic semiconductors by analyzing the trends in the electronic structure as the host is changed from GaN to GaP, GaAs and GaSb, keeping the transition metal impurity fixed. In contrast with earlier interpretations which depended on the host semiconductor, we found that a single mechanism is sufficient to explain the ferromagnetic stabilization energy for the entire series.

Other Condensed Matter
Port-Hamiltonian discontinuous Galerkin finite element methods

Nishant Kumar,J J W van der Vegt,H J Zwart

DOI: https://doi.org/10.1093/imanum/drae008

2024-04-27

IMA Journal of Numerical Analysis

Abstract:Abstract A port-Hamiltonian (pH) system formulation is a geometrical notion used to formulate conservation laws for various physical systems. The distributed parameter port-Hamiltonian formulation models infinite dimensional Hamiltonian dynamical systems that have a nonzero energy flow through the boundaries. In this paper, we propose a novel framework for discontinuous Galerkin (DG) discretizations of pH-systems. Linking DG methods with pH-systems gives rise to compatible structure preserving semidiscrete finite element discretizations along with flexibility in terms of geometry and function spaces of the variables involved. Moreover, the port-Hamiltonian formulation makes boundary ports explicit, which makes the choice of structure and power preserving numerical fluxes easier. We state the Discontinuous Finite Element Stokes–Dirac structure with a power preserving coupling between elements, which provides the mathematical framework for a large class of pH discontinuous Galerkin discretizations. We also provide an a priori error analysis for the port-Hamiltonian discontinuous Galerkin Finite Element Method (pH-DGFEM). The port-Hamiltonian discontinuous Galerkin finite element method is demonstrated for the scalar wave equation showing optimal rates of convergence.

mathematics, applied
About Dissipative and Pseudo Port-Hamiltonian Formulations of Irreversible Newtonian Compressible Flows

Luis A. Mora,Yann Le Gorrec,Denis Matignon,Hector Ramirez,Juan Yuz

DOI: https://doi.org/10.48550/arXiv.2003.11422

2020-03-25

Abstract:In this paper we consider the problem of obtaining a general port-Hamiltonian formulation of Newtonian fluids. We propose the port-Hamiltonian models to describe the energy flux of rotational three-dimensional isentropic and non-isentropic fluids, whose boundary flows and efforts can be used for control purposes or for power-preserving interconnection with other physical systems. In case of two-dimensional flows, we include the considerations about the operators associated with fluid vorticity, preserving the port-Hamiltonian structure of the models proposed.

Fluid Dynamics,Dynamical Systems
An Extended Discrete Gradient Formula for Oscillatory Hamiltonian Systems

Kai Liu,Wei Shi,Xinyuan Wu

DOI: https://doi.org/10.1088/1751-8113/46/16/165203

2015-01-01

Abstract:In this paper, incorporating the idea of the discrete gradient method into the extended Runge-Kutta-Nystrom integrator, we derive and analyze an extended discrete gradient formula for the oscillatory Hamiltonian system with the Hamiltonian H(p, q) = 1/2p(T) p+1/2q(T) Mq+U(q), where q : R -> R-d represents generalized positions, p : R -> R-d represents generalized momenta and M is an element of R (dxd) is a symmetric and positive semi-definite matrix. The solution of this system is a nonlinear oscillator. Basically, many nonlinear oscillatory mechanical systems with a partitioned Hamiltonian function lend themselves to this approach. The extended discrete gradient formula presented in this paper exactly preserves the energy H(p, q). We derive some properties of the new formula. The convergence is analyzed for the implicit schemes based on the discrete gradient formula, and it turns out that the convergence of the implicit schemes based on the extended discrete gradient formula is independent of parallel to M parallel to, which is a significant property for the oscillatory Hamiltonian system. Thus, it transpires that a larger step size can be chosen for the new energy-preserving schemes than that for the traditional discrete gradient methods when applied to the oscillatory Hamiltonian system. Illustrative examples show the competence and efficiency of the new schemes in comparison with the traditional discrete gradient methods in the scientific literature.
Port-Hamiltonian structures in infinite-dimensional optimal control: Primal-Dual gradient method and control-by-interconnection

Hannes Gernandt,Manuel Schaller

2024-06-03

Abstract:In this note, we consider port-Hamiltonian structures in numerical optimal control of ordinary differential equations. By introducing a novel class of nonlinear monotone port-Hamiltonian (pH) systems, we show that the primal-dual gradient method may be viewed as an infinite-dimensional nonlinear pH system. The monotonicity and the particular block structure arising in the optimality system is used to prove exponential stability of the dynamics towards its equilibrium, which is a critical point of the first order optimality conditions. Leveraging the port-based modeling, we propose an optimization-based controller in a suboptimal receding horizon control fashion. To this end, the primal-dual gradient based optimizer-dynamics is coupled to a pH plant dynamics in a power-preserving manner. We show that the resulting model is again monotone pH system and prove that the closed-loop exhibits local exponential convergence towards the equilibrium.

Optimization and Control
On implicit and explicit representations for 1D distributed port-Hamiltonian systems

Antoine Bendimerad-Hohl,Denis Matignon,Ghislain Haine,Laurent Lefèvre

2024-02-12

Abstract:First, two examples of 1D distributed port-Hamiltonian systems with dissipation, given in explicit (descriptor) form, are considered: the Dzekster model for the seepage of underground water and a nanorod model with non-local viscous damping. Implicit representations in Stokes-Lagrange subspaces are formulated. These formulations lead to modified Hamiltonian functions with spatial differential operators. The associated power balance equations are derived, together with the new boundary ports. Second, the port-Hamiltonian formulations for the Timoshenko and the Euler-Bernoulli beams are recalled, the latter being a flow-constrained version of the former. Implicit representations of these models in Stokes-Lagrange subspaces and corresponding power balance equations are derived. Bijective transformations are proposed between the explicit and implicit representations. It is proven these transformations commute with the flow-constraint projection operator.

Dynamical Systems,Analysis of PDEs
Discrete-time port-Hamiltonian systems: A definition based on symplectic integration

Paul Kotyczka,Laurent Lefèvre

DOI: https://doi.org/10.1016/j.sysconle.2019.104530

2019-11-01

Abstract:<p>We introduce a new definition of discrete-time port-Hamiltonian (PH) systems, which results from structure-preserving discretization of explicit PH systems in time. We discretize the underlying continuous-time Dirac structure with the collocation method and add discrete-time dynamics by the use of symplectic numerical integration schemes. The conservation of a structural discrete-time energy balance – expressed in terms of the discrete-time Dirac structure – extends the notion of symplecticity of geometric integration schemes to open systems. We discuss the energy approximation errors in the context of the presented definition and show that their order for linear PH systems is consistent with the order of the numerical integration scheme. Implicit Gauss–Legendre methods and Lobatto IIIA/IIIB pairs for partitioned systems are examples for integration schemes that are covered by our definition. The statements on the numerical energy errors are illustrated by elementary numerical experiments.</p>

automation & control systems,operations research & management science
Representing the dissipation of infinite-dimensional linear port-Hamiltonian systems

Friedrich M. Philipp

2024-05-28

Abstract:It is well known that linear and non-linear dissipative port-Hamiltonian systems in finite dimensions admit an energy balance, relating the energy increase in the system with the supplied energy and the dissipated energy. The integrand in the dissipation term is then a function of the state variable. In this note, we answer the question of when this is possible for linear port-Hamiltonian systems in infinite dimensions.

Analysis of PDEs
Existence of solutions to port-Hamiltonian systems: initial value problems and optimal control

Willem Esterhuizen,Bernhard Maschke,Till Preuster,Manuel Schaller,Karl Worthmann

2024-10-25

Abstract:We investigate the existence of solutions of reversible and irreversible port-Hamiltonian systems. To this end, we utilize the associated exergy, a function that is composed of the system's Hamiltonian and entropy, to prove global existence in time for bounded control functions. The results are then leveraged to prove existence of solutions of energy- and entropy-optimal control problems. Last, we explore model predictive control tailored to irreversible port-Hamiltonian systems by means of a numerical case study with a heat exchanger network.

Optimization and Control,Systems and Control
The port-Hamiltonian structure of continuum mechanics

Ramy Rashad,Stefano Stramigioli

2024-04-18

Abstract:In this paper we present a novel approach to the geometric formulation of solid and fluid mechanics within the port-Hamiltonian framework, which extends the standard Hamiltonian formulation to non-conservative and open dynamical systems. Leveraging Dirac structures, instead of symplectic or Poisson structures, this formalism allows the incorporation of energy exchange within the spatial domain or through its boundary, which allows for a more comprehensive description of continuum mechanics. Building upon our recent work in describing nonlinear elasticity using exterior calculus and bundle-valued differential forms, this paper focuses on the systematic derivation of port-Hamiltonian models for solid and fluid mechanics in the material, spatial, and convective representations using Hamiltonian reduction theory. This paper also discusses constitutive relations for stress within this framework including hyper-elasticity, for both finite- and infinite-strains, as well as viscous fluid flow governed by the Navier-Stokes equations.

Mathematical Physics,Differential Geometry,Dynamical Systems
A Hybrid Phase-Flow Method for Hamiltonian Systems with Discontinuous Hamiltonians

Shi Jin,Hao Wu,Zhongyi Huang

DOI: https://doi.org/10.1137/070709505

IF: 2.968

2008-01-01

SIAM Journal on Scientific Computing

Abstract:In this paper, we propose a new phase-flow method for Hamiltonian systems with discontinuous Hamiltonians. In the original phase-flow method introduced by Ying and Candès [J. Comput. Phys., 220 (2006), pp. 184-215], the phase map should be smooth to ensure the accuracy of the interpolation. Such an interpolation is inaccurate if the phase map is nonsmooth, for example, when the Hamiltonian is discontinuous. We modify the phase-flow method using a discontinuous Hamiltonian solver and establish the stability (for piecewise constant potentials) of such a solver. This extends the applicability of the highly efficient phase-flow method to singular Hamiltonian systems, with a mild increase of algorithm complexity. Such a particle method can be useful for the computation of high frequency waves through interfaces.

Discrete gradient methods for irreversible port-Hamiltonian systems

Energy-Preserving and Passivity-Consistent Numerical Discretization of Port-Hamiltonian Systems

Structure-Preserving Time Discretization of Port-Hamiltonian Systems via Discrete Gradient Pairs

Discrete nonlinear elastodynamics in a port-Hamiltonian framework

Port maps of Irreversible Port Hamiltonian Systems

A geometric framework for discrete time port-Hamiltonian systems

Conformal variational discretisation of infinite dimensional Hamiltonian systems with gradient flow dissipation

On discrete-time dissipative port-Hamiltonian (descriptor) systems

Discrete-Time Models for Implicit Port-Hamiltonian Systems

Trends in ferromagnetism, hole localization, and acceptor level depth for Mn substitution in

Port-Hamiltonian discontinuous Galerkin finite element methods

About Dissipative and Pseudo Port-Hamiltonian Formulations of Irreversible Newtonian Compressible Flows

An Extended Discrete Gradient Formula for Oscillatory Hamiltonian Systems

Port-Hamiltonian structures in infinite-dimensional optimal control: Primal-Dual gradient method and control-by-interconnection

On implicit and explicit representations for 1D distributed port-Hamiltonian systems

Discrete-time port-Hamiltonian systems: A definition based on symplectic integration

Representing the dissipation of infinite-dimensional linear port-Hamiltonian systems

Existence of solutions to port-Hamiltonian systems: initial value problems and optimal control

The port-Hamiltonian structure of continuum mechanics

A Hybrid Phase-Flow Method for Hamiltonian Systems with Discontinuous Hamiltonians