非定常大步长隐式TVD格式途径.pdfVIP

中国工程热物理学会 热机气动热力学 第十届年会 编号:012002 非定常大步长隐式TO格式途径 陈宏冀 中国科学院工程热物理研究所，北京1000080 电话:(010E-mail:chj@mail.etp.ac.en 摘要: 模拟非定常问题的时间精确计算需要高精度大步长TVD格式，然而用最常见的 二阶时间精度隐式差分方法来构造大步长TVD格式目前仍缺乏理论依据，本文找到了一 个新的途径，并证明了得出的格式为无条件TVD格式，有望满足对非定常计算的迫切需 求。 一、引 言 显式格式时间推进计算的步长由于稳定性的要求而受CFL条件限制;若是TO格 式，也要在类似的条件下是TO的川。采用隐式格式就是为了克服这样的限制。时间后 差Euler全隐式格式对任意步长稳定，用其构造的TV。格式也是无条件TV。的[z.”。尽 管它只有一阶时间精度。但很适用于定常问题的渐近解。对于非定常问题的时间精确 解，则需要高一些的时间积分精度。二阶精度的Crank-Yicolson格式是无条件稳定的， 但构造成TO格式要在小于2的Courant数下是TV。的“〔，，，在步长上得到的好处太小。 另一种达到二阶精度的办法是采用三层时间后差，但未见有人给出这样的TVD格式条 件。可见适用于时间精确解的隐式TVD格式的大步长推进还缺少理论支持。由于这个原 因，在国际文献中，与定常计算大量采用了隐式TVD格式的情况相比，非定常计算的情 况很不一样。Yee在文献仁31中对定常解采用隐式TVD格式，而对时间精确解采用了显 式的。Rai和Chakravarthy在文献[4]中用全隐式TO格式进行定常计算，其后Rai在 文献[5,6〕中计算透平中的非定常流动，提到当需要二阶时间精度时采用三层时间后 差，但没有用TVD格式。Yoon和Rwak在文献[7〕中将Euler全隐式LU-SGS算法与TVD 格式结合起来，在文献[8〕中用的也是这样的格式，然而他们虽然在不少论文中提到二 阶时间精度Crank-Nicolson差分，却没与TO格式结合起来，例如文献幻〔。在众多文 献中非定常流计算用TO 格式的多用显式，例如文献[10,11]:或用隐式方法而不用 TO格式，例如文献 1〔2];很少有用隐式TO格式的，尤其是在比较权威的和在理论上 比较慎重的工作中。另一方面，也有少数作者从实践的角度试用了隐式TVD格式，在较 大步长下对具体问题得到了可以接受的计算结果，其中有的还观察到了收敛上的异常现 象洲。 随着对非定常流动现象的观注程度的增大，大步长时间精确解隐式TVD格式的理论 基础迫切需要解决。本文就是作者致力于此的一个开始。 二、单参数隐式格式族 对于标量双曲型守恒律 231 au af(u) - 十 =0 (1) at ax 考虑它的带参数0的差分格式 u1.+mh+vz一h;-+i/z=ui一“(1一0)lhi+vz一h,-，二) (2) 其中A=At/Ax,h为数值通量。当0=0时得显式格式，0=1时得Euler全隐式格 式，都只有一阶时间差分精度。当0=1/2时则得到二阶时间精度的Crank-Nicolson 格式。将格式改写成 ul+1一A叫C二J+v/2zA△,”+v./2zu-卜-Cctj-,1//2zA△jj-一,，//2zuIn+i 一u,+A(1一0)(C;+i/2A;+vzu-Cj-vzAt-1/2u). (3) 的形式。已经证明显式格式为TVD格式的一个充分条件是G[1 兄件1/2-0