A quasi-Monte Carlo Method for an Optimal Control Problem Under Uncertainty

by   Philipp A. Guth, et al.

We study an optimal control problem under uncertainty, where the target function is the solution of an elliptic partial differential equation with random coefficients, steered by a control function. The robust formulation of the optimization problem is stated as a high-dimensional integration problem over the stochastic variables. It is well known that carrying out a high-dimensional numerical integration of this kind using a Monte Carlo method has a notoriously slow convergence rate; meanwhile, a faster rate of convergence can potentially be obtained by using sparse grid quadratures, but these lead to discretized systems that are non-convex due to the involvement of negative quadrature weights. In this paper, we analyze instead the application of a quasi-Monte Carlo method, which retains the desirable convexity structure of the system and has a faster convergence rate compared to ordinary Monte Carlo methods. In particular, we show that under moderate assumptions on the decay of the input random field, the error rate obtained by using a specially designed, randomly shifted rank-1 lattice quadrature rule is essentially inversely proportional to the number of quadrature nodes. The overall discretization error of the problem, consisting of the dimension truncation error, finite element discretization error and quasi-Monte Carlo quadrature error, is derived in detail. We assess the theoretical findings in numerical experiments.


Parabolic PDE-constrained optimal control under uncertainty with entropic risk measure using quasi-Monte Carlo integration

We study the application of a tailored quasi-Monte Carlo (QMC) method to...

The Vlasov-Fokker-Planck Equation with High Dimensional Parametric Forcing Term

We consider the Vlasov-Fokker-Planck equation with random electric field...

Speeding up Monte Carlo Integration: Control Neighbors for Optimal Convergence

A novel linear integration rule called control neighbors is proposed in ...

Random geometries and Quasi Monte Carlo methods for optimal control PDE problems based on fictitious domain FEMS and cut elements

This work investigates an elliptic optimal control problem defined on un...

Quasi-Monte Carlo finite element approximation of the Navier-Stokes equations with initial data modeled by log-normal random fields

In this paper, we analyze the numerical approximation of the Navier-Stok...

On dropping the first Sobol' point

Quasi-Monte Carlo (QMC) points are a substitute for plain Monte Carlo (M...

Constrained Ensemble Langevin Monte Carlo

The classical Langevin Monte Carlo method looks for i.i.d. samples from ...

Please sign up or login with your details

Forgot password? Click here to reset