In this work, we analyze Parametrized Advection-Dominated distributed Optimal Control Problems with random inputs in a Reduced Order Model (ROM) context. All the simulations are initially based on a finite element method (FEM) discretization; moreover, a space-time approach is considered when dealing with unsteady cases. To overcome numerical instabilities that can occur in the optimality system for high values of the Péclet number, we consider a Streamline Upwind Petrov-Galerkin technique applied in an optimize-then-discretize approach. We combine this method with the ROM framework in order to consider two possibilities of stabilization: Offline-Only stabilization and Offline-Online stabilization. Moreover we consider random parameters and we use a weighted Proper Orthogonal Decomposition algorithm in a partitioned approach to deal with the issue of uncertainty quantification. Several quadrature techniques are used to derive weighted ROMs: tensor rules, isotropic sparse grids, Monte-Carlo and quasi Monte-Carlo methods. We compare all the approaches analyzing relative errors between the FEM and ROM solutions and the computational efficiency based on the speedup-index.

In this paper we will consider distributed Linear-Quadratic Optimal Control Problems dealing with Advection-Diffusion PDEs for high values of the Péclet number. In this situation, computational instabilities occur, both for steady and unsteady cases. A Streamline Upwind Petrov-Galerkin technique is used in the optimality system to overcome these unpleasant effects. We will apply a finite element method discretization in a optimize-then-discretize approach. For the parabolic case, a space-time framework will be considered and stabilization will also occur in the bilinear forms involving time derivatives. Then we will build Reduced Order Models on this discretization procedure and two possible settings can be analyzed: whether or not stabilization is needed in the online phase, too. In order to build the reduced bases for state, control, and adjoint variables we will consider a Proper Orthogonal Decomposition algorithm in a partitioned approach. The discussion is supported by computational experiments, where relative errors between the FEM and ROM solutions are studied together with the respective computational times.