Quasi-Monte Carlo and discontinuous Galerkin

Abstract

In this study, we consider the development of tailored quasi-Monte Carlo (QMC) cubatures for non-conforming discontinuous Galerkin (DG) approximations of elliptic partial differential equations (PDEs) with random coefficients. We consider both the affine and uniform and the lognormal models for the input random field, and investigate the use of QMC cubatures to approximate the expected value of the PDE response subject to input uncertainty. In particular, we prove that the resulting QMC convergence rate for DG approximations behaves in the same way as if continuous finite elements were chosen. Notably, the parametric regularity bounds for DG, which are developed in this work, are also useful for other methods such as sparse grids. Numerical results underline our analytical findings.