CFD数值模拟原理-7讲稿.pptVIP

§7－5 SIMPLE算法的讨论 1. 速度修正方程中略去： 虽不影响最终解，但是加重了修正值 的负担 ——收敛速度慢 速度场收敛， 但是夸大了压力修正值，采用亚松弛 2.压力修正方程可以化成一个poisson的方程 联解，u，v，p， 可压缩流，必须用压缩流的 修正方程 3.流体流动是相对压力作用的结果，参考点的选择 在边界上可以相差任何一个常数—对流动量没有影响 ＝0 ——帕斯卡原理 4. 的求解 ：ADI 块修正 5. 收敛数据 a） b）连续方程余量的代数和 节点上余量最大绝对值 c）连续方程余量范数: 控制体上各内点求和 d）内点动量方程余量之和/入口动量 §7－6 SIMPLE算法的发展与改进 一.SIMPLER （1980） SIMPLE的 仅用于修正速度，而压力修正采用新法 提高收敛速度 思想方法：由已知速度场 直接求压力值 p 由已知的速度场求出（已知） 类似得到v速度的修正方程 上述代入连续性方程 求解上面方程，直接 p, 而非p= p*+ p′ ①速度边界条件 ②压力边界条件：易用SIMPLE方法 SIMPLER u0，v0 动量方程系数 求压力方程 u= u*+ u′,v= v* + v′ 输出 收敛 （ 连方 ） （ 动量方程 ） （ 修正压力 的求解，类似于SIMPLE 方法） （ 修正速度） 二.SIMPLEST （1981） 1.对流采用迎风格式 2.对流项全部放入源项 去 求解思想类似于SIMPLER 该方法有利于强化非线性问题的迭代过程收敛，防止扩散。 扩散项 对流项 修正时假设此两项假定不变 三.SIMPLESC SIMPLE方法中略去 太大， 两边同时减去 再采用类似SIMPLE方法 加速了修正速度，收敛性提高 （我的经验，易发散） 四.SIMPLE的Data 修正方案（1986） 原因： SIMPLE：设 修正方法 显然不完全合理 Data 思想： 预测——校正 预测阶段：同SIMPLE 校正阶段 ： 校正方程 ：动量方程的源项 中速度有关的 , 随速度 变化而变化 该法每一源项上几乎是SIMPLE计算量的两倍，但总的计算量下降。 §7－7 开口系统的流场计算 引言： 关键问题要注意出口截面流量的确定及截面上的法向速度的确定。 ①常见 ：充分发展： 但受到限制，仅用于有直段的流量。 计算方法：设邻近出口的 ②如果出口附近没有直段:采用局部单向化，如 （一） 取均匀的流场 人为按质量守恒，给出出口截面上的平均流速—粗糙，但有实用意义 （二） 从内点的速度分布来获得出口截面上的速度分布 原则：①质量守恒 ②出口截面上的每一点以其上游的一点或数点获得信息，构成已知点的初速度值。具体方法： a）假定出口截面上各点的法向速度的相对变化率为一常数，即对不同的j ． ． ． ． ． ． ． ． LI L2 j 由质量守恒，确定f 2.假定出口截面上各点的法向速度的一阶导数做为常数 由质量守恒 注意，出口截面上出现负流的问题，在计算中决不允许出现，否则必须调整出口截面的位置。 如果在出口截面上出现负速度，会影响到计算的区域内部，这也就意味出口截面离得太紧近。 处理方法 如图： 在速度面上找出最小的速度值 （最大的负速度） umin 令 umin y u (a) (b) 或 最终变成（b）所示的出口 突然扩大计算结果比较（290） Lr： 流线与大圆重接触点 一般出口截面取总长为回流区长度的4倍 Lr §十一章 网格生成技术 §11—1 引言 计算区域的不规范，复杂性，难以用比较规则的网格去划分与计算。 传统的处理方法，最普遍