English

Goal-Oriented Adaptive Finite Element Multilevel Quasi-{M}onte {C}arlo

Numerical Analysis 2025-08-06 v1 Numerical Analysis

Abstract

The efficient approximation of quantity of interest derived from PDEs with lognormal diffusivity is a central challenge in uncertainty quantification. In this study, we propose a multilevel quasi-Monte Carlo framework to approximate deterministic, real-valued, bounded linear functionals that depend on the solution of a linear elliptic PDE with a lognormal diffusivity coefficient {parameterized by a 49-dimensional Gaussian random vector} and deterministic geometric singularities in bounded domains of Rd\mathbb{R}^d. We analyze the parametric regularity and develop the multilevel implementation based on a sequence of adaptive meshes, developed in "Goal-oriented adaptive finite element multilevel Monte Carlo with convergence rates", \emph{CMAME}, 402 (2022), p. 115582. For further variance reduction, we incorporate importance sampling and introduce a level-0 control variate within the multilevel hierarchy. {Introducing such control variate can alter the optimal choice of initial mesh, further highlighting the advantages of adaptive meshes.} Numerical experiments demonstrate that our adaptive QMC algorithm achieves a prescribed accuracy at substantially lower computational cost than the standard multilevel Monte Carlo method.

Keywords

Cite

@article{arxiv.2508.02925,
  title  = {Goal-Oriented Adaptive Finite Element Multilevel Quasi-{M}onte {C}arlo},
  author = {Joakim Beck and Yang Liu and Erik von Schwerin and Raúl Tempone},
  journal= {arXiv preprint arXiv:2508.02925},
  year   = {2025}
}