NS方程激波计算的摄动有限差分方法研究.pdfVIP

NS方程激波计算的摄动有限差分方法 NS方程激波计算的摄动有限差分方法 申义庆高智 杨国伟 中国科学院力学研究高温气体动力学实验室，北京，100080 摘要：摄动有限差分(PFD)方法从～阶迎风差分格式出发，将差分系数展开为网格步长的幂级数，通过提 高修正微分方程的逼近精度来获得更高精度的差分格式。由于格式基于一阶迎风格式，因此具有迎风效应、 网格节点少等特点。本文首先通过对Burgers方程的摄动差分格式的推导，将摄动有限差分格式引入时间相 关法的计算，并构造了守恒形式的摄动有限差分格式，然后推广到一维Navier-Stokes方程组的计算。数值 比较研究表明：本文构造的Ns方程摄动有限差分格式具有比～阶迎风较高的精度和分辨率，而且保持了一阶 迎风格式的无振荡性质。 O引言 高精度、高分辨率差分格式是计算流体力学数值方法研究的目的之一。一般来说，差分 格式的精度越高，使用的网格点越多，而对间断解的高分辨率性质则通过采用附加条件的限 制来实现，如二十年来提出和发展的TVD、ENO、NND等格式u。30，但由于所涉及的网格点较多， 给边界的处理带来一定的困难。紧致格式H1利用较少的节点能获得较高的精度，但在计算激波 时会带来非物理振荡，而且由于格式中利用了节点上的导数值，对于复杂边界的情况也难以 给定。数值摄动方法基于一阶迎风差分或二阶中心差分格式，对差分方程的非微商项进行数 值摄动展开，通过提高修正微分方程的逼近精度来获得高精度的格式，由于摄动格式只利用 了三个节点，具有迎风性等特点，因此在计算流体力学中将具有广阔的应用前景。 数值摄动有限差分格式的高精度、高分辨能力已在众多的模型方程和不可压缩流动的计 算中得到了应用和验证㈣。，但对可压缩流动的计算尚需进一步的发展。本文首先通过对 Burgers方程的摄动差分格式的推导，将摄动有限差分格式引入时间相关法的计算，并构造了 守恒形式的摄动有限差分格式，然后推广到一维Navier—Stokes方程组的计算。数值比较研 究表明：本文构造的NS方程摄动有限差分格式比一阶迎风具有较高的精度和分辨率，而且保 持了一阶迎风格式的无振荡性质。 1模型万程搬动有限爱分万法 对Burgers方程 一十一=———_ ㈤1L， 要+u罢：去窘 m 孤Re舐2 其在定常状态时为 “耋：土垒 蟛’(2) “面2面萨 对方程(2)中一阶导数采用一阶迎风格式，二阶导数采用二阶中心格式，有如下差分方程 一l+ot“，竺[堕+坐z，，燮：土坚二靼(3)2 o 缸 2 。 缸 Re 缸2 其中，Q=sign(u，)。设对差商系数U，进行如下摄动展开 (4) 彦-1／J(1+么l△x+么2Ax2+彳3△x3+么4△x4+…) 51 第十一届全国激波与激波管学术会议 将(4)代入(3)后可得修正微分方程为 嘣詈卜善燕碧恤卜去r窘”善赢器归2Ⅲ5， 在J点对方程(2)进行局部点凝固系数，可得关系式 ≥沪(眦广1罢I， (6) 在(5)中对缸”利用待定系数法，则可求得系数4。 么。=丽1(口Reuj)”，刀=o，1，2，… 由此，将“归代替(3)中的差商系数”／，则可得n阶得摄动差分格式： 半“归警+半“印警=丙1笪等塑 如对方程(2)利用时间相关法求解，即求方程(1)的定常解，则n阶格式为 “』+l一2u』+材／一 詈”半“，訾+半“，警=南 Il 其中， 4 4缸 对时间项詈I／的求解可采用Runge—Kutta方