English

Quasi-Monte Carlo and discontinuous Galerkin

Numerical Analysis 2024-12-12 v6 Numerical Analysis

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.

Keywords

Cite

@article{arxiv.2207.07698,
  title  = {Quasi-Monte Carlo and discontinuous Galerkin},
  author = {Vesa Kaarnioja and Andreas Rupp},
  journal= {arXiv preprint arXiv:2207.07698},
  year   = {2024}
}