Abstract:In this paper, we apply quasi-Monte Carlo (QMC) methods with an initial preintegration step to estimate cumulative distribution functions and probability density functions in uncertainty quantification (UQ). The distribution and density functions correspond to a quantity of interest involving the solution to an elliptic partial differential equation (PDE) with a lognormally distributed coefficient and a normally distributed source term. There is extensive previous work on using QMC to compute expected values in UQ, which have proven very successful in tackling a range of different PDE problems. However, the use of QMC for density estimation applied to UQ problems will be explored here for the first time. Density estimation presents a more difficult challenge compared to computing the expected value due to discontinuities present in the integral formulations of both the distribution and density. Our strategy is to use preintegration to eliminate the discontinuity by integrating out a carefully selected random parameter, so that QMC can be used to approximate the remaining integral. First, we establish regularity results for the PDE quantity of interest that are required for smoothing by preintegration to be effective. We then show that an $N$-point lattice rule can be constructed for the integrands corresponding to the distribution and density, such that after preintegration the QMC error is of order $\mathcal{O}(N^{-1+\epsilon})$ for arbitrarily small $\epsilon>0$. This is the same rate achieved for computing the expected value of the quantity of interest. Numerical results are presented to reaffirm our theory.

Uncertainty quantification for random domains using periodic random variables

Quantifying Domain Uncertainty in Linear Elasticity

Stochastic basis adaptation and spatial domain decomposition for PDEs with random coefficients

Quasi-Monte Carlo for partial differential equations with generalized Gaussian input uncertainty

Fast approximation by periodic kernel-based lattice-point interpolation with application in uncertainty quantification

Spectral Convergence of a Semi-discretized Numerical System for the Spatially Homogeneous Boltzmann Equation with Uncertainties

Quasi-Monte Carlo and discontinuous Galerkin

Spline-Based Stochastic Collocation Methods for Uncertainty Quantification in Nonlinear Hyperbolic PDEs

Uncertainty Quantification by Random Measures and Fields

Random time step probabilistic methods for uncertainty quantification in chaotic and geometric numerical integration

A polynomial chaos approach for uncertainty quantification of Monte Carlo transport codes

Density estimation for elliptic PDE with random input by preintegration and quasi-Monte Carlo methods

Efficient Distribution Estimation and Uncertainty Quantification for Elliptic Problems on Domains with Stochastic Boundaries

Total Uncertainty Quantification in Inverse PDE Solutions Obtained with Reduced-Order Deep Learning Surrogate Models

Uncertainty Quantification of RF Circuits Using Stochastic Collocation Techniques

Uncertainty quantification for electrical impedance tomography using quasi-Monte Carlo methods

Uncertainty Quantification for the Homogeneous Landau-Fokker-Planck Equation via Deterministic Particle Galerkin methods

Adaptive Uncertainty Quantification for Stochastic Hyperbolic Conservation Laws

Error Estimate of a Quasi-Monte Carlo Time-Splitting Pseudospectral Method for Nonlinear Schrödinger Equation with Random Potentials

An introduction to uncertainty quantification for kinetic equations and related problems