Inverse uncertainty quantification of TRACE physical model parameters using sparse gird stochastic collocation surrogate model

被引:24
作者
Wu, Xu [1 ]
Mui, Travis [1 ]
Hu, Guojun [1 ]
Meidani, Hadi [2 ]
Kozlowski, Tomasz [1 ]
机构
[1] Univ Illinois, Dept Nucl Plasma & Radiol Engn, Urbana, IL 61801 USA
[2] Univ Illinois, Dept Civil & Environm Engn, Urbana, IL 61801 USA
关键词
POLYNOMIAL CHAOS EXPANSION; SENSITIVITY-ANALYSIS; BAYESIAN FRAMEWORK; SIMULATIONS; CALIBRATION; VALIDATION; INDEXES; GRIDS;
D O I
10.1016/j.nucengdes.2017.05.011
中图分类号
TL [原子能技术]; O571 [原子核物理学];
学科分类号
0827 ; 082701 ;
摘要
Within the BEPU (Best Estimate plus Uncertainty) methodology uncertainties must be quantified in order to prove that the investigated design remains within acceptance criteria. For best-estimate system thermal-hydraulics codes like TRACE and RELAP5, significant uncertainties come from the closure laws which are used to describe transfer terms in the balance equations. The accuracy and uncertainty information of these correlations are usually unknown to the code users, which results in the user simply ignoring or describing them using expert opinion or personal judgment during uncertainty and sensitivity analysis. The purpose of this paper is to replace such ad-hoc expert judgment of the uncertainty information of TRACE physical model parameters with inverse Uncertainty Quantification (UQ) based on OECD/NRC BWR Full-size Fine-Mesh Bundle Tests (BFBT) benchmark steady-state void fraction data. Inverse UQ seeks statistical descriptions of the physical model random input parameters that are consistent with the experimental data. Inverse UQ always captures the uncertainty of its estimates rather than merely determining point estimates of the best-fit input parameters. Bayesian analysis is used to establish the inverse UQ problems based on experimental data, with systematic and rigorously derived surrogate models based on Sparse Gird Stochastic Collocation (SGSC). Global sensitivity analysis including Sobol' indices and correlation coefficients are used to identify the important TRACE input parameters. Several adaptive Markov Chain Monte Carlo (MCMC) sampling techniques are investigated and implemented to explore the posterior probability density functions. This research solves the problem of lack of uncertainty information for TRACE physical model parameters for the closure relations. The quantified uncertainties are necessary for future uncertainty and sensitivity study of TRACE code in nuclear reactor system design and safety analysis. (C) 2017 Elsevier B.V. All rights reserved.
引用
收藏
页码:185 / 200
页数:16
相关论文
共 43 条
[1]   A tutorial on adaptive MCMC [J].
Andrieu, Christophe ;
Thoms, Johannes .
STATISTICS AND COMPUTING, 2008, 18 (04) :343-373
[2]  
[Anonymous], 2008, GLOBAL SENSITIVITY A
[3]  
[Anonymous], 2274 AIAA
[4]  
[Anonymous], 2010, Verification and Validation in Scientific Computing, DOI DOI 10.1017/CBO9780511760396
[5]  
[Anonymous], SAND20144633
[6]   High dimensional polynomial interpolation on sparse grids [J].
Barthelmann, V ;
Novak, E ;
Ritter, K .
ADVANCES IN COMPUTATIONAL MATHEMATICS, 2000, 12 (04) :273-288
[7]   Statistical calibration of computer simulations [J].
Campbell, Katherine .
RELIABILITY ENGINEERING & SYSTEM SAFETY, 2006, 91 (10-11) :1358-1363
[8]   On the Lebesgue constant of Leja sequences for the complex unit disk and of their real projection [J].
Chkifa, Moulay Abdellah .
JOURNAL OF APPROXIMATION THEORY, 2013, 166 :176-200
[9]   Polynomial chaos expansion for sensitivity analysis [J].
Crestaux, Thierry ;
Le Maitre, Olivier ;
Martinez, Jean-Marc .
RELIABILITY ENGINEERING & SYSTEM SAFETY, 2009, 94 (07) :1161-1172
[10]  
De Crecy A., 2001, M C C SALT LAK CIT U