Source author record

Martin Mönnigmann

Martin Mönnigmann 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

15works
5topics
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

15 published item(s)

preprint2025arXiv

Efficient Adjoint Petrov-Galerkin Reduced Order Models for fluid flows governed by the incompressible Navier-Stokes equations

This research paper investigates the Adjoint Petrov-Galerkin (APG) method for reduced order models (ROM) and fluid dynamics governed by the incompressible Navier-Stokes equations. The Adjoint Petrov-Galerkin ROM, derived using the Mori-Zwanzig formalism, demonstrates superior accuracy and stability compared to standard Galerkin ROMs. However, challenges arise due to the time invariance of the test basis vectors, resulting in high computational requirements. To address this, we introduce a new efficient Adjoint Petrov-Galerkin (eAPG) ROM formulation, extending its application to the incompressible Navier-Stokes equations by exploiting the polynomial structure inherent in these equations. The offline and online phases partition eliminates the need for repeated test basis vector evaluations. This improves computational efficiency in comparison to the general Adjoint Petrov-Galerkin ROM formulation. A novel approach to augmenting the memory length, a critical factor influencing the stability and accuracy of the APG-ROM, is introduced, employing a data-driven optimization. Numerical results for the 3D turbulent flow around a circular cylinder demonstrate the efficacy of the proposed approach. Error measures and computational cost evaluations, considering metrics such as floating point operations and simulation time, provide a comprehensive analysis.

preprint2020arXiv

Accelerated Nonlinear Model Predictive Control by Exploiting Saturation

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.

preprint2020arXiv

Accelerating MPC by online detection of state space sets with common optimal feedback laws

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.

preprint2020arXiv

An observer for partially obstructed wood particles in industrial drying processes

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.

preprint2020arXiv

Controllability analysis and optimal control of biomass drying with reduced order models

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.

preprint2020arXiv

Regional predictive control with suboptimally extended regions of validity

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.

preprint2020arXiv

State space sets with common optimal feedback laws for nonlinear MPC

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.

preprint2020arXiv

Supervisory model predictive control for PV battery and heat pump system with phase change slurry thermal storage

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.

preprint2019arXiv

A dynamic programming approach to solving constrained linear-quadratic optimal control problems

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.

preprint2019arXiv

On the structure of the set of active sets in constrained linear quadratic regulation

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.

preprint2015arXiv

Improved automatic computation of Hessian matrix spectral bounds

This paper presents a fast and powerful method for the computation of eigenvalue bounds for Hessian matrices $\nabla^2 φ(x) $ of nonlinear functions $φ: 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(φ)$, where $N(φ)$ denotes the number of operations needed to evaluate $φ$ at a point in its domain. We improve the previous method by exploiting sparsity, which naturally arises in the underlying codelists.

preprint2015arXiv

Simulation studies on online constraint removal with a Lyapunov function

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.

preprint2013arXiv

Robust optimization of periodically operated nonlinear uncertain processes

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.

preprint2012arXiv

Efficient Computation of Spectral Bounds for Hessian Matrices on Hyperrectangles for Global Optimization

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.