高等计算流体力学-03.pptxVIP

求解流体力学方程组的有限差分方法;控制方程;贴体网格 H型 C型 O型;一、贴体坐标中的基本方程;导数的变换 ;二阶导数;二阶导数;度量系数及其计算方法 ;;;;;任意曲线坐标系中流体力学方程组的守恒形式;;;;二、守恒型Euler方程;二维守恒型Euler方程可以改写为拟线性的形式： ;其中， 为绝热指数，对于完全气体， ， 和 的定义为：;三、守恒型的Navier－Stokes方程;有限差分方法;有限差分方法;;有限差分方法的主要特点是： 1）有限差分方法只需构造偏导数的离散方法，这使得它比较容易推广到高阶精度，对于多维问题也是如此。 2）有限差分方法对网格的光滑性有较高要求，不易在复杂形状求解域上实施。 3）在曲线坐标系中，有限差分方法要对几何量和物理量的确定组合进行差分离散，这样做的后果之一是有限差分方法可能产生所谓几何诱导误差（geometrically induced errors） 。例如，考虑一无 界均匀流场，易知，该流动应是定常的。即 。在由一般曲 线坐标构成的网格上采用有限差分方法计算此流动，可能有 （三维 以及动网格中的有限差分方法更容易出现这种情况）。而采用有限体积方 法，则恒有 。有限差分方法在某些情况下不能复现均匀流 场，是几何诱导误差的表现之一。;如何计算数值通量？;激波间断和广义解;;一、激波的形成;考虑Bergers方程;;二、广义解; 下面讨论一下广义解的意义。 如果已知 是光滑的，设 围成的 区域为 ，则由 式利用Green公式知;;间断关系的推导;以一维Euler方程为例，R-H关系为：;三、熵条件;另外，考虑;如何在不唯一的广义解中挑选所谓物理解呢？最直观的方法是认为物理解是 一个光滑的粘性问题的解在粘性趋于零时的极限。也就是说：;这种选择物理解的方法，虽然非常直观，但并不实用，因为粘性问题的解 一般也是未知的。对于标量非线性双曲型守恒律，Oleinik 提出了下面的熵 条件：;40;;含激波的流场的计算方法;;;激波捕捉方法的基本思想是：在计算包含激波的流场时，采用统一的计算 格式，??对激波进行任何特殊处理。 当计算格式满足一定要求时，可以自动计算出流场中的间断。 这种方法的优点是：编程计算比较简单，且可以适用于具有任意激波结构（包括激波的产生、消失、运动等各种情况）的流场计算。 其缺点是：此时激波不再是理想的间断，而是厚度为一至数个网格宽度的连续变化的结构；另外，由于激波内部流场的梯度较大，常规的计算格式或者会把激波“抹平”，或者在激波附近产生虚假数值振荡。 因此，为了高分辨率、无振荡的计算激波，对数值计算方法提出了很高的要求。 事实上，能有效计算激波的高分辨率格式，是最近30年以来，计算流体力学研究的热点问题之一。 随着高分辨率格式研究的进展，目前，激波捕捉方法已经成为计算包含激波的可压缩流场的主流方法。;一、守恒格式和Lax-Wendroff定理;守恒格式的特点是： 对于任意的 ，数值通量用相同的函数关系计算， 且函数涉及的自变量（离散点上的数值解）与 的相对位置关系不变。 这一特点使得相邻网格点，如 和 点上的计算格式在公共界面 处的数值通量是相同的。稍后我们将说明，正是这个特点使守恒格式具有自 动捕捉间断的能力。;守恒格式的相容性条件;;因此，;整理上式，可得：; 任取一与间断只有有限个交点并包含在;;;二、人工粘性和格式粘性; 在我们已经介绍过线性波动方程 的修正方程。修正方程 的偶数阶导数项影响数值解的振幅，当 等关系成立时，偶数阶 导数项起耗散作用，相当于计算格式中的数值粘性。 某些差分格式本身就带有数值粘性这部分数值粘性称为“格式粘性”。 某些格式缺乏必要的格式粘性（如中心差分格式）或者格式粘性不足以 抑制间断附近的非物理振荡时，常常需要在计算格式中添加显式 的数值粘性项 ，称为“人工粘性”。根据修正方程的性质，人工粘性项一般可以 写成 或者 的离散形式。 等称为人工粘性系数。;;考虑Euler显式格式;注意到： ;线性波动方