数值计算在热工中应用B第六章-1.pptVIP

所谓多重网格方法就是为了克服固定网格的缺点而发展起来的迭代解法。采用此法时可先在较细的网格上进行迭代，已把短波误差分量衰减掉，然后再在较粗一点的网格上进行迭代，以把次短波误差分量衰减掉。 如此逐步使网格变得越来越粗，以把各种波长的误差分量基本上都消去。到最后一层粗网格时，节点数已不多，可以采用直接解法。然后再由粗网格依次返回到各级细网格上进行计算，如此反复数次，最后在最细的网格上获得所需的解。 predominant : 6.3 线性代数方程组迭代收敛条件与加速方法 6.3.1 Jakob,G-S 迭代收敛充分条件及其分析 1. 充分条件－Scarborough准则 系数矩阵不可约且按行或按列对角占优-diagonal 对于各行成立，且至少对一行不等号成立。 2. 对流－扩散问题离散方程系数特性分析 1)系数矩阵不可约－所谓可约是指由矩阵下标元素组 成的集合W一定可以分为两个非空子集R与S，它们 满足W＝R＋S，而且从R和S中各取任一元素k及 l， 一定有：ak ,l ≡ 0 ；如果W不满足这样的条件则称不 可约－non-reducible. 分析：离散方程的系数代表了节点间的相互影响；对于 空间上为椭圆型的问题，区域中的任一点，必对其周 围的点有影响，而矩阵可约就相当于一定可以将计算 区域分割为两个互不影响的两块，这是不可能的。 矩阵不可约是流动 与热量传递过程总是互 相影响的这一物理事实 的数学反映。 2) 矩阵对角占优特性分析－本课程推荐的方法一定 满足对角占优，而且至少一行不等号成立！ (1) 非稳态全隐格式： (2) 稳态有源项： (3) 稳态无源项：对内点 从边界点必可找出至少一点满足不等号： 设Tw已知，则 TW 求解时变为： 故此点有： 对第三类边界条件的边界点，则有附件源项 如果边界点均为第二类边界条件，不可能！至少有一 点为第一类或者第三类；否则无确定解！ 本课程推荐的方法一定满足对角占优，而且至少一 行不等号成立。 6.3.2 影响迭代收敛速度的因素分析 1. 从边界条件影响传入内部的角度： 常物性无内热源物体稳态导热，控制方程为Laplace 方程 ；从均匀初场开始迭代；均匀场满足 Laplace方程，但不是问题的解，原因在于不满足边 界条件；因此边界条件影响传入的快慢必影响迭代收敛 速度。 2. 从满足守恒条件的角度： 对于第一类边界条件的问题，可以将全部边界条 件组织到迭代初场中，但是这样的场不满足守恒定律， 因此能加快守恒条件满足的方法必可加快收敛； 3. 从误差矢量衰减的角度： 迭代过程是使误差矢量逐渐衰减的过程，误差矢量 由各谐波分量组成，能使各种谐波分量协调地衰减的方 法能促使迭代过程地收敛； 4. 从增加直接求解分量的角度： 直接求解是能使各点同时满足守恒条件与边界条件 的最强方式，适当增加直接求解的分量有助于提高迭代 收敛速度。 6.2.3 加速边界条件影响传入计算区域的方法 Jakob迭代：每一轮迭代 边界条件的直接影响只能 传入一个网格，收敛速度 最慢。 G－S迭代：每一轮迭代 起始边的影响传到整个计 算区域内，收敛速度加快。 线迭代：每一轮迭代起始 边、上下边界条件的影响 全部传入计算区域，收敛 进一步加快。 ADI线迭代：每一轮迭全 部边界条件的影响均传入 计算区域，收敛进一步最 快。 ADI线迭代线迭代G-S迭代Jakob迭代 Jakob迭代收敛速度最慢，相邻两次迭代间的变化 最小，有时被用来克服外迭代不收敛的困难；SIMPLEST 算法中对流项采用Jakob迭代，在规定迭代次数的情况 有利于非线性强烈的问题外迭代收敛。 6.4.1 块修正技术的需要 6.4.2 块修正技术的基本思想 6.4 促进守恒条件满足的块修正技术 6.4.3 单块修正值的代数方程及其边界条件 6.4.4 应用块修正技术注意事项 6.4 促进守恒条件满足的块修正技术 6.4.1 块修正技术的需要 对下图示二维稳态无内热源的导热问题，采用常 规的ADI方法求解时，迭代收敛速度大大减慢：东西 方向影响最强，但离散方程系数很小；南北绝热边界 没有提供什么信息，但系数很大－因此需要一种能加 速收敛的方法。 而是一个块加一个共同的修正值， 记为 或者 要求 或者 能满足该块的守恒要求。 6.4.2 块修正技术的基本思想 从物理上看，所谓迭代收敛就是计算结果逐渐满 足守恒条件的过程；所谓迭代一轮，就是在原来不满 是代数方程求解过程。 足守恒特性的解 上加上一个修正值