哈尔滨工业大学计算传热学第六章求解椭圆型流动与换热问题的原始变量法分析报告.ppt

§6－1 控制方程求解中的关键问题 A ‘Checker-board’ pressure field §６―2交错网格及动量方程的离散 A u-control volume and its neighbouring velocity components §6－3 求解N－S方程的压力修正算法 The SIMPLE algorithm §6－4 SIMPLE算法讨论 §6－5 SIMPLE算法的发展与改进 The SIMPLER algorithm 三、SIMPLEC算法 原算法的问题： 略去 ，在计算收敛时 ，故不影响解的正确性。但在计算中， ，故等于令 ，造成等式不协调 而 即 等式右边 而左边不等于0。 当扩散略而不计时，动量方程为点迭代求解。虽然点迭代收敛较慢，但它有利于克服强烈非线性问题计算中的发散。 改进方法： 因 和 同一数量级，故略去 比略 去 的影响小得多。 与SIMPLEST计算步骤的区别： 不再亚松弛，即 。 1. 2. 此法收敛速度远胜于原法。 3. 四、Date修正方案 原方法的问题： 略去 的问题，即速度修正对于对流－扩散项的影响。 1. 略去了速度修正对源项的影响。 2. 即 改进方案：预测－校正法 用原方法算出 ，得到预测值： 1. 再用新方法计算 2. 上式代入连续方程的离散形式，可得 修正方程。 其中a的表达式与 方程的相同。 获得校正值 3. 每次修正均无需亚松弛 4. 讨论： 虽然在每层计算时，Date法是原方法的两倍，但由于收敛速度大增，故总计算量减少了许多。 五、PUMPIN算法 各种算法总结： 设速度场 算a、b系数。 0. 设或计算一个压力场 1. 解动量方程得 。 2. 解压力修正方程（Date需2次）得 。修正速度 ，其满足连续方程。 3. 更新压力值。 4. （1）亚松弛 ；（2）C、D直接修正 ；（3）R解压力泊松方程。 重复1～4直至收敛。 5. 改进方案： 因为到第三步已得到满足连续方程得u、v，则可以采用一种省时的方法计算压力更新值。 只要在一个合适的参考点给定压力值，则其他各点的压力可逐一算出。（用多重路径积分获得压力更新值算法的思路）。具体计算时，从参考点出发，沿多种可能路径向待求点积分，取平均值为该点压力。 参考点的选择对稳定和收敛性有影响，需要一定的经验和技巧。 6-7 同位网格法 引入交错网格: 消除压力速度间的失藕, 但网格系统复杂化。 压力信号的传播：一种扩散过程 采用同位网格：关键在于界面上物理量的插值方式和压力修正 方程如何建立，并保证速度压力不失藕。 * 第六章 求解椭圆型流动与换热问题的原始变量法 上一章的通用方程适用于没有压力梯度项的能量、质量方程，或压力分布已知，将其并入源项的动量方程。 本章介绍速度、压力为未知量时，动量方程即N－S 方程的求解。 为了形式简洁，且与上一章形式一致，采用求和符号规则，写出控制方程 或 求解中存在的问题 一、压力场检测问题 以一维问题为例： 用中心差分 或 i–1 w i e i+1 x 当流场均匀 。压力场应为 ，但一个锯齿或台阶压力场也满足上述差分格式，因为离散方程中不包含 点，造成检测压力场的能力不强。 二、压力、速度耦合求解问题 直接求解 的代数方程组占大量内存，费机时。目前多采用分离式（顺序）求解。 但压力无控制方程，它含在动量方程的源项中，其与速度的关系隐含在连续方程里。 因此构造一个计算压力修正值的方程，使速度场满足连续方程。亦即使各变量同步修正，以提高迭代收敛速度。 一、交错网格 功用： 存放除速度 以外的其它参数， 。 a 作为除动量方程外其他方程的控制容积，以稳态能量 方程为例： b 主控制容积网格。 1. N W P E S NE u控制容积网格。 2. 以 的位置为中心，与以P为中心的主控