基于UQLab的COBRA不确定性量化分析.docx

?

?

基于UQLab的COBRA不确定性量化分析

?

?

姜强刘天才杨宏伟

摘?要：本文以1/8堆芯子通道模型作为COBRA计算实例，以通道冷却剂入口温度、通道平均入口质量流速、系统压力等参数作为输入不确定性源，选择燃料最大温度作为模型响应值。使用UQLab，得到了不同抽样样本数下模型响应值的概率密度分布、可能的范围;同时应用不同的敏感性分析方法，定量描述了模型响应值受每个输入变量不确定性影响的大小。

关键词：UQLab?COBRA?不确定性量化

：TL333????????????????：A???????????：1674-098X（2020）02（b）-0076-04

Abstract：Thispaperusesa1/8coresubchannelmodelasaCOBRAcalculationexample.Channelinletcoolanttemperature，averageinletmassflowrate，systempressure，andotherparametersareusedasthesourcesofinputuncertainty.Themaximumfueltemperatureistakenasthemodelresponse.TheprobabilitydensitydistributionandthepossiblerangeofthemodelresponseunderdifferentsampleswereobtainedbasedonUQLab.Atthesametime，differentsensitivityanalysismethodswereappliedtoquantitativelydescribetheinfluenceofeachinputuncertaintyonmodelresponse.

KeyWords：UQLab;COBRA;UncertaintyQuantification

核能开发与应用领域拥有大量不同规模的数值模拟软件，并被广泛应用于模拟核电站、反应堆等安全攸关系统，具体包括反应堆设计、正常工况模拟、事故序列预测等方面。评估软件计算的可靠性、模型响应的稳健性、可能的结果范围，对于反应堆工程设计至关重要，需要有科学、完善的不确定性量化技术。本文以COBRA软件作为UQLab的分析算例，有效地量化了COBRA感兴趣输出参数的不确定性。

1?不确定性量化分析概述

主流不确定性量化分析的技术路线分为两种：输入不确定性传播方法和输出不确定性的推断方法。前者主要思路是使用少量的不确定输入参数，将主观概率范围和概率分布分配给这些参数，通过核心模型传播不确定性来确定输出参数的统计特性;后者的主要思路是使用一组相关实验数据为大量输入参数建立不确定性数据的数据库，再创建表征各种条件、瞬态等物理参数的超立方体，通过比较导出输出参数的误差带。

通常，输出不确定性的推断方法要求有一个大型而全面的实验数据库可用，轻水堆和沸水反应堆的不确定性研究可以使用这种方法。但是，在高温反应堆或者快堆领域，只有非常有限的实验和运行数据存在，输入不确定性传播方法往往是唯一可行的方法。

输入不确定性传播方法由于其较好的普适性，在国内外得到了较为广泛的应用：德国GRS（GesellschaftfürAnlagen-undReaktorsicherheit）基于输入不确定性传播方法，开发了SUSA（SoftwareforUncertaintyandSensitivityAnalyses）不確定性分析软件，广泛地应用于核反应堆不确定性分析计算中。美国SANDIA国家实验室开发的多级并行面向对象框架DAKOTA可用于设计优化、参数估计、不确定性量化和灵敏度分析等领域，其可嵌入的特性在SNAP、TRACE等软件中得到了广泛的应用。

2?UQLab简介

UQLab是苏黎世联邦理工学院（瑞士）开发的基于MATLAB的通用不确定性量化框架。由开源科学模块组成，基于输入不确定性传播方法，通过蒙特卡罗采样、敏感性分析、可靠性分析、构建替代模型等进行不确定性量化。

3?COBRA简介

COBRA是欧洲研究人员于开发的子通道分析程序，并成功应用到轻水堆计算程序NORMA、QUARK以及NORMA-FP的热工水力计算部分。程序能够进行堆芯分析和子通道分析，在反应堆工程热工水力的设计中，COBRA使用极为广泛。

COBRA程序采用固定格式的