6.4 原始变量顺序求解流场的压力修正方法.pptVIP

压力修正中的亚松弛 引入亚松弛： 松弛因子 Patankar Demirdzic Peric 修正速度方程不需要进行松弛 4．SIMPLE算法中的迭代 三重迭代含义 1. SIMPLE算法中动量方程、压力修正方程、速度修正方程之间的迭代 2. 非线性方程线化之后求解的迭代 3. 线化之后的离散动量方程、离散压力修正方程的迭代求解 外迭代 内迭代 5 迭代收敛判据 1) 停止内迭代的判据 内迭代: 采用迭代方法求解系数和非齐次项暂被固定冻结而得以线化的代数方程组(临时的线化方程组) 可以“适时地”停止 停止内迭代往往以压力修正值p’方程为依据 (1) 如果采用交替方向线迭代与块修正，则规定运算轮数即可 (2) 规定方程余量范数小数 范数：函数空间中某函数大小的量度。Euclid范数：平方求和再开方 判据： (3) 规定余量范数缩小的比例 经过k次迭代后的余量范数 初始的余量范数 2）停止外迭代的判据 动量方程、压力修正方程之间的迭代 非线性方程系数冻结线化的迭代 压力修正值p’方程 其中非齐次源项 bp’ 代表 u*，v* 不符合 连续方程的衡量： (1) bp’绝对值之和 或 绝对值最大者 某参考质量流量 或者： (2) bp’ 的范数小于某个参考流量 (3) 动量方程余量绝对值或范数小于某个参考动量 或： (4) 令某特征物理量在连续数个层次迭代的相对偏差小于允许值 Nusselt 数 6.4 原始变量顺序求解流场的压力修正方法 胡 茂 彬 /~humaobin/ humaobin@ 6.4.1 压力修正方法的基本思想 压力修正方法是为了解决流场中压力P没有单独的控制方程，无法求解的问题 压力修正方法发挥了连续方程的作用，通过连续方程推导得到了压力修正值方程 分 析 控 制 方 程 有了压力场，才可以计算速度场 然而，连续方程没事可干 假设压力场p* 压力修正方法 流 程 图 求解动量方程 得到u*，v* 利用连续方程 改进压力场 得到p*+p’,u*+u’, v*+v’ 进入下一时层求解 解收敛否？ Yes No 以新压力 代替老值 关 键 要求使改进的压力场 p*+p’ 对应的改进速度场 u*+u’, v*+v’ 能够满足连续方程 如何构造压力修正值 p’ 方程? 如何得到修正速度 u’, v’? 6.4.2 速度修正值的简化近似计算 u*, v*, p* 必满足离散的动量方程： 期望：修正后的速度也能满足动量方程 速度修正值u’, v’方程 邻点的耦合项 压力修正值梯度项 大胆的假设 修正的速度 6.4.3 将连续方程离散式转化为压力修正值p’方程 连续方程： 全隐式离散： 将修正的速度方程代入 压力修正值p’方程 其中： 比较 连续方程离散形式： 为0，则u*, v* 已经满足连续方程，即方程组已经收敛 的绝对值反映了靠近收敛的距离，迭代的过程就是不断减少这个值 6.4.4 p’方程的边界条件 1、给定边界压力分布，速度待求： 2、给定边界速度 6.5 SIMPLE算法 压力关联方程的半隐式方法 (Semi-Implicit Method for Pressure Linked Equations) Patankar,Spalding (1972) 半 隐 Semi-Implicit 设定初始速度场u0,v0,计算系数 流 程 图 （非稳态情况） 设定初始压力场p* 求解动量方程， 得到u*, v* 求解压力修正值方程, 得到p’, u’, v’ 计算其它物理量 进入下一时层求解 解收敛否？ Yes No 以新值 替代老值 SIMPLE算法的三个简化假设 1半隐假设：速度修正值计算没有计及邻点速度修正值的影响 2 冻结系数法线化代数方程组： 3 初始速度场 u0, v0 和初始压力场 p* 的设定各自独立，一般不可能匹配 6.5.2 SIMPLE算法若干问题讨论 1．算法的简化假设不影响流场最终收敛解，但影响收敛的速度 三个简化假设 1 半隐假设：构造能够迭代求解的方法 2 冻结系数法线化代数方程组：处理非线性问题的基本做法 3 u0, v0 和 p* 不匹配：解决迭代起步问题 2．P’方程的数学特征 压力修正方程 如果p’是方程的解，p’+C 也符合方程 P’ 方程的多值特性 边界条件不能唯一确定 p’ 值 给定边界压力或给定边界法向速度，都令其方程中与边界相应的系数为零，结果p’还是不能唯一确定 P’ 方程的多值特性 实际上，由该方程确定 p’ 的绝对值意义不大，重要的是 p’ 的相对值 一般做法 在每个迭代层次求解 p’ 的方程之前，令整个初场的 p’ =0 ，这样做可