Source author record

Arjan van der Schaft

Arjan van der Schaft appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

33works
15topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

33 published item(s)

preprint2026arXiv

Hopfield neural networks as port-Hamiltonian and gradient systems

The structure of continuous Hopfield networks is revisited from a system-theoretic point of view. After adopting a novel electrical network interpretation involving nonlinear capacitors, it is shown that Hopfield networks admit a port-Hamiltonian formulation provided an extra passivity condition is satisfied. Subsequently it is shown that any Hopfield network can be represented as a gradient system, with Riemannian metric given by the inverse of the Hessian matrix of the total energy stored in the nonlinear capacitors. On the other hand, the well-known 'energy' function employed by Hopfield turns out to be the dissipation potential of the gradient system, and this potential is shown to satisfy a dissipation inequality that can be used for analysis and interconnection.

preprint2023arXiv

Reciprocity of nonlinear systems

One of the key contributions of the 1972 seminal paper by Willems was the analysis of symmetry (also called reciprocity) of input-state-output systems, both from an external (input-output) and internal (state) point of view. The developed theory also included the combination of reciprocity with passivity, and the consideration of relaxation systems, which are passive reciprocal systems without any oscillatory behavior. The paper was motivated from a fundamental system-theoretic point of view (how is external structure reflected into internal structure), as well as by a wide range of application areas, including electrical network synthesis, thermodynamics, and viscoelastic materials. On the other hand, the obtained results are for linear systems, and the extension to the nonlinear case, even for subclasses of nonlinear systems, is far from trivial. The present paper aims at taking some steps into this direction.

preprint2022arXiv

Indirect Adaptive Control of Nonlinearly Parameterized Nonlinear Dissipative Systems

In this note we address the problem of indirect adaptive (regulation or tracking) control of nonlinear, input affine dissipative systems. It is assumed that the supply rate, the storage and the internal dissipation functions may be expressed as nonlinearly parameterized regression equations where the mappings (depending on the unknown parameters) satisfy a monotonicity condition -- this encompasses a large class of physical systems, including passive systems. We propose to estimate the system parameters using the "power-balance" equation, which is the differential version of the classical dissipation inequality, with a new estimator that ensures global, exponential, parameter convergence under the very weak assumption of interval excitation of the power-balance equation regressor. To design the indirect adaptive controller we make the standard assumption of existence of an asymptotically stabilizing controller that depends -- possibly nonlinearly -- on the unknown plant parameters, and apply a certainty-equivalent control law. The benefits of the proposed approach, with respect to other existing solutions, are illustrated with examples.

preprint2022arXiv

Port-Hamiltonian systems and monotonicity

The relationships between port-Hamiltonian systems modeling and the notion of monotonicity are explored. The earlier introduced notion of incrementally port-Hamiltonian systems is extended to maximal cyclically monotone relations, together with their generating functions. This gives rise to new classes of incrementally port-Hamiltonian systems, with examples stemming from physical systems modeling as well as from convex optimization. An in-depth treatment is given of the composition of maximal monotone and maximal cyclically monotone relations, where in the latter case the resulting maximal cyclically monotone relation is shown to be computable through the use of generating functions. Furthermore, connections are discussed with incremental versions of passivity, and it is shown how incrementally port-Hamiltonian systems with strictly convex Hamiltonians are (maximal) equilibrium independent passive. Finally, the results on compositionality of monotone relations are employed for a convex optimization approach to the computation of the equilibrium of interconnected incrementally port-Hamiltonian systems.

preprint2020arXiv

A family of virtual contraction based controllers for tracking of flexible-joints port-Hamiltonian robots: theory and experiments

In this work we present a constructive method to design a family of virtual contraction based controllers that solve the standard trajectory tracking problem of flexible-joint robots (FJRs) in the port-Hamiltonian (pH) framework. The proposed design method, called virtual contraction based control (v-CBC), combines the concepts of virtual control systems and contraction analysis. It is shown that under potential energy matching conditions, the closed-loop virtual system is contractive and exponential convergence to a predefined trajectory is guaranteed. Moreover, the closed-loop virtual system exhibits properties such as structure preservation, differential passivity and the existence of (incrementally) passive maps.

preprint2020arXiv

Limits to Energy Conversion

A consequence of the Second Law of thermodynamics is that no thermodynamic system with a single heat source at constant temperature can convert heat into mechanical work in a recurrent manner. First we note that this is equivalent to cyclo-passivity at the mechanical port of the thermodynamic system, while the temperature at the thermal port of the system is kept constant. This leads to the general system-theoretic question which systems with two power ports have similar behavior: when is a system cyclo-passive at one of its ports, while the output variable at the other port (such as the temperature in the thermodynamic case) is kept constant? This property is called `one-port cyclo-passivity', and entails, whenever it holds, a fundamental limitation to energy transfer from one port (where the output is kept constant) to the other port. Sufficient conditions for one-port cyclo-passivity are derived for systems formulated in general port-Hamiltonian form. This is illustrated by a variety of examples from different {(multi-)}physical domains; from coupled inductors and capacitor microphones to synchronous machines.

preprint2017arXiv

Converse passivity theorems

Passivity is an imperative concept and a widely utilized tool in the analysis and control of interconnected systems. It naturally arises in the modelling of physical systems involving passive elements and dynamics. While many theorems on passivity are known in the theory of robust control, very few converse passivity results exist. This paper establishes various versions of converse passivity theorems for nonlinear feedback systems. In particular, open-loop passivity is shown to be necessary to ensure closed-loop passivity from an input-output perspective. Moreover, the stability of the feedback interconnection of a specific system with an arbitrary passive system is shown to imply passivity of the system itself.

preprint2016arXiv

A unifying energy-based approach to optimal frequency and market regulation in power grids

In this paper we provide a unifying energy-based approach to the modeling, analysis and control of power systems and markets, which is based on the port-Hamiltonian framework. Using a primal-dual gradient method applied to the social welfare problem, a distributed dynamic pricing algorithm in port-Hamiltonian form is obtained. By interconnection with the physical model a closed-loop port-Hamiltonian system is obtained, whose properties are exploited to prove asymptotic stability to the optimal points. This result is extended such that also general nodal power constraints are included into the social welfare problem. Additionally, the cases of line congestion and power transmission cost in (a)cyclic networks considering (non)linear power flow models are covered as well. Finally, we provide port-Hamiltonian descriptions and analysis of the well studied distributed averaging proportional integral (DAPI) and certain internal-model-based controllers, which solve an optimal frequency regulation problem.

preprint2016arXiv

A unifying energy-based approach to stability of power grids with market dynamics

In this paper a unifying energy-based approach is provided to the modeling and stability analysis of power systems coupled with market dynamics. We consider a standard model of the power network with a third-order model for the synchronous generators involving voltage dynamics. By applying the primal-dual gradient method to a social welfare optimization, a distributed dynamic pricing algorithm is obtained, which can be naturally formulated in port-Hamiltonian form. By interconnection with the physical model a closed-loop port-Hamiltonian system is obtained, whose properties are exploited to prove asymptotic stability to the set of optimal points. This result is extended to the case that also general nodal power constraints are included into the social welfare problem. Additionally, the case of line congestion and power transmission costs in acyclic networks is covered. Finally, a dynamic pricing algorithm is proposed that does not require knowledge about the power supply and demand.

preprint2016arXiv

Bisimulation equivalence of DAE systems

In this paper the notion of bisimulation relation for linear input-state-output systems is extended to general linear differential-algebraic (DAE) systems. Geometric control theory is used to derive a linear-algebraic characterization of bisimulation relations, and an algorithm for computing the maximal bisimulation relation between two linear DAE systems. The general definition is specialized to the case where the matrix pencil $sE - A$ is regular. Furthermore, by developing a one-sided version of bisimulation, characterizations of simulation and abstraction are obtained.

preprint2016arXiv

Interconnections of input-output Hamiltonian systems with dissipation

Recently, negative imaginary and counter-clockwise systems have attracted attention as an interesting class of systems, which is well-motivated by applications. In this paper first the formulation and extension of negative imaginary and counter-clockwise systems as (nonlinear) input-output Hamiltonian systems with dissipation is summarized. Next it is shown how by considering the time-derivative of the outputs a port-Hamiltonian system is obtained, and how this leads to the consideration of alternate passive outputs for port-Hamiltonian systems. Furthermore, a converse result to positive feedback interconnection of input-output Hamiltonian systems with dissipation is obtained, stating that the positive feedback interconnection of two linear systems is an input-output Hamiltonian system with dissipation if and only if the systems themselves are input-output Hamiltonian systems with dissipation. This implies that the Poisson and resistive structure matrices can be redefined in such a way that the interaction between the two systems only takes place via the coupling term in the Hamiltonian of the interconnected system. Subsequently, it is shown how network interconnection of (possibly nonlinear) input-output Hamiltonian systems with dissipation results in another input-output Hamiltonian system with dissipation, and how this leads to a stability analysis of the interconnected system in terms of the Hamiltonians and output mappings of the systems associated to the vertices, as well as the topology of the network.

preprint2016arXiv

Nonlinear Analysis of an Improved Swing Equation

In this paper, we investigate the properties of an improved swing equation model for synchronous generators. This model is derived by omitting the main simplifying assumption of the conventional swing equation, and requires a novel analysis for the stability and frequency regulation. We consider two scenarios. First we study the case that a synchronous generator is connected to a constant load. Second, we inspect the case of the single machine connected to an infinite bus. Simulations verify the results.

preprint2016arXiv

Optimal power dispatch in networks of high-dimensional models of synchronous machines

This paper investigates the problem of optimal frequency regulation of multi-machine power networks where each synchronous machine is described by a sixth order model. By analyzing the physical energy stored in the network and the generators, a port-Hamiltonian representation of the multi-machine system is obtained. Moreover, it is shown that the open-loop system is passive with respect to its steady states which implies that passive controllers can be used to control the multi-machine network. As a special case, a distributed consensus based controller is designed that regulates the frequency and minimizes a global quadratic generation cost in the presence of a constant unknown demand. In addition, the proposed controller allows freedom in choosing any desired connected undirected weighted communication graph.

preprint2016arXiv

Stability Analysis of Networked Systems Containing Damped and Undamped Nodes

This paper answers the question if a qualitatively heterogeneous passive networked system containing damped and undamped nodes shows consensus in the output of the nodes in the long run. While a standard Lyapunov analysis shows that the damped nodes will always converge to a steady-state value, the convergence of the undamped nodes is much more delicate and depends on the parameter values of the network as well as on the topology of the graph. A complete stability analysis is presented based on an eigenvector analysis involving the mass values and the topology of both the original graph and the reduced graph obtained by a Kron reduction that eliminates the damped nodes.

preprint2015arXiv

A geometric approach to differential Hamiltonian systems and differential Riccati equations

Motivated by research on contraction analysis and incremental stability/stabilizability the study of 'differential properties' has attracted increasing attention lately. Previously lifts of functions and vector fields to the tangent bundle of the state space manifold have been employed for a geometric approach to differential passivity and dissipativity. In the same vein, the present paper aims at a geometric underpinning and elucidation of recent work on 'control contraction metrics' and 'generalized differential Riccati equations'.

preprint2015arXiv

A network dynamics approach to chemical reaction networks

A crisp survey is given of chemical reaction networks from the perspective of general nonlinear network dynamics, in particular of consensus dynamics. It is shown how by starting from the complex-balanced assumption the reaction dynamics governed by mass action kinetics can be rewritten into a form which allows for a very simple derivation of a number of key results in chemical reaction network theory, and which directly relates to the thermodynamics of the system. Central in this formulation is the definition of a balanced Laplacian matrix on the graph of chemical complexes together with a resulting fundamental inequality. This directly leads to the characterization of the set of equilibria and their stability. Both the form of the dynamics and the deduced dynamical behavior are very similar to consensus dynamics. The assumption of complex-balancedness is revisited from the point of view of Kirchhoff's Matrix Tree theorem, providing a new perspective. Finally, using the classical idea of extending the graph of chemical complexes by an extra 'zero' complex, a complete steady-state stability analysis of mass action kinetics reaction networks with constant inflows and mass action outflows is given.

preprint2015arXiv

A Port-Hamiltonian Approach to Optimal Frequency Regulation in Power Grids

This paper studies the problem of frequency regulation in power grids, while maximizing the social welfare. Two price-based controllers are proposed; the first one an internal-model-based controller and the second one based on a continuous gradient method for optimization. Both controllers can be implemented in a fully distributed fashion, with freedom in choosing a controller communication network. As a result, two real-time dynamic pricing models described by port- Hamiltonian systems are obtained. By coupling with the port- Hamiltonian description of the physical network we obtain a closed-loop port-Hamiltonian system, whose properties are ex- ploited to prove asymptotic stability of the set of optimal points. Numerical results show the performance of both controllers in a simple case study.

preprint2015arXiv

Complex and detailed balancing of chemical reaction networks revisited

The characterization of the notions of complex and detailed balancing for mass action kinetics chemical reaction networks is revisited from the perspective of algebraic graph theory, in particular Kirchhoff's Matrix Tree theorem for directed weighted graphs. This yields an elucidation of previously obtained results, in particular with respect to the Wegscheider conditions, and a new necessary and sufficient condition for complex balancing, which can be verified constructively.

preprint2014arXiv

Stability of dynamical distribution networks with arbitrary flow constraints and unknown in/outflows

A basic model of a dynamical distribution network is considered, modeled as a directed graph with storage variables corresponding to every vertex and flow inputs corresponding to every edge, subject to unknown but constant inflows and outflows. We analyze the dynamics of the system in closed-loop with a distributed proportional-integral controller structure, where the flow inputs are constrained to take value in closed intervals. Results from our previous work are extended to general flow constraint intervals, and conditions for asymptotic load balancing are derived that rely on the structure of the graph and its flow constraints.

preprint2014arXiv

Structure-preserving model reduction of physical network systems by clustering

In this paper, we establish a method for model order reduction of a certain class of physical network systems. The proposed method is based on clustering of the vertices of the underlying graph, and yields a reduced order model within the same class. To capture the physical properties of the network, we allow for weights associated to both the edges as well as the vertices of the graph. We extend the notion of almost equitable partitions to this class of graphs. Consequently, an explicit model reduction error expression in the sense of H2-norm is provided for clustering arising from almost equitable partitions. Finally the method is extended to second-order systems.

preprint2013arXiv

Explicit Simplicial Discretization of Distributed-Parameter Port-Hamiltonian Systems

Simplicial Dirac structures as finite analogues of the canonical Stokes-Dirac structure, capturing the topological laws of the system, are defined on simplicial manifolds in terms of primal and dual cochains related by the coboundary operators. These finite-dimensional Dirac structures offer a framework for the formulation of standard input-output finite-dimensional port-Hamiltonian systems that emulate the behavior of distributed-parameter port-Hamiltonian systems. This paper elaborates on the matrix representations of simplicial Dirac structures and the resulting port-Hamiltonian systems on simplicial manifolds. Employing these representations, we consider the existence of structural invariants and demonstrate how they pertain to the energy shaping of port-Hamiltonian systems on simplicial manifolds.

preprint2013arXiv

On differential passivity of physical systems

Differential passivity is a property that allows to check with a pointwise criterion that a system is incrementally passive, a property that is relevant to study interconnected systems in the context of regulation, synchronization, and estimation. The paper investigates how restrictive is the property, focusing on a class of open gradient systems encountered in the coenergy modeling framework of physical systems, in particular the Brayton-Moser formalism for nonlinear electrical circuits.

preprint2012arXiv

A Graph-Theoretical Approach for the Analysis and Model Reduction of Complex-Balanced Chemical Reaction Networks

In this paper we derive a compact mathematical formulation describing the dynamics of chemical reaction networks that are complex-balanced and are governed by mass action kinetics. The formulation is based on the graph of (substrate and product) complexes and the stoichiometric information of these complexes, and crucially uses a balanced weighted Laplacian matrix. It is shown that this formulation leads to elegant methods for characterizing the space of all equilibria for complex-balanced networks and for deriving stability properties of such networks. We propose a method for model reduction of complex-balanced networks, which is similar to the Kron reduction method for electrical networks and involves the computation of Schur complements of the balanced weighted Laplacian matrix.

preprint2012arXiv

Discrete Exterior Geometry Approach to Structure-Preserving Discretization of Distributed-Parameter Port-Hamiltonian Systems

This paper addresses the issue of structure-preserving discretization of open distributed-parameter systems with Hamiltonian dynamics. Employing the formalism of discrete exterior calculus, we introduce a simplicial Dirac structure as a discrete analogue of the Stokes-Dirac structure and demonstrate that it provides a natural framework for deriving finite-dimensional port-Hamiltonian systems that emulate their infinite-dimensional counterparts. The spatial domain, in the continuous theory represented by a finite-dimensional smooth manifold with boundary, is replaced by a homological manifold-like simplicial complex and its augmented circumcentric dual. The smooth differential forms, in discrete setting, are mirrored by cochains on the primal and dual complexes, while the discrete exterior derivative is defined to be the coboundary operator. This approach of discrete differential geometry, rather than discretizing the partial differential equations, allows to first discretize the underlying Stokes-Dirac structure and then to impose the corresponding finite-dimensional port-Hamiltonian dynamics. In this manner, a number of important intrinsically topological and geometrical properties of the system are preserved.

preprint2012arXiv

Hamiltonian Perspective on Compartmental Reaction-Diffusion Networks

Inspired by the recent developments in modeling and analysis of reaction networks, we provide a geometric formulation of the reversible reaction networks under the influence of diffusion. Using the graph knowledge of the underlying reaction network, the obtained reaction-diffusion system is a distributed-parameter port-Hamiltonian system on a compact spatial domain. Motivated by the need for computer based design, we offer a spatially consistent discretization of the PDE system and, in a systematic manner, recover a compartmental ODE model on a simplicial triangulation of the spatial domain. Exploring the properties of a balanced weighted Laplacian matrix of the reaction network and the Laplacian of the simplicial complex, we characterize the space of equilibrium points and provide a simple stability analysis on the state space modulo the space of equilibrium points. The paper rules out the possibility of the persistence of spatial patterns for the compartmental balanced reaction-diffusion networks.

preprint2012arXiv

Model-order reduction of biochemical reaction networks

In this paper we propose a model-order reduction method for chemical reaction networks governed by general enzyme kinetics, including the mass-action and Michaelis-Menten kinetics. The model-order reduction method is based on the Kron reduction of the weighted Laplacian matrix which describes the graph structure of complexes in the chemical reaction network. We apply our method to a yeast glycolysis model, where the simulation result shows that the transient behaviour of a number of key metabolites of the reduced-order model is in good agreement with those of the full-order model.

preprint2012arXiv

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.

preprint2012arXiv

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.

preprint2011arXiv

On the Mathematical Structure of Balanced Chemical Reaction Networks Governed by Mass Action Kinetics

Motivated by recent progress on the interplay between graph theory, dynamics, and systems theory, we revisit the analysis of chemical reaction networks described by mass action kinetics. For reaction networks possessing a thermodynamic equilibrium we derive a compact formulation exhibiting at the same time the structure of the complex graph and the stoichiometry of the network, and which admits a direct thermodynamical interpretation. This formulation allows us to easily characterize the set of equilibria and their stability properties. Furthermore, we develop a framework for interconnection of chemical reaction networks. Finally we discuss how the established framework leads to a new approach for model reduction.

preprint2011arXiv

Structure-preserving tangential interpolation for model reduction of port-Hamiltonian Systems

Port-Hamiltonian systems result from port-based network modeling of physical systems and are an important example of passive state-space systems. In this paper, we develop the framework for model reduction of large-scale multi-input/multi-output port-Hamiltonian systems via tangential rational interpolation. The resulting reduced-order model not only is a rational tangential interpolant but also retains the port-Hamiltonian structure; hence is passive. This reduction methodology is described in both energy and co-energy system coordinates. We also introduce an $\mathcal{H}_2$-inspired algorithm for effectively choosing the interpolation points and tangential directions. The algorithm leads a reduced port-Hamiltonian model that satisfies a subset of $\mathcal{H}_2$-optimality conditions. We present several numerical examples that illustrate the effectiveness of the proposed method showing that it outperforms other existing techniques in both quality and numerical efficiency.