Robust Optimization of PDEs with Random Coefficients Using a Multilevel Monte Carlo Method

11/07/2017
by   Andreas Van Barel, et al.
0

This paper addresses optimization problems constrained by partial differential equations with uncertain coefficients. In particular, the robust control problem and the average control problem are considered for a tracking type cost functional with an additional penalty on the variance of the state. The expressions for the gradient and Hessian corresponding to either problem contain expected value operators. Due to the large number of uncertainties considered in our model, we suggest to evaluate these expectations using a multilevel Monte Carlo (MLMC) method. Under mild assumptions, it is shown that this results in the gradient and Hessian corresponding to the MLMC estimator of the original cost functional. Furthermore, we show that the use of certain correlated samples yields a reduction in the total number of samples required. Two optimization methods are investigated: the nonlinear conjugate gradient method and the Newton method. For both, a specific algorithm is provided that dynamically decides which and how many samples should be taken in each iteration. The cost of the optimization up to some specified tolerance τ is shown to be proportional to the cost of a gradient evaluation with requested root mean square error τ. The algorithms are tested on a model elliptic diffusion problem with lognormal diffusion coefficient. An additional nonlinear term is also considered.

READ FULL TEXT

page 5

page 23

research
09/29/2021

Multilevel Quasi-Monte Carlo for Optimization under Uncertainty

This paper considers the problem of optimizing the average tracking erro...
research
05/29/2020

MG/OPT and MLMC for Robust Optimization of PDEs

An algorithm is proposed to solve robust control problems constrained by...
research
08/31/2020

Algorithms for Solving High Dimensional PDEs: From Nonlinear Monte Carlo to Machine Learning

In recent years, tremendous progress has been made on numerical algorith...
research
03/03/2021

Full history recursive multilevel Picard approximations for ordinary differential equations with expectations

We consider ordinary differential equations (ODEs) which involve expecta...
research
02/01/2023

Multilevel Monte Carlo FEM for Elliptic PDEs with Besov Random Tree Priors

We develop a multilevel Monte Carlo (MLMC)-FEM algorithm for linear, ell...
research
11/14/2019

Space-time multilevel Monte Carlo methods and their application to cardiac electrophysiology

We present a novel approach aimed at high-performance uncertainty quanti...
research
07/29/2019

Enhanced Multi-Index Monte Carlo by means of Multiple Semi-Coarsened Multigrid for Anisotropic Diffusion Problems

In many models used in engineering and science, material properties are ...

Please sign up or login with your details

Forgot password? Click here to reset