CFD数值模拟原理-5

§5-3迭代法收敛性及收敛速度的讨论

1、收敛性

Jacobi与G－S迭代法收敛的一个充分必要条件：系数矩阵不可约，而且按行或按列对角占优

对各方程成立，且其中至少有一行是不等号成立。不可约：线性无关。物理意义，可约就意味着可以把流动物理区域分成两个互不影响的区域—这是不可能的。

采用本书推荐各种离散方法，弱对角占优的条件一定满足PwnesNESWΔyΔxδxwδxeδynδys由S分解为两项由此可见，采用上述方法的离散方程，迭代必收敛。2、收敛速度收敛速度：其中，使初始误差缩小成为ρ倍，n为达到ρ倍所需的步数。不同的迭代方法，收敛速度如表3、影响因素R受扫描方向影响与不同的边界条件对区域内的影响程度有关。第一类边界条件，规定了边点上值，具有最强烈的影响。（最强）第二类边界条件，（最弱）迭代方式点迭代线迭代Jacobih2/2h2G-Sh22h2SOR2hh<1h－均分网格的步长（重点格式——正方形区域）分析结果4、提高R的方法①增加迭代法中直接解法的分量，从点—>线—>ADI②适当地选择扫描的起始边界条件，最好是第一类边界条件5、加快收敛方法①亚松弛迭代：使相邻两层次间未知量的变化不太大②在稳态问题中的离散方程中加入拟非稳态项加大了系统系数矩阵的对角优势，有利于收敛例：正方形中的传热问题，达到同样的结果Jacobi：17步G－S：11步线迭代：7步分母中加入拟非稳态项§5－4加速迭代解法收敛速度的块修正技术基本思想：采用质量守恒或能量守恒进行物理参数的中间修正。如能量守恒，点守恒，块也要守恒，如图：在迭代过程的不一定完全满足能量守恒关系（对各个控制体而言）而引入修正，[如对i=I的整块（红斜影部分）]要求在整个块中能量守恒关系得到满足（积分守恒）计算区域中共有(L1-2)个修正值i=2,3,……,L1-2(L1-1)个块(L1-2)个修正值，组成了一维问题，可用TDMA求解解出i=12i=IL2L1j=1M1M2使之在块内完全满足能量方程 利用列平均值或行平均值来促进代数方程迭代收敛速度的方法称为块修正法。注意：1.块修正技术不是一个独立的代数方程解出，应与其他迭代方法结合使用。2.当离散方程求解的变量有十分明确的取值范围，不适宜。∵加后实际是使解产生整体平移，有可能无物理意义，如组分（0－1），加上后会使组分超出[0,1]的范围，无定义。§5－5TDMA算法的扩展一、双三对角阵算法（DTDMA）解两个相互耦合问题的方法时，为提高收敛性而采用的方法。例如速度与温度耦合（平面）其中，当行定后(i=const)(A)wj,Фj组成了相互耦合的一维三对角方程，具体求解方法=>(边界条件处理后)综合处理,消元法（消去Ф

j）此步的目的是使方程中仅有j与j+1Pj,Rj,Qj,Uj,：由下面的方法获得对j-1(B)将此代入(A)方程(B)

=

，，

由此的表达式

比较上述式与B方程：其中采用迭代递推。利用（B）方程，可以逐一回代，求下一步各j点上的，同一列上采用TDMA方法————双三对角算法。循环三对角阵算法（CTDMA）

周期问题

圆柱现问题

双块修正技术采用块修正时，由于块较大，有可能造成块向个别控制体（网格）上守恒关系得到破坏，反而不利于收敛，甚至发散。可以分成更小的块区，平衡修正——双块（多块）修正的思想（对列或行而言）

§5－7多重网格方法

CFX中一套网格的求解时，初始收敛快，后期收敛慢，为什么？以下面方程为例

采用G-S迭代法

设n次迭代后，与真解的差为，

——差分方程

采用Von，Neumamm方法分析前后两次迭代之间的误差矢量的变动情况

—波数波数×波长＝2与步长有关高频——短波误差低频——长波误差代入

注意：以i点为中心，在i点上：＝0以典型的之值计算r可见，为较大的量，连续3次迭代，衰减快而小的分量,仅下降到：＝0.764倍，衰减慢。可见，如果求解中，不变

采用不变的迭代，高频（短波）衰减快，低频（长波）衰减慢。

多重网格思路：（1）在较细的网格上进行迭代解，将短波误差分量衰减掉；（2）在较粗一点的网格上进行迭代，将次短波误差分量衰减掉；（3）以此逐步使网格变得越来越粗；（a）(b)(c)（4）到最后一层网格，节点数少，可直接解；（5）然后再由粗网格返到细网格；（6）反复进行。注意：多重网格使用上有灵活性，复杂区：多网格；