Publication

Bayesian Uncertainty Quantification of Physical Models in Thermal-Hydraulics System Codes

Damar Canggih Wicaksono
2018
EPFL thesis
Abstract

Nuclear thermal-hydraulics (TH) system codes use several parametrized physical or empirical models to describe complex two-phase flow phenomena. The reliability of their predictions is as such primarily affected by the uncertainty associated with the parameters of the models. Because these model parameters often cannot be measured, nor have inherent physical meanings, their uncertainties are mostly based on expert judgment.

The present doctoral research aims to quantify the uncertainty of physical model parameters implemented in a TH system code based on experimental data. Specifically, this thesis develops a methodology to use experimental data to inform these uncertainties in a more objective manner. The methodology is based on a probabilistic framework and consists of three steps adapted from recent developments in applied statistics: global sensitivity analysis (GSA), metamodeling, and Bayesian calibration.

The methodology is applied to reflood experiments from the FEBA separate effect test facility (SETF), which are modeled with the TH system code TRACE. Reflood is chosen as a relevant phenomenon for the safety analysis of light water reactors (LWRs) and three typical time-dependent outputs are investigated: clad temperature, pressure drops and liquid carryover.

First, GSA allows screening out input parameters that have a low impact on the reflood transient. Functional data analysis (FDA) is then used to reduce the dimensionality of the time-dependent code outputs, while preserving their interpretability. The resulting quantities can be used once more with GSA to investigate, quantitatively, the effect of the input parameters on the overall time-dependent outputs.

Second, a Gaussian process (GP) metamodel is developed and validated as a surrogate for the TRACE model. The average prediction error of the metamodel is sufficiently low to predict all considered outputs, and its computational cost is less than 5 [s] as compared to 6−15 [min] per TRACE run.

Third and finally, the a posteriori model parameter uncertainties are quantified by calibration on a selected test from the FEBA experiments. Several posterior probability density functions (PDFs) corresponding to different calibration schemes – with and without model bias term and for different types of output – are formulated and directly sampled using a Markov Chain Monte Carlo (MCMC) ensemble sampler and the GP metamodel. The posterior samples are then propagated in a set of FEBA experiments to check the validity of the posterior model parameter values and uncertainties.

The calibration is performed on different types of output to inform model parameters that would have otherwise remained non-identifiable. The calibration scheme with model bias term is able to constrain the prior uncertainties of the model parameters while keeping the nominal TRACE parameters values within the posterior uncertainty interval. That is in contrast with the results of the calibration without model bias term, in which the posterior uncertainties are concentrated on either side of the prior range, and at times do not include the nominal TRACE parameters values.

The methodology was shown to successfully inform the uncertainty of the model parameters involved in a reflood transient. In the future, the methodology can be applied to model parameters involved in other TH phenomena using data from SETFs and contributes to achieve the goal of quantifying uncertainties in the safety assessment of LWRs.

About this result
This page is automatically generated and may contain information that is not correct, complete, up-to-date, or relevant to your search query. The same applies to every other page on this website. Please make sure to verify the information with EPFL's official sources.

Graph Chatbot

Chat with Graph Search

Ask any question about EPFL courses, lectures, exercises, research, news, etc. or try the example questions below.

DISCLAIMER: The Graph Chatbot is not programmed to provide explicit or categorical answers to your questions. Rather, it transforms your questions into API requests that are distributed across the various IT services officially administered by EPFL. Its purpose is solely to collect and recommend relevant references to content that you can explore to help you answer your questions.