Special Session 79: 

Quasi-Monte Carlo methods in Uncertainty Quantification

Dirk Nuyens
KU Leuven
Belgium
Co-Author(s):    
Abstract:
The approximation of very high-dimensional integrals is a common challenge in many uncertainty quantification problems. These integrals pop up as expectations of functionals of the solution of a PDE with random coefficients, or they appear when calculating the Bayesian posterior mean of a quantity of interest in an inverse problem. Sometimes the uncertainties can be modelled by very smooth random fields and then polynomial chaos expansions (PCE) methods can be used to approximate the uncertainty with very few terms (say 10), however, when the needed dimensionality is really high (and we are talking thousands or millions), e.g., for a rough random field with possibly high variance, the PCE method will have problems. In such a context quasi-Monte Carlo (QMC) methods can deliver a solution. Sometimes giving dimension-independent convergence of $1/N$ or even higher algebraic rates.