Sie haben Javascript deaktiviert!
Sie haben versucht eine Funktion zu nutzen, die nur mit Javascript möglich ist. Um sämtliche Funktionalitäten unserer Internetseite zu nutzen, aktivieren Sie bitte Javascript in Ihrem Browser.

Show image information

Journal Articles, Conference Proceedings & Preprints

Open list in Research Information System


Nüske, F., Peitz, S., Philipp, F., Schaller, M., & Worthmann, K. (n.d.). Finite-data error bounds for Koopman-based prediction and control. Journal of Nonlinear Science, 33, Article 14.

The Koopman operator has become an essential tool for data-driven approximation of dynamical (control) systems in recent years, e.g., via extended dynamic mode decomposition. Despite its popularity, convergence results and, in particular, error bounds are still quite scarce. In this paper, we derive probabilistic bounds for the approximation error and the prediction error depending on the number of training data points; for both ordinary and stochastic differential equations. Moreover, we extend our analysis to nonlinear control-affine systems using either ergodic trajectories or i.i.d. samples. Here, we exploit the linearity of the Koopman generator to obtain a bilinear system and, thus, circumvent the curse of dimensionality since we do not autonomize the system by augmenting the state by the control inputs. To the best of our knowledge, this is the first finite-data error analysis in the stochastic and/or control setting. Finally, we demonstrate the effectiveness of the proposed approach by comparing it with state-of-the-art techniques showing its superiority whenever state and control are coupled.


Gebken, B. (2022). Computation and analysis of Pareto critical sets in smooth and nonsmooth multiobjective optimization.

Multiobjective optimization is concerned with the simultaneous optimization of multiple scalar-valued functions. In this case, a point is optimal if there is no other point that is at least as good in all objectives and better in at least one objective. A necessary condition for optimality can be derived based on first-order information of the objectives. The set of points that satisfy this necessary condition is called the Pareto critical set. This thesis presents new results about Pareto critical sets for smooth and nonsmooth multiobjective optimization problems, both in terms of their efficient computation and structural properties. In the smooth case they are computed via a continuation method and in the nonsmooth case via a descent method. Afterwards, the structure of the boundary of the Pareto critical set is analyzed, which consists of Pareto critical sets of smaller subproblems. Finally, inverse problems are considered, where a data set is given and an objective vector is sought for which the data points are critical.

Dellnitz, M., Hüllermeier, E., Lücke, M., Ober-Blöbaum, S., Offen, C., Peitz, S., & Pfannschmidt, K. (n.d.). Efficient time stepping for numerical integration using reinforcement  learning. In SIAM Journal on Scientific Computing.

Many problems in science and engineering require an efficient numerical approximation of integrals or solutions to differential equations. For systems with rapidly changing dynamics, an equidistant discretization is often inadvisable as it either results in prohibitively large errors or computational effort. To this end, adaptive schemes, such as solvers based on Runge–Kutta pairs, have been developed which adapt the step size based on local error estimations at each step. While the classical schemes apply very generally and are highly efficient on regular systems, they can behave sub-optimal when an inefficient step rejection mechanism is triggered by structurally complex systems such as chaotic systems. To overcome these issues, we propose a method to tailor numerical schemes to the problem class at hand. This is achieved by combining simple, classical quadrature rules or ODE solvers with data-driven time-stepping controllers. Compared with learning solution operators to ODEs directly, it generalises better to unseen initial data as our approach employs classical numerical schemes as base methods. At the same time it can make use of identified structures of a problem class and, therefore, outperforms state-of-the-art adaptive schemes. Several examples demonstrate superior efficiency. Source code is available at

Peitz, S., Dellnitz, M., & Bannenberg, S. (2022). Efficient Virtual Design and Testing of Autonomous Vehicles. In H. G. Bock, K.-H. Küfer, P. Maas, A. Milde, & V. Schulz (Eds.), German Success Stories in Industrial Mathematics (Vol. 35). Springer International Publishing.

With the ever increasing capabilities of sensors and controllers, autonomous driving is quickly becoming a reality. This disruptive change in the automotive industry poses major challenges for manufacturers as well as suppliers as entirely new design and testing strategies have to be developed to remain competitive. Most importantly, the complexity of autonomously driving vehicles in a complex, uncertain, and safety-critical environment requires new testing procedures to cover the almost infinite range of potential scenarios.

Sonntag, K., & Peitz, S. (2022). Fast Multiobjective Gradient Methods with Nesterov Acceleration via  Inertial Gradient-like Systems. In arXiv:2207.12707.

We derive efficient algorithms to compute weakly Pareto optimal solutions for smooth, convex and unconstrained multiobjective optimization problems in general Hilbert spaces. To this end, we define a novel inertial gradient-like dynamical system in the multiobjective setting, whose trajectories converge weakly to Pareto optimal solutions. Discretization of this system yields an inertial multiobjective algorithm which generates sequences that converge weakly to Pareto optimal solutions. We employ Nesterov acceleration to define an algorithm with an improved convergence rate compared to the plain multiobjective steepest descent method (Algorithm 1). A further improvement in terms of efficiency is achieved by avoiding the solution of a quadratic subproblem to compute a common step direction for all objective functions, which is usually required in first order methods. Using a different discretization of our inertial gradient-like dynamical system, we obtain an accelerated multiobjective gradient method that does not require the solution of a subproblem in each step (Algorithm 2). While this algorithm does not converge in general, it yields good results on test problems while being faster than standard steepest descent by two to three orders of magnitude.

Klus, S., Nüske, F., & Peitz, S. (2022). Koopman analysis of quantum systems. Journal of Physics A: Mathematical and Theoretical, 55(31), 314002.

Koopman operator theory has been successfully applied to problems from various research areas such as fluid dynamics, molecular dynamics, climate science, engineering, and biology. Applications include detecting metastable or coherent sets, coarse-graining, system identification, and control. There is an intricate connection between dynamical systems driven by stochastic differential equations and quantum mechanics. In this paper, we compare the ground-state transformation and Nelson's stochastic mechanics and demonstrate how data-driven methods developed for the approximation of the Koopman operator can be used to analyze quantum physics problems. Moreover, we exploit the relationship between Schrödinger operators and stochastic control problems to show that modern data-driven methods for stochastic control can be used to solve the stationary or imaginary-time Schrödinger equation. Our findings open up a new avenue towards solving Schrödinger's equation using recently developed tools from data science.

Berkemeier, M. B., & Peitz, S. (2022). Multi-Objective Trust-Region Filter Method for Nonlinear Constraints using Inexact Gradients. In arXiv:2208.12094.

In this article, we build on previous work to present an optimization algorithm for nonlinearly constrained multi-objective optimization problems. The algorithm combines a surrogate-assisted derivative-free trust-region approach with the filter method known from single-objective optimization. Instead of the true objective and constraint functions, so-called fully linear models are employed and we show how to deal with the gradient inexactness in the composite step setting, adapted from single-objective optimization as well. Under standard assumptions, we prove convergence of a subset of iterates to a quasi-stationary point and if constraint qualifications hold, then the limit point is also a KKT-point of the multi-objective problem.

Gebken, B., Bieker, K., & Peitz, S. (2022). On the structure of regularization paths for piecewise differentiable regularization terms. Journal of Global Optimization.

Regularization is used in many different areas of optimization when solutions are sought which not only minimize a given function, but also possess a certain degree of regularity. Popular applications are image denoising, sparse regression and machine learning. Since the choice of the regularization parameter is crucial but often difficult, path-following methods are used to approximate the entire regularization path, i.e., the set of all possible solutions for all regularization parameters. Due to their nature, the development of these methods requires structural results about the regularization path. The goal of this article is to derive these results for the case of a smooth objective function which is penalized by a piecewise differentiable regularization term. We do this by treating regularization as a multiobjective optimization problem. Our results suggest that even in this general case, the regularization path is piecewise smooth. Moreover, our theory allows for a classification of the nonsmooth features that occur in between smooth parts. This is demonstrated in two applications, namely support-vector machines and exact penalty methods.

Bieker, K., Gebken, B., & Peitz, S. (2022). On the Treatment of Optimization Problems with L1 Penalty Terms via Multiobjective Continuation. IEEE Transactions on Pattern Analysis and Machine Intelligence, 44(11), 7797–7808.

We present a novel algorithm that allows us to gain detailed insight into the effects of sparsity in linear and nonlinear optimization, which is of great importance in many scientific areas such as image and signal processing, medical imaging, compressed sensing, and machine learning (e.g., for the training of neural networks). Sparsity is an important feature to ensure robustness against noisy data, but also to find models that are interpretable and easy to analyze due to the small number of relevant terms. It is common practice to enforce sparsity by adding the ℓ1-norm as a weighted penalty term. In order to gain a better understanding and to allow for an informed model selection, we directly solve the corresponding multiobjective optimization problem (MOP) that arises when we minimize the main objective and the ℓ1-norm simultaneously. As this MOP is in general non-convex for nonlinear objectives, the weighting method will fail to provide all optimal compromises. To avoid this issue, we present a continuation method which is specifically tailored to MOPs with two objective functions one of which is the ℓ1-norm. Our method can be seen as a generalization of well-known homotopy methods for linear regression problems to the nonlinear case. Several numerical examples - including neural network training - demonstrate our theoretical findings and the additional insight that can be gained by this multiobjective approach.

Peitz, S., & Bieker, K. (n.d.). On the Universal Transformation of Data-Driven Models to Control Systems. Automatica.

As in almost every other branch of science, the major advances in data science and machine learning have also resulted in significant improvements regarding the modeling and simulation of nonlinear dynamical systems. It is nowadays possible to make accurate medium to long-term predictions of highly complex systems such as the weather, the dynamics within a nuclear fusion reactor, of disease models or the stock market in a very efficient manner. In many cases, predictive methods are advertised to ultimately be useful for control, as the control of high-dimensional nonlinear systems is an engineering grand challenge with huge potential in areas such as clean and efficient energy production, or the development of advanced medical devices. However, the question of how to use a predictive model for control is often left unanswered due to the associated challenges, namely a significantly higher system complexity, the requirement of much larger data sets and an increased and often problem-specific modeling effort. To solve these issues, we present a universal framework (which we call QuaSiModO: Quantization-Simulation-Modeling-Optimization) to transform arbitrary predictive models into control systems and use them for feedback control. The advantages of our approach are a linear increase in data requirements with respect to the control dimension, performance guarantees that rely exclusively on the accuracy of the predictive model, and only little prior knowledge requirements in control theory to solve complex control problems. In particular the latter point is of key importance to enable a large number of researchers and practitioners to exploit the ever increasing capabilities of predictive models for control in a straight-forward and systematic fashion.

Banholzer, S., Gebken, B., Dellnitz, M., Peitz, S., & Volkwein, S. (2022). ROM-Based Multiobjective Optimization of Elliptic PDEs via Numerical Continuation. In H. Michael, H. Roland, K. Christian, U. Michael, & U. Stefan (Eds.), Non-Smooth and Complementarity-Based Distributed Parameter Systems (pp. 43–76). Springer.

Multiobjective optimization plays an increasingly important role in modern applications, where several objectives are often of equal importance. The task in multiobjective optimization and multiobjective optimal control is therefore to compute the set of optimal compromises (the Pareto set) between the conflicting objectives. Since the Pareto set generally consists of an infinite number of solutions, the computational effort can quickly become challenging which is particularly problematic when the objectives are costly to evaluate as is the case for models governed by partial differential equations (PDEs). To decrease the numerical effort to an affordable amount, surrogate models can be used to replace the expensive PDE evaluations. Existing multiobjective optimization methods using model reduction are limited either to low parameter dimensions or to few (ideally two) objectives. In this article, we present a combination of the reduced basis model reduction method with a continuation approach using inexact gradients. The resulting approach can handle an arbitrary number of objectives while yielding a significant reduction in computing time.


Gebken, B., & Peitz, S. (2021). An efficient descent method for locally Lipschitz multiobjective optimization problems. Journal of Optimization Theory and Applications, 188, 696–723.

In this article, we present an efficient descent method for locally Lipschitz continuous multiobjective optimization problems (MOPs). The method is realized by combining a theoretical result regarding the computation of descent directions for nonsmooth MOPs with a practical method to approximate the subdifferentials of the objective functions. We show convergence to points which satisfy a necessary condition for Pareto optimality. Using a set of test problems, we compare our method to the multiobjective proximal bundle method by M\"akel\"a. The results indicate that our method is competitive while being easier to implement. While the number of objective function evaluations is larger, the overall number of subgradient evaluations is lower. Finally, we show that our method can be combined with a subdivision algorithm to compute entire Pareto sets of nonsmooth MOPs.

Berkemeier, M. B., & Peitz, S. (2021). Derivative-Free Multiobjective Trust Region Descent Method Using Radial  Basis Function Surrogate Models. Mathematical and Computational Applications, 26(2).

We present a flexible trust region descend algorithm for unconstrained and convexly constrained multiobjective optimization problems. It is targeted at heterogeneous and expensive problems, i.e., problems that have at least one objective function that is computationally expensive. The method is derivative-free in the sense that neither need derivative information be available for the expensive objectives nor are gradients approximated using repeated function evaluations as is the case in finite-difference methods. Instead, a multiobjective trust region approach is used that works similarly to its well-known scalar pendants. Local surrogate models constructed from evaluation data of the true objective functions are employed to compute possible descent directions. In contrast to existing multiobjective trust region algorithms, these surrogates are not polynomial but carefully constructed radial basis function networks. This has the important advantage that the number of data points scales linearly with the parameter space dimension. The local models qualify as fully linear and the corresponding general scalar framework is adapted for problems with multiple objectives. Convergence to Pareto critical points is proven and numerical examples illustrate our findings.

Goelz, C., Mora, K., Stroehlein, J. K., Haase, F. K., Dellnitz, M., Reinsberger, C., & Vieluf, S. (2021). Electrophysiological signatures of dedifferentiation differ between fit and less fit older adults. Cognitive Neurodynamics.

Ober-Blöbaum, S., & Peitz, S. (2021). Explicit multiobjective model predictive control for nonlinear systems  with symmetries. International Journal of Robust and Nonlinear Control, 31(2), 380–403.

Model predictive control is a prominent approach to construct a feedback control loop for dynamical systems. Due to real-time constraints, the major challenge in MPC is to solve model-based optimal control problems in a very short amount of time. For linear-quadratic problems, Bemporad et al. have proposed an explicit formulation where the underlying optimization problems are solved a priori in an offline phase. In this article, we present an extension of this concept in two significant ways. We consider nonlinear problems and - more importantly - problems with multiple conflicting objective functions. In the offline phase, we build a library of Pareto optimal solutions from which we then obtain a valid compromise solution in the online phase according to a decision maker's preference. Since the standard multi-parametric programming approach is no longer valid in this situation, we instead use interpolation between different entries of the library. To reduce the number of problems that have to be solved in the offline phase, we exploit symmetries in the dynamical system and the corresponding multiobjective optimal control problem. The results are verified using two different examples from autonomous driving.

Gebken, B., & Peitz, S. (2021). Inverse multiobjective optimization: Inferring decision criteria from data. Journal of Global Optimization, 80, 3–29.

It is a challenging task to identify the objectives on which a certain decision was based, in particular if several, potentially conflicting criteria are equally important and a continuous set of optimal compromise decisions exists. This task can be understood as the inverse problem of multiobjective optimization, where the goal is to find the objective function vector of a given Pareto set. To this end, we present a method to construct the objective function vector of an unconstrained multiobjective optimization problem (MOP) such that the Pareto critical set contains a given set of data points with prescribed KKT multipliers. If such an MOP can not be found, then the method instead produces an MOP whose Pareto critical set is at least close to the data points. The key idea is to consider the objective function vector in the multiobjective KKT conditions as variable and then search for the objectives that minimize the Euclidean norm of the resulting system of equations. By expressing the objectives in a finite-dimensional basis, we transform this problem into a homogeneous, linear system of equations that can be solved efficiently. Potential applications of this approach include the identification of objectives (both from clean and noisy data) and the construction of surrogate models for expensive MOPs.

Nüske, F., Koltai, P., Boninsegna, L., & Clementi, C. (2021). Spectral Properties of Effective Dynamics from Conditional Expectations. Entropy.

<jats:p>The reduction of high-dimensional systems to effective models on a smaller set of variables is an essential task in many areas of science. For stochastic dynamics governed by diffusion processes, a general procedure to find effective equations is the conditioning approach. In this paper, we are interested in the spectrum of the generator of the resulting effective dynamics, and how it compares to the spectrum of the full generator. We prove a new relative error bound in terms of the eigenfunction approximation error for reversible systems. We also present numerical examples indicating that, if Kramers–Moyal (KM) type approximations are used to compute the spectrum of the reduced generator, it seems largely insensitive to the time window used for the KM estimators. We analyze the implications of these observations for systems driven by underdamped Langevin dynamics, and show how meaningful effective dynamics can be defined in this setting.</jats:p>

Klus, S., Gelß, P., Nüske, F., & Noé, F. (2021). Symmetric and antisymmetric kernels for machine learning problems in quantum physics and chemistry. Machine Learning: Science and Technology, Article 045016.

Nüske, F., Gelß, P., Klus, S., & Clementi, C. (2021). Tensor-based computation of metastable and coherent sets. Physica D: Nonlinear Phenomena, Article 133018.

Gerlach, R. (2021). The Computation and Analysis of Invariant Sets of Infinite-Dimensional Systems.

One central aspect in the study of dynamical systems is the analysis of its invariant sets such as the global attractor and (un)stable manifolds. In particular, when the underlying system depends on a parameter it is crucial to efficiently track those set with respect to this parameter. For the computation of invariant sets we rely on numerical algorithms for their approximation but typically those tools can only be applied to finite-dimensional dynamical systems. Thus, in thesis we present a numerical framework for the global dynamical analysis of infinite-dimensional systems. We will use embedding techniques for the definition of the core dynamical system (CDS) which is a dynamically equivalent finite-dimensional system. The CDS is then used for the approximation of related embedded invariant sets, i.e, one-to-one images, by set-oriented numerical methods. For the construction of the CDS it is crucial to choose an appropriate observation map and to design its corresponding inverse. Therefore, we will present suitable numerical realizations of the CDS for DDEs and PDEs. For a subsequent geometric analysis of the embedded invariant set we will consider a manifold learning technique called diffusion maps which reveals its intrinsic geometry and estimates its dimension. Finally, we apply our develop numerical tools on some well-known infinite-dimensional dynamical systems such as the Mackey-Glass equation, the Kuramoto-Sivashinsky equation and the Navier-Stokes equation.


Gerlach, R., Ziessler, A., Eckhardt, B., & Dellnitz, M. (2020). A Set-Oriented Path Following Method for the Approximation of Parameter Dependent Attractors. SIAM Journal on Applied Dynamical Systems, 705–723.

In this work we present a set-oriented path following method for the computation of relative global attractors of parameter-dependent dynamical systems. We start with an initial approximation of the relative global attractor for a fixed parameter λ0 computed by a set-oriented subdivision method. By using previously obtained approximations of the parameter-dependent relative global attractor we can track it with respect to a one-dimensional parameter λ > λ0 without restarting the whole subdivision procedure. We illustrate the feasibility of the set-oriented path following method by exploring the dynamics in low-dimensional models for shear flows during the transition to turbulence and of large-scale atmospheric regime changes .

Klus, S., Nüske, F., Peitz, S., Niemann, J.-H., Clementi, C., & Schütte, C. (2020). Data-driven approximation of the Koopman generator: Model reduction, system identification, and control. Physica D: Nonlinear Phenomena, 406.

We derive a data-driven method for the approximation of the Koopman generator called gEDMD, which can be regarded as a straightforward extension of EDMD (extended dynamic mode decomposition). This approach is applicable to deterministic and stochastic dynamical systems. It can be used for computing eigenvalues, eigenfunctions, and modes of the generator and for system identification. In addition to learning the governing equations of deterministic systems, which then reduces to SINDy (sparse identification of nonlinear dynamics), it is possible to identify the drift and diffusion terms of stochastic differential equations from data. Moreover, we apply gEDMD to derive coarse-grained models of high-dimensional systems, and also to determine efficient model predictive control strategies. We highlight relationships with other methods and demonstrate the efficacy of the proposed methods using several guiding examples and prototypical molecular dynamics problems.

Peitz, S., Otto, S. E., & Rowley, C. W. (2020). Data-Driven Model Predictive Control using Interpolated Koopman  Generators. SIAM Journal on Applied Dynamical Systems, 19(3), 2162–2193.

In recent years, the success of the Koopman operator in dynamical systems analysis has also fueled the development of Koopman operator-based control frameworks. In order to preserve the relatively low data requirements for an approximation via Dynamic Mode Decomposition, a quantization approach was recently proposed in [Peitz & Klus, Automatica 106, 2019]. This way, control of nonlinear dynamical systems can be realized by means of switched systems techniques, using only a finite set of autonomous Koopman operator-based reduced models. These individual systems can be approximated very efficiently from data. The main idea is to transform a control system into a set of autonomous systems for which the optimal switching sequence has to be computed. In this article, we extend these results to continuous control inputs using relaxation. This way, we combine the advantages of the data efficiency of approximating a finite set of autonomous systems with continuous controls. We show that when using the Koopman generator, this relaxation --- realized by linear interpolation between two operators --- does not introduce any error for control affine systems. This allows us to control high-dimensional nonlinear systems using bilinear, low-dimensional surrogate models. The efficiency of the proposed approach is demonstrated using several examples with increasing complexity, from the Duffing oscillator to the chaotic fluidic pinball.

Bieker, K., Peitz, S., Brunton, S. L., Kutz, J. N., & Dellnitz, M. (2020). Deep model predictive flow control with limited sensor data and online learning. Theoretical and Computational Fluid Dynamics, 34, 577–591.

The control of complex systems is of critical importance in many branches of science, engineering, and industry, many of which are governed by nonlinear partial differential equations. Controlling an unsteady fluid flow is particularly important, as flow control is a key enabler for technologies in energy (e.g., wind, tidal, and combustion), transportation (e.g., planes, trains, and automobiles), security (e.g., tracking airborne contamination), and health (e.g., artificial hearts and artificial respiration). However, the high-dimensional, nonlinear, and multi-scale dynamics make real-time feedback control infeasible. Fortunately, these high- dimensional systems exhibit dominant, low-dimensional patterns of activity that can be exploited for effective control in the sense that knowledge of the entire state of a system is not required. Advances in machine learning have the potential to revolutionize flow control given its ability to extract principled, low-rank feature spaces characterizing such complex systems.We present a novel deep learning modelpredictive control framework that exploits low-rank features of the flow in order to achieve considerable improvements to control performance. Instead of predicting the entire fluid state, we use a recurrent neural network (RNN) to accurately predict the control relevant quantities of the system, which are then embedded into an MPC framework to construct a feedback loop. In order to lower the data requirements and to improve the prediction accuracy and thus the control performance, incoming sensor data are used to update the RNN online. The results are validated using varying fluid flow examples of increasing complexity.

Hernández Castellanos, C. I., Ober-Blöbaum, S., & Peitz, S. (2020). Explicit Multi-objective Model Predictive Control for Nonlinear Systems  Under Uncertainty. International Journal of Robust and Nonlinear Control, 30(17), 7593–7618.

In real-world problems, uncertainties (e.g., errors in the measurement, precision errors) often lead to poor performance of numerical algorithms when not explicitly taken into account. This is also the case for control problems, where optimal solutions can degrade in quality or even become infeasible. Thus, there is the need to design methods that can handle uncertainty. In this work, we consider nonlinear multi-objective optimal control problems with uncertainty on the initial conditions, and in particular their incorporation into a feedback loop via model predictive control (MPC). In multi-objective optimal control, an optimal compromise between multiple conflicting criteria has to be found. For such problems, not much has been reported in terms of uncertainties. To address this problem class, we design an offline/online framework to compute an approximation of efficient control strategies. This approach is closely related to explicit MPC for nonlinear systems, where the potentially expensive optimization problem is solved in an offline phase in order to enable fast solutions in the online phase. In order to reduce the numerical cost of the offline phase, we exploit symmetries in the control problems. Furthermore, in order to ensure optimality of the solutions, we include an additional online optimization step, which is considerably cheaper than the original multi-objective optimization problem. We test our framework on a car maneuvering problem where safety and speed are the objectives. The multi-objective framework allows for online adaptations of the desired objective. Alternatively, an automatic scalarizing procedure yields very efficient feedback controls. Our results show that the method is capable of designing driving strategies that deal better with uncertainties in the initial conditions, which translates into potentially safer and faster driving strategies.

Peitz, S., & Klus, S. (2020). Feedback Control of Nonlinear PDEs Using Data-Efficient Reduced Order Models Based on the Koopman Operator. In Lecture Notes in Control and Information Sciences (Vol. 484, pp. 257–282). Cham: Springer.

In the development of model predictive controllers for PDE-constrained problems, the use of reduced order models is essential to enable real-time applicability. Besides local linearization approaches, proper orthogonal decomposition (POD) has been most widely used in the past in order to derive such models. Due to the huge advances concerning both theory as well as the numerical approximation, a very promising alternative based on the Koopman operator has recently emerged. In this chapter, we present two control strategies for model predictive control of nonlinear PDEs using data-efficient approximations of the Koopman operator. In the first one, the dynamic control system is replaced by a small number of autonomous systems with different yet constant inputs. The control problem is consequently transformed into a switching problem. In the second approach, a bilinear surrogate model is obtained via a convex combination of these autonomous systems. Using a recent convergence result for extended dynamic mode decomposition (EDMD), convergence of the reduced objective function can be shown. We study the properties of these two strategies with respect to solution quality, data requirements, and complexity of the resulting optimization problem using the 1-dimensional Burgers equation and the 2-dimensional Navier–Stokes equations as examples. Finally, an extension for online adaptivity is presented.

Klus, S., Nüske, F., & Hamzi, B. (2020). Kernel-Based Approximation of the Koopman Generator and Schrödinger Operator. Entropy.

<jats:p>Many dimensionality and model reduction techniques rely on estimating dominant eigenfunctions of associated dynamical operators from data. Important examples include the Koopman operator and its generator, but also the Schrödinger operator. We propose a kernel-based method for the approximation of differential operators in reproducing kernel Hilbert spaces and show how eigenfunctions can be estimated by solving auxiliary matrix eigenvalue problems. The resulting algorithms are applied to molecular dynamics and quantum chemistry examples. Furthermore, we exploit that, under certain conditions, the Schrödinger operator can be transformed into a Kolmogorov backward operator corresponding to a drift-diffusion process and vice versa. This allows us to apply methods developed for the analysis of high-dimensional stochastic differential equations to quantum mechanical systems.</jats:p>

Schütze, O., Cuate, O., Martín, A., Peitz, S., & Dellnitz, M. (2020). Pareto Explorer: a global/local exploration tool for many-objective optimization problems. Engineering Optimization, 52(5), 832–855.

Multi-objective optimization is an active field of research that has many applications. Owing to its success and because decision-making processes are becoming more and more complex, there is a recent trend for incorporating many objectives into such problems. The challenge with such problems, however, is that the dimensions of the solution sets—the so-called Pareto sets and fronts—grow with the number of objectives. It is thus no longer possible to compute or to approximate the entire solution set of a given problem that contains many (e.g. more than three) objectives. On the other hand, the computation of single solutions (e.g. via scalarization methods) leads to unsatisfying results in many cases, even if user preferences are incorporated. In this article, the Pareto Explorer tool is presented—a global/local exploration tool for the treatment of many-objective optimization problems (MaOPs). In the first step, a solution of the problem is computed via a global search algorithm that ideally already includes user preferences. In the second step, a local search along the Pareto set/front of the given MaOP is performed in user specified directions. For this, several continuation-like procedures are proposed that can incorporate preferences defined in decision, objective, or in weight space. The applicability and usefulness of Pareto Explorer is demonstrated on benchmark problems as well as on an application from industrial laundry design.

Flaßkamp, K., Ober-Blöbaum, S., & Peitz, S. (2020). Symmetry in Optimal Control: A Multiobjective Model Predictive Control Approach. In O. Junge, O. Schütze, G. Froyland, S. Ober-Blöbaum, & K. Padberg-Gehle (Eds.), Advances in Dynamics, Optimization and Computation. Cham: Springer.

Many dynamical systems possess symmetries, e.g. rotational and translational invariances of mechanical systems. These can be beneficially exploited in the design of numerical optimal control methods. We present a model predictive control scheme which is based on a library of precomputed motion primitives. The primitives are equivalence classes w.r.t. the symmetry of the optimal control problems. Trim primitives as relative equilibria w.r.t. this symmetry, play a crucial role in the algorithm. The approach is illustrated using an academic mobile robot example.

Gerlach, R., & Ziessler, A. (2020). The Approximation of Invariant Sets in Infinite Dimensional Dynamical Systems. In O. Junge, O. Schütze, S. Ober-Blöbaum, & K. Padberg-Gehle (Eds.), Advances in Dynamics, Optimization and Computation (Vol. 304, pp. 55–85). Springer.

In this work we review the novel framework for the computation of finite dimensional invariant sets of infinite dimensional dynamical systems developed in [6] and [36]. By utilizing results on embedding techniques for infinite dimensional systems we extend a classical subdivision scheme [8] as well as a continuation algorithm [7] for the computation of attractors and invariant manifolds of finite dimensional systems to the infinite dimensional case. We show how to implement this approach for the analysis of delay differential equations and partial differential equations and illustrate the feasibility of our implementation by computing the attractor of the Mackey-Glass equation and the unstable manifold of the one-dimensional Kuramoto-Sivashinsky equation.


Nüske, F., Boninsegna, L., & Clementi, C. (2019). Coarse-graining molecular systems by spectral matching. The Journal of Chemical Physics.

Sahai, T., Ziessler, A., Klus, S., & Dellnitz, M. (2019). Continuous relaxations for the traveling salesman problem. Nonlinear Dynamics.

Hanke, S., Peitz, S., Wallscheid, O., Böcker, J., & Dellnitz, M. (2019). Finite-Control-Set Model Predictive Control for a Permanent Magnet Synchronous Motor Application with Online Least Squares System Identification. In 2019 IEEE International Symposium on Predictive Control of Electrical Drives and Power Electronics (PRECEDE).

In comparison to classical control approaches in the field of electrical drives like the field-oriented control (FOC), model predictive control (MPC) approaches are able to provide a higher control performance. This refers to shorter settling times, lower overshoots, and a better decoupling of control variables in case of multi-variable controls. However, this can only be achieved if the used prediction model covers the actual behavior of the plant sufficiently well. In case of model deviations, the performance utilizing MPC remains below its potential. This results in effects like increased current ripple or steady state setpoint deviations. In order to achieve a high control performance, it is therefore necessary to adapt the model to the real plant behavior. When using an online system identification, a less accurate model is sufficient for commissioning of the drive system. In this paper, the combination of a finite-control-set MPC (FCS-MPC) with a system identification is proposed. The method does not require high-frequency signal injection, but uses the measured values already required for the FCS-MPC. An evaluation of the least squares-based identification on a laboratory test bench showed that the model accuracy and thus the control performance could be improved by an online update of the prediction models.

Peitz, S., & Klus, S. (2019). Koopman operator-based model reduction for switched-system control of PDEs. Automatica, 106, 184–191.

We present a new framework for optimal and feedback control of PDEs using Koopman operator-based reduced order models (K-ROMs). The Koopman operator is a linear but infinite-dimensional operator which describes the dynamics of observables. A numerical approximation of the Koopman operator therefore yields a linear system for the observation of an autonomous dynamical system. In our approach, by introducing a finite number of constant controls, the dynamic control system is transformed into a set of autonomous systems and the corresponding optimal control problem into a switching time optimization problem. This allows us to replace each of these systems by a K-ROM which can be solved orders of magnitude faster. By this approach, a nonlinear infinite-dimensional control problem is transformed into a low-dimensional linear problem. Using a recent convergence result for the numerical approximation via Extended Dynamic Mode Decomposition (EDMD), we show that the value of the K-ROM based objective function converges in measure to the value of the full objective function. To illustrate the results, we consider the 1D Burgers equation and the 2D Navier–Stokes equations. The numerical experiments show remarkable performance concerning both solution times and accuracy.

Dellnitz, M., Gebken, B., Gerlach, R., & Klus, S. (2019). On the equivariance properties of self-adjoint matrices. Dynamical Systems, 1–19.

We investigate self-adjoint matrices A∈Rn,n with respect to their equivariance properties. We show in particular that a matrix is self-adjoint if and only if it is equivariant with respect to the action of a group Γ2(A)⊂O(n) which is isomorphic to ⊗nk=1Z2. If the self-adjoint matrix possesses multiple eigenvalues – this may, for instance, be induced by symmetry properties of an underlying dynamical system – then A is even equivariant with respect to the action of a group Γ(A)≃∏ki=1O(mi) where m1,…,mk are the multiplicities of the eigenvalues λ1,…,λk of A. We discuss implications of this result for equivariant bifurcation problems, and we briefly address further applications for the Procrustes problem, graph symmetries and Taylor expansions.

Gebken, B., Peitz, S., & Dellnitz, M. (2019). On the hierarchical structure of Pareto critical sets. Journal of Global Optimization, 73(4), 891–913.

In this article we show that the boundary of the Pareto critical set of an unconstrained multiobjective optimization problem (MOP) consists of Pareto critical points of subproblems where only a subset of the set of objective functions is taken into account. If the Pareto critical set is completely described by its boundary (e.g., if we have more objective functions than dimensions in decision space), then this can be used to efficiently solve the MOP by solving a number of MOPs with fewer objective functions. If this is not the case, the results can still give insight into the structure of the Pareto critical set.

Gerlach, R., Koltai, P., & Dellnitz, M. (2019). Revealing the intrinsic geometry of finite dimensional invariant sets of  infinite dimensional dynamical systems. In arXiv:1902.08824.

Embedding techniques allow the approximations of finite dimensional attractors and manifolds of infinite dimensional dynamical systems via subdivision and continuation methods. These approximations give a topological one-to-one image of the original set. In order to additionally reveal their geometry we use diffusion mapst o find intrinsic coordinates. We illustrate our results on the unstable manifold of the one-dimensional Kuramoto--Sivashinsky equation, as well as for the attractor of the Mackey-Glass delay differential equation.

Ziessler, A., Dellnitz, M., & Gerlach, R. (2019). The Numerical Computation of Unstable Manifolds for Infinite Dimensional Dynamical Systems by Embedding Techniques. SIAM Journal on Applied Dynamical Systems, 1265–1292.

In this work we extend the novel framework developed by Dellnitz, Hessel-von Molo, and Ziessler to the computation of finite dimensional unstable manifolds of infinite dimensional dynamical systems. To this end, we adapt a set-oriented continuation technique developed by Dellnitz and Hohmann for the computation of such objects of finite dimensional systems with the results obtained in the work of Dellnitz, Hessel-von Molo, and Ziessler. We show how to implement this approach for the analysis of partial differential equations and illustrate its feasibility by computing unstable manifolds of the one-dimensional Kuramoto--Sivashinsky equation as well as for the Mackey--Glass delay differential equation.


Gebken, B., Peitz, S., & Dellnitz, M. (2018). A Descent Method for Equality and Inequality Constrained Multiobjective Optimization Problems. In Numerical and Evolutionary Optimization – NEO 2017. Cham.

In this article we propose a descent method for equality and inequality constrained multiobjective optimization problems (MOPs) which generalizes the steepest descent method for unconstrained MOPs by Fliege and Svaiter to constrained problems by using two active set strategies. Under some regularity assumptions on the problem, we show that accumulation points of our descent method satisfy a necessary condition for local Pareto optimality. Finally, we show the typical behavior of our method in a numerical example.

Peitz, S., & Dellnitz, M. (2018). A Survey of Recent Trends in Multiobjective Optimal Control—Surrogate Models, Feedback Control and Objective Reduction. Mathematical and Computational Applications, 23(2).

Multiobjective optimization plays an increasingly important role in modern applications, where several criteria are often of equal importance. The task in multiobjective optimization and multiobjective optimal control is therefore to compute the set of optimal compromises (the Pareto set) between the conflicting objectives. The advances in algorithms and the increasing interest in Pareto-optimal solutions have led to a wide range of new applications related to optimal and feedback control, which results in new challenges such as expensive models or real-time applicability. Since the Pareto set generally consists of an infinite number of solutions, the computational effort can quickly become challenging, which is particularly problematic when the objectives are costly to evaluate or when a solution has to be presented very quickly. This article gives an overview of recent developments in accelerating multiobjective optimal control for complex problems where either PDE constraints are present or where a feedback behavior has to be achieved. In the first case, surrogate models yield significant speed-ups. Besides classical meta-modeling techniques for multiobjective optimization, a promising alternative for control problems is to introduce a surrogate model for the system dynamics. In the case of real-time requirements, various promising model predictive control approaches have been proposed, using either fast online solvers or offline-online decomposition. We also briefly comment on dimension reduction in many-objective optimization problems as another technique for reducing the numerical effort.

Vieluf, S., Mora, K., Gölz, C., Reuter, E.-M., Godde, B., Dellnitz, M., … Voelcker-Rehage, C. (2018). Age- and Expertise-Related Differences of Sensorimotor Network Dynamics during Force Control. Neuroscience, 203–213.

Klus, S., Peitz, S., & Schuster, I. (2018). Analyzing high-dimensional time-series data using kernel transfer  operator eigenfunctions. ArXiv:1805.10118.

Kernel transfer operators, which can be regarded as approximations of transfer operators such as the Perron-Frobenius or Koopman operator in reproducing kernel Hilbert spaces, are defined in terms of covariance and cross-covariance operators and have been shown to be closely related to the conditional mean embedding framework developed by the machine learning community. The goal of this paper is to show how the dominant eigenfunctions of these operators in combination with gradient-based optimization techniques can be used to detect long-lived coherent patterns in high-dimensional time-series data. The results will be illustrated using video data and a fluid flow example.

Peitz, S. (2018). Controlling nonlinear PDEs using low-dimensional bilinear approximations  obtained from data. ArXiv:1801.06419.

In a recent article, we presented a framework to control nonlinear partial differential equations (PDEs) by means of Koopman operator based reduced models and concepts from switched systems. The main idea was to transform a control system into a set of autonomous systems for which the optimal switching sequence has to be computed. These individual systems can be approximated very efficiently by reduced order models obtained from data, and one can guarantee equality of the full and the reduced objective function under certain assumptions. In this article, we extend these results to continuous control inputs using convex combinations of multiple Koopman operators corresponding to constant controls, which results in a bilinear control system. Although equality of the objectives can be carried over when the PDE depends linearly on the control, we show that this approach is also valid in other scenarios using several flow control examples of varying complexity.

Klus, S., Nüske, F., Koltai, P., Wu, H., Kevrekidis, I., Schütte, C., & Noé, F. (2018). Data-Driven Model Reduction and Transfer Operator Approximation. Journal of Nonlinear Science, 985–1010.

Gölz, C., Voelcker-Rehage, C., Mora, K., Reuter, E.-M., Godde, B., Dellnitz, M., … Vieluf, S. (2018). Improved Neural Control of Movements Manifests in Expertise-Related Differences in Force Output and Brain Network Dynamics. Frontiers in Physiology.

Hanke, S., Peitz, S., Wallscheid, O., Klus, S., Böcker, J., & Dellnitz, M. (2018). Koopman Operator-Based Finite-Control-Set Model Predictive Control for  Electrical Drives. ArXiv:1804.00854.

Predictive control of power electronic systems always requires a suitable model of the plant. Using typical physics-based white box models, a trade-off between model complexity (i.e. accuracy) and computational burden has to be made. This is a challenging task with a lot of constraints, since the model order is directly linked to the number of system states. Even though white-box models show suitable performance in most cases, parasitic real-world effects often cannot be modeled satisfactorily with an expedient computational load. Hence, a Koopman operator-based model reduction technique is presented which directly links the control action to the system's outputs in a black-box fashion. The Koopman operator is a linear but infinite-dimensional operator describing the dynamics of observables of nonlinear autonomous dynamical systems which can be nicely applied to the switching principle of power electronic devices. Following this data-driven approach, the model order and the number of system states are decoupled which allows us to consider more complex systems. Extensive experimental tests with an automotive-type permanent magnet synchronous motor fed by an IGBT 2-level inverter prove the feasibility of the proposed modeling technique in a finite-set model predictive control application.

Peitz, S., Ober-Blöbaum, S., & Dellnitz, M. (2018). Multiobjective Optimal Control Methods for the Navier-Stokes Equations Using Reduced Order Modeling. Acta Applicandae Mathematicae, 161(1), 171–199.

In a wide range of applications it is desirable to optimally control a dynamical system with respect to concurrent, potentially competing goals. This gives rise to a multiobjective optimal control problem where, instead of computing a single optimal solution, the set of optimal compromises, the so-called Pareto set, has to be approximated. When the problem under consideration is described by a partial differential equation (PDE), as is the case for fluid flow, the computational cost rapidly increases and makes its direct treatment infeasible. Reduced order modeling is a very popular method to reduce the computational cost, in particular in a multi query context such as uncertainty quantification, parameter estimation or optimization. In this article, we show how to combine reduced order modeling and multiobjective optimal control techniques in order to efficiently solve multiobjective optimal control problems constrained by PDEs. We consider a global, derivative free optimization method as well as a local, gradient-based approach for which the optimality system is derived in two different ways. The methods are compared with regard to the solution quality as well as the computational effort and they are illustrated using the example of the flow around a cylinder and a backward-facing-step channel flow.

Beermann, D., Dellnitz, M., Peitz, S., & Volkwein, S. (2018). POD-based multiobjective optimal control of PDEs with non-smooth objectives. In PAMM (pp. 51–54).

A framework for set‐oriented multiobjective optimal control of partial differential equations using reduced order modeling has recently been developed [1]. Following concepts from localized reduced bases methods, error estimators for the reduced cost functionals are utilized to construct a library of locally valid reduced order models. This way, a superset of the Pareto set can efficiently be computed while maintaining a prescribed error bound. In this article, this algorithm is applied to a problem with non‐smooth objective functionals. Using an academic example, we show that the extension to non‐smooth problems can be realized in a straightforward manner. We then discuss the implications on the numerical results.

Hruska, E., Abella, J. R., Nüske, F., Kavraki, L. E., & Clementi, C. (2018). Quantitative comparison of adaptive sampling methods for protein dynamics. The Journal of Chemical Physics.

Litzinger, F., Boninsegna, L., Wu, H., Nüske, F., Patel, R., Baraniuk, R., … Clementi, C. (2018). Rapid Calculation of Molecular Kinetics Using Compressed Sensing. Journal of Chemical Theory and Computation, 2771–2783.

Beermann, D., Dellnitz, M., Peitz, S., & Volkwein, S. (2018). Set-Oriented Multiobjective Optimal Control of PDEs Using Proper Orthogonal Decomposition. In Reduced-Order Modeling (ROM) for Simulation and Optimization (pp. 47–72). Cham.

In this chapter, we combine a global, derivative-free subdivision algorithm for multiobjective optimization problems with a posteriori error estimates for reduced-order models based on Proper Orthogonal Decomposition in order to efficiently solve multiobjective optimization problems governed by partial differential equations. An error bound for a semilinear heat equation is developed in such a way that the errors in the conflicting objectives can be estimated individually. The resulting algorithm constructs a library of locally valid reduced-order models online using a Greedy (worst-first) search. Using this approach, the number of evaluations of the full-order model can be reduced by a factor of more than 1000.

Boninsegna, L., Nüske, F., & Clementi, C. (2018). Sparse learning of stochastic dynamical equations. The Journal of Chemical Physics.

Klus, S., Gelß, P., Peitz, S., & Schütte, C. (2018). Tensor-based dynamic mode decomposition. Nonlinearity, 31(7), 3359–3380.

Dynamic mode decomposition (DMD) is a recently developed tool for the analysis of the behavior of complex dynamical systems. In this paper, we will propose an extension of DMD that exploits low-rank tensor decompositions of potentially high-dimensional data sets to compute the corresponding DMD modes and eigenvalues. The goal is to reduce the computational complexity and also the amount of memory required to store the data in order to mitigate the curse of dimensionality. The efficiency of these tensor-based methods will be illustrated with the aid of several different fluid dynamics problems such as the von Kármán vortex street and the simulation of two merging vortices.


Peitz, S. (2017).   Exploiting structure in multiobjective optimization and optimal control.

Multiobjective optimization plays an increasingly important role in modern applications, where several criteria are often of equal importance. The task in multiobjective optimization and multiobjective optimal control is therefore to compute the set of optimal compromises (the Pareto set) between the conflicting objectives. Since – in contrast to the solution of a single objective optimization problem – the Pareto set generally consists of an infinite number of solutions, the computational effort can quickly become challenging. This is even more the case when many problems have to be solved, when the number of objectives is high, or when the objectives are costly to evaluate. Consequently, this thesis is devoted to the identification and exploitation of structure both in the Pareto set and the dynamics of the underlying model as well as to the development of efficient algorithms for solving problems with additional parameters, with a high number of objectives or with PDE-constraints. These three challenges are addressed in three respective parts. In the first part, predictor-corrector methods are extended to entire Pareto sets. When certain smoothness assumptions are satisfied, then the set of parameter dependent Pareto sets possesses additional structure, i.e. it is a manifold. The tangent space can be approximated numerically which yields a direction for the predictor step. In the corrector step, the predicted set converges to the Pareto set at a new parameter value. The resulting algorithm is applied to an example from autonomous driving. In the second part, the hierarchical structure of Pareto sets is investigated. When considering a subset of the objectives, the resulting solution is a subset of the Pareto set of the original problem. Under additional smoothness assumptions, the respective subsets are located on the boundary of the Pareto set of the full problem. This way, the “skeleton” of a Pareto set can be computed and due to the exponential increase in computing time with the number of objectives, the computations of these subsets are significantly faster which is demonstrated using an example from industrial laundries. In the third part, PDE-constrained multiobjective optimal control problems are addressed by reduced order modeling methods. Reduced order models exploit the structure in the system dynamics, for example by describing the dynamics of only the most energetic modes. The model reduction introduces an error in both the function values and their gradients, which has to be taken into account in the development of algorithms. Both scalarization and set-oriented approaches are coupled with reduced order modeling. Convergence results are presented and the numerical benefit is investigated. The algorithms are applied to semi-linear heat flow problems as well as to the Navier-Stokes equations.

Peitz, S., Schäfer, K., Ober-Blöbaum, S., Eckstein, J., Köhler, U., & Dellnitz, M. (2017). A multiobjective MPC approach for autonomously driven electric vehicles. Proceedings of the 20th World Congress of the International Federation of Automatic Control (IFAC), 50(1), 8674–8679.

We present a new algorithm for model predictive control of non-linear systems with respect to multiple, conflicting objectives. The idea is to provide a possibility to change the objective in real-time, e.g. as a reaction to changes in the environment or the system state itself. The algorithm utilises elements from various well-established concepts, namely multiobjective optimal control, economic as well as explicit model predictive control and motion planning with motion primitives. In order to realise real-time applicability, we split the computation into an online and an offline phase and we utilise symmetries in the open-loop optimal control problem to reduce the number of multiobjective optimal control problems that need to be solved in the offline phase. The results are illustrated using the example of an electric vehicle where the longitudinal dynamics are controlled with respect to the concurrent objectives arrival time and energy consumption.

Peitz, S., Schäfer, K., Ober-Blöbaum, S., Eckstein, J., Köhler, U., & Dellnitz, M. (2017). A Multiobjective MPC Approach for Autonomously Driven Electric Vehicles * *This research was funded by the German Federal Ministry of Education and Research (BMBF) within the Leading-Edge Cluster Intelligent Technical Systems OstWestfalenLippe (it’s OWL). IFAC-PapersOnLine, 8674–8679.

Dellnitz, M., Klus, S., & Ziessler, A. (2017). A Set-Oriented Numerical Approach for Dynamical Systems with Parameter Uncertainty. SIAM Journal on Applied Dynamical Systems, 120–138.

Peitz, S., & Dellnitz, M. (2017). Gradient-Based Multiobjective Optimization with Uncertainties. In NEO 2016 (pp. 159–182). Cham.

In this article we develop a gradient-based algorithm for the solution of multiobjective optimization problems with uncertainties. To this end, an additional condition is derived for the descent direction in order to account for inaccuracies in the gradients and then incorporated into a subdivision algorithm for the computation of global solutions to multiobjective optimization problems. Convergence to a superset of the Pareto set is proved and an upper bound for the maximal distance to the set of substationary points is given. Besides the applicability to problems with uncertainties, the algorithm is developed with the intention to use it in combination with model order reduction techniques in order to efficiently solve PDE-constrained multiobjective optimization problems.

Nüske, F., Wu, H., Prinz, J.-H., Wehmeyer, C., Clementi, C., & Noé, F. (2017). Markov state models from short non-equilibrium simulations—Analysis and correction of estimation bias. The Journal of Chemical Physics.

Dellnitz, M., & Klus, S. (2017). Sensing and control in symmetric networks. Dynamical Systems, 61–79.


Eckstein, J., Peitz, S., Schäfer, K., Friedel, P., Köhler, U., Hessel von Molo, M., Ober-Blöbaum, S., & Dellnitz, M. (2016). A comparison of two predictive approaches to control the longitudinal dynamics of electric vehicles. Procedia Technology, 3rd International Conference on System-Integrated Intelligence: New Challenges for Product and Production Engineering, 26, 465–472.

In this contribution we compare two different approaches to the implementation of a Model Predictive Controller in an electric vehicle with respect to the quality of the solution and real-time applicability. The goal is to develop an intelligent cruise control in order to extend the vehicle range, i.e. to minimize energy consumption, by computing the optimal torque profile for a given track. On the one hand, a path-based linear model with strong simplifications regarding the vehicle dynamics is used. On the other hand, a nonlinear model is employed in which the dynamics of the mechanical and electrical subsystem are modeled.

Peitz, S., Gräler, M., Henke, C., Molo, M. H., Dellnitz, M., & Trächtler, A. (2016). Multiobjective Model Predictive Control of an Industrial Laundry. In Procedia Technology (pp. 483–490).

In a wide range of applications, it is desirable to optimally control a system with respect to concurrent, potentially competing goals. This gives rise to a multiobjective optimal control problem where, instead of computing a single optimal solution, the set of optimal compromises, the so-called Pareto set, has to be approximated. When it is not possible to compute the entire control trajectory in advance, for instance due to uncertainties or unforeseeable events, model predictive control methods can be applied to control the system during operation in real time. In this article, we present an algorithm for the solution of multiobjective model predictive control problems. In an offline scenario, it can be used to compute the entire set of optimal compromises whereas in a real time scenario, one optimal compromise is computed according to an operator's preference. The results are illustrated using the example of an industrial laundry. A logistics model of the laundry is developed and then utilized in the optimization routine. Results are presented for an offline as well as an online scenario.

Dellnitz, M., Eckstein, J., Flaßkamp, K., Friedel, P., Horenkamp, C., Köhler, U., … Tiemeyer, S. (2016). Multiobjective Optimal Control Methods for the Development of an Intelligent Cruise Control. In Progress in Industrial Mathematics at ECMI (Vol. 22, pp. 633–641). Cham: Springer International Publishing.

Dellnitz, M., Eckstein, J., Flaßkamp, K., Friedel, P., Horenkamp, C., Köhler, U., … Tiemeyer, S. (2016). Multiobjective Optimal Control Methods for the Development of an Intelligent Cruise Control. In Mathematics in Industry. Cham.

Ziessler, A., Molo, M. H.-V., & Dellnitz, M. (2016). On the computation of attractors for delay differential equations. Journal of Computational Dynamics, 5–5.


Demoures, F., Gay-Balmaz, F., Leyendecker, S., Ober-Blöbaum, S., Ratiu, T. S., & Weinand, Y. (2015). Discrete variational Lie group formulation of geometrically exact beam dynamics. Numerische Mathematik, 73–123.

Campos, C. M., Ober-Blöbaum, S., & Trélat, E. (2015). High order variational integrators in the optimal control of mechanical systems. Discrete & Continuous Dynamical Systems - A, 35(9), 4193–4223.

Ober-Blöbaum, S., & Padberg-Gehle, K. (2015). Multiobjective optimal control of fluid mixing. PAMM, 639–640.

Peitz, S., & Dellnitz, M. (2015). Multiobjective Optimization of the Flow Around a Cylinder Using Model Order Reduction. PAMM, 15(1), 613–614.

In this article an efficient numerical method to solve multiobjective optimization problems for fluid flow governed by the Navier Stokes equations is presented. In order to decrease the computational effort, a reduced order model is introduced using Proper Orthogonal Decomposition and a corresponding Galerkin Projection. A global, derivative free multiobjective optimization algorithm is applied to compute the Pareto set (i.e. the set of optimal compromises) for the concurrent objectives minimization of flow field fluctuations and control cost. The method is illustrated for a 2D flow around a cylinder at Re = 100.

Peitz, S., & Dellnitz, M. (2015). Multiobjective Optimization of the Flow Around a Cylinder Using Model Order Reduction. In PAMM (pp. 613–614).

n this article an efficient numerical method to solve multiobjective optimization problems for fluid flow governed by the Navier Stokes equations is presented. In order to decrease the computational effort, a reduced order model is introduced using Proper Orthogonal Decomposition and a corresponding Galerkin Projection. A global, derivative free multiobjective optimization algorithm is applied to compute the Pareto set (i.e. the set of optimal compromises) for the concurrent objectives minimization of flow field fluctuations and control cost. The method is illustrated for a 2D flow around a cylinder at Re = 100.

Keuck, L., Frohleke, N., Bocker, J., & Ziessler, A. (2015). PFC-control for improved inductor utilization. In 2015 17th European Conference on Power Electronics and Applications (EPE’15 ECCE-Europe).

Ringkamp, M., Ober-Blöbaum, S., & Leyendecker, S. (2015). Relaxing mixed integer optimal control problems using a time transformation. PAMM, 27–30.

Flasskamp, K., Hage-Packhäuser, S., & Ober-Blöbaum, S. (2015). Symmetry exploiting control of hybrid mechanical systems. Journal of Computational Dynamics, 2(1), 25–50.


Flaßkamp, K., Timmermann, J., Ober-Blöbaum, S., & Trächtler, A. (2014). Control strategies on stable manifolds for energy-efficient swing-ups of double pendula. International Journal of Control, 1886–1905.

Anacker, H., Dellnitz, M., Flaßkamp, K., Groesbrink, S., Hartmann, P., Heinzemann, C., … Ziegert, S. (2014). Methods for the Design and Development. In Lecture Notes in Mechanical Engineering. Berlin, Heidelberg.

Gail, T., Leyendecker, S., & Ober-Blöbaum, S. (2014). On the role of quadrature rules and system dimensions in variational multirateintegrators. In The 3rdJoint International Conference on Multibody System Dynamics.

Flasskamp, K., Ober-Blobaum, S., Schneider, T., & Bocker, J. (2014). Optimal control of a switched reluctance drive by a direct method using a discrete variational principle. In 52nd IEEE Conference on Decision and Control.

Baier, R., Dellnitz, M., Hessel-von Molo, M., Sertl, S., & Kevrekidis, I. G. (2014). The computation of convex invariant sets via Newton’s method. Journal of Computational Dynamics, 1(1), 39–69.


Ober-Blöbaum, S., & Seifried, A. (2013). A multiobjective optimization approach for optimal control problems of mechanical systems with uncertainties. In 2013 European Control Conference (ECC).

Witting, K., Ober-Blöbaum, S., & Dellnitz, M. (2013). A variational approach to define robustness for parametric multiobjective optimization problems. Journal of Global Optimization, 331–345.

Noé, F., & Nüske, F. (2013). A Variational Approach to Modeling Slow Processes in Stochastic Dynamical Systems. Multiscale Modeling & Simulation, 635–655.

Timmermann, R., & Dellnitz, M. (2013). Analysis of team and player performance using recorded trajectory data. In Performance Analysis of Sport IX, Part 8, Routledge.

Demoures, F., Gay-Balmaz, F., Leitz, T., Leyendecker, S., Ober-Blöbaum, S., & Ratiu, T. S. (2013). Asynchronous variational Lie group integration for geometrically exact beam dynamics. PAMM, 45–46.

Gail, T., Leyendecker, S., & Ober-Blöbaum, S. (2013). Computing time investigations for variational multirate integration. PAMM, 43–44.

Specht, A., Ober-Blobaum, S., Wallscheid, O., Romaus, C., & Bocker, J. (2013). Discrete-time model of an IPMSM based on variational integrators. In 2013 International Electric Machines & Drives Conference.

Flasskamp, K., Murphey, T., & Ober-Blobaum, S. (2013). Discretized switching time optimization problems. In 2013 European Control Conference (ECC).

Froyland, G. (2013). Markov modelling for random dynamical systems. In Equadiff 99.

Flaßkamp, K., Murphey, T., & Ober-Blöbaum, S. (2013). Optimization for discretized switched systems. PAMM, 401–402.

Jakobsmeyer, R., Schnittker, R., Timmermann, R., Zorn, R., Rückert, U., & Baumeister, J. (2013). Running performance analysis in basketball using recorded trajectory data. In Performance Analysis of Sport IX, Part 8, Routledge.

Schütze, O., Witting, K., Ober-Blöbaum, S., & Dellnitz, M. (2013). Set Oriented Methods for the Numerical Treatment of Multiobjective Optimization Problems. In EVOLVE- A Bridge between Probability, Set Oriented Numerics and Evolutionary Computation. Berlin, Heidelberg.


Meyer, A. (2012). Discontinuity Induced Bifurcations in Timed Continuous Petri Nets. IFAC Proceedings Volumes, 28–33.

Flaskamp, K., & Ober-Blobaum, S. (2012). Energy efficient control for mechanical systems based on inherent dynamical structures. In 2012 American Control Conference (ACC).

Ringkamp, M., Ober-Blöbaum, S., Dellnitz, M., & Schütze, O. (2012). Handling high-dimensional problems with multi-objective continuation methods via successive approximation of the tangent space. Engineering Optimization, 1117–1146.

Zanzottera, A., Mingotti, G., Castelli, R., & Dellnitz, M. (2012). Intersecting invariant manifolds in spatial restricted three-body problems: Design and optimization of Earth-to-halo transfers in the Sun–Earth–Moon scenario. Communications in Nonlinear Science and Numerical Simulation, 832–843.

Dellnitz, M., Dignath, F., Flaßkamp, K., Molo, M. H., Krüger, M., Timmermann, R., & Zheng, Q. (2012). Modelling and Analysis of the Nonlinear Dynamics of the Transrapid and Its Guideway. In Mathematics in Industry. Berlin, Heidelberg.

Dellnitz, M., Dignath, F., Flaßkamp, K., Molo, M. H., Krüger, M., Timmermann, R., & Zheng, Q. (2012). Modelling and Analysis of the Nonlinear Dynamics of the Transrapid and Its Guideway. In Mathematics in Industry. Berlin, Heidelberg.

Dellnitz, M., & Schütze, O. (2012). Multilevel Subdivision Techniques for Scalar Optimization Problems. In Global Analysis of Nonlinear Dynamics. New York, NY.

Baier, R., & Molo, M. H. (2012). Newton’s Method and Secant Method for Set-Valued Mappings. In Large-Scale Scientific Computing. Berlin, Heidelberg.

Flaßkamp, K., Timmermann, J., Ober-Blöbaum, S., Dellnitz, M., & Trächtler, A. (2012). Optimal Control on Stable Manifolds for a Double Pendulum. PAMM, 723–724.

Flaßkamp, K., Ober-Blöbaum, S., & Kobilarov, M. (2012). Solving Optimal Control Problems by Exploiting Inherent Dynamical Systems Structures. Journal of Nonlinear Science, 599–629.

Flasskamp, K., Murphey, T., & Ober-Blobaum, S. (2012). Switching time optimization in discretized hybrid dynamical systems. In 2012 IEEE 51st IEEE Conference on Decision and Control (CDC).

Meyer, A., & Silva, M. (2012). Symmetry Reductions in Timed Continuous Petri Nets Under Infinite Server Semantics. IFAC Proceedings Volumes, 153–159.


Klus, S., Sahai, T., Liu, C., & Dellnitz, M. (2011). An efficient algorithm for the parallel solution of high-dimensional differential equations. Journal of Computational and Applied Mathematics, 3053–3062.

Krüger, M., Witting, K., Trächtler, A., & Dellnitz, M. (2011). Parametric Model-Order Reduction in Hierarchical Multiobjective Optimization of Mechatronic Systems*. IFAC Proceedings Volumes, 12611–12619.

Hage-Packhäuser, S., & Dellnitz, M. (2011). Stabilization via symmetry switching in hybrid dynamical systems. Discrete & Continuous Dynamical Systems - B, 16(1), 239–263.

Meyer, A., Dellnitz, M., & Hessel-von Molo, M. (2011). Symmetries in timed continuous Petri nets. Nonlinear Analysis: Hybrid Systems, 125–135.


Dellnitz, M., Padberg, K., Preis, R., & Thiere, B. (2010). Continuous and Discrete Concepts for Detecting Transport Barriers in the Planar Circular Restricted Three Body Problem. In Nonlinear Science and Complexity. Dordrecht.

Thiere, B., Ober-Blöbaum, S., & Pergola, P. (2010). Detecting Initial Guesses for Trajectories in the (P)CRTBP. In AIAA/AAS Astrodynamics Specialist Conference.

Flaßkamp, K., Ober-Blöbaum, S., & Kobilarov, M. (2010). Solving optimal control problems by using inherent dynamical properties. PAMM, 577–578.

Tantardini, M., Fantino, E., Ren, Y., Pergola, P., Gómez, G., & Masdemont, J. J. (2010). Spacecraft trajectories to the L3 point of the Sun–Earth three-body problem. Celestial Mechanics and Dynamical Astronomy, 215–232.


Biasco, L., & Di Gregorio, L. (2009). A Birkhoff–Lewis Type Theorem for the Nonlinear Wave Equation. Archive for Rational Mechanics and Analysis, 303–362.

Dellnitz, M., Ober-Blöbaum, S., Post, M., Schütze, O., & Thiere, B. (2009). A multi-objective approach to the design of low thrust space trajectories using optimal control. Celestial Mechanics and Dynamical Astronomy, 33–59.

Dellnitz, M., & Witting, K. (2009). Computation of robust Pareto points. International Journal of Computing Science and Mathematics.

Schütze, O., Vasile, M., Junge, O., Dellnitz, M., & Izzo, D. (2009). Designing optimal low-thrust gravity-assist trajectories using space pruning and a multi-objective approach. Engineering Optimization, 155–181.

Schneider, T., Schulz, B., Henke, C., Witting, K., Steenken, D., & Bocker, J. (2009). Energy transfer via linear doubly-fed motor in different operating modes. In 2009 IEEE International Electric Machines and Drives Conference.

Padberg, K., Thiere, B., Preis, R., & Dellnitz, M. (2009). Local expansion concepts for detecting transport barriers in dynamical systems. Communications in Nonlinear Science and Numerical Simulation, 4176–4190.

Wang, F., & Dellnitz, M. (2009). Multi-objective shape optimization for piezoceramics. In 2008 Symposium on Piezoelectricity, Acoustic Waves, and Device Applications.

Blesken, M., Ruckert, U., Steenken, D., Witting, K., & Dellnitz, M. (2009). Multiobjective optimization for transistor sizing of CMOS logic standard cells using set-oriented numerical techniques. In 2009 NORCHIP.

Dellnitz, M., Froyland, G., Horenkamp, C., & Padberg, K. (2009). On the Approximation of Transport Phenomena - a Dynamical Systems Approach. GAMM-Mitteilungen, 47–60.

Romaus, C., Bocker, J., Witting, K., Seifried, A., & Znamenshchykov, O. (2009). Optimal energy management for a hybrid energy storage system combining batteries and double layer capacitors. In 2009 IEEE Energy Conversion Congress and Exposition.

Dellnitz, M., Froyland, G., Horenkamp, C., Padberg-Gehle, K., & Sen Gupta, A. (2009). Seasonal variability of the subpolar gyres in the Southern Ocean: a numerical investigation based on transfer operators. Nonlinear Processes in Geophysics, 655–663.

Meyer, A., & Dellnitz, M. (2009). Symmetries in Timed Continuous Petri Nets. IFAC Proceedings Volumes, 68–73.


Witting, K., Schulz, B., Dellnitz, M., Böcker, J., & Fröhleke, N. (2008). A new approach for online multiobjective optimization of mechatronic systems. International Journal on Software Tools for Technology Transfer, 223–231.

Schütze, O., Laumanns, M., Coello Coello, C. A., Dellnitz, M., & Talbi, E.-G. (2008). Convergence of stochastic search algorithms to finite size pareto set approximations. Journal of Global Optimization, 559–577.

Klöpper, B., Podlogar, H., Gausemeier, J., & Witting, K. (2008). Domain Spanning Search for the Identification of Solution Patterns for the Conceptual Design of Self-Optimizing Systems. In 2008 19th International Conference on Database and Expert Systems Applications.

Schütze, O., Coello Coello, C. A., Mostaghim, S., Talbi, E.-G., & Dellnitz, M. (2008). Hybridizing evolutionary strategies with continuation methods for solving multi-objective problems. Engineering Optimization, 383–402.

Geisler, M. S. J., Witting, D.-M. K., Trächtler, A., & Dellnitz, M. (2008). Multiobjective Optimization of Control Trajectories for the Guidance of a Rail-bound Vehicle. IFAC Proceedings Volumes, 4380–4386.

Geisler, J., Witting, K., Trächtler, A., & Dellnitz, M. (2008). Self-Optimization of the Guidance Module of a Rail-bound Vehicle. In 7th International Heinz Nixdorf Symposium: Self-optimizing Mechatronic Systems: Designing the Future.


Schütze, O., Talbi, E., Pulido, G. T., Coello, C. C., & Santana-Quintero, L. V. (2007). A Memetic PSO Algorithm for Scalar Optimization Problems. In 2007 IEEE Swarm Intelligence Symposium.

Preis, R., Monien, B., & Schamberger, S. (2007). Approximation Algorithms for Multilevel Graph Partitioning. In Handbook of Approximation Algorithms and Metaheuristics.

Froyland, G., Padberg, K., England, M. H., & Treguier, A. M. (2007). Detection of Coherent Oceanic Structures via Transfer Operators. Physical Review Letters.

Leyendecker, S., Ober-Blöbaum, S., Marsden, J. E., & Ortiz, M. (2007). Discrete Mechanics and Optimal Control for Constrained Multibody Dynamics. In Volume 5: 6th International Conference on Multibody Systems, Nonlinear Dynamics, and Control, Parts A, B, and C.

<jats:p>This paper formulates the dynamical equations of mechanics subject to holonomic constraints in terms of the states and controls using a constrained version of the Lagrange-d’Alembert principle. Based on a discrete version of this principle, a structure preserving time-stepping scheme is derived. It is shown that this respect for the mechanical structure (such as a reliable computation of the energy and momentum budget, without numerical dissipation) is retained when the system is reduced to its minimal dimension by the discrete null space method. Together with initial and final conditions on the configuration and conjugate momentum, the reduced time-stepping equations serve as nonlinear equality constraints for the minimisation of a given cost functional. The algorithm yields a sequence of discrete configurations together with a sequence of actuating forces, optimally guiding the system from the initial to the desired final state. The resulting discrete optimal control algorithm is shown to have excellent energy and momentum properties, which are illustrated by two specific examples, namely reorientation and repositioning of a rigid body subject to external forces and the reorientation of a rigid body with internal momentum wheels.</jats:p>

Grüne, L., & Junge, O. (2007). Global Optimal Control of Perturbed Systems. Journal of Optimization Theory and Applications, 411–429.

Padberg, K., Hauff, T., Jenko, F., & Junge, O. (2007). Lagrangian structures and transport in turbulent magnetized plasmas. New Journal of Physics, 400–400.

Dellnitz, M., Padberg, K., Post, M., & Thiere, B. (2007). Set Oriented Approximation of Invariant Manifolds: Review of Concepts for Astrodynamical Problems. In AIP Conference Proceedings.

Vasile, M., Schütze, O., Junge, O., Radice, G., & Dellnitz, M. (2007). Spiral Trajectories in Global Optimisation of Interplanetary and Orbital Transfers. European Space Agency.


Knoke, T., Romaus, C., Bocker, J., Dell’Aere, A., & Witting, K. (2006). Energy Management for an Onboard Storage System Based on Multi-Objective Optimization. In IECON 2006 - 32nd Annual Conference on IEEE Industrial Electronics.

Sertl, S., & Dellnitz, M. (2006). Global Optimization using a Dynamical Systems Approach. Journal of Global Optimization, 569–587.

Dellnitz, M., Molo, M. H., Metzner, P., Preis, R., & Schütte, C. (2006). Graph Algorithms for Dynamical Systems. In Analysis, Modeling and Simulation of Multiscale Problems.

Dell’Aere, A. (2006). Multi-Objective Optimization in Self-Optimizing Systems. In IECON 2006 - 32nd Annual Conference on IEEE Industrial Electronics.

Dellnitz, M., Junge, O., Post, M., & Thiere, B. (2006). On target for Venus – set oriented computation of energy efficient low thrust trajectories. Celestial Mechanics and Dynamical Astronomy, 357–370.

Junge, O., & Ober-Blobaum, S. (2006). Optimal Reconfiguration of Formation Flying Satellites. In Proceedings of the 44th IEEE Conference on Decision and Control.

Junge, O., Marsden, J. E., & Ober-Blobaum, S. (2006). Optimal Reconfiguration of Formation Flying Spacecraft ---a Decentralized Approach. In Proceedings of the 45th IEEE Conference on Decision and Control.

Thiere, B. (2006). Return Time Dynamics as a Tool for Finding Almost Invariant Sets. Annals of the New York Academy of Sciences, 44–54.

Dellnitz, M., & Junge, O. (2006). Set Oriented Numerical Methods in Space Mission Design. In Elsevier Astrodynamics Series.

Mehta, P. G., Hessel-von Molo, M., & Dellnitz, M. (2006). Symmetry of attractors and the Perron-Frobenius operator. Journal of Difference Equations and Applications, 1147–1178.


Grüne, L., & Junge, O. (2005). A set oriented approach to optimal feedback stabilization. Systems & Control Letters, 169–180.

Krauskopf, B., Osinga, H. M., Doedel, E. J., Henderson, M. E., Guckenheimer, J., Vladimirsky, A., … Junge, O. (2005). A Survey of Methods for Computing (un)stable Manifolds of Vector Fields. International Journal of Bifurcation and Chaos, 763–791.

<jats:p> The computation of global invariant manifolds has seen renewed interest in recent years. We survey different approaches for computing a global stable or unstable manifold of a vector field, where we concentrate on the case of a two-dimensional manifold. All methods are illustrated with the same example — the two-dimensional stable manifold of the origin in the Lorenz system. </jats:p>

Dellnitz, M., Sch�tze, O., & Hestermeyer, T. (2005). Covering Pareto Sets by Multilevel Subdivision Techniques. Journal of Optimization Theory and Applications, 113–136.

Li, R., Pottharst, A., Frohieke, N., Becker, J., Witting, K., Dellnitz, M., … Feldmann, R. (2005). Design and implementation of a hybrid energy supply system for railway vehicles. In Twentieth Annual IEEE Applied Power Electronics Conference and Exposition, 2005. APEC 2005.

Krishnamurthy, A., & Preis, R. (2005). Satellite Formation, a Mobile Sensor Network in Space. In 19th IEEE International Parallel and Distributed Processing Symposium.


Day, S., Junge, O., & Mischaikow, K. (2004). A Rigorous Numerical Method for the Global Analysis of Infinite-Dimensional Discrete Dynamical Systems. SIAM Journal on Applied Dynamical Systems, 117–160.

Junge, O., & Osinga, H. M. (2004). A set oriented approach to global optimal control. ESAIM: Control, Optimisation and Calculus of Variations, 259–270.

Chaudhuri, I., Sertl, S., Hajnal, Z., Dellnitz, M., & Frauenheim, T. (2004). Global optimization of silicon nanoclusters. Applied Surface Science, 108–113.

Bezrukov, S., Elsässer, R., Monien, B., Preis, R., & Tillich, J.-P. (2004). New spectral lower bounds on the bisection width of graphs. Theoretical Computer Science, 155–174.

Elsässer, R., Monien, B., Preis, R., & Frommer, A. (2004). Optimal Diffusion Schemes and Load Balancing on Product Graphs. Parallel Processing Letters, 61–73.

<jats:p> We discuss nearest neighbor load balancing schemes on processor networks which are represented by a cartesian product of graphs and present a new optimal diffusion scheme for general graphs. In the first part of the paper, we introduce the Alternating-Direction load balancing scheme, which reduces the number of load balance iterations by a factor of 2 for cartesian products of graphs. The resulting flow is theoretically analyzed and can be very high for certain cases. Therefore, we further present the Mixed-Direction scheme which needs the same number of iterations but computes in most cases a much smaller flow. In the second part of the paper, we present a simple optimal diffusion scheme for general graphs, calculating a balancing flow which is minimal in the l<jats:sub>2</jats:sub> norm. It is based on the spectra of the graph representing the network and needs only m-1 iterations to balance the load with m being the number of distinct eigenvalues. Known optimal diffusion schemes have the same performance, however the optimal scheme presented in this paper can be implemented in a very simple manner. The number of iterations of optimal diffusion schemes is independent of the load scenario and, thus, they are practical for networks which represent graphs with known spectra. Finally, our experiments exhibit that the new optimal scheme can successfully be combined with the Alternating-Direction and Mixed-Direction schemes for efficient load balancing on product graphs. </jats:p>


Schütze, O. (2003). A New Data Structure for the Nondominance Problem in Multi-objective Optimization. In Lecture Notes in Computer Science. Berlin, Heidelberg.

Aston, P. J., & Dellnitz, M. (2003). Computation of the dominant Lyapunov exponent via spatial integration using matrix norms. Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences, 2933–2955.

Dellnitz, M., & Preis, R. (2003). Congestion and Almost Invariant Sets in Dynamical Systems. In Lecture Notes in Computer Science. Berlin, Heidelberg.

Schütze, O., Mostaghim, S., Dellnitz, M., & Teich, J. (2003). Covering Pareto Sets by Multilevel Evolutionary Subdivision Techniques. In Lecture Notes in Computer Science. Berlin, Heidelberg.

Froyland, G., & Dellnitz, M. (2003). Detecting and Locating Near-Optimal Almost-Invariant Sets and Cycles. SIAM Journal on Scientific Computing, 1839–1863.


Elsässer, R., Monien, B., & Preis, R. (2002). Diffusion Schemes for Load Balancing on Heterogeneous Networks. Theory of Computing Systems, 305–320.

Dellnitz, M. (2002). Finding zeros by multilevel subdivision techniques. IMA Journal of Numerical Analysis, 167–185.

Dellnitz, M., & Junge, O. (2002). Set Oriented Numerical Methods for Dynamical Systems. In Handbook of Dynamical Systems.


Junge, O. (2001). An adaptive subdivision technique for the approximation of attractors and invariant measures: proof of convergence. Dynamical Systems, 213–222.

Froyland, G. (2001). Extracting Dynamical Behavior via Markov Models. In Nonlinear Dynamics and Statistics. Boston, MA.

Froyland, G., Junge, O., & Ochs, G. (2001). Rigorous computation of topological entropy with respect to a finite partition. Physica D: Nonlinear Phenomena, 68–84.

Dellnitz, M., Froyland, G., & Junge, O. (2001). The Algorithms Behind GAIO — Set Oriented Numerical Methods for Dynamical Systems. In Ergodic Theory, Analysis, and Efficient Simulation of Dynamical Systems. Berlin, Heidelberg.


Dellnitz, M., Froyland, G., & Sertl, S. (2000). A Conjecture on the Existence of Isolated Eigenvalues of the Perron-Frobenius Operator. In Equadiff 99.

Dellnitz, M., Froyland, G., & Sertl, S. (2000). On the isolated spectrum of the Perron-Frobenius operator. Nonlinearity, 1171–1188.

Junge, O. (2000). Rigorous discretization of subdivision techniques. In Equadiff 99.

Dellnitz, M., Junge, O., Rumpf, M., & Strzodka, R. (2000). The computation of an unstable invariant set inside a cylinder containing a knotted flow. In Equadiff 99 (pp. 1053--1059).

Aston, P. J., & Dellnitz, M. (2000). The Computation of Lyapunov Exponents via Spatial Integration Using Vector Norms. In Equadiff 99.


Deuflhard, P., Dellnitz, M., Junge, O., & Schütte, C. (1999). Computation of Essential Molecular Dynamics by Subdivision Techniques. In Computational Molecular Dynamics: Challenges, Methods, Ideas. Berlin, Heidelberg.

Dellnitz, M., & Junge, O. (1999). On the Approximation of Complicated Dynamical Behavior. SIAM Journal on Numerical Analysis, 491–515.

Aston, P. J., & Dellnitz, M. (1999). The computation of lyapunov exponents via spatial integration with application to blowout bifurcations. Computer Methods in Applied Mechanics and Engineering, 223–237.


Dellnitz, M., & Junge, O. (1998). An adaptive subdivision technique for the approximation of attractors and invariant measures. Computing and Visualization in Science, 63–68.


Dellnitz, M., & Hohmann, A. (1997). A subdivision algorithm for the computation of unstable manifolds and global attractors. Numerische Mathematik, 75, 293–317.

Dellnitz, M., & Junge, O. (1997). Almost Invariant Sets in Chua’s Circuit. International Journal of Bifurcation and Chaos, 2475–2485.

<jats:p> Recently multilevel subdivision techniques have been introduced in the numerical investigation of complicated dynamical behavior. We illustrate the applicability and efficiency of these methods by a detailed numerical study of Chua's circuit. In particular we will show that there exist two regions in phase space which are almost invariant in the sense that typical trajectories stay inside each of these sets on average for quite a long time. </jats:p>

Guder, R., Dellnitz, M., & Kreuzer, E. (1997). An adaptive method for the approximation of the generalized cell mapping. Chaos, Solitons & Fractals, 525–534.

Dellnitz, M., Hohmann, A., Junge, O., & Rumpf, M. (1997). Exploring invariant sets and invariant measures. Chaos: An Interdisciplinary Journal of Nonlinear Science, 221–228.


Dellnitz, M., & Hohmann, A. (1996). The Computation of Unstable Manifolds Using Subdivision and Continuation. In Nonlinear Dynamical Systems and Chaos. Basel.


Dellnitz, M., & Melbourne, I. (1995). A note on the shadowing lemma and symmetric periodic points. Nonlinearity, 1067–1075.

Dellnitz, M., & Heinrich, C. (1995). Admissible symmetry increasing bifurcations. Nonlinearity, 1039–1066.

Dellnitz, M., Field, M., Golubitsky, M., Ma, J., & Hohmann, A. (1995). Cycling Chaos. International Journal of Bifurcation and Chaos, 1243–1247.

Dellnitz, M., Golubitsky, M., Hohmann, A., & Stewart, I. (1995). Spirals in Scalar Reaction–Diffusion Equations. International Journal of Bifurcation and Chaos, 1487–1501.

<jats:p> Spiral patterns have been observed experimentally, numerically, and theoretically in a variety of systems. It is often believed that these spiral wave patterns can occur only in systems of reaction–diffusion equations. We show, both theoretically (using Hopf bifurcation techniques) and numerically (using both direct simulation and continuation of rotating waves) that spiral wave patterns can appear in a single reaction–diffusion equation [ in u(x, t)] on a disk, if one assumes "spiral" boundary conditions (u<jats:sub>r</jats:sub> = mu<jats:sub>θ</jats:sub>). Spiral boundary conditions are motivated by assuming that a solution is infinitesimally an Archimedian spiral near the boundary. It follows from a bifurcation analysis that for this form of spirals there are no singularities in the spiral pattern (technically there is no spiral tip) and that at bifurcation there is a steep gradient between the "red" and "blue" arms of the spiral. </jats:p>

Aston, P. J., & Dellnitz, M. (1995). Symmetry Breaking Bifurcations of Chaotic Attractors. International Journal of Bifurcation and Chaos, 1643–1676.

<jats:p> In an array of coupled oscillators, synchronous chaos may occur in the sense that all the oscillators behave identically although the corresponding motion is chaotic. When a parameter is varied this fully symmetric dynamical state can lose its stability, and the main purpose of this paper is to investigate which type of dynamical behavior is expected to be observed once the loss of stability has occurred. The essential tool is a classification of Lyapunov exponents based on the symmetry of the underlying problem. This classification is crucial in the derivation of the analytical results but it also allows an efficient computation of the dominant Lyapunov exponent associated with each symmetry type. We show how these dominant exponents determine the stability of invariant sets possessing various instantaneous symmetries, and this leads to the idea of symmetry breaking bifurcations of chaotic attractors. Finally, the results and ideas are illustrated for several systems of coupled oscillators. </jats:p>

Golubitsky, M., Marsden, J., Stewart, I., & Dellnitz, M. (1995). The constrained Liapunov-Schmidt procedure and periodic orbits. In Normal Forms and Homoclinic Chaos. Providence, Rhode Island.


Dellnitz, M. (1994). Collisions of chaotic attractors. Schlaglichter Der Forschung: Zum 75. Jahrestag Der Universität Hamburg, 411–428.

Dellnitz, M., & Scheurle, J. (1994). Eigenvalue Movement for a Class of Reversible Hamiltonian Systems with Three Degrees of Freedom. In Dynamics, Bifurcation and Symmetry. Dordrecht.

Dellnitz, M., & Melbourne, I. (1994). Generic movement of eigenvalues for equivariant self-adjoint matrices. Journal of Computational and Applied Mathematics, 249–259.

Dellnitz, M., Golubitsky, M., & Nicol, M. (1994). Symmetry of Attractors and the Karhunen-Loève Decomposition. In Trends and Perspectives in Applied Mathematics. New York, NY.


Barany, E., Dellnitz, M., & Golubitsky, M. (1993). Detecting the symmetry of attractors. Physica D: Nonlinear Phenomena, 66–87.

Melbourne, I., & Dellnitz, M. (1993). Normal forms for linear Hamiltonian vector fields commuting with the action of a compact Lie group. Mathematical Proceedings of the Cambridge Philosophical Society, 235–268.

<jats:title>Abstract</jats:title><jats:p>We obtain normal forms for infinitesimally symplectic matrices (or linear Hamiltonian vector fields) that commute with the symplectic action of a compact Lie group of symmetries. In doing so we extend Williamson's theorem on normal forms when there is no symmetry present.</jats:p><jats:p>Using standard representation-theoretic results the symmetry can be factored out and we reduce to finding normal forms over a real division ring. There are three real division rings consisting of the real, complex and quaternionic numbers. Of these, only the real case is covered in Williamson's original work.</jats:p>

Dellnitz, M. (1993). The equivariant Darboux theorem. Lectures in Applied Mathematics, 29, 163–169.


Dellnitz, M. (1992). Computational bifurcation of periodic solutions in systems with symmetry. IMA Journal of Numerical Analysis, 12(3), 429–455.

Dellnitz, M., Melbourne, I., & Marsden, J. E. (1992). Generic bifurcation of Hamiltonian vector fields with symmetry. Nonlinearity, 979–996.

Dellnitz, M., Marsden, J. E., Melbourne, I., & Scheurle, J. (1992). Generic Bifurcations of Pendula. In Bifurcation and Symmetry. Basel.

Dellnitz, M., Golubitsky, M., & Melbourne, I. (1992). Mechanisms of Symmetry Creation. In Bifurcation and Symmetry. Basel.


Dellnitz, M., & Werner, B. (1989). Computational methods for bifurcation problems with symmetries—with special attention to steady state and Hopf bifurcation points. Journal of Computational and Applied Mathematics, 97–123.

Open list in Research Information System

Further information:

The University for the Information Society