Convergence and rate optimality of adaptive multilevel stochastic Galerkin FEM

Alex Bespalov, Dirk Praetorius, Michele Ruggeri

Research output: Contribution to journalArticlepeer-review

32 Downloads (Pure)


We analyze an adaptive algorithm for the numerical solution of parametric elliptic partial differential equations in two-dimensional physical domains, with coefficients and right-hand-side functions depending on infinitely many (stochastic) parameters. The algorithm generates multilevel stochastic Galerkin approximations; these are represented in terms of a sparse generalized polynomial chaos expansion with coefficients residing in finite element spaces associated with different locally refined meshes. Adaptivity is driven by a two-level a posteriori error estimator and employs a Dörfler-type marking on the joint set of spatial and parametric error indicators. We show that, under an appropriate saturation assumption, the proposed adaptive strategy yields optimal convergence rates with respect to the overall dimension of the underlying multilevel approximation spaces.
Original languageEnglish
JournalI M A Journal of Numerical Analysis
Early online date19 May 2021
Publication statusE-pub ahead of print - 19 May 2021


Dive into the research topics of 'Convergence and rate optimality of adaptive multilevel stochastic Galerkin FEM'. Together they form a unique fingerprint.

Cite this