# Optimization with Partial Differential Equations

### Open loop Control and Efficient Solution of Linear Saddle Point Systems

The first category is the so-called open loop control problem where the system input is computed purely on the basis of a-priori knowledge. Here, we focus on the efficient solution of the optimization problem with state-of-the-art numerical techniques. The task of solving optimal control problems is challenging as most realistic problems are posed in three dimensions and the discretization, often done with finite elements, results in many degrees of freedom. At the heart of the problem typically lies the solution of the first order optimality conditions or Karush–Kuhn–Tucker (KKT) conditions and the corresponding linear system in saddle point form. For realistic scenarios it is not possible to use techniques based on factorizations to solve this problem. Hence iterative Krylov subspace solvers have to be employed. These are usually only feasible if some form of preconditioning is applied; a technique that solves a modified, simpler problem. We derived preconditioners that take the structure of the KKT system and the underlying infinite-dimensional problem into account. For steady and unsteady problems we constructed efficient approximations to the Schur-complement of the saddle point problem, which involves the forward and the adjoint discretized PDE. In the case of unsteady problems we showed that multilevel techniques approximating only one matrix provide sufficient approximations for the Schur-complement and that independence with respect to the mesh-parameter is given. For a variety of problems, such as the (time-periodic) heat equation or the unsteady Stokes equation, we illustrated that only crude approximations to the discretized PDE suffice to obtain accurate results to the optimal control problem.

## Internal Projects

#### Optimal Control of Chemical Processes

#### Fast Solvers for Phase Field Models (with NDS group)

#### Low-rank Solver (with NDS group)

#### Preconditioning for PDE-constrained optimization (with NDS group)

## Concluded Projects

#### High Performance Computing and PDE-constrained optimization (with NDS group)

### Closed Loop Control and Solution of Large Scale Matrix Equations

The second category of control problems studied within the research group is the one of regulator problems. These problems are described by the minimization of an objective function over a possibly infinite time horizon incorporating a feedback loop, i.e., the incorporation of certain output (observed quantities) of the current system behavior into the control input. Again, one is interested in the minimization of a functional subject to unsteady, possibly nonlinear, PDEs. In general the PDEs are rewritten as an infinite dimensional evolution equation, or Cauchy problem. In the case of a linear PDE the optimization problem can then be posed as a linear quadratic regulator (LQR) problem that we solve via a system theoretic approach leading to algebraic or differential Riccati equations. For nonlinear PDEs this approach must be embedded in a model predictive control (MPC) scheme, where local linearizations around working points (or working trajectories) are employed. In any case after semi-discretization using the method of lines, a very large system of ordinary differential equations (ODEs) and related very large matrix equations are obtained. Our special focus in these approaches is to prove the numerical feasibility of such matrix equation based approaches in the PDE context.