预处理技术稳定性分析与改进.pdfVIP

中国下程热物理学会 热机气动热力学 第卜一届年会 编}}：052001 预处理技术稳定性分析及改进 李雪松 徐建中 4‘ 中国科学院工程热物理研究所 北京：100080 E—mail：lixs78@sina．COlll TEL：010 摘要： 依据Roe型格式，探讨r产生预处理计算的不稳定性的原因。通过对预处理后的特征矩阵的 分析．认为其中台有不稳定的结构，而这是小稳定产生的原医【之一。依据这一思路，改进了预处理Roe 格式的实现形式，减少，特征矩阵中不稳定的结构。计算实姨表明，这一改进，相当显著的提商了预 处理计算的稳定性。这一结果也反过来证实，对计算不稳定原因的分析，从而为进一步的改进，提供 r罂路。 关键词：预处理特机矩阵稳定性Roe型格式 1前言 采用预处理技术，可以将高速可压缩流中r。泛应_I|=】的各种激波捕获格式，如TvD格 式等．戍用到低速不可压流动中。这是因为，一方面，预处理技术克服了低速流条件r 系统刚性人的弱点，使收敛加速，无论流速有多低，收敛速度都与计算高速流时相似， 从而人人减少了计算量，流速越低，这一优势越明显；另一方面，通过预处理改造，修 止了激波捕获格式在低速流下过火的耗散，使得格式耗散与实际流动的耗散同尺度”。”， 这保证了激波捕获格式在低速不可压缩流条件下，能够计算出正确的结果。因此，预处 理技术，使得应用一套程序，可以处理从不可压缩流到超音速流的所有马赫数的流动， 具有极人的好处。这即可以使我们将已有的可压缩流代码，延伸应用到不可压缩流中， 也使得计算可压不可压混和流动成为可能。 然而，预处理技术也存在一些缺陷，其中主要的一点，是预处理计算稳定性欠佳， 为了使得计算收敛，不得不对一个关键的预处理因子0做截断处理，从而影响了计算的 精度。为了改进这一点，也有很多工作将0局部化，以期获得更好的精度，但同时也导 致计算稳定性F降“。 本文将探讨引起预处理计算稳定性的原因，并采取改进措施。 2预处理技术 预处理技术就是在N—s方程组时间项上乘以预处理矩阵，从而改变系统特征值 进而政变系统刚性与格式的数值耗散。不妨采用欧拉方稃，三维任意曲线坐标系(‘ q．‘)F，守恒形式的预处理方程如F： 一旦+堡+塑+塑：0 (1) a 粥 卸 凿 Q=JQ F=J({、F+E。O+乞H) 。 其中：一 、， ’、 G=J【”。F+11yG+11：H) H=J(‘。F+‘，G+‘：H) PuPVPwP 其中Q=[P E]7为守恒变量，F、G、tt为对流通量，其具体形式可以 参见一般的计算流体力学讲义。r为守恒形式下的预处理矩阵。 uvws]T卜．，取预处理矩阵形式为对角阵”1： 在以压力、熵定义的基本变量w=[P 1 1 l F。=diag[01]。这与文献”1所提出的预处理阵本质上是一致的。r与r。，的关 系为： r：盟r．型 (2) 8W”．!}O 这里的关键是0的定义，0应与马赫数M同量级。根据Turke]的建议“、…，0可定 义为如r形式： 0。。)，1)0m。=KM。f2 (3) e=min(max(M2 M“为参考马赫数，对于外流可取远场来流马赫数，对于内流可取流场内平均或最大 马赫数。 理想中，0。。应为一个小