CFD第五章 流动计算.pptVIP

第五章 问题一二三 交错网格法 第二节动量方程 动量方程-1 SIMPLE算法 Simple-1 连续性方程离散化 压力修正方程 实施步骤 边界条件 压力修正边界条件 压力的相对性 计算实例 交错网格的实施 速度分布 压力分布 压力三维 速度star 速度star-1 压力分布star * 不可压缩流体定常流动的动量方程 是否已经可解？ H S J 第五章 流动计算 第一节 困难与办法 连续性方程 问题一：压力如何与速度相关联 H S J p=10 p=1 问题二：压力梯度如何表达 P E W N S 问题三：速度梯度如何表达 H S J 解决办法之一：交错网格法 u v p u、v、p各在不同的点上计算 H S J 第二节 动量方程 令 ： 即得离散的动量方程 将压力项单独列出，便于计算： H S J Simi-Implicit Methd for Pressure-Linked Equations 先假定压力场求速度场，再用连续性方程修正压力 H S J 第三节 SIMPLE算法 一、基本思路 二、推导过程 H S J 将不精确的值表示为： u=u*+u? v=v*+v? 将修正关系表示为： 对修正值也有： e n P w s E N W S 化简 H S J 连续性方程的离散化（定常不可压） 整理，得压力修正方程 H S J 注意b的表达式！ 注意事项 对u、v和p进行计算时，其单元体是不同的，因此方程中各系数的表达式中的?x、 ? y和? x、 ?y可能是不同的。 有些值需要进行插值计算。 三、SIMPLE算法的实施步骤 1、估计压力场p*； 2、求解动量方程，得u*、v*、w*； 3、求解压力修正方程得p； 4、计算p＝p*＋p； 5、修正速度u＝u*＋u、v＝v*＋v、w＝w*＋w； 6、计算其它物理量（如温度、浓度、湍流度等）； 7、将p 视为新的p*，返回第二步。 H S J 四、边界条件 1、动量方程的边界条件 H S J 速度边界条件 固壁边界 速度进口边界 出口边界 单向坐标法 速度已知 压力边界条件 压力修正方程中，aE=0 H S J 2、压力修正方程的边界条件 p’=0 速度（流量）边界条件 ue ue 五、压力的相对性 绝对压力不确定 压力差计算的相对误差增大 当边界条件不包含压力的时候，给定一点的压力为零 H S J 六、计算实例：方腔驱动问题 H S J u＝1 v＝0 固 体 边 界 1×1 流体域 ?＝1.0 ?＝1／Re 交错网格的实施 H S J v2 u1 v1 u2 速度分布 H S J 压力分布 H S J 压力分布 H S J 速度分布：Star－CD（单元） H S J 速度分布：Star－CD（顶点） H S J 压力分布：Star（顶点） H S J *