混合不确定度量化方法及其在计算流体动力学迎风格式中的应用.pdfVIP

第36卷第4期 爆炸与冲击 V01．36，No．4 2016年7月 EXPLOSl0NANDSHOCKWAVES Jul．。2016 混合不确定度量化方法及其在计算流体动力学 迎风格式中的应用+ 梁 霄1’2，王瑞利1 (1．北京应用物理与计算数学研究所，北京100094； 2．山东科技大学数学学院，山东青岛266590) 摘要：针对流体力学数值求解间断问题时，初始状态含有偶然和认知混合型的不确定性，将认知不确定 chaos， 度作为外层，偶然不确定度作为内层，分别使用非嵌人多项式混沌方法(non—intrusivepolynomial NIPC)和概率盒(P-box)理论处理偶然不确定度和认知不确定度，发展了流体力学数值求解过程中，初始状态 含有混合不确定度传播量化的一种方法。以迎风格式和黎曼解法器求解Sod问题为例，评估了左状态密度 (偶然不确定度)和理想气体多方指数(认知不确定度)对模型输出结果的影响，有效验证了该方法的可行性。 关键词：爆炸力学；建模与模拟；非嵌入多项式混沌法；混合不确定度量化；概率盒 中图分类号：0385 国标学科代码：13035 文献标志码：A 在重大复杂工程领域，众多工程事故已经引发社会各界的关注。人们总是抱着一个良好的愿望，希 望能逼真建模，研制高可信度的数值模拟软件，精确模拟各类事故的发生过程，科学分析其产生的机理。 因此，数值模拟软件必须能真实反映实际系统的外部表征、内在特性及其动作行为等现象，其数值模拟 结果才能使人置信。但由于物理过程的复杂性及人们的认知缺陷，物理建模与数值模拟始终存在误差 和不确定性。为此，要利用数值模拟结果分析机理和模拟过程，必须量化数值模拟结果的不确定性，评 估建模与模拟的可信度。但长期以来，物理建模与数值模拟软件的研制者重点关注的是软件自身的研 制和对实际问题的顺利模拟，即算下去的问题，对软件自身的正确性、模拟结果的不确定性和建模与模 拟的可信度，即算多准的问题，一直没有给予足够的重视。尤其是对于理论、实验无法解释或难于解决 的科学问题，或做实验耗资巨大的工程设计，很想利用数值模拟真实、全过程，全系统、反复精密地进行 模拟的研究者持一种矛盾的心态：既想利用数值模拟这条快捷经济的重要研究途径，又对程序的正确性 和数值模拟提供的结果心存疑虑。为了给复杂装置设计提供高效可靠的数值模拟工具，必须展开计算 模型正确性和可信度的评估研究。 andvalidation，V 中的最基本模型，此计算模型的验证与确认(verification 8LV)受到高度重视[11]。W． 本概念、原理及步骤和发展过程。王瑞利等∞o综述了VV的相关概念、术语、规范、可信度评估方法 和应用等的研究进展，就MS的VV的几个关键问题进行了概括，构建了复杂工程M＆S的VV 的知识指南。确认是建立在定量准确评价的概念上，是评价模型描述实际物理过程／实验(试验)的准确 程度、模型形式与参数的适应范围，确认刻画实际问题／实验(试验)区域的形状，预测可信度评估等重要 活动。对于流体力学确认的核心是随机输入参数对模拟输出响应量不确定传播进行量化的方法。 不确定度量化(uncertainty *收稿151期：2014-12—24；修回日期：2015-03—27 山东科技大学人才引进项目(2013RCJJ027) 第一作者：梁霄(1984一)，男，博士，讲师，mathlx@163．corn。 5lO 爆 炸 与 冲 击 第36卷 认知缺陷尽可能给出分布或区间描述，建立建模和模拟中随机或认知不确定性因素和建模与模拟描述 系统响应量之间的关系，进而量化影响系统性能的不确定性，给出系统响应量性能分析的不确定度度量 的方法。对于不确定度量化有多种方法，其中敏感度分析方法、抽样方法、代理模型方法、随机微分方程