This paper introduces a simple and computationally tractable probabilistic framework for forward uncertainty quantification based on Gaussian process regression, also known as Kriging. The aim is to equip uncertainty measures in the propagation of input uncertainty to simulator outputs with predictive uncertainty and confidence bounds accounting for the limited accuracy of the surrogate model, which is mainly due to using a finite amount of training data. The additional uncertainty related to the estimation of some of the prior model parameters (namely, trend coefficients and kernel variance) is further accounted for. Two different scenarios are considered. In the first one, the Gaussian process surrogate is used to emulate the actual simulator and propagate input uncertainty in the framework of a Monte Carlo analysis, i.e., as computationally cheap replacement of the original code. In the second one, semianalytical estimates for the statistical moments of the output quantity are obtained directly based on their integral definition. The estimates for the first scenario are more general, more tractable, and they naturally extend to inputs of higher dimensions. The impact of noise on the target function is also discussed. Our findings are demonstrated based on a simple illustrative function and validated by means of several benchmark functions and a high-dimensional test case with more than one hundred uncertain variables.

PROBABILISTIC UNCERTAINTY PROPAGATION USING GAUSSIAN PROCESS SURROGATES / Manfredi, Paolo. - In: INTERNATIONAL JOURNAL FOR UNCERTAINTY QUANTIFICATION. - ISSN 2152-5099. - 14:6(2024), pp. 71-104. [10.1615/int.j.uncertaintyquantification.2024052162]

PROBABILISTIC UNCERTAINTY PROPAGATION USING GAUSSIAN PROCESS SURROGATES

Manfredi, Paolo
2024

Abstract

This paper introduces a simple and computationally tractable probabilistic framework for forward uncertainty quantification based on Gaussian process regression, also known as Kriging. The aim is to equip uncertainty measures in the propagation of input uncertainty to simulator outputs with predictive uncertainty and confidence bounds accounting for the limited accuracy of the surrogate model, which is mainly due to using a finite amount of training data. The additional uncertainty related to the estimation of some of the prior model parameters (namely, trend coefficients and kernel variance) is further accounted for. Two different scenarios are considered. In the first one, the Gaussian process surrogate is used to emulate the actual simulator and propagate input uncertainty in the framework of a Monte Carlo analysis, i.e., as computationally cheap replacement of the original code. In the second one, semianalytical estimates for the statistical moments of the output quantity are obtained directly based on their integral definition. The estimates for the first scenario are more general, more tractable, and they naturally extend to inputs of higher dimensions. The impact of noise on the target function is also discussed. Our findings are demonstrated based on a simple illustrative function and validated by means of several benchmark functions and a high-dimensional test case with more than one hundred uncertain variables.
File in questo prodotto:
File Dimensione Formato  
jnl-2024-IJUQ-GPR.pdf

accesso riservato

Tipologia: 2a Post-print versione editoriale / Version of Record
Licenza: Non Pubblico - Accesso privato/ristretto
Dimensione 1.34 MB
Formato Adobe PDF
1.34 MB Adobe PDF   Visualizza/Apri   Richiedi una copia
manfredi-IJ4UQ-GPR-UQ-final.pdf

embargo fino al 30/05/2025

Tipologia: 2. Post-print / Author's Accepted Manuscript
Licenza: Pubblico - Tutti i diritti riservati
Dimensione 1.26 MB
Formato Adobe PDF
1.26 MB Adobe PDF   Visualizza/Apri   Richiedi una copia
Pubblicazioni consigliate

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11583/2994561