Global Optimal Control Using Direct Multiple Shooting

The goal of this thesis is the development of a novel and efficient algorithm to determine the global optimum of an optimal control problem. In contrast to previous methods, the approach presented here is based on the direct multiple shooting method for discretizing the optimal control problem, which results in a significant increase of efficiency. To relax the discretized optimal control problems, the so-called alpha-branch-and-bound method in combination with validated integration is used. For the direct comparison of the direct single-shooting-based relaxations with the direct multipleshooting-based algorithm, several theoretical results are proven that build the basis for the efficiency increase of the novel method. A specialized branching strategy takes care that the additionally introduced variables due to the multiple shooting approach do not increase the size of the resulting branch-and-bound tree. An adaptive scaling technique of the commonly used Gershgorin method to estimate the eigenvalues of interval matrices leads to optimal relaxations and therefore leads to a general improvement of the alpha-branch-and-bound relaxations in a single shooting and a multiple shooting framework, as well as for the corresponding relaxations of non-dynamic nonlinear problems. To further improve the computational time, suggestions regarding the necessary second-order interval sensitivities are presented in this thesis, as well as a heuristic that relaxes on a subspace only. The novel algorithm, as well as the single-shooting-based alternative for a direct comparison, are implemented in a newly developed software package called GloOptCon. The new method is used to globally solve both a number of benchmark problems from the literature, and so far in the context of global optimal control still little considered applications. The additional problems pose new challenges either due to their size or due to having its origin in mixed integer optimal control with an integer-valued time-dependent control variable. The theoretically proven increase of efficiency is validated by the numerical results. Compared to the previous approach from the literature, the number of iterations for typical problems is more than halved, meanwhile the computation time is reduced by almost an order of magnitude. This in turn allows the global solution of significantly larger optimal control problems

Numerical Solution of Optimal Control Problems with Explicit and Implicit Switches

This dissertation deals with the efficient numerical solution of switched optimal control problems whose dynamics may coincidentally be affected by both explicit and implicit switches. A framework is being developed for this purpose, in which both problem classes are uniformly converted into a mixed–integer optimal control problem with combinatorial constraints. Recent research results relate this problem class to a continuous optimal control problem with vanishing constraints, which in turn represents a considerable subclass of an optimal control problem with equilibrium constraints. In this thesis, this connection forms the foundation for a numerical treatment. We employ numerical algorithms that are based on a direct collocation approach and require, in particular, a highly accurate determination of the switching structure of the original problem. Due to the fact that the switching structure is a priori unknown in general, our approach aims to identify it successively. During this process, a sequence of nonlinear programs, which are derived by applying discretization schemes to optimal control problems, is solved approximatively. After each iteration, the discretization grid is updated according to the currently estimated switching structure. Besides a precise determination of the switching structure, it is of central importance to estimate the global error that occurs when optimal control problems are solved numerically. Again, we focus on certain direct collocation discretization schemes and analyze error contributions of individual discretization intervals. For this purpose, we exploit a relationship between discrete adjoints and the Lagrange multipliers associated with those nonlinear programs that arise from the collocation transcription process. This relationship can be derived with the help of a functional analytic framework and by interrelating collocation methods and Petrov–Galerkin finite element methods. In analogy to the dual-weighted residual methodology for Galerkin methods, which is well–known in the partial differential equation community, we then derive goal–oriented global error estimators. Based on those error estimators, we present mesh refinement strategies that allow for an equilibration and an efficient reduction of the global error. In doing so we note that the grid adaption processes with respect to both switching structure detection and global error reduction get along with each other. This allows us to distill an iterative solution framework. Usually, individual state and control components have the same polynomial degree if they originate from a collocation discretization scheme. Due to the special role which some control components have in the proposed solution framework it is desirable to allow varying polynomial degrees. This results in implementation problems, which can be solved by means of clever structure exploitation techniques and a suitable permutation of variables and equations. The resulting algorithm was developed in parallel to this work and implemented in a software package. The presented methods are implemented and evaluated on the basis of several benchmark problems. Furthermore, their applicability and efficiency is demonstrated. With regard to a future embedding of the described methods in an online optimal control context and the associated real-time requirements, an extension of the well–known multi–level iteration schemes is proposed. This approach is based on the trapezoidal rule and, compared to a full evaluation of the involved Jacobians, it significantly reduces the computational costs in case of sparse data matrices