A weighted POD-reduction approach for parametrized PDE-constrained Optimal Control Problems with random inputs and applications to environmental sciences

Reduced basis approximations of Optimal Control Problems (OCPs) governed by steady partial differential equations (PDEs) with random parametric inputs are analyzed and constructed. Such approximations are based on a Reduced Order Model, which in this work is constructed using the method of weighted Proper Orthogonal Decomposition. This Reduced Order Model then is used to efficiently compute the reduced basis approximation for any outcome of the random parameter. We demonstrate that such OCPs are well-posed by applying the adjoint approach, which also works in the presence of admissibility constraints and in the case of non linear-quadratic OCPs, and thus is more general than the conventional Lagrangian approach. We also show that a step in the construction of these Reduced Order Models, known as the aggregation step, is not fundamental and can in principle be skipped for noncoercive problems, leading to a cheaper online phase. Numerical applications in three scenarios from environmental science are considered, in which the governing PDE is steady and the control is distributed. Various parameter distributions are taken, and several implementations of the weighted Proper Orthogonal Decomposition are compared by choosing different quadrature rules.

[1]  Anders Logg,et al.  Automated Solution of Differential Equations by the Finite Element Method: The FEniCS Book , 2012 .

[2]  Max Gunzburger,et al.  POD and CVT-based reduced-order modeling of Navier-Stokes flows , 2006 .

[3]  J. Hesthaven,et al.  Certified Reduced Basis Methods for Parametrized Partial Differential Equations , 2015 .

[4]  Gianluigi Rozza,et al.  RBniCS - reduced order modelling in FEniCS , 2015 .

[5]  Gianluigi Rozza,et al.  Reduced Basis Methods for Uncertainty Quantification , 2017, SIAM/ASA J. Uncertain. Quantification.

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

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

[8]  Gianluigi Rozza,et al.  Supremizer stabilization of POD–Galerkin approximation of parametrized steady incompressible Navier–Stokes equations , 2015 .

[9]  Christoph Schwab,et al.  Karhunen-Loève approximation of random fields by generalized fast multipole methods , 2006, J. Comput. Phys..

[10]  G. Rozza,et al.  On the stability of the reduced basis method for Stokes equations in parametrized domains , 2007 .

[11]  Traian Iliescu,et al.  B-spline based finite-element method for the stationary quasi-geostrophic equations of the ocean , 2015 .

[12]  Mark Kärcher,et al.  A Certified Reduced Basis Approach for Parametrized Optimal Control Problems , 2012 .

[13]  Gianluigi Rozza,et al.  Reduced basis approximation and a posteriori error estimation for Stokes flows in parametrized geometries: roles of the inf-sup stability constants , 2013, Numerische Mathematik.

[14]  Fabio Cavallini,et al.  Quasi-Geostrophic Theory of Oceans and Atmosphere: Topics in the Dynamics and Thermodynamics of the Fluid Earth , 2012 .

[15]  Gianluigi Rozza,et al.  Model Reduction for Parametrized Optimal Control Problems in Environmental Marine Sciences and Engineering , 2017, SIAM J. Sci. Comput..

[16]  Dominique Chapelle,et al.  A Galerkin strategy with Proper Orthogonal Decomposition for parameter-dependent problems – Analysis, assessments and applications to parameter estimation , 2013 .

[17]  D. Rovas,et al.  Reliable Real-Time Solution of Parametrized Partial Differential Equations: Reduced-Basis Output Bound Methods , 2002 .

[18]  Bernard Haasdonk,et al.  Chapter 2: Reduced Basis Methods for Parametrized PDEs—A Tutorial Introduction for Stationary and Instationary Problems , 2017 .

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

[20]  J. Lions.  Optimal Control of Systems Governed by Partial Differential Equations , 1971 .

[21]  Luca Dedè,et al.  Reduced Basis Method and A Posteriori Error Estimation for Parametrized Linear-Quadratic Optimal Control Problems , 2010, SIAM J. Sci. Comput..

[22]  Karen Veroy,et al.  Certified Reduced Basis Methods for Parametrized Elliptic Optimal Control Problems with Distributed Controls , 2017, Journal of Scientific Computing.

[23]  A. Patera,et al.  Reduced basis approximation and a posteriori error estimation for affinely parametrized elliptic coercive partial differential equations , 2007 .

[24]  Sören Bartels,et al.  Numerical Approximation of Partial Differential Equations , 2016 .

[25]  Università degli Studi di Trieste Dipartimento di Matematica e Geoscienze Corso di Laurea Magistrale in Matematica WEIGHTED REDUCED ORDER METHODS FOR PARAMETRIZED PDEs IN UNCERTAINTY QUANTIFICATION PROBLEMS , 2017 .

[26]  A. Patera,et al.  Reduced basis approximation and a posteriori error estimation for affinely parametrized elliptic coercive partial differential equations , 2007 .

[27]  Mark Kärcher,et al.  A certified reduced basis method for parametrized elliptic optimal control problems , 2014 .

[28]  N. Nguyen,et al.  An ‘empirical interpolation’ method: application to efficient reduced-basis discretization of partial differential equations , 2004 .

[29]  T. Sullivan.  Introduction to Uncertainty Quantification , 2015 .

[30]  A. Quarteroni,et al.  Numerical Approximation of Partial Differential Equations , 2008 .

[31]  Gianluigi Rozza,et al.  Stabilized Weighted Reduced Basis Methods for Parametrized Advection Dominated Problems with Random Inputs , 2017, SIAM/ASA J. Uncertain. Quantification.

[32]  Edriss S. Titi,et al.  Existence of solutions to the Stommel-Charney model of the Gulf Stream , 1988 .

[33]  M. Griebel,et al.  Singular value decomposition versus sparse grids: refined complexity estimates , 2018, IMA Journal of Numerical Analysis.

[34]  Gianluigi Rozza,et al.  A Weighted POD Method for Elliptic PDEs with Random Inputs , 2018, J. Sci. Comput..

[35]  G. Rozza,et al.  Stabilized Weighted Reduced Basis Methods for Parametrized Advection Dominated Problems with Random Inputs , 2017, SIAM/ASA J. Uncertain. Quantification.