第2章-计算流体动力学分析-CFD软件原理与应用.ppt

* * 4. 混合格式 将上述两种方法综合起来，就是混合格式hybrid scheme。当|Pe| 2时，使用二阶中心差分格式； |Pe| 2时，采用一阶迎风格式。 经过一系列化简，混合格式插值方程如下所示，各物理量的意义如前所述： 优点：不会出现震荡解，无条件稳定；与高阶离散格式相比，计算效率较高； 缺点：只具有一阶精度。 其中， 原始公式： * * 5. 指数格式 指数格式（exponential scheme）是利用方程的精确解建立起来的一种离散格式。综合考虑了传导和对流两方面的影响。 其中， 原始公式： 优点：在应用一维问题时，对于任何的P值、任意网格数量，都可以获得精确解； 缺点：指数运算耗费时间；只在一维且源项为0的情况下才可以实现。 * * 6. 乘方格式 乘方格式（power-law scheme）则是与上述指数格式非常接近的一种离散格式。当P 10时，扩散项D按0计算；当 0 P 10时，单位面积上的通量按照一个多项式来计算。 其中， 化简公式： 这种离散格式的计算精度与指数格式较为接近，但比指数格式要省时。它与混合格式有着类似的性质，常作为混合格式的替代格式。使用也较为普遍。 * * 7. 各种离散格式的汇总 离散格式 系数aW 系数aE 中心差分格式 一阶迎风格式 混合格式 指数格式 乘方格式 对于一维、稳态、无源项的对流-扩散方程，最终都生成了相同的离散方程： 且有 * * 8. 低阶格式中的假扩散与人工粘性 总结： 1.本节介绍的各种离散形式，都属于低阶离散。在对流-扩散方程中，对流项引起的计算误差称之为“假扩散”（false diffusion）。因为这种离散格式截差的首项包括一阶导数，使得数值计算过程中扩散的作用被人为的放大了。也就相当于引进了人工粘性（artificial viscosity）或数值粘性（numberical viscosity）； 2.假扩散的成因：非稳定项和对流项的一阶导数离散；流动方向和网格线成倾斜交叉；非常数的源项的影响； 3.为了消除或减轻数值计算中的假扩散，可以采用截断误差较高的离散格式；或者采用自适应网格，使得网格朝向和流动方向一致。 五、空间离散的高阶离散格式 1. 二阶迎风格式 * * 上述的各种离散格式，虽然可以保证计算的稳定性、满足流动方向的要求，但一阶精度会导致假扩散。因此，本章提出高阶离散格式，它的改变主要有两点：引入更多的相邻节点；考虑流动方向性的问题。高阶离散格式可以有效降低这种误差。 二阶迎风格式与一阶迎风格式类似。但区别在于，要求出P处的未知量φ前者需要用到上游两个节点的值（WW和W），后者只需要用到上游一个点（W）的值。 * * 1. 二阶迎风格式 二阶迎风格式的对流-扩散方程的离散方程： 化简公式： 其中， 其中，当流动沿着正方向时，有Fw 0 和Fe 0，α =1；当流动沿着负方向时，有Fw 0 和Fe 0，α =0. 二阶迎风格式是在一阶迎风格式的基础上，考虑了物理量在节点间分布曲线的曲率影响。实际上，对流项采用二阶迎风，扩散项仍是中心差分。显然，其截断误差是二阶的，可以减少假扩散。 * * 2. QUICK格式 对流运动的二次迎风插值（Quadratic Upwind Interpolation of Convective Kinematics） 2.1QUICK格式的数学描述 化简公式： 其中， 其中，当Fw 0时，有αw =1；当Fe 0时，有αe =1；当Fe 0时，有αw =0；当Fe 0时，有αe =0。 * * 2. QUICK格式 2.2QUICK格式的特点及其改进格式 优点：对于与流动方向对齐的结构网格而言，QUICK格式将可产生比二阶迎风格式更为精确的计算结果，例如六面体网格（三维问题）和四边形网格（二维问题）。对于其他类型的网格，一般使用二阶迎风格式。 缺点：QUICK格式并不是绝对稳定。因此，很多人提出了改进的办法。 其中， 化简公式： 其中，当Fw 0时，有αw =1；当Fe 0时，有αe =1；当Fe 0时，有αw =0；当Fe 0时，有αe =0。 但它总能得到稳定解。 * * 3. 对高阶格式的讨论 离散格式的三个要求： （1）良好的稳定性； （2）较高的精度；