松弛迭代实施方法对收敛的影响.pdfVIP

中国工程热物理学会 传热传质学 123341 学术会议论文 编号： 松弛迭代实施方法对收敛的影响 1,2 3 3 3 1* 刘豫龙 ，龚金海 ，刘德绪 ，莫琼 ，徐明海 1. 266555 （中国石油大学（华东）储运与建筑工程学院，山东青岛经济技术开发区 ，徐明海， minghai@ ） 2. 2570872 （胜利石油管理局胜利发电厂，山东东营 ） (3.中原石油勘探局勘察设计研究院，河南省濮阳市中原路332 号 457001) 摘 要：本文通过导热、方腔自然对流和管内入口段对流换热三个算例说明了两种松弛实施技巧对离 散方程收敛性的影响，结果表明把松弛融合到离散方程中的实施方法在亚松弛的情况下，大大降低了 收敛的速度，在超松弛的情况下可以加速收敛，存在最佳松弛因子，超过最佳松弛因子不但不能加速 收敛，还会导致迭代过程振荡或发散。而先解离散方程后松弛处理的方法在松弛因子的有效变化范围 内，也存在最佳松弛因子，但迭代次数与松弛因子变化不太敏感，且最佳松弛因子较前一种方法稍大 些。亚松弛的情况下应采用后松弛而超松弛时应采用先松弛后解方程组的实施方法。一维稳态导热的 算例表明，松弛因子为0.2 时后松弛所用的迭代次数不到先松弛迭代次数的六分之一，而方腔自然对 25 流换热问题计算需要的时间在小松弛因子的条件下差别可达 倍，管内对流换热相差更大。松弛因子 1 接近于 时两种方法的迭代次数趋于相同，而超松弛时先松弛所用的迭代次数小于后松弛，达到最佳 松弛因子效果最好。 关键词：对流扩散方程；松弛迭代；收敛性 0 引言 流体流动与传热的控制方程可用对流扩散方程统一描述，通过改变扩散系数、源项 和密度的实际代表意义，可以表述不同的实际对流扩散过程。对流扩散方程由于其非线 性和各变量耦合特点，一般采用迭代法求解。迭代法中往往采用松弛技术提高收敛速度 （超松弛）或者提高算法稳定性（亚松弛）。松弛技术的具体实施目前有两种方法，其 一为Patankar 教授[1]提出的把松弛因子融合到离散方程中，求解离散方程得到松弛后的 解。另一种实施方法先求解离散方程获得离散方程的解，然后将解进行亚松弛处理，例 如，求解不可压流动与传热的经典分离算法SIMPLE，就是速度和其它标量采用先把松 弛结合到离散方程中去，而压力则采用解完压力方程再松弛的技术。这两种松弛处理方 法在计算方法教科书[2,3,4]或者Patankar 的著作中被认为是一致的，然而这个结论并没有 进行严格的数学证明，也没有具体的数值试验报道。本文以稳态无源导热和流体流动与 传热问题的离散方程求解过程两个方面的实例说明了这两种松弛实施的方法并不一致。 1 两种松弛处理实施 二维直角坐标系下的稳态对流扩散方程一般表示为： ur div (ρVφ ) div (Γgradφ ) + Sφ (1) 其中包括对流项、扩散项和源项。 N n w