---
title: Port-Hamiltonian Systems
url: https://www.emergentmind.com/topics/port-hamiltonian-system
type: topic
---

# Port-Hamiltonian Systems

A port-Hamiltonian system (PHS) is a geometric, energy-based framework for the modeling, analysis, and control of open nonlinear, linear, finite- or infinite-dimensional, and multi-physical dynamical systems. PHS generalize classical Hamiltonian systems to allow for energy dissipation, constraints, boundary interactions, port-based interconnection, and coupling across physical domains. The structure is grounded in the use of Dirac structures for expressing energy-conserving interconnection laws and Lagrangian or resistive structures for constitutive and dissipative relations. At its core, a PHS encodes the dynamics and energetics of complex systems through the careful assignment of energy (Hamiltonian function), skew-adjoint interconnection, and dissipative mappings, often in state space or exterior-algebraic/generalized geometric representations. The theory unifies diverse domains—mechanics, electromagnetics, circuits, fluids, thermodynamics—under a compositional paradigm that is robust to modularity, discretization, and feedback interconnection.

## 1. Mathematical Framework and Dirac Structures

The port-Hamiltonian paradigm is built on three key components:

1. **Hamiltonian Function**: $H: X \to \mathbb{R}$, the total stored energy (generalized potential plus kinetic, electromagnetic, thermal, etc.), typically smooth and bounded below.

2. **Interconnection via Dirac Structures**: An energy-conserving relation $D \subseteq TX \oplus T^*X$ (in finite dimensions, $J^T=-J$), enforcing power-conservation $\langle e, \dot{x} \rangle + \langle e^p, f^p \rangle = 0$ and allowing for ports that interface subsystems or the environment [2405.01241, 2412.19673, 1107.2006, 2301.02024]. In networked or discretized settings, this extends naturally to open graphs or discrete manifolds using incidence matrices or discrete exterior calculus [1107.2006, 1201.5764].

3. **Dissipative or Resistive Structure**: $R(x)=R(x)^T \succeq 0$ captures irreversible (dissipative) effects, with positive semi-definite quadratic form or generally monotone relations [2412.19673, 2211.06676, 2206.09139]. The combination of Dirac (lossless interconnection) and resistive coupling yields a maximally monotone structure, crucial for Lyapunov stability and passivity [2211.06676].

The generic (finite-dimensional) state-space representation,
\[
\dot{x} = [J(x) - R(x)]\nabla H(x) + G(x)u,\qquad y = G(x)^T \nabla H(x),
\]
captures energy flow, with $J(x)$ skew-symmetric, $R(x)$ positive semi-definite, and $G(x)$ the port matrix for input $u$ and output $y$ [2301.02024].

### Energy Balance

The fundamental dissipation (power-balance) inequality holds:
\[
\dot{H}(x) = -\nabla H(x)^T R(x) \nabla H(x) + y^T u \leq y^T u,
\]
demonstrating the system is (input–output) passive, with $H$ as a storage or Lyapunov function [2412.19673, 2301.02024, 2211.06676].

## 2. Extensions: DAEs, Boundary Ports, and Infinite-Dimensional Systems

PHS generalize to systems with algebraic constraints, infinite-dimensional PDEs, and systems on time-varying domains:

- **PH-DAEs**: Descriptor forms account for algebraic constraints (e.g., circuits, constrained mechanics) by considering singular or rectangular mass matrices and algebraic constraint equations. The port-Hamiltonian-DAE structure, $E\dot{x} = Jz(x)-r(z(x)) + Bu$, unifies modeling for general networked, possibly index-2, systems [2211.06676, 1903.10451, 2004.12951].

- **Boundary/Distributed Parameter Systems**: Distributed-parameter PHS are governed by skew-adjoint differential (Stokes–Dirac) structures encoding spatial/temporal energy exchanges and include boundary ports for energy flow at domain boundaries. In PDEs, power-balance is ensured by integrating over the spatial domain and accounting for boundary fluxes [2404.12078, 2501.14930]. The Dirac formalism also applies to moving boundaries, allowing dynamic meshing and energy-stable discretizations [2501.14930].

- **Discrete Manifolds and Graphs**: The structure extends naturally to graph-based and discretized systems via discrete Dirac structures. For example, the coupling of flow and effort variables across simplicial complexes captures topological constraints and energy balances exactly in finite-dimensional settings [1107.2006, 1201.5764].

## 3. Modular Interconnection, Control, and Coupling

### Modular Coupling

PHS are closed under power-conserving interconnection: the composition of Dirac structures remains a Dirac structure, whether interconnecting subsystems, implementing boundary feedback, or assembling large-scale multiphysics models [2511.20150, 2301.02024, 2004.12951]. This guarantees that properties such as passivity, energy conservation/dissipation, and structural integrity are globally preserved.

### Control by Interconnection

Compositional control strategies exploit PHS structure:

- **Set-point stabilization and shaping**: Passivity-based feedback, Casimir function shaping, and interconnection with controller-structure PHS yield simple, physically interpretable, and robust control laws [2412.19673].
- **Energy/Power Ports**: Newer extensions distinguish between power ports (effort × flow = power) and energy ports (conjugate variables entering directly in the Hamiltonian), facilitating direct energy shaping [2405.01241, 2412.19673].
- **Distributed and Decoupled Simulation**: Structure-preserving couplings (linear/skew-symmetric interconnection of subsystems) and operator splitting for distributed simulation retain PHS properties at all scales [2511.20150].

## 4. Discretization, Numerical Methods, and Stochastic PHS

### Structure-Preserving Discretization

Time discretization using symplectic Runge–Kutta methods (collocation, Gauss–Legendre, discrete gradients) preserves discrete Dirac structures and yields discrete-time PHS with exact or high-order discrete-assured energy balances [1811.07852, 1903.10451, 2307.01351]. The passage from continuous to discrete PHS is systematically achieved by discretizing the underlying Dirac relation and enforcing energy contracts at each time step.

### Energy-Stable and Exergetic Formulations

Recent generalizations include energy-stable PHS (es-pH), which blend port-Hamiltonian and energy-stable system formalisms, and exergetic PHS, in which the Hamiltonian expresses exergy (maximum available work) in thermodynamic settings [2506.06471, 2008.04091]. These formulations support model reduction, optimal control, and extend to nonequilibrium thermodynamics (GENERIC structure correspondence).

### Stochastic Port-Hamiltonian Systems

PHS with stochastic ports enable power-conserving modeling under random perturbations. Each port can embed a noise process (semimartingale), and the stochastic Dirac structure ensures energy-balance and interconnection properties are maintained in expectation. The resulting S(PHS) encompass noise-driven systems, stochastic optimal control, and uncertainty quantification in physical networks [1910.01901].

## 5. Theoretical Properties: Passivity, Monotonicity, and Generic Controllability

### Passivity and Lyapunov Stability

Passivity—the property that the system cannot deliver more energy than it receives—follows from the structural conditions $R(x)\geq 0$, $J=-J^T$, and under mild convexity of $H$ ensures Lyapunov or asymptotic stability for $u=0$ configurations [2412.19673, 2211.06676, 2206.09139].

### Monotonicity and Incremental Passivity

PHS can be cast as monotone or maximally monotone operators. Incrementally port-Hamiltonian systems admit relations that are cyclically or maximally monotone, enabling convex optimization approaches to equilibrium computation and facilitating the analysis of incremental and differential passivity [2206.09139].

### Generic Controllability

In the linear case, controllability is a generic property of the port-Hamiltonian class: the set of controllable linear PHS is a relative generic subset of all such models, meaning that for almost all choices of parameters, the system is controllable [2104.02111]. Exceptional uncontrollable cases form a set of measure zero.

## 6. Physical and Engineering Applications

PHS are foundational in modeling:

- **Electrical and electronic circuits**: Including nonlinear networks, interconnecting components via Kirchhoff-Dirac structures, and subsuming MNA equations into an energy-structured DAE model [2004.10821, 2301.02024].
- **Mechanical and multi-body systems**: Mass–spring–damper networks, constrained multibody dynamics.
- **Distributed-parameter systems**: Continuum mechanics, Navier–Stokes fluids, Maxwell's equations, and coupled fluid–structure or electromagnetics–circuit models [2404.12078, 2501.14930, 2301.02024].
- **Thermodynamics and exergy-based systems**: Incorporating exergy as Hamiltonian, bond-graph methods, and facilitating thermodynamic optimization [2008.04091].

Beyond direct modeling, PHS provide systematic methodologies for structure-preserving model reduction, energy-based optimal control, structure-preserving numerical schemes, modular system construction, and robust fault-tolerant control for large multi-physics systems [2511.20150, 2506.06471, 1811.07852, 2301.02024].

## 7. Generalizations and Recent Developments

Research continues to expand the PHS framework:

- **Singular and Beyond-Passivity Systems**: Including formulations where the vector field admits singularities, yielding cyclo-dissipative systems capable of active power injection, relevant for power electronics and microgrid control [2602.10855].
- **Descriptor and Under/Overdetermined Systems**: Generalizations to under- and over-determined systems, higher-index DAEs, and arbitrary differentiable Hamiltonians, with Dirac structures accommodating algebraic constraints and invariance under nonlinear coordinate transformations [1903.10451, 2211.06676].
- **Discrete- and Time-Varying Domains**: Extension to moving boundary problems, discrete manifolds, and adaptive mesh settings, leveraging geometric structure for robust spatial-temporal discretization [2501.14930, 1201.5764].

A unifying thread across this body of work is that the port-Hamiltonian structure provides a geometric, modular, and physically interpretable framework that transcends disciplinary boundaries, ensuring that core energetic properties—conservation, dissipation, interconnection passivity—are maintained at every level of modeling, discretization, reduction, and control [2405.01241, 2412.19673, 2511.20150].

Source: https://www.emergentmind.com/topics/port-hamiltonian-system