Generalized dimension truncation error analysis for high-dimensional numerical integration: lognormal setting and beyond

by   Philipp A. Guth, et al.

Partial differential equations (PDEs) with uncertain or random inputs have been considered in many studies of uncertainty quantification. In forward uncertainty quantification, one is interested in analyzing the stochastic response of the PDE subject to input uncertainty, which usually involves solving high-dimensional integrals of the PDE output over a sequence of stochastic variables. In practical computations, one typically needs to discretize the problem in several ways: approximating an infinite-dimensional input random field with a finite-dimensional random field, spatial discretization of the PDE using, e.g., finite elements, and approximating high-dimensional integrals using cubatures such as quasi-Monte Carlo methods. In this paper, we focus on the error resulting from dimension truncation of an input random field. We show how Taylor series can be used to derive theoretical dimension truncation rates for a wide class of problems and we provide a simple checklist of conditions that a parametric mathematical model needs to satisfy in order for our dimension truncation error bound to hold. Some of the novel features of our approach include that our results are applicable to non-affine parametric operator equations, dimensionally-truncated conforming finite element discretized solutions of parametric PDEs, and even compositions of PDE solutions with smooth nonlinear quantities of interest. As a specific application of our method, we derive an improved dimension truncation error bound for elliptic PDEs with lognormally parameterized diffusion coefficients. Numerical examples support our theoretical findings.


page 1

page 2

page 3

page 4


Application of dimension truncation error analysis to high-dimensional function approximation

Parametric mathematical models such as partial differential equations wi...

Analyticity and sparsity in uncertainty quantification for PDEs with Gaussian random field inputs

We establish summability results for coefficient sequences of Wiener-Her...

A Dimension-adaptive Combination Technique for Uncertainty Quantification

We present an adaptive algorithm for the computation of quantities of in...

A survey of unsupervised learning methods for high-dimensional uncertainty quantification in black-box-type problems

Constructing surrogate models for uncertainty quantification (UQ) on com...

Wavenumber-explicit parametric holomorphy of Helmholtz solutions in the context of uncertainty quantification

A crucial role in the theory of uncertainty quantification (UQ) of PDEs ...

Quasi-Monte Carlo and discontinuous Galerkin

In this study, we design and develop Quasi-Monte Carlo (QMC) cubatures f...

Please sign up or login with your details

Forgot password? Click here to reset