Efficient Adaptive Stochastic Galerkin Methods for Parametric Operator Equations

Abstract

This paper is concerned with the design and implementation of efficient solution algorithms for elliptic PDE problems with correlated random data. The energy orthogonality that is built into stochastic Galerkin approximations is cleverly exploited to give an innovative energy error estimation strategy that utilizes the tensor product structure of the approximation space. An associated error estimator is constructed and shown theoretically and numerically to be an effective mechanism for driving an adaptive refinement process. The codes used in the numerical studies are available online.