Polynomial chaos expansions for dependent random variables

Abstract

Polynomial chaos expansions (PCE) are well-suited to quantifying uncertainty in models parameterized by independent random variables. The assumption of independence leads to simple strategies for evaluating PCE coefficients. In contrast, the application of PCE to models of dependent variables is much more challenging. Three approaches can be used. The first approach uses mapping methods where measure transformations, such as the Nataf and Rosenblatt transformation, can be used to map dependent random variables to independent ones; however we show that this can significantly degrade performance since the Jacobian of the map must be approximated. A second strategy is the class of dominating support methods which build PCE using independent random variables whose distributional support dominates the support of the true dependent joint density; we provide evidence that this approach appears to produce approximations with suboptimal accuracy. A third approach, the novel method proposed here, uses Gram-Schmidt orthogonalization (GSO) to numerically compute orthonormal polynomials for the dependent random variables. This approach has been used successfully when solving differential equations using the intrusive stochastic Galerkin method, and in this paper we use GSO to build PCE using a non-intrusive stochastic collocation method. The stochastic collocation method treats the model as a black box and builds approximations of model output from a set of samples. Building PCE from samples can introduce ill-conditioning which does not plague stochastic Galerkin methods. To mitigate this ill-conditioning we generate weighted Leja sequences, which are nested sample sets, to build accurate polynomial interpolants. We show that our proposed approach produces PCE which are orders of magnitude more accurate than PCE constructed using mapping or dominating support methods.

27 pages, 10 figures

References in corpus (1)

• Compressive Sampling of Polynomial Chaos Expansions: Convergence Analysis and Sampling Strategies

Cited by in corpus (14)

• Adaptive Multi-index Collocation for Uncertainty Quantification and Sensitivity Analysis
• Emulation of stochastic simulators using generalized lambda models
• Manifold learning-based polynomial chaos expansions for high-dimensional surrogate models
• Multifidelity uncertainty quantification with models based on dissimilar parameters
• Multivariate sensitivity-adaptive polynomial chaos expansion for high-dimensional surrogate modeling and uncertainty quantification
• An $hp$-adaptive multi-element stochastic collocation method for surrogate modeling with information re-use
• Global Sensitivity Analysis on Numerical Solver Parameters of Particle-In-Cell Models in Particle Accelerator Systems
• Adaptive Sparse Polynomial Chaos Expansions via Leja Interpolation
• Global sensitivity analysis with limited data via sparsity-promoting D-MORPH regression: Application to char combustion
• Polynomial Chaos Expansions on Principal Geodesic Grassmannian Submanifolds for Surrogate Modeling and Uncertainty Quantification
• High-Dimensional Uncertainty Quantification of Electronic and Photonic IC with Non-Gaussian Correlated Process Variations
• A Resolution Independent Neural Operator
• Neural Chaos: A Spectral Stochastic Neural Operator
• Embedded Model Form Uncertainty Quantification with Measurement Noise for Bayesian Model Calibration