热物理过程的数值模拟-计算传热学3汇总

四、非线笥问题迭代式解法的收敛性

每一层次上满足迭代法求解的收敛条件+相邻次间代数方程的系数变化不太大（亦即未知量

的变化不太大J多数情形下非线性问题迭代式解法是可以收敛的）。

使相邻两层次间未知量变化不太大的措施：

1欠松弛迭代常用逐次欠弛线迭法（SLUR）：一组临时系数下逐线迭代求解+对所得的解

施以欠松弛，再用欠松弛后的解去计算新的系数，常数，以进入下一层次的迭代。

实施：常把欠松弛处理纳入迭代过程，而不是在一个层次迭代完成后再行欠松弛。

.（n 1）川）.'a n bt n b

t p =t p （t p ）

a p

（先）t p n1） = 7a n b t n b b （1一•）屯t p n）

co o

a'p t p n 9 、a n bt n b b'

a'p -a^ ■, b' = b （^ ）（a p ）t p n），用交替方向线迭代法求解这一方程，就实现了SLUR

的迭代求解。为一般化起见，上式中t n b上没有标以迭代层次的符号（J, GS时不相同）。

2、采用拟非稳态法

前面已指出，稳态问题的迭代解法与非稳态问题的步进法十分相似。对于非线性稳态问题，

从代数方程的一组临时系数进入到另一组临时系数亦好象非稳态问题前进了一个时间层，非稳态问题的物理特性：系数热惯性越大（a； = PM v/也I ），温度变化越慢，仿此，对稳态非线性

问题，可在离散方程中加入拟非稳态项，以减小未知量托两个层次间的变化，即

由

（=a n b -S p：V）t p n。= Ua n bt n b b=（3a n b - S p：V a；）t p n。=二a n bt n b b a p tf

Za n bt n b - b - a；t p n）

（n 1）

t p o

Ea n b -S p心V +a p

一直进行到t p,t n b收敛，虚拟时间步的大小通过计算实践确定。

3、采用Jacobi点迭代法

中止迭代的判据（该层次迭代）除前述变化率判据外，还可以规定迭代的轮数，例如规定进

行4-6次ADI线迭代就结束该层次上的计算。此时，用收敛速度低的丁迭代也就起到了欠松弛的作用。

五、迭代法的收敛速度

1收敛速度

对给定的代数方程组（包括是临时系数的情形），采用不同的迭代方法求解时，使一定的初始误差缩小成：•倍所需要的迭代轮数K是不相的

e" 32 e (0)

因为G 是对称的，所以 G 2 = . T (G T

G)=匸(G 2) = "G) 所以 |門2 斗G k ||2|e (0)||2 = P(G))k |e (0)|2

a =|e ⑹|2/|e (

T 2 兰(P (G))k ，E 兰kl.P(G)

即 k _1八 /嶄 MG) - 一1「/(—嶄 ”G))

:- , ?(G)均 <1

对于给定的：•，所需迭代轮数k 与一 l n T(G)成反比，规定用

R = _l n 「(G)

表示迭代法的收敛速度，则

k - -l n ： /R or R _ -l n ： /k

即所需迭代轮数与收敛速度成反比，收敛速度又与谱半径成反比，收敛速度愈快，迭代轮数 愈少。 注意：不同的迭代方法每进行一轮迭代所需的运算次数不同，最终所需的计算时间的多少取 决于迭代轮数及每一轮迭代所需的时间。

2、收敛性的定性分析

为什么不同的迭代方法的收敛速度不同，亦即为达到满足一定精度要求所需的迭代轮数不 同？以二维常物性、无内热源、稳态导热问题来进行讨论。

1B 、C a p t p 二a E t E a w t w a N t ” a $t s

迭代法需要假定一个初场，例如假定一个均场， 从微分方程

为看，均匀是其一个解，但却不是所研究 问题的解，为什么？因为它虽然满足内部节点上的离 散方程，却不满足与边界有关的节点的离散方程(图 中红点)，即不满足边界条件。 所以迭代法的实质是要 通过迭代，尽快建立起与边界条件相适应的

••变量场,

关键：必须使B 、C 的影响迅速传入计算区域内部， 以改进节点••变量值，尽快与B 、C 相应，B 、

C 的影响传入愈快，逼近真解就愈快，收敛就越快！

B 、

C 的影响传入计算区域内部的快慢与哪些因素有关？ (1)

与迭代方法有关

J 迭代：节点温度的更新均用上轮迭代所得的“旧”值来计算，所以完成一轮迭代后，

2 ■ 2

y

y

X

t BB

t RB

的影响只能传入与边界相邻的一批节点上，即仅可传入一个网格，且扫描方向与收敛快慢无关。

要在以后各轮迭代中，B、C的影响才由这些节点逐步向内渗透，所以收敛慢。

GS迭代：假设从左向右扫描，则每做完一轮迭代，左边界和下边界的影响传遍全区域，而右边界的影响只能传入一个网格，且收敛速度受迭代扫描方向的影响。

线迭代：GS线迭代。自左T右扫描，完成一轮迭代不仅左边界的影响逐步传入，而且在每一列的直接求解中，上、下边界的影响全部传入到该列的各节点上，即一轮迭代使左、上、下边界影响传入全区域，但右边界影响仍仅传入一个网格。

ADI :—轮迭代包括一次逐行、一次逐列的扫描；所以在每一轮迭代后所有边界的影响均传入计算区域内部，从而加快了收敛速度。

收敛速度的比较，正方形区域，1B、C, Laplace方程五点格式，均匀网格步长为h。

迭代方法点迭代线迭代

Jacobi h2/2h2

Gouss-Seidel h22h2

SOR2h 2 J2 h

(2)与边界条件的性质有关

x XX

1B、C规定了边界节点的温度，影响直接传入计算区域内部；

t

3B、C规定了环境温度及定向点位置，t - t -,对边界温度的限定程度比1B、C时

l cX

弱，所以对内部的影响也较弱；或将t f视为外部温度，其对计算区域内部的影响被外部换热热阻

削弱，而1B、C可视为「，一•-■或外部热阻一；0的极限情况，故3B、C的影响比1B、C时弱。

2B、C仅规定了壁面的钭率，壁温完全不确定，对内部节点温度值的改进提供的信息最少，收敛最慢。

可见，为了提高代数方程迭代解法的收敛速度，应力求使边界条件的影响迅速传入计算区域内部，措施：

①增加迭代解法中直接解法的成分，从点迭代T线迭代T ADI ;

②适当选择扫描的始边，多以1B、C或3B、C的边界为始边，少以2B、C (尤其是绝热边

界)为扫描始边。

5-4不规则区域的处理一网格生成技术

如何对不规则区域进行有效的处理，以便于进行传热与流动过程的数值模拟，是近年来计算传热研究中的一个重要课题。