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 | Dimensione | Formato | |
---|---|---|---|
jnl-2024-IJUQ-GPR.pdf
non disponibili
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.
https://hdl.handle.net/11583/2994561