Library Subscription: Guest
Begell Digital Portal Begell Digital Library eBooks Journals References & Proceedings Research Collections
International Journal for Uncertainty Quantification
IF: 0.967 5-Year IF: 1.301

ISSN Print: 2152-5080
ISSN Online: 2152-5099

Open Access

International Journal for Uncertainty Quantification

DOI: 10.1615/Int.J.UncertaintyQuantification.2013006821
pages 63-93

DIMENSIONALITY REDUCTION FOR COMPLEX MODELS VIA BAYESIAN COMPRESSIVE SENSING

Khachik Sargsyan
P.O.Box 969, MS 9051; Sandia National Laboratories, Livermore, California 94551, USA
Cosmin Safta
P.O.Box 969, MS 9051; Sandia National Laboratories, Livermore, California 94551, USA
Habib N. Najm
Sandia National Laboratories, Livermore, CA, 94551
Bert J. Debusschere
P.O.Box 969, MS 9051; Sandia National Laboratories, Livermore, California 94551, USA
Daniel Ricciuto
Oak Ridge National Laboratory, 1 Bethel Valley Road, MS 6301, Oak Ridge, Tennessee 37831, USA
Peter Thornton
Oak Ridge National Laboratory, 1 Bethel Valley Road, MS 6301, Oak Ridge, Tennessee 37831, USA

ABSTRACT

Uncertainty quantification in complex physical models is often challenged by the computational expense of these models. One often needs to operate under the assumption of sparsely available model simulations. This issue is even more critical when models include a large number of input parameters. This "curse of dimensionality", in particular, leads to a prohibitively large number of basis terms in spectral methods for uncertainty quantification, such as polynomial chaos (PC) methods. In this work, we implement a PC-based surrogate model construction that "learns" and retains only the most relevant basis terms of the PC expansion, using sparse Bayesian learning. This dramatically reduces the dimensionality of the problem, making it more amenable to further analysis such as sensitivity or calibration studies. The model of interest is the community land model with about 80 input parameters, which also exhibits nonsmooth input-output behavior. We enhanced the methodology by a clustering and classifying procedure that leads to a piecewise-PC surrogate thereby dealing with nonlinearity. We then obtain global sensitivity information for five outputs with respect to all input parameters using less than 10,000 model simulations−a very small number for an 80-dimensional input parameter space.