Sebastian Peitz

Tensor-based dynamic mode decomposition

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 Karman vortex street and the simulation of two merging vortices.

Multiobjective Optimal Control Methods for the Development of an Intelligent Cruise Control

During the last years, alternative drive technologies, for example electrically powered vehicles (EV), have gained more and more attention, mainly caused by an increasing awareness of the impact of CO2 emissions on climate change and by the limitation of fossil fuels. However, these technologies currently come with new challenges due to limited lithium ion battery storage density and high battery costs which lead to a considerably reduced range in comparison to conventional internal combustion engine powered vehicles. For this reason, it is desirable to increase the vehicle range without enlarging the battery. When the route and the road slope are known in advance, it is possible to vary the vehicles velocity within certain limits in order to reduce the overall drivetrain energy consumption. This may either result in an increased range or, alternatively, in larger energy reserves for comfort functions such as air conditioning. In this presentation, we formulate the challenge of range extension as a multiobjective optimal control problem. We then apply different numerical methods to calculate the so-called Pareto set of optimal compromises for the drivetrain power profile with respect to the two concurrent objectives battery state of charge and mean velocity. In order to numerically solve the optimal control problem by means of a direct method, a time discretization of the drivetrain power profile is necessary. In combination with a vehicle dynamics simulation model, the optimal control problem is transformed into a high dimensional nonlinear optimization problem. For the approximation of the Pareto set, two different optimization algorithms implemented in the software package GAIO are used. The first one yields a global optimal solution by applying a set-oriented subdivision technique to parameter space. By construction, this technique is limited to coarse discretizations of the drivetrain power profile. In contrast, the second technique, which is based on an image space continuation method, is more suitable when the number of parameters is large while the number of objectives is less than five. We compare the solutions of the two algorithms and study the influence of different discretizations on the quality of the solutions. A MATLAB/Simulink model is used to describe the dynamics of an EV. It is based on a drivetrain efficiency map and considers vehicle properties such as rolling friction and air drag, as well as environmental conditions like slope and ambient temperature. The vehicle model takes into account the traction battery too, enabling an exact prediction of the batterys response to power requests of drivetrain and auxiliary loads, including state of charge.

Set-Oriented Multiobjective Optimal Control of PDEs using Proper Orthogonal Decomposition

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.

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).

Abstract 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.

Gradient-Based Multiobjective Optimization with Uncertainties

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.

Multiobjective Optimal Control Methods for the Navier-Stokes Equations Using Reduced Order Modeling

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.

A Descent Method for Equality and Inequality Constrained Multiobjective Optimization Problems

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.