计算流体CFD2004-6.pptVIP

第六章 可压缩流动的有限差分法 6.1 位势方程 6.2 Euler方程 6.2.5 高分辨率的二阶迎风格式（TVD格式） 6.3 N-S方程组 动量预测 PISO步骤 6.5边界条件 * 华中科技大学能源与动力工程学院 在可压缩粘性流动中，激波紊流边界层的相互影响就构成了计算流动动力学中最重要的物理现象之一。空气在大于100m/s（相应的马赫数近似为0.3）又小于1700m/s（马赫数近似为5）的范围内（即0.3≤M≤5），流动可以认为是可压缩非粘性的。 为了考虑高速流动的可压缩性和密度的变化，应用密度相关控制方程的守恒形式。对于可压缩流动，也会遇到一些流动区域是低马赫数或不可压缩的（如近壁面）。在这种情况下，密度相关的方程会不适合于计算，并造成结果收敛很慢。 用位势方程比较好分析，但是一些重要的物理现象如旋涡，非等熵，非等温的影响在控制方程中未考虑。因此，绝大多数分析可压缩非粘性流都采用Euler方程。 可压缩非粘性流包括旋涡，非等熵，非等温的影响，需要同时求解连续性方程，动量方程和能量方程。除跨音速流外，位势方程中的小扰动或非线性化处理在这种方法中是不允许的。因此，在激波间断的跨音速流中遇到的困难需要特殊的计算格式求解。 Euler方程求解的最基本要求是：求解格式能提供适当的人工粘性以使解迅速收敛于精确解。因此，与对流项相关的特征值和相容性在激波和膨胀波的求解中是重要的因数。 Euler方程的求解格式可以分为三种主要形式：（1）中差格式（2）一阶迎风格式（3）二阶迎风格式和基本不振荡格式。 CFD最重要的发展也许就是高分辨率的二阶迎风格式，它被称为总变差减少(TVD)格式 熵条件——不允许膨胀激波中存在熵的减少。 单调性条件——为了防止数值格式中的震荡现象。 总变差减小（TVD）——任何物理上可容解的总变差不允许增加。 （1）高分辨率格式的定义 （2）具有限制器的TVD格式 （3）TVD的时间积分方法 6.3.3 可压缩流体的PISO格式 6.3.1 显格式 6.3.2 隐含格式 方程组联合求解 方程组独立求解 在处理可压缩流体时，可以按照与不可压缩流体的PISO格式类似的步骤，但要引入一个附加的修正，这是由于动量，能量，压力（也就是连续性）方程的耦合包含了密度和温度。 SIMPLE(Semi-Implicit Method for Pressure-Linked Equations) PISO“Pressure implicit with Split--Operator”(分裂算子的压力隐式) δ 动量第一次修正 ρ ρ （a）动量的预测 （b）动量的第一次修正 (c) 能量的预测 (d) 动量的第二次修正 （e）能量修正 （f）动量的第三次修正 最后求得确定， 能量预测 c- c0 c+ c+ c0 c- 进口：三个物理边界条件 出口：三个数值边界条件 原始变量的Dirchlet边界条件的数目等于雅可比矩阵的正特征值的数目，即物理边界条件的数目。负特征值代表着从流域中外推出来的数字边界条件。 图6.5.1中在进口点，特征C0和C+的斜率分别为u和u+a ，通常它的值在x轴的正方向都是正的。第三个特征C-的斜率取决于进口马赫数的大小。在超音速进口处，是正值，而在亚音速流中它是负值。因此，在亚音速进口处，不需要给定和C-相关的边界条件。同理，在出口不需要为C+和C0指定边界条件，而对C-而言，必须对亚音速出口提供边界条件（而不为超音速出口提供边界条件）。 * 华中科技大学能源与动力工程学院