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 条
[31]   Inverse uncertainty quantification of input model parameters for thermal-hydraulics simulations using expectation-maximization under Bayesian framework [J].
Shrestha, Rijan ;
Kozlowski, Tomasz .
JOURNAL OF APPLIED STATISTICS, 2016, 43 (06) :1011-1026
[32]  
Smolyak S. A., 1963, Sov. Math. Dokl, V4, P240
[33]  
Sobol I.M., 1990, Matematicheskoe Modelirovanie, V2, P112
[34]   Global sensitivity indices for nonlinear mathematical models and their Monte Carlo estimates [J].
Sobol, IM .
MATHEMATICS AND COMPUTERS IN SIMULATION, 2001, 55 (1-3) :271-280
[35]   A dynamically adaptive sparse grids method for quasi-optimal interpolation of multidimensional functions [J].
Stoyanov, Miroslav K. ;
Webster, Clayton G. .
COMPUTERS & MATHEMATICS WITH APPLICATIONS, 2016, 71 (11) :2449-2465
[36]  
Stoyanov MK, 2013, TECHNICAL REPORT
[37]   Global sensitivity analysis using polynomial chaos expansions [J].
Sudret, Bruno .
RELIABILITY ENGINEERING & SYSTEM SAFETY, 2008, 93 (07) :964-979
[38]   Calibration, validation, and sensitivity analysis: What's what [J].
Trucano, T. G. ;
Swiler, L. P. ;
Igusa, T. ;
Oberkampf, W. L. ;
Pilch, M. .
RELIABILITY ENGINEERING & SYSTEM SAFETY, 2006, 91 (10-11) :1331-1357
[39]  
USNRC, 2014, TRAC RELAP ADV COMP, V1
[40]   Historical insights in the development of Best Estimate Plus Uncertainty safety analysis [J].
Wilson, Gary E. .
ANNALS OF NUCLEAR ENERGY, 2013, 52 :2-9