Multilevel Monte Carlo Analysis for Optimal Control of Elliptic PDEs with Random Coefficients

This work is motivated by the need to study the impact of data uncertainties and material imperfections on the solution to optimal control problems constrained by partial differential equations. We consider a pathwise optimal control problem constrained by a diffusion equation with random coefficient together with box constraints for the control. For each realization of the diffusion coefficient we solve an optimal control problem using the variational discretization [M. Hinze, Comput. Optim. Appl., 30 (2005), pp. 45-61]. Our framework allows for lognormal coefficients whose realizations are not uniformly bounded away from zero and infinity. We establish finite element error bounds for the pathwise optimal controls. This analysis is nontrivial due to the limited spatial regularity and the lack of uniform ellipticity and boundedness of the diffusion operator. We apply the error bounds to prove convergence of a multilevel Monte Carlo estimator for the expected value of the pathwise optimal controls. In addition we analyze the computational complexity of the multilevel estimator. We perform numerical experiments in 2D space to confirm the convergence result and the complexity bound.

[1]  Jangwoon Lee,et al.  A Stochastic Galerkin Method for Stochastic Control Problems , 2012 .

[2]  L. S. Hou,et al.  Finite element approximations of stochastic optimal control problems constrained by stochastic elliptic PDEs , 2011 .

[3]  Matthias Heinkenschloss,et al.  Inexact Objective Function Evaluations in a Trust-Region Algorithm for PDE-Constrained Optimization under Uncertainty , 2014, SIAM J. Sci. Comput..

[4]  Gianluigi Rozza,et al.  Multilevel and weighted reduced basis method for stochastic optimal control problems constrained by Stokes equations , 2016, Numerische Mathematik.

[5]  J. Zabczyk,et al.  Stochastic Equations in Infinite Dimensions , 2008 .

[6]  Alfio Borzì.  Multigrid and sparse-grid schemes for elliptic control problems with random coefficients , 2010, Comput. Vis. Sci..

[7]  Elisabeth Ullmann,et al.  Mixed finite element analysis of lognormal diffusion and multilevel Monte Carlo methods , 2013, 1312.6047.

[8]  Hehu Xie,et al.  A Multilevel Correction Method for Optimal Controls of Elliptic Equations , 2014, SIAM J. Sci. Comput..

[9]  Zdzisław Denkowski,et al.  Set-Valued Analysis , 2021 .

[10]  Gianluigi Rozza,et al.  Stochastic Optimal Robin Boundary Control Problems of Advection-Dominated Elliptic Equations , 2013, SIAM J. Numer. Anal..

[11]  Aaas News,et al.  Book Reviews , 1893, Buffalo Medical and Surgical Journal.

[12]  Stefan Ulbrich,et al.  Optimization with PDE Constraints , 2008, Mathematical modelling.

[13]  R. Ghanem,et al.  Stochastic Finite Elements: A Spectral Approach , 1990 .

[14]  Michael Hinze,et al.  A Globalized Semi-smooth Newton Method for Variational Discretization of Control Constrained Elliptic Optimal Control Problems , 2012, Constrained Optimization and Optimal Control for Partial Differential Equations.

[15]  Angela Kunoth,et al.  Analytic Regularity and GPC Approximation for Control Problems Constrained by Linear Parametric Elliptic and Parabolic PDEs , 2013, SIAM J. Control. Optim..

[16]  K. A. Cliffe,et al.  Multilevel Monte Carlo methods and applications to elliptic PDEs with random coefficients , 2011, Comput. Vis. Sci..

[17]  Andrea Barth,et al.  Multi-level Monte Carlo Finite Element method for elliptic PDEs with stochastic coefficients , 2011, Numerische Mathematik.

[18]  Max D. Gunzburger,et al.  Optimal Control of Stochastic Flow over a Backward-Facing Step Using Reduced-Order Modeling , 2011, SIAM J. Sci. Comput..

[19]  Catherine E. Powell,et al.  An Introduction to Computational Stochastic PDEs , 2014 .

[20]  Drew P. Kouri,et al.  A Multilevel Stochastic Collocation Algorithm for Optimization of PDEs with Uncertain Coefficients , 2014, SIAM/ASA J. Uncertain. Quantification.

[21]  Michael Ulbrich,et al.  Constrained Optimization with Low-Rank Tensors and Applications to Parametric Problems with PDEs , 2017, SIAM J. Sci. Comput..

[22]  Wolfgang Hackbusch,et al.  Elliptic Differential Equations: Theory and Numerical Treatment , 2017 .

[23]  Jangwoon Lee,et al.  Error Estimates of Stochastic Optimal Neumann Boundary Control Problems , 2011, SIAM J. Numer. Anal..

[24]  Juan A. Lazzús.  Hybrid Particle Swarm-Ant Colony Algorithm to Describe the Phase Equilibrium of Systems Containing Supercritical Fluids with Ionic Liquids , 2013 .

[25]  Enrique Zuazua,et al.  Averaged control and observation of parameter-depending wave equations , 2014 .

[26]  Robert Scheichl,et al.  Finite Element Error Analysis of Elliptic PDEs with Random Coefficients and Its Application to Multilevel Monte Carlo Methods , 2013, SIAM J. Numer. Anal..

[27]  Fabio Nobile,et al.  A Stochastic Collocation Method for Elliptic Partial Differential Equations with Random Input Data , 2007, SIAM Rev..

[28]  Alfio Borzì,et al.  On the treatment of distributed uncertainties in PDE‐constrained optimization , 2010 .

[29]  Alfio Borzì,et al.  A POD framework to determine robust controls in PDE optimization , 2011, Comput. Vis. Sci..

[30]  T. J. Dodwell,et al.  A Hierarchical Multilevel Markov Chain Monte Carlo Algorithm with Applications to Uncertainty Quantification in Subsurface Flow , 2013, SIAM/ASA J. Uncertain. Quantification.

[31]  Peter Benner,et al.  Low-rank solvers for unsteady Stokes–Brinkman optimal control problem with random data , 2016, Computer Methods in Applied Mechanics and Engineering.

[32]  Michael Hinze,et al.  A Variational Discretization Concept in Control Constrained Optimization: The Linear-Quadratic Case , 2005, Comput. Optim. Appl..

[33]  Arnd Rösch,et al.  Error estimates for linear-quadratic control problems with control constraints , 2006, Optim. Methods Softw..

[34]  Julia Charrier,et al.  Strong and Weak Error Estimates for Elliptic Partial Differential Equations with Random Coefficients , 2012, SIAM J. Numer. Anal..

[35]  Alfio Borzì,et al.  Multigrid Methods and Sparse-Grid Collocation Techniques for Parabolic Optimal Control Problems with Random Coefficients , 2009, SIAM J. Sci. Comput..

[36]  Bart G. van Bloemen Waanders,et al.  A Trust-Region Algorithm with Adaptive Stochastic Collocation for PDE Optimization under Uncertainty , 2012, SIAM J. Sci. Comput..

[37]  Peter Benner,et al.  Block-Diagonal Preconditioning for Optimal Control Problems Constrained by PDEs with Uncertain Inputs , 2016, SIAM J. Matrix Anal. Appl..

[38]  Elisabeth Ullmann,et al.  Further analysis of multilevel Monte Carlo methods for elliptic PDEs with random coefficients , 2012, Numerische Mathematik.

[39]  Gianluigi Rozza,et al.  Reduced basis approximation of parametrized optimal flow control problems for the Stokes equations , 2015, Comput. Math. Appl..

[40]  Klaus Ritter,et al.  Multilevel Monte Carlo Approximation of Distribution Functions and Densities , 2015, SIAM/ASA J. Uncertain. Quantification.

[41]  Ralf Kornhuber,et al.  Multilevel Monte Carlo Finite Element Methods for Stochastic Elliptic Variational Inequalities , 2014, SIAM J. Numer. Anal..

[42]  Garth N. Wells,et al.  Optimal control with stochastic PDE constraints and uncertain controls , 2011, ArXiv.

[43]  Robert Michael Kirby,et al.  Stochastic Collocation for Optimal Control Problems with Stochastic PDE Constraints , 2012, SIAM J. Control. Optim..

[44]  Alexey Chernov,et al.  Convergence analysis of multilevel Monte Carlo variance estimators and application for random obstacle problems , 2015, Numerische Mathematik.

[45]  Alfio Quarteroni,et al.  Weighted Reduced Basis Method for Stochastic Optimal Control Problems with Elliptic PDE Constraint , 2014, SIAM/ASA J. Uncertain. Quantification.

[46]  Gianluigi Rozza,et al.  Reduced Basis Method for Parametrized Elliptic Optimal Control Problems , 2013, SIAM J. Sci. Comput..

[47]  Guannan Zhang,et al.  Stochastic finite element methods for partial differential equations with random input data* , 2014, Acta Numerica.

[48]  Enrique Zuazua,et al.  Averaged control , 2014, Autom..