INVERSE UNCERTAINTY QUANTIFICATION OF CTF PHYSICAL MODEL PARAMETERS USING BAYESIAN INFERENCE

Junhong Liu,Shanfang Huang,Xiaoyu Guo,Jiageng Wang,Kan Wang
DOI: https://doi.org/10.1299/jsmeicone.2019.27.1435
2019-01-01
The Proceedings of the International Conference on Nuclear Engineering (ICONE)
Abstract:This work illustrates a quantitative uncertainty analysis of CTF, one thermal-hydraulics sub-channel code for nuclear engineering applications. Like many other thermal-hydraulics codes, large uncertainties of CTF stem from the inaccuracy of parameters in the empirical correlations. However, many previous uncertainty analyses of CTF are based on the inexact parameter distributions described by the “expert opinions”, which are generally overestimations. To resolve these issues, inverse uncertainty quantification (IUQ) is conducted to determine the parameter uncertainties. IUQ is a process of quantifying the uncertainties of physical model parameters given relevant experimental measurements. In this study, the OECD/NRC BWR Full-size Fine-Mesh Bundle Tests (BFBT) benchmark steady-state void fraction data is used as the prior knowledge and Bayesian inference is implemented to formulate the IUQ problem. Two Markov Chain Monte Carlo (MCMC) sampling techniques, the adaptive Metropolis-Hasting with global scaling algorithm and the adaptive Metropolis-Hasting with component-wise scaling algorithm, are implemented and investigated to explore the posterior probability density functions. Once the framework is established, approximate 100000 runs of CTF are needed which is a big computational burden. Gaussian process regression (GPR) is therefore introduced to construct the surrogate model for CTF, which can significantly reduce the simulation time and make the Bayesian inference possible. The outcomes contain the posterior distributions of 7 physical model parameters that are crucial to the BFBT experiment. Results of IUQ can be used for the forward uncertainty propagation and validation analysis. The results also show that two MCMC algorithms take the advantages of self-adaptive method, capable of inferring the posterior distribution of model parameters with arbitrary form of prior information. These two MCMC algorithms give consistent results and the posterior accepted samples of each model parameter are then fitted to the classic distribution.
What problem does this paper attempt to address?