Investigation of Identification and Uncertainty Quantification Methods for Constitutive Model
-
摘要: 典型基于输入的最佳估算加不确定性(BEPU)分析方法需要全面考虑各种来源的不确定性,但作为重要不确定性来源之一的最佳估算程序内部本构模型不确定性却很难得到准确量化。针对传统BEPU方法的不足,采用函数型数据分析加次序敏感性分析方法识别重要模型,贝叶斯校准方法加高斯过程代理模型量化模型参数不确定性,得到了一种准确的关键本构模型识别与不确定性量化方法。使用该方法量化了FLECHT-SEASET实验中重要模型的不确定性,并将量化的不确定性抽样传播至包壳温度。另外对该方法与传统的关键本构模型识别方法进行对比。结果表明,该方法可以准确识别瞬态过程中的关键本构模型,传播计算结果能够很好地包络实验值。Abstract: Typical input-based best estimate plus uncertainty (BEPU) analysis methods require comprehensive consideration of uncertainty from all sources. Uncertainty of the input parameters is then propagated to the output by many sampling calculations through the best estimate (BE) code. In the past, the uncertainty of the code constitutive model, which was one of the most important sources of uncertainty, was often not properly analyzed. For the purpose of quantifying the uncertainty of code constitutive models and including it into the BEPU analysis, a methodology for the identification and uncertainty quantification of key constitutive models was developed in this paper. First of all, the figure of merit (FOM) of the sensitivity analysis needed to be defined. Time series data were chosen as the target in order to better characterize the effect of the parameters on the whole transient process, and a dimension reduction was performed for the time series data using functional data analysis (FDA) method to obtain the FOMs. Then, sensitivity analysis was performed on several FOMs to get the importance ranking of the constitutive models. The sensitivity analysis approach was a two-step process, where the input model parameters were firstly screened by the Morris method, and then the screened model parameters were subjected to a quantitative sensitivity analysis by the Sobol method to obtain a parameter importance ranking. Finally, uncertainty quantification was performed for the high ranking models selected based on the sensitivity analyses. The posterior probability distributions of the parameters were calculated using the Markov chain Monte Carlo (MCMC) method based on Bayesian calibration approach. In order to validate the applicability of the above method, it was used to quantify the uncertainties of the important models based on the FLECHT-SEASET experiment, and to propagate the quantified uncertainties to the cladding temperature. In addition, the result of the sensitivity analysis of the FOM obtained by FDA analysis was compared with two traditional FOMs, which were the integral of the difference between the sampling calculation results and the experimental values, and the peak cladding temperature (PCT), respectively. The results show that it is difficult to obtain an accurate ranking of model importance when using PCT as FOM. This approach can accurately identify the key constitutive models in transient processes, and the propagation calculation results can well envelope the experimental values.
-
-
[1] NRC. Emergency core cooling systems, revisions to acceptance criteria[R]. US:NRC, 1988.
[2] NRC. Acceptance criteria for emergency core cooling systems for light-water nuclear power reactors[R]. US:NRC, 1988.
[3] IAEA. Best estimate safety analysis for nuclear power plants:Uncertainty evaluation[R]. Vienna:IAEA, 2008.
[4] 熊青文,苟军利,杜鹏,等. 本构模型不确定性评价结构化方法开发及应用[J]. 核动力工程,2022,43(4):147-153. XIONG Qingwen, GOU Junli, DU Peng, et al. Development and application of structural method for uncertainty evaluation of constitutive models[J]. Nuclear Power Engineering, 2022, 43(4):147-153(in Chinese). [5] FREIXA J, de ALFONSO E, REVENTóS F. Testing methodologies for quantifying physical models uncertainties:A comparative exercise using CIRCE and IPREM (FFTBM)[J]. Nuclear Engineering and Design, 2016, 305:653-665.
[6] KOVTONYUK A, LUTSANYCH S, MORETTI F, et al. Development and assessment of a method for evaluating uncertainty of input parameters[J]. Nuclear Engineering and Design, 2017, 321:219-229.
[7] REVENTÓS F, DE ALFONSO A, ZHANG J, et al. PREMIUM, a benchmark on the quantification of the uncertainty of the physical models in the system thermal-hydraulic codes:Methodologies and data review[R]. Brussel:OECD, 2016.
[8] WU X, XIE Z, ALSAFADI F, et al. A comprehensive survey of inverse uncertainty quantification of physical model parameters in nuclear system thermal-hydraulics codes[J]. Nuclear Engineering and Design, 2021, 384:111460.
[9] WICAKSONO D C, ZERKAK O, PAUTZ A. Bayesian calibration of thermal-hydraulics model with time-dependent output[C]//Presented at the Proceedings of the 11th International Topical Meeting on Nuclear Reactor Thermal Hydraulics, Operation and Safety (NUTHOS-11). Korea:[s. n.], 2016.
[10] WU X, KOZLOWSKI T, MEIDANI H, et al. Inverse uncertainty quantification using the modular Bayesian approach based on Gaussian process, Part 1:Theory[J]. Nuclear Engineering and Design, 2018, 335:339-355.
[11] WU X, KOZLOWSKI T, MEIDANI H, et al. Inverse uncertainty quantification using the modular Bayesian approach based on Gaussian process, Part 2:Application to TRACE[J]. Nuclear Engineering and Design, 2018, 335:417-431.
[12] LI D, LIU X, YANG Y. Investigation of uncertainty quantification method for BE models using MCMC approach and application to assessment with FEBA data[J]. Annals of Nuclear Energy, 2017, 107:62-70.
[13] XIONG Q, GOU J, CHEN W, et al. Investigation of uncertainty quantification methods for constitutive models and the application to LOFT LBLOCA[J]. Annals of Nuclear Energy, 2019, 132:119-133.
[14] MCRAE G J, TILDEN J W, SEINFELD J H. Global sensitivity analysis:A computational implementation of the Fourier amplitude sensitivity test (FAST)[J]. Computers & Chemical Engineering, 1982, 6(1):15-25.
[15] 李睿. Sobol灵敏度分析方法在结构动态特性分析中的应用研究[D]. 长沙:湖南大学,2003. [16] MORRIS M D. Factorial sampling plans for preliminary computational experiments[J]. Technometrics, 1991, 33(2):161-174.
[17] 熊青文,苟军利,杜鹏,等. 适用于核电厂BEPU的高效全局敏感性分析方法开发及应用[J]. 原子能科学技术,2022,56(7):1321-1327. XIONG Qingwen, GOU Junli, DU Peng, et al. Development of efficient global sensitivity analysis method for BEPU in nuclear power plant[J]. Atomic Energy Science and Technology, 2022, 56(7):1321-1327(in Chinese). [18] XIONG Q, GOU J, WEN Y, et al. Global sensitivity analysis of LOFT large break loss of coolant accident with optimized moment-independent method[J]. Annals of Nuclear Energy, 2020, 139:107289.
[19] SALTELLI A, RATTO M, TARANTOLA S, et al. Sensitivity analysis practices:Strategies for model-based inference[J]. Reliability Engineering & System Safety, 2006, 91(10-11):1109-1125.
[20] RAMSAY J O, SILVERMAN B W. Fitting differential equations to functional data:Principal differential analysis[M]. Berlin:Springer, 2005.
[21] CAMPBELL K, MCKAY M D, WILLIAMS B J. Sensitivity analysis when model outputs are functions[J]. Reliability Engineering & System Safety, 2006, 91(10-11):1468-1472.
[22] LAMBONI M, MONOD H, MAKOWSKI D. Multivariate sensitivity analysis to measure global contribution of input factors in dynamic models[J]. Reliability Engineering & System Safety, 2011, 96(4):450-459.
[23] CHAO F, SHAN J, GOU J, et al. Study on the algorithm for solving two-fluid seven-equation two-pressure model[J]. Annals of Nuclear Energy, 2018, 111:379-390.
[24] SHAN J. Development of multi-scale high-fidelity system analysis code NUSOL-SYS[C]//Proceedings of the CDNET-2018. Chengdu, China:[s. n.], 2018.
[25] HOCHREITER L. FLECHT SEASET program:Final report[R]. USA:Westinghouse Electric Corp., 1985.
[26] VRUGT J A, TER BRAAK C, DIKS C, et al. Accelerating Markov chain Monte Carlo simulation by differential evolution with self-adaptive randomized subspace sampling[J]. International Journal of Nonlinear Sciences and Numerical Simulation, 2009, 10(3):273-290.
[27] GELMAN A, RUBIN D B. Inference from iterative simulation using multiple sequences[J]. Statistical Science, 1992, 7(4):457-472.
计量
- 文章访问数: 58
- HTML全文浏览量: 6
- PDF下载量: 21