数值传热学第四章-数值计算要点.ppt

第四章 热传导 主要包含以下内容： §4.1本章的对象 §4.2 一维稳态热传导 §4.3 不稳态一维热传导 §4.4 二维与三维问题 §4.5超松弛与欠松弛 §4.6某些几何上的考虑 §4.1 本章的对象 一、本章研究对象 本章以导热问题为代表，介绍扩散方程的数值求解法。将通用微分方程中的对流项略去，整个方法的介绍将在第五章完成。 二、以导热问题的数值解作为学习起点的原因 热传导作为物理过程易于理解，而且在数学上的复杂性最小，计算方法也比较成熟； 工程流动与换热过程中的不少现象，其控制方程类似于热传导方程。如二维位势流动；常物性流体在直管内的充分发展对流换热；质扩散过程；轴承的润滑流动；某些通过多孔介质的流动。 导热问题数值解过程中所采用的一些方法与技巧对于对流问题的数值解也适用。如边界条件的处理、源项的线性化及代数方程组的求解方法等。 §4.4-3 代数方程的求解 可以采用任何一种合适的求解方法，选择方法时两者无相互影响； 在计算机程序中，可以很方便地把两个步骤分成两个独立的阶段，其中任何一段都可以单独进行修改。 这个准则是充分条件，并非必要条件。即有时虽然违背了这个条件，但迭代仍可能收敛。 对于满足正系数法则的方程，如非稳态全隐格式 , 由于各系数均大于零，且 3. ?值的选取 不存在什么选取最佳?值的一般法则，一般需要通过计算实践来确定。通常根据实际经验以及对所给定的问题作试探性的计算求得一个合适的?值。 在整个计算期间都保持相同的?值是没有必要的。 线迭代：自左向右扫描(G-S法)，不仅左边界的影响逐步传入，且在每列的直接求解中，上下端点的影响全部传入到该列的各个点上。因此，每完成一轮迭代，左、上、下边界的影响全部传入，但右边界的影响则仅前进一个网格。 交替方向隐式扫描法（ADI法）：每一轮的迭代包括了逐行与逐列的扫描，因而在每一轮迭代中，所有边界的影响均已传入区域内部，从而加快了收敛速度。 §4. 7 结 语 §4.5 超松弛与欠松弛 1. 定义： 在代数方程迭代求解的过程中，往往希望加快或是减慢前后两次迭代间因变量的变化。加快变化为超松弛，减慢变化为欠松弛。 超松弛常用于和G-S法相结合，这种组合起来的模式是所谓的持续超松弛。超松弛很少与逐行或逐列迭代法结合使用。 在强烈的非线性方程组的迭代求解过程中，采用欠松弛避免发散。 2. 方程式的写法 本次迭代所产生的TP变化 ?为松弛因子。 ?1为逐次超松弛(SOR)， ?1为逐次亚松弛(SUR)。显然，当相邻两轮的迭代值之差永远具有相同的正负号时，采用超松弛迭代可以加速收敛过程。 对于一个很小的?值，TP的变化变得很慢。 4. 通过惯量进行松弛 相应的迭代方程 i 为惯量，i 0欠松弛， i 0 超松弛。 对i 取值与相应松弛的理解 对TP产生的变化加以修正 ?1, 超松弛。左侧系数比原来小了，右侧多了一个负值项；在惯性松弛中， i 0 的结果与上述相对应，故也为超松弛。 ?1，亚松弛。左侧系数比原来大了，右侧增加了一个正项；在惯性松弛中， i 0 的结果与上述相对应，故也为亚松弛。 5. i 值的选取 不存在选取最佳i 值的一般法则，对特定的问题，可根据实际经验确定。 i 值应当与aP值不相上下，同时i 的模值越大，松弛作用也愈剧烈。 6. 稳定问题与非稳定问题求解方法上的类似 将稳定问题作为非稳定问题来求解是一种有效的方法。这时问题的初始条件就是迭代过程的初始试探值，每一时刻的解就相当于问题的一次迭代解，非稳态问题中的 相应表示上次迭代值 ， 与 i 起着相同的作用。因 恒正，相当于一种特殊的欠松弛。 ?? 越小， 越大，欠松弛作用愈强； 当 ??很大时， 很小，不再有欠松弛作用，而直接转化为稳态问题的求解。 7. 加速边界条件的影响传入计算区域的方法 Jacobi点迭代：完成一轮后，边界的影响只能传入到与边界相邻的一批节点上，即传入一个网格，且扫描的方向与收敛快慢无关。 G-S点迭代：从左到右完成一轮扫描后，左边界的影响就传到整个区域。但右、上、下边界的影响也只能传入一个网格，迭代速度比Jacobi法快，且受迭代扫描方向的影响。 §4.2-7 线性代数方程的求解 1. 求解方法 线性代数方程组的解法通常有迭代法和直接解法两种。 由于一维导热数值求解的离散方程中，待求温度仅与左右两个节点的温度有关，这样形成的代数方程组的系数矩阵将是三对角矩阵，采用追赶法（TDMA）. 2. TDMA算法要点 TDMA是一种简单、方便、