Sparse Approximation of Data-Driven Polynomial Chaos Expansions: An Induced Sampling Approach

Abstract

One of the open problems in the field of forward uncertainty quantification (UQ) is the ability to form accurate assessments of uncertainty having only incomplete information about the distribution of random inputs. Another challenge is to efficiently make use of limited training data for UQ predictions of complex engineering problems, particularly with high dimensional random parameters. We address these challenges by combining data-driven polynomial chaos expansions with a recently developed preconditioned sparse approximation approach for UQ problems. The first task in this two-step process is to employ the procedure developed in (Ahlfeld et al. 2016, [1]) to construct an "arbitrary" polynomial chaos expansion basis using a finite number of statistical moments of the random inputs. The second step is a novel procedure to effect sparse approximation via 1 minimization in order to quantify the forward uncertainty. To enhance the performance of the preconditioned 1 minimization problem, we sample from the so-called induced distribution, instead of using Monte Carlo (MC) sampling from the original, unknown probability measure. We demonstrate on test problems that induced sampling is a competitive and often better choice compared with sampling from asymptotically optimal measures (such as the equilibrium measure) when we have incomplete information about the distribution. We demonstrate the capacity of the proposed induced sampling algorithm via sparse representation with limited data on test functions, and on a Kirchoff plating bending problem with random Young's modulus.

Optimal design for kernel interpolation: Applications to uncertainty quantification

“…"Non-intrusive" effectively means that existing black-box tools can be used in their current form. Most existing work concerns the collocation based polynomial approximations, this include the interpolation technic based on sparse grid [9,21,23,35], the least-squares projection onto polynomial spaces [3,11,22,13], the compressive sampling method with 1 minimization [5,16,12,37]. Another approach to construct the surrogate using D N is the so called Gaussian process (GP) regression [2,1,26,31].…”