传热学4 数值解法.pptVIP

2 代数方程的求解方法 2 ） 迭代法： 先对要计算的场作出假设（设 定初场），在迭代计算中不断予以改进，直 到计算前的假定值与计算结果相差小于允许 值为止的方法，称迭代计算收敛。 1 ） 直接解法： 通过有限次运算获得精确 解的方法，如：矩阵求解，高斯消元法。 迭代法目前应用较多的是： 1 ）高斯 —— 赛德尔迭代法： 每次迭代计算， 均是使用节点温度的最新值。 2 ）用雅可比迭代法： 每次迭代计算，均用 上一次迭代计算出的值。 设有一三元方程组 ： 11 1 12 2 13 3 1 21 1 22 2 23 3 2 31 1 32 2 33 3 3 a t a t a t b a t a t a t b a t a t a t b ? ? ? ? ? ? ? ? ? 其中 （ i=1,2,3 ； j=1,2,3 ）及 是 已知的系数（均不为零）及常数。 , i j a i b 采用高斯 —— 赛德尔迭代法的步骤： （ 1 ）将三元方程变形为迭式方程： 1 1 12 2 13 3 11 2 2 21 1 23 3 22 3 3 31 1 32 2 33 1 ( ) 1 ( ) 1 ( ) t b a t a t a t b a t a t a t b a t a t a ? ? ? ? ? ? ? ? ? （ 2 ）假设一组解（迭代初场），记为 ： 并代入迭代方程求得第一 次解 每次计算均用最新值代入。 (0) (0) (0) 1 2 3 t t t 、 、 (1) (1) (1) 1 2 3 t t t 、 、 （ 3 ）以新的初场重复计算，直到相邻两 次迭代值之差小于允许值，则称迭代收敛， 计算终止。 ) ( 1 ) 0 ( 3 13 ) 0 ( 2 12 1 11 ) 1 ( 1 t a t a b a t ? ? ? ) ( 1 ) 0 ( 3 23 ) 1 ( 1 21 2 22 ) 1 ( 2 t a t a b a t ? ? ? (1) (1) (1) 3 3 31 1 32 2 33 1 ( ) t b a t a t a ? ? ? 判断迭代是否收敛的准则： ? ? ? ? ? ? ? ? ? ? ? ? ) ( max ) ( ) 1 ( ) ( ) ( ) 1 ( ) ( ) 1 ( max max max k k i k i k i k i k i k i k i t t t t t t t t k 及 k+1 表示迭代次数； — 第 k 次迭代得到的最大值； (k) max t 当计算区域中有接近于零的 t 时，第三个较好。 3 6 10 ~ 10 E ? ? ? — 允许的偏差； 相对偏差 值一般 取 E 说明： 1 ）对于一个代数方程组，若选用的迭代方 式不合适，有可能导致发散，即称 迭代过程 发散 ； 2 ）对于常物性导热问题，组成的差分方程 组，迭代公式的选择应使一个迭代变量的系 数总是大于或等于该式中其他变量系数绝对 值的代数和，此时，结果一定收敛。 第四章 导热问题数值解基础 1 、重点内容： ① 掌握导热问题数值解法的基本思路； ② 利用热平衡法和泰勒级数展开法建立 节点的离散方程。 2 、掌握内容： 数值解法的实质。 3 、了解内容： 了解非稳态导热问题的两 种差分格式及其稳定性。 求解导热问题实际上就是对导热微分方程在 定解条件下的积分求解，从而获得分析解。随着 计算机技术的迅速发展，对物理问题进行离散求 解的数值方法发展得十分迅速，这些数值解法主 要有以下几种： （ 1 ）有限差分法 （ 2 ）有限元方法 （ 3 ）边界元方法 基本思想是把连续的定解区域用有限个离散点构成的网格来代替， 这些离散点 称作网格的节点；把连续定解区域上的连续变量的函数用在网格上定义的离散 变量函数来近似；把原方程和定解条件中的微商用差商来近似， 积分用积分和 来近似，于是原微分方程和定解条件就近似地代之以代数方程组，即有限差分 方程组 ， 解此方程组就可以得到原问题在离散点上的近似解。然后再利用插值 方法便可以从离散解得到定解问题在整个区域上的近似解。 （ 1 ）有限差分法 有限元方法的基础是变分原理和加权余量法，其基本求解思想是把计算域划分 为有限个互不重叠的单元，在每个单元内，选择一些合适的节点作为求解函数 的插值点，将微分方程中的变量改写成由各变量或其导数的节点值与所选用的 插值函数组成的线性表达式 ，借助于变分原理或加权余量法，将微分方程离 散求解。 （ 2 ）有限元法 边界元法是在有限元法之后发展起来的一种较精确有效的工程数值分析方法 。 又称边界积分方程 - 边界元法。它以定义在边界上的边界积分方程为控制方程， 通过对边界分元插值离散，化为代数方程组求解。它与基于偏