(7.2.1)--6.4-关于L-W和MacComark方法应用范围的讨论讲义.pdfVIP

《计算流体力学基础》讲义 第六章 简单的计算流体力学技术 6.4 关于L-W 和MacComark 方法应用范围的讨论 本节主要讨论LaxWendroff 方法和MacCormack 方法的应用范围，包括黏性流动、守恒 形式方程、空间推进方法等内容。前面在假设无黏流动的基础上，使用Euler 方程的非守恒 形式，我们已经描述了LaxWendroff 方法和MacCormack 方法，并且讨论了时间推进方法。 实际上，这些方法都能应用到黏性流动计算中，下面依次进行讨论。 对于黏性流动的计算来说，定常黏性流动的控制方程 Navier-Stokes 方程具有部分椭圆 型的数学特性，而LaxWendroff 方法和MacCormack 方法并不适用于求解椭圆型偏微分方程。 然而，非定常 N-S 方程在包含时间轴的任何平面上具有双曲型和抛物型偏微分方程的混合 性质，因此LaxWendroff 方法和MacCormack 方法都是适用的。事实上，MacCormack 方法 曾广泛应用于非定常Navier-Stokes 方程的时间推进求解，求解时对流项的空间导数项在预 估步和校正步依次由前向差分和后向差分来代替，黏性项均使用中心差分代替。 如果流动控制方程写成守恒形式，那么LaxWendroff 方法和MacCormack 方法也适用的。 考虑式(6.23) 所示的守恒型二维Euler 方程： U F G =− − +J (6.23) t x y 式中列向量 为通量变量， 和 为通量项， 为源项。注意的是：守恒形式方程的 U F G J 因变量不是原始变量，而是通量变量，即列向量 中的变量，所以在每个时间步后原始变 U 量必须通过求解(2.100) ~ (2.104)等式方可获得。   (2.100) u u (2.101)  v v (2.102)  (e + V 2 / 2 ) u2 + v2 e =− (2.104)  2 下面介绍空间推进方法。为了阐明空间推进的思路，将 MacCormack 方法应用于如图 6.3 所示的一维流动中，其流动方向为从左到右。 - 1 -