Marburg University

Almost all real-world applications involve a degree of uncertainty. This may be the result of noisy measurements, restrictions on observability, or simply unforeseen events. Since many models in both engineering and the natural sciences make use of partial differential equations (PDEs), it is natural to consider PDEs with random inputs. In this context, passing from modelling and simulation to optimization or control results in stochastic PDE-constrained optimization problems. This leads to a number of theoretical, algorithmic, and numerical challenges.

 From a mathematical standpoint, the solution of the underlying PDE is a random field, which in turn makes the quantity of interest or the objective function an implicitly defined random variable. In order to minimize this distributed objective, one can use, e.g., stochastic order constraints, a distributionally robust approach, or risk measures. In this talk, we will make use of risk measures.

After motivating the approach via a model for the mitigation of an airborne pollutant, we build up an analytical framework and introduce some useful risk measures. This allows us to prove the existence of solutions and derive optimality conditions. We then present several approximation schemes for handling non-smooth risk measures in order to leverage existing numerical methods from PDE-constrained optimization. Finally, we discuss solutions techniques and illustrate our results with numerical examples.

The notion quantization originates from information theory, where it refers to the approximation of a continuous signal on a discrete set. Our research on quantization is mainly motivated by applications in quadrature problems. In that context, one aims at finding for a given probability measure $\mu$ on a metric space a discrete approximation that is supported on a finite number of points, say $N$, and is close to $\mu$ in a Wasserstein metric.

In general it is a hard problem to find close to optimal quantizations, if  $N$ is large and/or  $\mu$ is given implicitly, e.g. being the marginal distribution of a stochastic differential equation. In this talk we analyse the efficiency of empirical measures in the constructive quantization problem. That means the random approximating measure is the uniform distribution on $N$ independent $\mu$-distributed elements.

We show that this approach is order order optimal in many cases. Further, we give fine asymptotic estimates for the quantization error that involve moments of the density of the absolutely continuous part of $\mu$, so called high resolution formulas. The talk ends with an outlook on possible applications and open problems.

The talk is based on joint work with Michael Scheutzow (TU Berlin) and Reik Schottstedt (U Marburg).