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

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