Volume 60, pp. 589-617, 2024.

Quasi-Monte Carlo and Discontinuous Galerkin

Vesa Kaarnioja and Andreas Rupp

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.

Full Text (PDF) [443 KB], BibTeX

Key words

diffusion equation, discontinuous Galerkin, quasi-Monte Carlo, random coefficient

AMS subject classifications

65C05, 65N30

< Back