网格之间流场信息传递的高精度方法和验证-空气动力学学报.PDF

第28卷第1期 空气动力学学报 V01．28，No．1 2010年2月 ACTAAERoDYNAMICASINICA Feb．．2010 文章编号：0258—1825(2010)01—0001—06 网格之间流场信息传递的高精度方法和验证 刘 君1，白晓征1，王 巍2 摘要：理论分析了插值方法引入误差的原因．提出一种利用时间推进方法和移动网格技术实现网格之间流场信 息传递的思想。并构造计算方法。在介绍新的信息传递算法原理的基础上，根据三维情况下传递过程需要满足的 时间推进步长稳定性条件，推导出最小推进步数的要求。最后，通过对激波诱导物体运动的非定常流动模拟。比较 了信息传递算法和线性插值、加权最小二乘插值方法．计算结果表明，与插值方法相比．信息传递过程没有引入误 差。能有效抑制非物理振荡，保持流场计算精度。 关键词：多块网格；插值方法，计算流体力学；非定常流场模拟 中图分类号：V211．3 文献标识码：A 方法，给出了传递过程中时间推进步长的稳定性条 0 引 言 件。 本文进一步分析了插值方法引入误差的原因，然 网格分区技术是结构网格求解复杂外形的基础， 后在介绍新的信息传递算法原理基础上，推导了三维 广泛应用的主要有重叠(overlapping)网格技术和对 情况下传递过程时间推进步长的稳定性条件，最后给 接(patched)网格技术两类，计算过程中，首先采用流 出了三维算例。 体力学解算器分别求解各子区域网格上物理量，然后 通过边界上信息传递，从而完成整个流场的计算。实 1插值方法引入误差的原因分析 际问题中相邻子区域的边界网格或重叠部分网格很 少在位置上正好一一对应，因此需要采用一定的算法 采用有限差分法求解波动方程L(乱)一甍+口笔 传递相邻子区域的信息。对于包含有激波等间断的 一0，等价于求解其修正方程： 流场模拟，目前常用线性插值或其他高精度插值方 法，插值表达式通常为多项式，用网格自身及其周围 L(“)=∑竹努 (1) ”山 Z一2 单元的信息来确定多项式的系数；线性插值优点是满 由于系数扎与网格空间分布间距△z相关，计算得到 足单调性、守恒性、形式简单、效率高，但是只有一阶 的流场参数与网格分布相关，在其他条件相同的情况 精度，数值耗散大；高精度插值提高精度，但是涉及到 下，不同网格分布在同一物理位置的数值解也不相 网格内物理量导数计算，而且在应用非结构动网格技 同。 术和插值方法求解分离问题时发现，网格重构以后流 设有两套网格，在同一物理位置z，处，根据差 场出现的非物理振荡会导致气动力明显的波动[1，。]。 分方程得到第二套网格数值解为“P，根据第一套网 在文献[1]中给出一种新的信息传递算法，充分 格结点变量插值得到的该点值为U丁。插值方法是相 利用时间推进的特点，采用移动网格的方法，随着时 同时刻不同空间位置物理量之间的关联，可以写为二 间推进完成旧网格信息到新网格的传递。这种新的 者之间距离的函数：