计算流体力学数值方法研讨.ppt

第三章 数值方法 主要内容 空间离散技术 标量输运方程 动量方程 时间离散技术 1 标量输运方程 有限体积法 一维对流扩散方程 扩散项的离散 源项的离散 代数方程的组装 二维和三维问题 对流项离散基础 离散特性 对流项高级离散方法 高阶对流方法的实现 曲线网格 边界条件 代数方程的求解 小结 对于任意控制体内一个守恒物理量: 控制体内物理量的变化率+经过控制体边界面的净通量＝控制体积内的源(汇) 经过控制体边界面的总的通量由对流(随流动的迁移)和扩散(由随机分子运动或湍流运动导致的净输运)两部分组成。如果用 表示单位质量流体的守恒量，则通用标量输运方程(或称对流－扩散方程)可表示为： 1，有限体积法直接对上式进行离散 2，本章只考虑稳态问题，即上式左边第一项为零 2 动量方程 动量的特征 压力－速度耦合 压力校正法 小结 每一个流体流动问题都满足质量和动量守恒方程。 而每个动量分量满足标量输运方程，对于单个控制体积，标量输运方程为： 当输运量为动量时： 3 时间相关的方法 时间相关的标量输运方程 单步法 多步法 CFD中时间步进法的使用 小结 交错网格的优点： [动量方程]:现在压力正 好储存在需要计算压力的 节点上： [质量方程]:现在速度分量正好储存在需要计算质量通量的点上： 可以看出，在交错网格技术中，pi的贡献不再被抵消，所以，也就不会出现奇偶解耦的问题。而且，得到的压力方程与前面的标量输运方程具有相同的形式 ，所以可以使用相同的子程序。 交错网格的缺点： 由于有多个不同的节点和 控制体积单元组，增加了几何 复杂性； 对于非笛卡儿网格，速度 节点会不在压力节点之间，如 右图所示。 (d) Rhie-Chow 速度插值 (Rhie 和 Chow,1983) 另外一种方法就是仍然 采用速度和压力同位储存， 但对用来计算质量通量的对 流速度采用不同的插值方法。 动量方程使得单元中心的速度与单元面上的压力联系起来： 即， Rhie-Chow的思路就是对 压力项和非压力项分别插值。 式中的上划线表示平均节点值得到上的值(严格来说是线性插值)。等价地(或者更方便地)加减中心压力差 基本上只是加减了一个与压力梯度成比例的量： 除了方程右端项第一项采用单元中心离散值的平均外，第二项采用经过单元面离散值的中心差分。采用上述插值技术后，质量方程为： 可以看出： 中心压力值不再被抵消，因而没有奇偶节点解耦问题； 这是一个压力方程，而且与标量输运方程的离散方程形式相同： 压力－速度耦合是纳维司托克斯方程的重要特征，对于笛卡儿网格，交错网格技术是最有效的办法，但是，对于非笛卡儿网格(曲线网格)同位存储网格是最常用的，同时也是几乎所有通用CFD软件采用的方法。 压力校正法 什么是压力校正法： (1) 它是一种求解不可压缩流的数值迭代方法 (2) 用于导出同时满足质量和动量守恒的速度场和压力场 (3) 包括下列速度场和压力场的更新： 首先，固定压力，通过动量方程求解速度场。 然后，计算压力校正值， 并用于校正速度值使之满足质量守恒方程 (4) 有两种常用的校正方法:SIMPLE和PISO 压力的变化可用于强制满足质量守恒： 净流入质量通量 净流出质量通量 使单元压力增加 使单元压力减小 1, SIMPLE法(Patankar,Spalding,1972) SIMPLE – Semi-Implicit Method for Pressure-Linked Equations 对于稳态情况： 第一步：用当前压力值 求解动量方程： 这样求得的速度一般是步满足 质量守恒方程的。 第二步：导出压力校正公式： (i)使速度的变化与压力的 变化联系起来：