王瑞利，梁 霄，林 忠

基于误差马尾图量化爆轰数值模拟结果的置信度＊

王瑞利1，梁 霄2，林 忠1

（1.北京应用物理与计算数学研究所，北京100094；2.山东科技大学数学学院，山东 青岛266590）

针对爆轰流体力学数值模拟过程中输入参数的不确定性，通过抽样技术，形成确定性爆轰流体力学程序的各种输入和数值求解，建立输入参数与输出响应量的样本，再通过概率框架下的误差累积分布函数与马尾图，给出了爆轰数值模拟过程中输入参数不确定度对模拟结果影响的置信度量化方法。通过一维黎曼问题、平面爆轰问题计算了误差马尾图，给出了二维爆轰拉氏自适应流体动力学LAD2D程序计算网格与模拟结果置信度的关系，对多物理爆轰过程发展高置信度数值模拟软件有很好的借鉴作用。

累积分布函数；误差马尾图；置信度量化；爆轰数值模拟

LAD2D程序［1］是一个非结构网格拉氏自适应的二维爆轰流体动力学软件，此程序已通过软件质量保证（SQA）和大量测试模型的考核，验证了程序的正确性［2－4］，但由于爆轰流体力学物理过程的复杂性和人们认识的缺陷，在建模过程中含有抽象、简化和近似，逼真建模很难，有时只能唯象建模和逐渐逼近真实情况，建模过程含有不确定性。加之描述其过程的数学物理模型是高度非线性的偏微分方程组，很难解析求解［5－6］。在数值求解过程中，由于连续到离散，存在计算模型误差、离散误差、计算机舍入误差和分析误差等，数值模拟过程始终是一种近似，含有不确定性。……