三维可压缩流欧拉方程组隐式时间推进法.pdfVIP

维普资讯 第 8卷 第 1期 空 气 动 力 学 学 报 Vo1．8．No． ]990年 3月 ACrA AERODYNAMICA SINI(A M ． 1990 三维可压缩流欧拉方程组的隐式时间推进法 权正五 (吉林工业大学) 摘要 本文提出了求解三维定常无粘可压缩泷场问题的一种隐武时间推进浩 在计 掉空间以造变速度分量为束知盎量的欧拉方程组为控制方程，从而简化壁面也界的处理。在 Beam—W arming 的 AF一方法基础上，用对角化厦矢通量分裂等方法提高 了计 算速 度 殛解的精度和稳定性。作为算侧，对伴有擞渡的方形收放喷管踌声速流动进行了计算，以 显示方法的有效性。 关键词 可压缩流，时闻推进法，逆变速度分量。 一 、 引 言 定常无粘跨声速流动的控制方程是混合型的，这给数值求解带来很大困难。针对这 一 困难已经提出了多种隐式时间推进法。但这些方法在处理无粘流的壁面条件时一般不 能把边界条件直接嵌入基本方程 中，而通常在边界上建立新的关系式，和内点交互进行 求解。因而往往 出现由于边界而影响解的收敛速度，当解收敛之后在边界上仍留有相当 量的剩余等 问题。 1986年大官司一山本提出了一种新的时间推进法…。它是在计算空闻以逆变速度分 量为未知变量进行求解。其优点是，可以把壁面条件简捷地直接嵌人基本方程 中，边界 点和内点可同时进行计算，从而克服了上述在边界上的困难。 本文把选一方法推广应用于三维可压缩跨声速流动 问题中，并改进其边界处理的具 体作法，使之更为合理。因为以逆变速度分量为未知变量的欧拉方程组为控制方程，实 际用来求解的控制方程和原来的方程相比形式上将大不相厢，把它写成完垒守恒型是困 难的。但从控制方程 中适 当地分离出附加项，那么那些为提高计算速度和解的精度及稳 定性而行之有效的，例如线性化方法、AF一方法 、对角化 方 法 以及 矢 通 量 分裂 法…等照样能应用，而且最终要求解的式子也不比以往的方法更复杂。本文把所 有这些 方法融合于本解法中，使解法更为有效适用。最后作为算例，用本方法计算了三维方形 喷管的踌声速流动。 本文于 19~．8年 3月 2日收到，7月 18日收到任改稿 。 维普资讯 空 气 动 力 擘 学 报 (j990年)第 8卷 二、控制方程 在计算空间，完垒气体无粘可压缩流动的控制方程可写成 “ 6一童+ 占=6+兽 一。 (1) 式 中 p PU i P “‘ puU，4-(a雪i／ax。)p 6=， P “ E ：= √ pu：u；4-(a |／Ox ) p “ pu。Ui+ ( i／ax ) P (P+p)Ui 其中，p为密度， “．、 、毗分别为物理空间w、 、 方 向的速度分量}e为单位体 积流体滞止内能，p为压力， ，一茕器 ，u．=(甍) s 以往的方法大都是以(1)式为计算空间内的控制方程。为简化物面边界处理，本文 将采用如下的以逆变速度分量为未知变量的欧拉方程为控制方程 石+啻{+一+c+#0t+等 +=o (8) 式中 p ] f