Goal-Oriented Adaptive Finite Element Multilevel Quasi-{M}onte {C}arlo
By: Joakim Beck , Yang Liu , Erik von Schwerin and more
Potential Business Impact:
Makes computer models of science problems faster.
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 $\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.
Similar Papers
Exploiting Inexact Computations in Multilevel Sampling Methods
Numerical Analysis
Makes computer simulations run much faster.
Multilevel lattice-based kernel approximation for elliptic PDEs with random coefficients
Numerical Analysis
Speeds up computer math for science problems.
Sufficient conditions for QMC analysis of finite elements for parametric differential equations
Numerical Analysis
Makes computer simulations of science problems faster.