第12讲-不可压流动.pptVIP

Copyright by Li Xinliang Copyright by Li Xinliang 计算流体力学讲义 第十二讲 不可压缩Navier-Stokes方程的求解 李新亮 lixl@ ；力学所主楼219； 知识点： * 讲义、课件上传至 (流体中文网） - “流体论坛” -“ CFD基础理论 ” 讲课录像及讲义上传至网盘 /browse.aspx/.Public Copyright by Li Xinliang 拟压缩性方法 求解压力Poisson方程法 涡流函数法 Simple方法 Copyright by Li Xinliang * 知识回顾 一、 代数方程组的求解 直接法 Gauss消元法 LU分解法 追赶法： Copyright by Li Xinliang * 迭代法 Jacobi迭代 n+1 n n n n Gauss-Seidel迭代 n+1 n n+1 n+1 * n+1 n+1 n n+1 n+1 n+1 n n n LU-ADI LU-SGS 知识回顾 Copyright by Li Xinliang * 知识回顾 二、 网格生成 A B C D E F A’ B’ C’ D’ E’ F’ 物理空间 计算空间 1. 代数网格生成法 2. 解椭圆型方程网格生成法 Copyright by Li Xinliang * 12. 1 不可压缩Navier-Stokes方程的特点 密度为常数的不可压缩Navier-Stokes方程组： 特点： 动量方程与能量方程解耦 压力属于约束变量而不是发展变量 温度对密度的影响可忽略不计 压力不能时间推进求解 可压缩 不可压缩 易处 压力可推进求解，易于使用显格式 压力方程具有椭圆性，无法推进求解。压力方程收敛性差 难处 可能出现间断 不会出现间断 研究重点 激波捕捉 压力处理 Copyright by Li Xinliang * 概念澄清： 压力—— 动力学压力及热力学压力 动力学压力—— 应力的中各向同性部分 连续介质微元体的受力平衡： 应力的概念 静止流体或无粘流体中力的平衡—— 动力学压力的概念 热力学压力—— 由分子动力学性质决定 —— 状态方程 热力学压力： 分子对固壁的碰撞，产生压力 完全气体： 动力学压力 热力学压力 可压缩N-S方程： 动力学与热力学耦合；动力学压力= 热力学压力 不可压缩N-S方程： 动力学与热力学解耦 由不可压缩条件确定压力 （纯动力学概念） 1） 压力的处理原则 Copyright by Li Xinliang * 奇偶失联与交错网格 压力项，通常采用中心差分离散 极端情况： 棋盘式压力场 高压 低压 特点： 高压-低压点间隔分布 采用中心差分格式计算出： 流场竟然“保持稳定” “奇偶失联” Copyright by Li Xinliang * 常用措施： 交错网格 压力p 速度 u 速度v 交错网格示意图 在u的网格点上离散 在v的网格点上离散 注： 对流项通常采用迎风格式离散 后差 前差 Copyright by Li Xinliang * 2) 对流项的处理原则 关系式1： 关系式2： 兰姆-葛罗米柯等式 总压 表达式 优点 不足 普通型 简单，易于迎风 有混淆误差 守恒型 简单，守恒 有混淆误差 旋度型 混淆误差小 计算量大 螺旋-对称型 混淆误差小 计算量略大 普通型 守恒型 Copyright by Li Xinliang * 12. 1 人工压缩性方法（求解定常方程） 人工压缩性因子 达到定常态 流动压缩时 （ ），压力升高 流动膨胀时 （ ），压力降低 增大b 可令压力收敛加快，但会增加方程的刚性（降低时间步长）。 人工压缩性因子 相当于 Copyright by Li Xinliang * 对于定常问题，需要迭代到收敛 对于非定常问题，需要内迭代 （效率较低） Step 1 : 得到n 时间步的值 Step 2: 进行如下内迭代直至收敛 Step 3: 收敛后的V即为 内迭代收敛慢，效率较低； 通常不使用人工压缩方法解非定常问题。 Copyright by Li Xinliang * 12.2 求解压力Poisson方法 （投影法） 1） 压力的控制方程 对动量方程求散度 Poisson方程——压力的控制方程 无法时间推进 需联立求解，通常采用时间分裂法 Cop