流体力学拉氏程序收敛性及数值计算不确定度初探.pdfVIP

第30卷第3期 计 算 物 理 V01．30，No．3 CHINESE OFCOMPUTATIONALPHYSICS 2013年5月 JOURNAL May，2013 文章编号：1001-246X(2013)03-0346-07 流体力学拉氏程序收敛性及数值 计算不确定度初探 刘 全1， 王瑞利1， 林 忠1， 刘希强2 (1．北京应用物理与计算数学研究所，北京100094；2．山东省聊城大学数学科学学院，山东252059) 阶、渐近收敛域、网格收敛指标(Gridconvergenceindex，GCI)等．给出GCI方法刻画数值计算不确定度的具体步骤． 并应用于N—R格式，验证了方法的有效性． 关键词：流体力学拉氏计算；渐近收敛性分析；网格收敛指标 中图分类号：0241 文献标识码：A 0 引言 流体力学欧氏方法(EulerianMethod)数值模拟中，经典收敛性分析描述在网格尺度趋于0的过程中，数 值解趋于精确解的程度和快慢．当得不到精确解时，需扩充收敛性的概念和方法．渐近收敛阶、渐近收敛域和 GCI等概念和方法就是在此基础上发展和提出的，已应用到数值模拟的验证和确认中¨11． 流体力学数值方法的研究一方面需要大力发展计算格式；另一方面需要评估实际应用格式的数值行为 (稳定性、收敛性、相容性和不确定度等)．收敛性分析是数值模拟不确定度研究的主要手段．这里的不确定 度是指当前计算值与计算渐近值之间的度量，不是数值解与精确解之间的度量．因为很多问题，不能确定其 精确解．对于实际问题，评估计算解的不确定度是复杂工程设计中很重要的一个环节．本文基于欧氏 (Eulerian)计算验证分析中的有关概念和方法并结合拉氏计算的特点，给出拉氏计算有关收敛性的概念和术 语，同时引进网格收敛指标、相对误差等表达式，用于分析数值计算本身的不确定度． 1 拉氏计算不确定度分析一些基本概念和方法 流体力学的数值方法按网格是否随流体运动，总体上可分为拉氏方法(Lagrangianmethod)、欧氏方法 (Eulerianmethod)和任意拉氏欧氏方法(Arbitrary Lagrangian 到高维和方程组．由于欧拉方法网格固定，微分方程的解析理论容易推广到离散情形，因此格式的稳定性、收 敛性分析相对比较容易．但在拉氏计算中，由于网格随流体运动，且拉氏格式的构造是直接从方程组出发，同 时流体力学方程组的解析理论远不成熟，因此离散格式的数值分析十分困难．现阶段人们更多用数值模拟的 方式来分析收敛性问题．首先须明确CFD计算结果的误差来源，一般可分为四类”1：①物理模型误差，②舍 入误差，③离散误差，④程序设计误差．此外对隐式求解问题，还有迭代误差．这里假定物理模型误差可以忽 略不计，只讨论数值离散和程序设计本身引起的误差．在交错型网格拉氏计算中，由于激波和网格非物理畸 变而引进的冲击波人工粘性和沙漏粘性(统称人为粘性)，对收敛性分析亦引进了相当的难度，计算过程中 网格的长宽比和均匀性对计算结果误差也会产生影响． 文献[1—2，5]等针对欧氏计算的特点，引进渐近收敛阶、渐近收敛域、网格收敛指标等基本概念和方法 收稿日期：2012—07—16；修回日期：2012—11—02 作者简介：刘全(1972一)，男，博士，副研究员，从事流体力学计算和程序设计，E-mail：liuquan@iapcm．ac．cn 万方数据 第3期 刘 全等：流体力学拉氏程序收敛性及数值计算不确定度初探 347 用于不确定度的分析．下面将这些概念和方法引进到拉氏计算中． 拉氏流体力学数值模拟，尽管网格大小随时间非均匀变化，但每个网格的质量保持不变(等于初始时刻 网格划分时的质量)．首先引进误差的定义，给定网格质量Ami和时间步长△￡，计算得到的任一物理量手，点 形式的误差估计是精确解和计算解之间的质量和时间步长的函数 亭：一f；…‘=占。+A(Ami