-
Minimal motifs for habituating systems
Authors:
Matthew Smart,
Stanislav Y. Shvartsman,
Martin Mönnigmann
Abstract:
Habituation - a phenomenon in which a dynamical system exhibits a diminishing response to repeated stimulations that eventually recovers when the stimulus is withheld - is universally observed in living systems from animals to unicellular organisms. Despite its prevalence, generic mechanisms for this fundamental form of learning remain poorly defined. Drawing inspiration from prior work on systems…
▽ More
Habituation - a phenomenon in which a dynamical system exhibits a diminishing response to repeated stimulations that eventually recovers when the stimulus is withheld - is universally observed in living systems from animals to unicellular organisms. Despite its prevalence, generic mechanisms for this fundamental form of learning remain poorly defined. Drawing inspiration from prior work on systems that respond adaptively to step inputs, we study habituation from a nonlinear dynamics perspective. This approach enables us to formalize classical hallmarks of habituation that have been experimentally identified in diverse organisms and stimulus scenarios. We use this framework to investigate distinct dynamical circuits capable of habituation. In particular, we show that driven linear dynamics of a memory variable with static nonlinearities acting at the input and output can implement numerous hallmarks in a mathematically interpretable manner. This work establishes a foundation for understanding the dynamical substrates of this primitive learning behavior and offers a blueprint for the identification of habituating circuits in biological systems.
△ Less
Submitted 25 July, 2024;
originally announced July 2024.
-
A minimal dynamical system and analog circuit for non-associative learning
Authors:
Matthew Smart,
Stanislav Y. Shvartsman,
Martin Mönnigmann
Abstract:
Learning in living organisms is typically associated with networks of neurons. The use of large numbers of adjustable units has also been a crucial factor in the continued success of artificial neural networks. In light of the complexity of both living and artificial neural networks, it is surprising to see that very simple organisms -- even unicellular organisms that do not possess a nervous syst…
▽ More
Learning in living organisms is typically associated with networks of neurons. The use of large numbers of adjustable units has also been a crucial factor in the continued success of artificial neural networks. In light of the complexity of both living and artificial neural networks, it is surprising to see that very simple organisms -- even unicellular organisms that do not possess a nervous system -- are capable of certain forms of learning. Since in these cases learning may be implemented with much simpler structures than neural networks, it is natural to ask how simple the building blocks required for basic forms of learning may be. The purpose of this study is to discuss the simplest dynamical systems that model a fundamental form of non-associative learning, habituation, and to elucidate technical implementations of such systems, which may be used to implement non-associative learning in neuromorphic computing and related applications.
△ Less
Submitted 9 May, 2024;
originally announced May 2024.
-
A reduced model for particle calcination for use in DEM/CFD simulations
Authors:
Lucas Reineking,
Torben Bergold,
Enric Illana,
Viktor Scherer,
Martin Mönnigmann
Abstract:
We treat the accurate simulation of the calcination reaction in particles, where the particles are large and, thus, the inner-particle processes must be resolved. Because these processes need to be described with coupled partial differential equations that must be solved numerically, the computation times for a single particle are too high for use in simulations that involve many particles. Simula…
▽ More
We treat the accurate simulation of the calcination reaction in particles, where the particles are large and, thus, the inner-particle processes must be resolved. Because these processes need to be described with coupled partial differential equations that must be solved numerically, the computation times for a single particle are too high for use in simulations that involve many particles. Simulations of this type arise when the Discrete Element Method (DEM) is combined with Computational Fluid Dynamics (CFD) to investigate industrial systems such as quick lime production in lime shaft kilns. We show that, based on proper orthogonal de-composition and Galerkin projection, reduced models can be derived for single particles that provide the same spatial and temporal resolution as the original PDE models at a considerably reduced computational cost. Replacing the finite-volume particle models with the reduced models results in an overall reduction of the reactor simulation time by about 60% for the simple example treated here.
△ Less
Submitted 26 November, 2023;
originally announced November 2023.
-
Parametric 3D Convolutional Autoencoder for the Prediction of Flow Fields in a Bed Configuration of Hot Particles
Authors:
Ali Mjalled,
Reza Namdar,
Lucas Reineking,
Mohammad Norouzi,
Fathollah Varnik,
Martin Mönnigmann
Abstract:
The use of deep learning methods for modeling fluid flow has drawn a lot of attention in the past few years. In situations where conventional numerical approaches can be computationally expensive, these techniques have shown promise in offering accurate, rapid, and practical solutions for modeling complex fluid flow problems. The success of deep learning is often due to its ability to extract hidd…
▽ More
The use of deep learning methods for modeling fluid flow has drawn a lot of attention in the past few years. In situations where conventional numerical approaches can be computationally expensive, these techniques have shown promise in offering accurate, rapid, and practical solutions for modeling complex fluid flow problems. The success of deep learning is often due to its ability to extract hidden patterns and features from the data, enabling the creation of data-driven reduced models that can capture the underlying physics of the domain. We present a data-driven reduced model for predicting flow fields in a bed configuration of hot particles. The reduced model consists of a parametric 3D convolutional autoencoder. The first part resolves the spatial and temporal dependencies present in the input sequence, while the second part of the architecture is responsible for predicting the solution at the subsequent timestep based on the information gathered from the preceding part. We also propose the utilization of a post-processing non-trainable output layer following the decoding path to incorporate the physical knowledge, e.g., no-slip condition, into the prediction. The evaluation of the reduced model for a bed configuration with variable particle temperature showed accurate results at a fraction of the computational cost required by traditional numerical simulation methods.
△ Less
Submitted 12 February, 2024; v1 submitted 6 July, 2023;
originally announced July 2023.
-
Exploiting symmetries in active set enumeration for constrained linear-quadratic optimal control
Authors:
Ruth Mitze,
Michal Kvasnica,
Martin Mönnigmann
Abstract:
This paper studies symmetric constrained linear-quadratic optimal control problems and their parametric solutions. The parametric solution of such a problem is a piecewise-affine feedback law that can be equivalently expressed as a set of active sets. We show symmetries of the optimal control problem entail symmetries of the active sets, which can be used to simplify finding the set of active sets…
▽ More
This paper studies symmetric constrained linear-quadratic optimal control problems and their parametric solutions. The parametric solution of such a problem is a piecewise-affine feedback law that can be equivalently expressed as a set of active sets. We show symmetries of the optimal control problem entail symmetries of the active sets, which can be used to simplify finding the set of active sets considerably. Specifically, we improve a recently proposed method for the dynamic-programming-based enumeration of all active sets. The achieved reduction of the computational effort is illustrated with an example.
△ Less
Submitted 21 March, 2023;
originally announced March 2023.
-
Estimating flow fields with Reduced Order Models
Authors:
Kamil David Sommer,
Lucas Reineking,
Yogesh Parry Ravichandran,
Romuald Skoda,
Martin Mönnigmann
Abstract:
The estimation of fluid flows inside a centrifugal pump in realtime is a challenging task that cannot be achieved with long-established methods like CFD due to their computational demands. We use a projection-based reduced order model (ROM) instead. Based on this ROM, a realtime observer can be devised that estimates the temporally and spatially resolved velocity and pressure fields inside the pum…
▽ More
The estimation of fluid flows inside a centrifugal pump in realtime is a challenging task that cannot be achieved with long-established methods like CFD due to their computational demands. We use a projection-based reduced order model (ROM) instead. Based on this ROM, a realtime observer can be devised that estimates the temporally and spatially resolved velocity and pressure fields inside the pump. The entire fluid-solid domain is treated as a fluid in order to be able to consider moving rigid bodies in the reduction method. A greedy algorithm is introduced for finding suitable and as few measurement locations as possible. Robust observability is ensured with an extended Kalman filter, which is based on a time-variant observability matrix obtained from the nonlinear velocity ROM. We present the results of the velocity and pressure ROMs based on a unsteady Reynolds-averaged Navier-Stokes CFD simulation of a 2D centrifugal pump, as well as the results for the extended Kalman filter.
△ Less
Submitted 7 July, 2023; v1 submitted 11 February, 2022;
originally announced February 2022.
-
Teaching MPC: Which Way to the Promised Land?
Authors:
Timm Faulwasser,
Sergio Lucia,
Moritz Schulze Darup,
Martin Mönnigmann
Abstract:
Since the earliest conceptualizations by Lee and Markus, and Propoi in the 1960s, Model Predictive Control (MPC) has become a major success story of systems and control with respect to industrial impact and with respect to continued and wide-spread research interest. The field has evolved from conceptually simple linear-quadratic (convex) settings in discrete and continuous time to nonlinear and d…
▽ More
Since the earliest conceptualizations by Lee and Markus, and Propoi in the 1960s, Model Predictive Control (MPC) has become a major success story of systems and control with respect to industrial impact and with respect to continued and wide-spread research interest. The field has evolved from conceptually simple linear-quadratic (convex) settings in discrete and continuous time to nonlinear and distributed settings including hybrid, stochastic, and infinite-dimensional systems. Put differently, essentially the entire spectrum of dynamic systems can be considered in the MPC framework with respect to both -- system theoretic analysis and tailored numerics. Moreover, recent developments in machine learning also leverage MPC concepts and learning-based and data-driven MPC have become highly active research areas.
However, this evident and continued success renders it increasingly complex to live up to industrial expectations while enabling graduate students for state-of-the-art research in teaching MPC. Hence, this position paper attempts to trigger a discussion on teaching MPC. To lay the basis for a fruitful debate, we subsequently investigate the prospect of covering MPC in undergraduate courses; we comment on teaching textbooks; and we discuss the increasing complexity of research-oriented graduate teaching of~MPC.
△ Less
Submitted 2 June, 2021;
originally announced June 2021.
-
Modeling pressure pulsation and backflow in progressing cavity pumps with deformable stator
Authors:
Jens Müller,
Yashar Kouhi,
Sebastian Leonow,
Martin Mönnigmann
Abstract:
This contribution studies the impact of the rotor-stator interaction in a single-stage progressing cavity pump on the flow rate and pressure. Specifically, we investigate the effect of the rotor movement on the sealings formed with deformable stators for various speeds and pressures. Sealings are reconstructed with the help of a geometric 3D model. We analyze the tangential and radial deviation of…
▽ More
This contribution studies the impact of the rotor-stator interaction in a single-stage progressing cavity pump on the flow rate and pressure. Specifically, we investigate the effect of the rotor movement on the sealings formed with deformable stators for various speeds and pressures. Sealings are reconstructed with the help of a geometric 3D model. We analyze the tangential and radial deviation of the rotor from its reference path and show that the radial deviation affects the flow rate, whereas the tangential deviation affects the pressure dynamics. The conjectures are confirmed with a laboratory test setup.
△ Less
Submitted 14 April, 2021;
originally announced April 2021.
-
Estimating load points of a motor-pump system using pressure and inverter drive data
Authors:
Yashar Kouhi,
Jens Müller,
Sebastian Leonow,
Martin Mönnigmann
Abstract:
We propose a novel method for the estimation of rotor position, speed, and torque of a motor-pump system consisting of a progressive cavity pump (PCP) driven by an induction motor which operates under V/f open-loop control. We compute the speed and rotor position of the PCP by applying a phase locked loop to the pressure signal at the pressure side of the pump. An extended Kalman filter is used to…
▽ More
We propose a novel method for the estimation of rotor position, speed, and torque of a motor-pump system consisting of a progressive cavity pump (PCP) driven by an induction motor which operates under V/f open-loop control. We compute the speed and rotor position of the PCP by applying a phase locked loop to the pressure signal at the pressure side of the pump. An extended Kalman filter is used to estimate the torque of the PCP based on the speed, effective value of the stator current of the induction motor and a nonlinear motor model. Furthermore, we derive a tractable condition under which the convergence of the observer is guaranteed. We use a laboratory experiment to verify our results.
△ Less
Submitted 8 October, 2020;
originally announced October 2020.
-
Accelerating MPC by online detection of state space sets with common optimal feedback laws
Authors:
Kai König,
Martin Mönnigmann
Abstract:
Model predictive control (MPC) samples a generally unknown and complicated feedback law point by point. The solution for the current state $x$ contains, however, more information than only the optimal signal $u$ for this particular state. In fact, it provides an optimal affine feedback law $x\rightarrow u(x)$ on a polytope $Π\subset \mathbb{R}^n$, i.e., on a full-dimensional state space set. It is…
▽ More
Model predictive control (MPC) samples a generally unknown and complicated feedback law point by point. The solution for the current state $x$ contains, however, more information than only the optimal signal $u$ for this particular state. In fact, it provides an optimal affine feedback law $x\rightarrow u(x)$ on a polytope $Π\subset \mathbb{R}^n$, i.e., on a full-dimensional state space set. It is an obvious idea to reuse this affine feedback law as long as possible. Reusing it on its polytope $Π$ is too conservative, however, because any $Π$ is a state space set with a common affine law $x\rightarrow (u_0^\prime (x), \dots, u_{N-1}^\prime (x))^\prime\in\mathbb{R}^{Nm}$ for the entire horizon $N$. We show a simple criterion exists for identifying the polytopes that have a common $x\rightarrow u_0(x)$, but may differ with respect to $u_1(x), \dots, u_{N-1}(x)$. Because this criterion is too computationally expensive for an online use, we introduce a simple heuristics for the fast construction of a subset of the polytopes of interest. Computational examples show (i) a considerable fraction of QPs can be avoided (10% to 40%) and (ii) the heuristics results in a reduction very close to the maximum one that could be achieved if the explicit solution was available. We stress the proposed approach is intended for use in online MPC and it does not require the explicit solution.
△ Less
Submitted 18 September, 2020;
originally announced September 2020.
-
Regional predictive control with suboptimally extended regions of validity
Authors:
Kai König,
Martin Mönnigmann
Abstract:
Model predictive control (MPC) is based on perpetually solving optimization problems. The solution of the optimization is usually interpreted as the optimal input for the current state. However, the solution of the optimization does not just provide an optimal input, but an entire optimal affine feedback law and a polytope on which this law is optimal. We recently proposed to use this feedback law…
▽ More
Model predictive control (MPC) is based on perpetually solving optimization problems. The solution of the optimization is usually interpreted as the optimal input for the current state. However, the solution of the optimization does not just provide an optimal input, but an entire optimal affine feedback law and a polytope on which this law is optimal. We recently proposed to use this feedback law as long as the system remains in its polytope. This can be interpreted as an event-based approach, where leaving the current polytope is the event that triggers the next optimization. This approach is especially appropriate for a networked control setting since the feedback laws and their polytopes can be evaluated with a low computational effort on lean local nodes. In this article the region of validity for a feedback law is extended. More precisely, the optimal polytopes are extended to nonlinearly bounded regions of validity resulting from the intersection of a feasibility and stability region. As a result, fewer quadratic programs need to be solved compared to the optimal approach. The new validity regions are still suitable for the evaluation on a lean local node. Moreover, the regions can be adjusted for a desired closed-loop performance.
△ Less
Submitted 18 September, 2020; v1 submitted 17 September, 2020;
originally announced September 2020.
-
Simulation studies on regional predictive control
Authors:
Kai König,
Martin Mönnigmann
Abstract:
We apply recently proposed regional predictive control approaches to several examples to evaluate their efficiency. We present numerical results as well as hardware-in-the-loop (HIL) results.
We apply recently proposed regional predictive control approaches to several examples to evaluate their efficiency. We present numerical results as well as hardware-in-the-loop (HIL) results.
△ Less
Submitted 17 September, 2020;
originally announced September 2020.
-
Reducing the computational effort of min-max model predictive control with regional feedback laws
Authors:
Kai König,
Martin Mönnigmann
Abstract:
Recently, a regional MPC approach has been proposed that exploits the piecewise affine structure of the optimal solution (without computing the entire explicit solution before). Here, regional refers to the idea of using the affine feedback law that is optimal in a vicinity of the current state of operation, and therefore provides the optimal input signal without requiring to solve a QP. In the pr…
▽ More
Recently, a regional MPC approach has been proposed that exploits the piecewise affine structure of the optimal solution (without computing the entire explicit solution before). Here, regional refers to the idea of using the affine feedback law that is optimal in a vicinity of the current state of operation, and therefore provides the optimal input signal without requiring to solve a QP. In the present paper, we apply the idea of regional MPC to min-max MPC problems. We show that the new robust approach can significantly reduce the number of QPs to be solved within min-max MPC resulting in a reduced overall computational effort. Moreover, we compare the performance of the new approach to an existing robust regional MPC approach using a numerical example with varying horizon. Finally, we provide a rule for choosing a suitable robust regional MPC approach based on the horizon.
△ Less
Submitted 12 May, 2021; v1 submitted 17 September, 2020;
originally announced September 2020.
-
Accelerated Nonlinear Model Predictive Control by Exploiting Saturation
Authors:
Raphael Dyrska,
Ruth Mitze,
Martin Mönnigmann
Abstract:
We present an approach for accelerating nonlinear model predictive control. If the current optimal input signal is saturated, also the optimal signals in subsequent time steps often are. We propose to use the open-loop optimal input signals whenever the first and some subsequent input signals are saturated. We only solve the next optimal control problem, when a non-saturated signal is encountered,…
▽ More
We present an approach for accelerating nonlinear model predictive control. If the current optimal input signal is saturated, also the optimal signals in subsequent time steps often are. We propose to use the open-loop optimal input signals whenever the first and some subsequent input signals are saturated. We only solve the next optimal control problem, when a non-saturated signal is encountered, or the end of the horizon is reached. In this way, we can save a significant number of NLPs to be solved while on the other hand keep the performance loss small. Furthermore, the NMPC is reactivated in time when it comes to controlling the system safely to its reference.
△ Less
Submitted 14 September, 2020;
originally announced September 2020.
-
State space sets with common optimal feedback laws for nonlinear MPC
Authors:
Ruth Mitze,
Raphael Dyrska,
Kai König,
Martin Mönnigmann
Abstract:
In model predictive control (MPC), an optimal control problem (OCP) is solved for the current state and the first input of the solution, the optimal feedback law, is applied to the system. This procedure requires to solve the OCP in every time step. Recently, a new approach was suggested for linear MPC. The parametric solution of a linear quadratic OCP is a piecewise-affine feedback law. The solut…
▽ More
In model predictive control (MPC), an optimal control problem (OCP) is solved for the current state and the first input of the solution, the optimal feedback law, is applied to the system. This procedure requires to solve the OCP in every time step. Recently, a new approach was suggested for linear MPC. The parametric solution of a linear quadratic OCP is a piecewise-affine feedback law. The solution at a point in state space provides an optimal feedback law and a domain on which this law is the optimal solution. As long as the system remains in the domain, the law can be reused and the calculation of an OCP is avoided. In some domains the optimal feedback laws are identical. By uniting the corresponding domains, bigger domains are achieved and the optimal feedback law can be reused more often. In the present paper, we investigate in how far this approach can be extended from linear to nonlinear MPC, we propose an algorithm and we illustrate the achieved savings with an example.
△ Less
Submitted 9 September, 2020;
originally announced September 2020.
-
Supervisory model predictive control for PV battery and heat pump system with phase change slurry thermal storage
Authors:
Yannik Löhr,
Daniel Wolf,
Clemens Pollerberg,
Alexander Hörsting,
Martin Mönnigmann
Abstract:
We present the design, implementation and experimental validation of a supervisory predictive control approach for an electrical heating system featuring a phase change slurry as heat storage and transfer medium. The controller optimizes the energy flows that are used as set points for the heat generation and energy distribution components. The optimization handles the thermal and electrical subsy…
▽ More
We present the design, implementation and experimental validation of a supervisory predictive control approach for an electrical heating system featuring a phase change slurry as heat storage and transfer medium. The controller optimizes the energy flows that are used as set points for the heat generation and energy distribution components. The optimization handles the thermal and electrical subsystems simultaneously and is able to switch between different objectives. We show the control can be implemented on low-cost embedded hardware and validate it with an experimental test bed comprising an installation of the complete heating system, including all hydraulic and all electrical components. Experimental results demonstrate the feasibility of both, a heat pump heating system with a phase change slurry, and the optimal control approach. The main control objectives, i.e., thermal comfort and maximum self-consumption of solar energy, can be met. In addition, the system and its controller provide a load shifting potential.
△ Less
Submitted 27 March, 2020; v1 submitted 26 March, 2020;
originally announced March 2020.
-
An observer for partially obstructed wood particles in industrial drying processes
Authors:
Marc Oliver Berner,
Viktor Scherer,
Martin Mönnigmann
Abstract:
In order for biomass drying processes to be efficient, it is crucial to achieve the target residual water content within a close margin, since more conservative drying would result in a waste of energy. A method for a reliable estimation of the water content is therefore of obvious importance. Ideally, such a method does not require any expensive sensors. We show reduced order models and extended…
▽ More
In order for biomass drying processes to be efficient, it is crucial to achieve the target residual water content within a close margin, since more conservative drying would result in a waste of energy. A method for a reliable estimation of the water content is therefore of obvious importance. Ideally, such a method does not require any expensive sensors. We show reduced order models and extended Kalman filters can be combined to reliably determine the water content and temperature of wood particles based on only surface temperature measurements. The proposed observer works reliably if measurements are only available for parts of a particle face. It can therefore still be applied if particle surfaces are partially obstructed, which is a prerequisite for use in industrial processes and units, such as rotary dryers. The extended Kalman filter uses a reduced order model that is obtained by applying proper orthogonal decomposition and Galerkin projection to coupled PDEs that model heat conduction and water diffusion in anisotropic particles. In contrast to the original PDE simulation model, the reduced model and the filter based on it are suitable for real time computations and monitoring.
△ Less
Submitted 27 March, 2020; v1 submitted 12 November, 2019;
originally announced November 2019.
-
Controllability analysis and optimal control of biomass drying with reduced order models
Authors:
Marc Oliver Berner,
Viktor Scherer,
Martin Mönnigmann
Abstract:
Complex industrial processes such as the drying of combustible biomass can be modeled with computational fluid dynamics simulations. Due to their complexity, it is not straightforward to use these models for the analysis of system properties or for solving optimal control problems. We show reduced order models can be derived and used for these purposes for industrial drying processes.
Complex industrial processes such as the drying of combustible biomass can be modeled with computational fluid dynamics simulations. Due to their complexity, it is not straightforward to use these models for the analysis of system properties or for solving optimal control problems. We show reduced order models can be derived and used for these purposes for industrial drying processes.
△ Less
Submitted 27 March, 2020; v1 submitted 11 November, 2019;
originally announced November 2019.
-
A dynamic programming approach to solving constrained linear-quadratic optimal control problems
Authors:
Ruth Mitze,
Martin Mönnigmann
Abstract:
The solution of a constrained linear-quadratic regulator problem is determined by the set of its optimal active sets. We propose an algorithm that constructs this set of active sets for a desired horizon N from that for horizon N-1. While it is not obvious how to extend the optimal feedback law itself for horizon N-1 to horizon N, a simple relation between the optimal active sets for two successiv…
▽ More
The solution of a constrained linear-quadratic regulator problem is determined by the set of its optimal active sets. We propose an algorithm that constructs this set of active sets for a desired horizon N from that for horizon N-1. While it is not obvious how to extend the optimal feedback law itself for horizon N-1 to horizon N, a simple relation between the optimal active sets for two successive horizon lengths exists. Specifically, every optimal active set for horizon N is a superset of an optimal active set for horizon N-1 if the constraints are ordered stage by stage. The stagewise treatment results in a favorable computational effort. In addition, it is easy to detect the solution of the current horizon is equal to the infinite-horizon solution, if such a finite horizon exists, with the proposed algorithm.
△ Less
Submitted 24 October, 2019;
originally announced October 2019.
-
A Complexity Analysis of Event-Triggered Model Predictive Control on Industrial Hardware
Authors:
Patrik Simon Berner,
Martin Mönnigmann
Abstract:
We implement a recently proposed event-triggered networked MPC approach on industrial hardware to analyze its practical relevance. There exist several alternatives for such an implementation that differ with respect to the distribution of computational load between local and central nodes, and with respect to network bandwidth requirements. These alternatives have been analyzed theoretically befor…
▽ More
We implement a recently proposed event-triggered networked MPC approach on industrial hardware to analyze its practical relevance. There exist several alternatives for such an implementation that differ with respect to the distribution of computational load between local and central nodes, and with respect to network bandwidth requirements. These alternatives have been analyzed theoretically before, but when implemented it becomes evident that their usefulness cannot be predicted based on theoretical considerations alone. It is the purpose of the present paper to account for both practical and theoretical aspects in determining which alternative is most appropriate for an implementation on industrial hardware. The smallest possible bandwidth is known to result for a variant in which only the active set of constraints is transmitted from the central to the local nodes. Since local nodes must determine the control law from the active set in this case, which requires matrix inversions, an unattractive computational cost results at first sight. Somewhat surprisingly, the computational cost scales practically linearly in the problem size when implemented. We confirm this result with a more detailed theoretical complexity analysis than given in previous papers. All results are illustrated with data obtained with an implementation on industrial hardware components.
△ Less
Submitted 10 May, 2019;
originally announced May 2019.
-
Normal Vectors on Modified Hopf Manifolds of Delay Differential Equations
Authors:
Jonas Otten,
Martin Mönnigmann
Abstract:
This document states the normal vector system for modified Hopf boundaries of delay differential systems with state and parameter dependent delays. Specifically, it states the proof for Proposition 1 in the paper entitled "Robust optimization of delay differential equations with state and parameter dependent delays" by the same authors [1].
This document states the normal vector system for modified Hopf boundaries of delay differential systems with state and parameter dependent delays. Specifically, it states the proof for Proposition 1 in the paper entitled "Robust optimization of delay differential equations with state and parameter dependent delays" by the same authors [1].
△ Less
Submitted 13 March, 2019;
originally announced March 2019.
-
A method for the optimization of nonlinear systems with delays that guarantees stability and robustness
Authors:
Jonas Otten,
Martin Mönnigmann
Abstract:
We present a method for the steady state optimization of nonlinear delay differential equations. The method ensures stability and robustness, where a system is called robust if it remains stable despite uncertain parameters. Essentially, we ensure stability of all steady states of the nonlinear system on the steady state manifold that results from the variation of the uncertain parameters. The unc…
▽ More
We present a method for the steady state optimization of nonlinear delay differential equations. The method ensures stability and robustness, where a system is called robust if it remains stable despite uncertain parameters. Essentially, we ensure stability of all steady states of the nonlinear system on the steady state manifold that results from the variation of the uncertain parameters. The uncertain parameters are characterized by finite intervals, which may be interpreted as error bars and therefore are of immediate practical relevance. Stability despite uncertain parameters can be guaranteed by enforcing a lower bound on the distance of the optimal steady state to submanifolds of saddle-node and Hopf bifurcations on the steady state manifold. We derive constraints that ensure this distance. The proposed method differs from previous ones in that stability and robustness are guaranteed with constraints instead of with the cost function. Because the cost function is not required to enforce stability and robustness, it can be used to state economic or similar goals, which is natural in applications. We illustrate the proposed method by optimizing a laser diode. The optimization finds a steady state of maximum intensity while guaranteeing asymptotic or exponential stability despite uncertain model parameters.
△ Less
Submitted 13 March, 2019;
originally announced March 2019.
-
On the structure of the set of active sets in constrained linear quadratic regulation
Authors:
Martin Mönnigmann
Abstract:
The constrained linear quadratic regulation problem is solved by a continuous piecewise affine function on a set of state space polytopes. It is an obvious question whether this solution can be built up iteratively by increasing the horizon, i.e., by extending the classical backward dynamic programming solution for the unconstrained case to the constrained case. Unfortunately, however, the piecewi…
▽ More
The constrained linear quadratic regulation problem is solved by a continuous piecewise affine function on a set of state space polytopes. It is an obvious question whether this solution can be built up iteratively by increasing the horizon, i.e., by extending the classical backward dynamic programming solution for the unconstrained case to the constrained case. Unfortunately, however, the piecewise affine solution for horizon N is in general not contained in the piecewise affine law for horizon N + 1. We show that a backward dynamic programming does, in contrast, result in a useful structure for the set of the active sets that defines the solution. Essentially, every active set for the problem with horizon N + 1 results from extending an active set for horizon N , if the constraints are ordered stage by stage. Consequently, the set for horizon N + 1 can be found by only considering the constraints of the additional stage. Furthermore, it is easy to detect which polytopes and affine pieces are invariant to increasing the horizon, and therefore persist in the limit N to infinity. Several other aspects of the structure of the set of active sets become evident if the active sets are represented by bit tuples. There exists, for example, a subset of special active sets that generates a positive invariant and persistent (i.e., horizon invariant) set around the origin. It is very simple to identify these special active sets, and the positive invariant and persistent region can be found without solving optimal control or auxiliary optimization problems. The paper briefly discusses the use of these results in model predictive control. Some opportunities for uses in computational methods are also briefly summarized.
△ Less
Submitted 12 April, 2019; v1 submitted 10 August, 2018;
originally announced August 2018.
-
Improved automatic computation of Hessian matrix spectral bounds
Authors:
Moritz Schulze Darup,
Martin Mönnigmann
Abstract:
This paper presents a fast and powerful method for the computation of eigenvalue bounds for Hessian matrices $\nabla^2 \varphi(x) $ of nonlinear functions $\varphi: U \subseteq R^n\rightarrow R$ on hyperrectangles $B \subset U$. The method is based on a recently proposed procedure for an efficient computation of spectral bounds using extended codelists. Both the previous approach and the one prese…
▽ More
This paper presents a fast and powerful method for the computation of eigenvalue bounds for Hessian matrices $\nabla^2 \varphi(x) $ of nonlinear functions $\varphi: U \subseteq R^n\rightarrow R$ on hyperrectangles $B \subset U$. The method is based on a recently proposed procedure for an efficient computation of spectral bounds using extended codelists. Both the previous approach and the one presented here substantially differ from established methods in that they do deliberately not use any interval matrices and thus result in a favorable numerical complexity of order $O(n)\,N(\varphi)$, where $N(\varphi)$ denotes the number of operations needed to evaluate $\varphi$ at a point in its domain. We improve the previous method by exploiting sparsity, which naturally arises in the underlying codelists.
△ Less
Submitted 22 July, 2015;
originally announced July 2015.
-
Simulation studies on online constraint removal with a Lyapunov function
Authors:
Michael Jost,
Gabriele Pannocchia,
Martin Mönnigmann
Abstract:
We apply a recently proposed method for the acceleration of model predictive control (MPC) to 36 MPC implementations, which result from combining six sample receding horizon control problems with six quadratic programming solvers. We implement each of the 36 system-solver-combinations both with and without constraint removal and compare computational times for statistically relevant numbers of run…
▽ More
We apply a recently proposed method for the acceleration of model predictive control (MPC) to 36 MPC implementations, which result from combining six sample receding horizon control problems with six quadratic programming solvers. We implement each of the 36 system-solver-combinations both with and without constraint removal and compare computational times for statistically relevant numbers of runs.
△ Less
Submitted 26 May, 2015; v1 submitted 4 November, 2014;
originally announced November 2014.
-
Robust optimization of periodically operated nonlinear uncertain processes
Authors:
Darya Kastsian,
Martin Mönnigmann
Abstract:
We present a method for determining optimal modes of operation for autonomously oscillating systems with uncertain parameters. In a typical application of the method, a nonlinear dynamical system is optimized with respect to an economic objective function with nonlinear programming methods, and stability is guaranteed for all points in a robustness region around the optimal point. The stability co…
▽ More
We present a method for determining optimal modes of operation for autonomously oscillating systems with uncertain parameters. In a typical application of the method, a nonlinear dynamical system is optimized with respect to an economic objective function with nonlinear programming methods, and stability is guaranteed for all points in a robustness region around the optimal point. The stability constraints are implemented by imposing a lower bound on the distance between the optimal point and all stability boundaries in its vicinity, where stability boundaries are described with notions from bifurcation theory. We derive the required constraints for a general class of periodically operated processes and show how these bounds can be integrated into standard nonlinear programming methods. We present results of the optimization of two chemical reaction systems for illustration.
△ Less
Submitted 17 August, 2013;
originally announced August 2013.
-
Efficient Computation of Spectral Bounds for Hessian Matrices on Hyperrectangles for Global Optimization
Authors:
Moritz Schulze Darup,
Martin Kastsian,
Stefan Mross,
Martin Mönnigmann
Abstract:
We compare two established and a new method for the calculation of spectral bounds for Hessian matrices on hyperrectangles by applying them to a large collection of 1522 objective and constraint functions extracted from benchmark global optimization problems. Both the tightness of the spectral bounds and the computational effort are assessed. Specifically, we compare eigenvalue bounds obtained wit…
▽ More
We compare two established and a new method for the calculation of spectral bounds for Hessian matrices on hyperrectangles by applying them to a large collection of 1522 objective and constraint functions extracted from benchmark global optimization problems. Both the tightness of the spectral bounds and the computational effort are assessed. Specifically, we compare eigenvalue bounds obtained with the interval variant of Gershgorin's circle criterion [2,6], Hertz and Rohn's [7,16] method for tight bounds of interval matrices, and a recently proposed Hessian matrix eigenvalue arithmetic [12], which deliberately avoids the computation of interval Hessians.
△ Less
Submitted 1 June, 2012;
originally announced June 2012.
-
Revised version of "Five-loop additive renormalization in the phi^4 theory and amplitude functions of the minimally renormalized specific heat in three dimensions"
Authors:
S. A. Larin,
M. Moennigmann,
M. Stroesser,
V. Dohm
Abstract:
We present an analytic five-loop calculation for the additive renormalization constant A(u,ε) and the associated renormalization-group function B(u) of the specific heat of the O(n) symmetric φ^4 theory within the minimal subtraction scheme. We show that this calculation does not require new five-loop integrations but can be performed on the basis of the previous five-loop calculation of the fou…
▽ More
We present an analytic five-loop calculation for the additive renormalization constant A(u,ε) and the associated renormalization-group function B(u) of the specific heat of the O(n) symmetric φ^4 theory within the minimal subtraction scheme. We show that this calculation does not require new five-loop integrations but can be performed on the basis of the previous five-loop calculation of the four-point vertex function combined with an appropriate identification of symmetry factors of vacuum diagrams. We also determine the amplitude function F+(u) of the specific heat in three dimensions for n=1,2,3 above T_c and F-(u) for n=1 below T_c up to five-loop order, without using the ε=4-d expansion. Accurate results are obtained from Borel resummations of B(u) for n=1,2,3 and of the amplitude functions for n=1. Previous conjectures regarding the smallness of the resummed higher-order contributions are confirmed. Combining our results for B(u) and F+(u) for n=1,2,3 with those of a recent three-loop calculation of F-(u) for general n in d=3 dimensions we calculate Borel resummed universal amplitude ratios A+/A- for n=1,2,3. Our result for A+/A- = 1.056 +/- 0.004 for n=2 is significantly more accurate than the previous result obtained from the εexpansion up to O(ε^2) and agrees well with the high-precision experimental result A+/A- = 1.054 +/- 0.001 for He(4) near the superfluid transition obtained from a recent experiment in space.
△ Less
Submitted 4 May, 1998;
originally announced May 1998.
-
Five-loop additive renormalization in the phi^4 theory and amplitude functions of the minimally renormalized specific heat in three dimensions
Authors:
S. A. Larin,
M. Moennigmann,
M. Stroesser,
V. Dohm
Abstract:
We present an analytic five-loop calculation for the additive renormalization constant A(u,epsilon) and the associated renormalization-group function B(u) of the specific heat of the O(n) symmetric phi^4 theory within the minimal subtraction scheme. We show that this calculation does not require new five-loop integrations but can be performed on the basis of the previous five-loop calculation of…
▽ More
We present an analytic five-loop calculation for the additive renormalization constant A(u,epsilon) and the associated renormalization-group function B(u) of the specific heat of the O(n) symmetric phi^4 theory within the minimal subtraction scheme. We show that this calculation does not require new five-loop integrations but can be performed on the basis of the previous five-loop calculation of the four-point vertex function combined with an appropriate identification of symmetry factors of vacuum diagrams. We also determine the amplitude functions of the specific heat in three dimensions for n=1,2,3 above T_c and for n=1 below T_c up to five-loop order. Accurate results are obtained from Borel resummations of B(u) for n=1,2,3 and of the amplitude functions for n=1. Previous conjectures regarding the smallness of the resummed higher-order contributions are confirmed. Borel resummed universal amplitude ratios A^+/A^- and a_c^+/a_c^- are calculated for n=1.
△ Less
Submitted 7 November, 1997;
originally announced November 1997.