CFD2020-第11讲-不可压流动.pptVIP

Copyright by Li Xinliang Copyright by Li Xinliang * 计算流体力学讲义2020 第十一讲 不可压缩Navier-Stokes方程的求解 李新亮 知识点： * Copyright by Li Xinliang 拟压缩性方法 求解压力Poisson方程法 涡流函数法 Simple方法 Copyright by Li Xinliang * 11. 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 * 11. 2 人工压缩性方法（求解定常方程） 人工压缩性因子 达到定常态 流动压缩时 （ ），压力升高 流动膨胀时 （ ），压力降低 增大b 可令压力收敛加快，但会增加方程的刚性（降低时间步长）。 人工压缩性因子 相当于 Copyright by Li Xinliang * 对于定常问题，需要迭代到收敛 对于非定常问题，需要内迭代 （效率较低） Step 1 : 得到n 时间步的值 Step 2: 进行如下内迭代直至收敛 Step 3: 收敛后的V即为 内迭代收敛慢，效率较低； 通常不使用人工压缩方法解非定常问题。 Copyright by Li Xinliang * 11.3 求解压力Poisson方法 （投影法） 1） 压力的控制方程 对动量方程求散度 Poisson方程——压力的控制方程 无法时间推进 需联立求解，通常采用时间分裂法 Copyright by Li Xinliang * 2） 投影法——求解微分型压力Poisson方程 原理： 将时间推进分成三个子步， 中间步解出压力 可时间推进 不能时间推进 Step 1: 预算步 Step 2: 压力修正步 求解，得到压力p Step 3: 最终步 得到n+1时刻的V 以 1阶精度时间推进方法为例，实际上可采用更高阶精度时间推进方法： Karniadakis GE, Israeli M, Orszag SA. 1991 High-order splitting methods for the incompressible Navier-stokes equations. J. Comp. Phys. 97:414-443. Copyright by Li Xinliang * 3） 投影法——求解离散型压力Poisson方程 压力修正步： 将离散的动量方