传热学4-1-屈治国.ppt

4.2.1 Taylor展开法 假定满足连续性条件，都可作Taylor展开 §4-2 内节点离散方程的建立方法 */41 两式左右两边相加，得到二阶导数的中心差分表达式 截断误差 未明确写出的级数余项中 的Δx的最低阶数为2 如果Δx= Δy */41 特点 数值解是一种近似解 各阶导数的差分表达式分子各项系数代数和为零 优点：简单、截断误差清晰 缺点（不适用）：网格不均匀，导热系数变化 */41 4.2.2 热平衡法 2.依据定律： 能量守恒定律；Fourier导热定律 1.基本思想： 对每个节点所代表的控制体列能量守恒方程式，从而得出该点与其它节点的关系式 (m,n) (m+1,n) (m,n+1) (m-1,n) (m,n-1) n e s w */41 3.具体推导 (m,n) (m+1,n) (m,n+1) (m-1,n) (m,n-1) n e s w 规定热量进入为正 */41 §4-3 边界节点离散方程建立及代数方程求解 原因： 第一类边界 ： 方程封闭 第二、三类边界：边界相邻内节点方程建立时， 引入边界未知温度，方程不封闭。 第二、三类边界： 4.3.1 边界节点离散方程的建立 边界条件表示（导入为正）： */41 存在内热源 Δx及Δy均匀 （1）平直边界的节点 x y qw (m, n) (m,n+1) (m,n-1) (m-1,n) */41 （2）外角点 x y qw (m, n) (m,n-1) (m-1,n) */41 （3）内角点 x y qw (m, n) (m,n+1) (m,n-1) (m-1,n) (m+1,n) */41 */41 （1）绝热边界， 平直边界 外部角点 内部角点 热流密度的讨论 （2）恒热流边界 （3）对流边界 阶梯逼近法 4.3.2 不规则区域的阶梯逼近法 */41 写出所有内节点和边界节点的温度差分方程 n个未知节点温度，n个代数方程式： 代数方程组的求解方法：直接解法、迭代解法 4.3.3 代数方程的求解 带状稀疏矩阵 */41 1. 直接接解法 （1） Cramer法则 */41 求解n阶线性方程组，需计算n+1个行列式， 如果取n=20，N=9.7×1020 每个行列式是n！个乘积的和， 每个乘积是n个数相乘，需要做n-1次乘法， 共需乘法次数(n+1)n!(n-1)， 1亿次乘法/s的计算机，需计算30.8万年 */41 (2) Gauss消元法 共需乘法次数n3/3+n2-n/3， 如果取n=20，N=3060 适用变量个数小于50 直接解法缺点 所需内存较大 方程数目多时不便 不适用于非线性问题 */41 假设初场 不断更新 收敛 关键：构造迭代格式 第k个迭代层次 3. 迭代解法 */41 迭代方程 */41 Jacobi迭代：迭代代入值总为上轮得到的值 */41 Gauss-Seidel迭代：将本轮最新值代入 计算机内存更省，收敛速度快 */41 是否合理？ 4. 迭代是否收敛判据 分母加小量 */41 主对角占优原则 迭代公式的选择应使每一个迭代变量的系数总是大于或等于该式中其他变量系数绝对值的代数和，此时，用迭代法求解代数方程，一定收敛 用每个方程采用导出该方程的中心节点的温度作为迭代变量，符合对角占优原则。 4. 迭代能否收敛判据 */41 */41 迭代次数 t1 t2 t3 0 1 0 3.625 0 5.675 0 3.769 2 1.735 4.545 4.996 3 1.864 4.038 5.058 4 1.983 3.980 5.013 5 2.003 3.994 5.000 6 2.0001 4.000 5.000 7 2.000 4.000 5.000 利用高斯赛德尔求解方程（例4.1） */41 迭代次数 t 0 1 2 3 4 t1 0 32 522 8722 143522 t2 0 -36 -396 -3996 -3996 t3 0 -155 -3355 -61755 -1068075 迭代发散 */41 习题4.23 （1）第一类边界条件： （2）第三类边界条件： */41 第一类边界条件 第三类边界条件 外壁面换热计算 思考：内外壁面的换热量是否相等？ * * * * * * * * * * * */35 传热学 Heat Transfer */35 传热学 Heat Transfer */35 传热学 Heat Transfer */35 传热学 He