差分格式的时间耗散效应对气动方程定常解精度的影响.pdfVIP

差分格式的时间耗散效应对气动方程定常解精度的影响 刘 君 (冒尊特技大孝t天擅求未． 鼍妒市410073J 囊要，本文疆过囊擅实验和理论分析，研究了应用时问相关法求■定常舟■过卷中。差分格式 的时阃麓度、时阔分曩法和边界处理对定，嘲的嚣响，认为在差分格式的构造和应用中，不但要注意 空问耗散，还应谖注童时间托散等因素． 关t词：●t分方程敦值■理论，差分括式．时润耗散 1 弓l言 对这一目题的兴趣是数值模拟壁面热流时产生的“’，当时莰们除了发现并从理论上证 明不同的差分格式对网捂雷祷鼓的要求不同外“；还发现一些差分格式得到澈波分辨率和 物面压力分布十分接近，但是热藏相差悬殊这样的现象．在调研后觉得现有的值微分方程 数值解法理论对这些问题阐述不够清楚，因此进行了较为深入的研究． 2 时间高精度格式的时问耗散 从复杂的非线性流体力学方程直接进行理论研究几乎是不可能的，所以只能采用非常 替Eulet方j莹，来进行为数值分析和差分格式构造与验证．可以看出这些方程是非定常 的，但从目前计算空气动力学发晨状掘来看，应用中求解的大部分问题是定常的，因此为 了把非定常差分格式用于定常问题，一般采用时间相关法。即将定常边再条件卞非定常解 长时间的渐近极限作为问题的定常解．不过从下面的简单例证可发现。非定常差分格式在 求解定常解的过程中其特征，如耗鼓、弥数和糟度等会发生变化，从非定常模型方}到寻出 的结论可能不再适用． 琏+4‘如=r·k (2．1) 咋+a·咚=0 (2．2) 5R 1] 上取△茸=o．05均匀只格，采用如下三种差分格式的时问相关法求其定常解．收敛标准 为击秘s+l_l／4I=o，这里。是指蝴计算机的机器零· 格式一： 掣+口警=，，华 械二：譬+口警一譬．学=，，学 格式兰：ltscCurmiek格式． 取不同的时问推避步长计算后，数值解和糟确解之爰△冀=矿一％的分布如图1．可 以发现：格式一的定常解不依棘△f，但是At对格式二和格式三的定常教值解有明显影 响；而且△f越小。△“的分布遗趋向于格式一的解． 这种差分格式的定常解依赣于时间步长的现象很早兢被发现鼬，但没有得到重视，直 翼绕藏时发现。定解与时阊推进步长有关…”；文献[s]研究1vD格式的人T耗散时发现； 结论与本文一致．。小的时闻步长得到的解更接近解析解”． 本文把上述定常解依族子时间步长的现象差分格式的时间耗敬效应：通过下面的理论 分析可以证明是构造格式过程中时间方向采用高精度离鼓是原因之一． 构造显式的时问二阶差分格式基础是如下Taylor展开式： ∥列一《尝y+笙2(窘丁 cz．。，…” L街J I掰) 根据(2．1)或(2．2)式将时问导数酱换为空向导致，并进行离散，得封差分格式 ∥=茸；一At‘％+o．钮2·％ Q．4) 式中‘为(2．3)式中时间一阶导数替换为空间导数后的离散算于，厶为(2．3)式中时闻二 ● 阶导数替换为空问导数后的离散算予．用于定常问题计算收敛时有：≈，1=≈n．，即 Ll，=O．5△l·L2』