A Trust-Region Algorithm with Adaptive Stochastic Collocation for PDE Optimization under Uncertainty

D. P. Kouri
Mathematics and Computer Science Division
Argonne National Laboratory

Matthias Heinkenschloss
Department of Computational and Applied Mathematics
Rice University

D. Ridzal
Optimization and Uncertainty Quantification
Sandia National Laboratories

B. G van Bloemen Waanders
Numerical Analysis and Applications
Sandia National Laboratories


SIAM Journal on Scientific Computing, 2013, Vol. 35, No. 4, pp. A1847-A1879.

Abstract

The numerical solution of optimization problems governed by partial differential equations (PDEs) with random coefficients is computationally challenging because of the large number of deterministic PDE solves required at each optimization iteration. This paper introduces an efficient algorithm for solving such problems based on a combination of adaptive sparse-grid collocation for the discretization of the PDE in the stochastic space and a trust-region framework for optimization and fidelity management of the stochastic discretization. The overall algorithm adapts the collocation points based on the progress of the optimization algorithm and the impact of the random variables on the solution of the optimization problem. It frequently uses few collocation points initially and increases the number of collocation points only as necessary, thereby keeping the number of deterministic PDE solves low while guaranteeing convergence. Currently an error indicator is used to estimate gradient errors due to adaptive stochastic collocation. The algorithm is applied to three examples, and the numerical results demonstrate a significant reduction in the total number of PDE solves required to obtain an optimal solution when compared with a Newton conjugate gradient algorithm applied to a fixed high-fidelity discretization of the optimization problem.

Keywords. PDE optimization, uncertainty, stochastic collocation, trust regions, sparse grids, adaptivity