CFD数值模拟原理-5.pptVIP

§5-3 迭代法收敛性及收敛速度的讨论 1、收敛性 Jacobi与G－S迭代法收敛的一个充分必要条件：系数矩阵不可约，而且按行或按列对角占优 对各方程成立，且其中至少有一行是不等号成立。 不可约：线性无关。物理意义，可约就意味着可以把流动物理区域分成两个互不影响的区域—这是不可能的。 采用本书推荐各种离散方法，弱对角占优的条件一定满足 2、收敛速度 收敛速度： 其中，使初始误差缩小成为ρ倍，n为达到ρ倍所需的步数。 不同的迭代方法，收敛速度如表 3、影响因素 R受扫描方向影响与不同的边界条件对区域内的影响程度有关。 第一类边界条件，规定了边点上值，具有最强烈的影响。（最强） 第二类边界条件，（最弱） 4、提高R的方法 ①增加迭代法中直接解法的分量，从点—线—ADI ②适当地选择扫描的起始边界条件，最好是第一类边界条件 5、加快收敛方法 ①亚松弛迭代：使相邻两层次间未知量的变化不太大 ②在稳态问题中的离散方程中加入拟非稳态项 加大了系统系数矩阵的对角优势，有利于收敛 例：正方形中的传热问题，达到同样的结果 Jacobi：17步 G－S：11步 线迭代：7步 §5－4 加速迭代解法收敛速度的块修正技术 基本思想：采用质量守恒或能量守恒进行物理参数的中间修正。 如能量守恒，点守恒，块也要守恒，如图： 利用列平均值或行平均值来促进代数方程迭代收敛速度的方法称为块修正法。 注意： 1.块修正技术不是一个独立的代数方程解出，应与其他迭代方法结合使用。 2.当离散方程求解的变量有十分明确的取值范围，不适宜。 ∵加 后实际是使解产生整体平移，有可能无物理意义，如组分（0－1），加上后会使组分超出[0,1]的范围，无定义。 §5－5 TDMA算法的扩展 一、双三对角阵算法（DTDMA） 解两个相互耦合问题的方法时，为提高收敛性而采用的方法。 例如速度与温度耦合（平面） 其中，当行定后(i=const) 此步的目的是使方程中仅有j与j+1 Pj,Rj,Qj,Uj,：由下面的方法获得 对j-1 其中 采用迭代递推。 利用（B）方程，可以逐一回代，求下一步各j点上的，同一列上采用TDMA方法————双三对角算法。 循环三对角阵算法（CTDMA） 周期问题 圆柱现问题 双块修正技术 采用块修正时，由于块较大，有可能造成块向个别控制体（网格）上守恒关系得到破坏，反而不利于收敛，甚至发散。可以分成更小的块区，平衡修正——双块（多块）修正的思想（对列或行而言） §5－7 多重网格方法 CFX中一套网格的求解时，初始收敛快，后期收敛慢，为什么？以下面方程为例 采用 G-S 迭代法 设n次迭代后， 与真解 的差为 ， 采用Von，Neumamm方法分析前后两次迭代之间的误差矢量的变动情况 —波数 波数×波长＝2 与步长有关 高频——短波误差 低频——长波误差 代入 以典型的 之值计算r 可见， 为较大的量，连续3次迭代， 衰减快 而 小的分量,仅下降到： ＝0.764 倍，衰减慢。 可见，如果求解中， 不变 采用 不变的迭代，高频（短波）衰减快，低频（长波）衰减慢。 多重网格思路： （1）在较细的网格上进行迭代解，将短波误差分量衰减掉； （2）在较粗一点的网格上进行迭代，将次短波误差分量衰减掉； （3）以此逐步使网格变得越来越粗； （a） (b) (c) （4）到最后一层网格，节点数少，可直接解； （5）然后再由粗网格返到细网格； （6）反复进行。 注意：多重网格使用上有灵活性，复杂区：多网格； 简单区：少网格。保证精度，节省时间。 * * P w n e s N E S W Δy Δx δxw δxe δyn δys 由S分解为两项 由此可见，采用上述方法的离散方程，迭代必收敛。 2h SOR 2h2 h2 G-S h2 h2/2 Jacobi 线迭代 点迭代 迭代方式 h1 h－均分网格的步长 （重点格式 ——正方形区域） 分析