Estimation of probability Density Functions for model input parameters using inverse uncertainty quantification with bias terms

The documentation of most nuclear thermal-hydraulics codes does not provide sufficient information on uncertainty of physical models (e.g. interfacial heat transfer coefficients). These models were derived based on experimental data and implemented as empirical correlations in the computational code. The uncertainty quantification for the relevant output quantity (e.g. Peak Cladding Temperature) requires estimation of the Probability Density Functions (PDFs) of the code inputs, such as physical models. In this paper, we investigate the effect of boundary conditions (outlet pressure, inlet liquid temperature, and inlet flow rate) on the uncertainty of two physical models (the interfacial friction coefficient and the wall to liquid friction coefficient). The boundary conditions effect was accounted for by adding a bias term to the mathematical framework of two existing methods for Inverse Uncertainty Quantification (IUQ): the Maximum Likelihood Estimation (MLE) method and the Maximum A Posterior (MAP) method. The two methods were demonstrated using the BFBT benchmark, experimental data was compared to code predictions of the RSTART thermal-hydraulics code for two different cases: without and with bias term. The results show an evident improvement in code prediction when the bias term is used. Finally, a validation set of experimental data was used to investigate the possibility of data overfitting, and the proposed methodology showed absence of overfitting when bias terms are used. (C) 2019 Elsevier Ltd. All rights reserved.